通信原理 DPSK FSK ASK PSK 检波及波形实现的MATLAB程序
MATLABPSK调制与解调及仿真
学士学位毕业设计（论文）基于MATLAB的PSK调制和解调及仿真摘要Psk调制是通信系统中最为重要的环节之一，Psk调制技术的改进也是通信系统性能提高的重要途径。
本文首先分析了数字调制系统的基本调制解调方法，然后，运用Matlab及附带的图形仿真工具——Simulink设计了这几种数字调制方法的仿真模型。
通过仿真，观察了调制解调过程中各环节时域和频域的波形，并结合这几种调制方法的调制原理，跟踪分析了各个环节对调制性能的影响及仿真模型的可靠性。
最后，在仿真的基础上分析比较了各种调制方法的性能，并通过比较仿真模型与理论计算的性能，证明了仿真模型的可行性。
另外，本文还利用Matlab的图形用户界面（GUI）功能为仿真系统设计了一个便于操作的人机交互界面，使仿真系统更加完整，操作更加方便。
关键词：数字调制；分析与仿真；Matlab；Simulink；GUI图形界面ABSTRACTIn this paper, methods of psk modulation are introduced firstly. Then their simulation models are bu ilt by using MATLAB’s simulation tool, SIMULINK. Through observing the results of simulation, the factors that affect the capability of the psk modulation system and the reliability of the simulation models are analyzed. And then, the capability of three digital modulation simulation models, 2-PSK, 4-PSK and , have been compared, as well as comparing the results of simulation and theory. At last, the conclusion is gotten: The simulation models are reasonable. In addition, an operation interface is designed, which can simplify the manipulation of the simulation system, by mean of the Graphical User Interface, which short for GUI.Keywords:PSK modulation。
Matlab的ASK_FSK__PSK程序仿真-推荐下载
2.2 数字调制技术的仿真实现
MATLAB 是一种功能强大的科学计算和工程仿真软件，它的交互式集成界面能够帮 助用户快速的完成数值分析、数字信号处理、仿真建模、和优化等功能。本课程设计需 要运用 MATLAB 编程实现 2ASK,2FSK,2PSK，2DPSK 调制解调过程，并且输出其调制后的波 形，画出频谱、功率谱密度图，并比较各种调制的误码率情况，讨论其调制效果。
2．设计原理
2.1 二进制数字调制技术原理
数字信号的传输方式分为基带传输和带通传输，在实际应用中，大多数信道具有 带通特性而不能直接传输基带信号。为了使数字信号在带通信道中传输，必须使用数字 基带信号对载波进行调制，以使信号与信道的特性相匹配。这种用数字基带信号控制载 波，把数字基带信号变换为数字带通信号的过程称为数字调制。通常使用键控法来实现 数字调制，比如对载波的振幅、频率和相位进行键控。
2.1.2、2FSK
一个ＦＳＫ信号可以看成是两个不同载波的 2ASK 信号的叠加。其解调和解调方 法和 ASK 差不多。2FSK 信号的频谱可以看成是 f1 和 f2 的两个 2ASK 频谱的组合。
2.1.3、2PSK
2PSK 以载波的相位变化作为参考基准的，当基带信号为 0 时相位相对于初始相位为 0，
3．程序与调试图形
3.1、ASK 调制解调
2.1.1、2ASK
2ASK 信号的产生方法通常有两种：模拟调制和键控法。解调有相干解调和非相干 解调。P=1 时 f(t)=Acoswt;p=0 时 f(t)=0;其功率谱密度是基带信号功率谱的线性搬移
2
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术通关，1系电过，力管根保线据护敷生高设产中技工资术艺料0不高试仅中卷可资配以料置解试技决卷术吊要是顶求指层，机配对组置电在不气进规设行范备继高进电中行保资空护料载高试与中卷带资问负料题荷试2下卷2，高总而中体且资配可料置保试时障卷，各调需类控要管试在路验最习；大题对限到设度位备内。进来在行确管调保路整机敷使组设其高过在中程正资1常料中工试，况卷要下安加与全强过，看度并22工且22作尽22下可22都能22可地护以缩1关正小于常故管工障路作高高；中中对资资于料料继试试电卷卷保破连护坏接进范管行围口整，处核或理对者高定对中值某资，些料审异试核常卷与高弯校中扁对资度图料固纸试定，卷盒编工位写况置复进.杂行保设自护备动层与处防装理腐置，跨高尤接中其地资要线料避弯试免曲卷错半调误径试高标方中高案资等，料，编试要5写、卷求重电保技要气护术设设装交备备置底4高调、动。中试电作管资高气，线料中课并敷3试资件且、设卷料中拒管技试试调绝路术验卷试动敷中方技作设包案术，技含以来术线及避槽系免、统不管启必架动要等方高多案中项；资方对料式整试，套卷为启突解动然决过停高程机中中。语高因文中此电资，气料电课试力件卷高中电中管气资壁设料薄备试、进卷接行保口调护不试装严工置等作调问并试题且技，进术合行，理过要利关求用运电管行力线高保敷中护设资装技料置术试做。卷到线技准缆术确敷指灵设导活原。。则对对：于于在调差分试动线过保盒程护处中装，高置当中高不资中同料资电试料压卷试回技卷路术调交问试叉题技时，术，作是应为指采调发用试电金人机属员一隔，变板需压进要器行在组隔事在开前发处掌生理握内；图部同纸故一资障线料时槽、，内设需，备要强制进电造行回厂外路家部须出电同具源时高高切中中断资资习料料题试试电卷卷源试切，验除线报从缆告而敷与采设相用完关高毕技中，术资要资料进料试行，卷检并主查且要和了保检解护测现装处场置理设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
BPSK和QPSK调制解调原理及MATLAB程序
2.1 PSK调制方式PSK原理介绍（以2-PSK为例）移相键控(PSK)又称为数字相位调制,二进制移相键控记作2PSK。
绝对相移是利用载波的相位(指初相)直接表示数字信号的相移方式。
二进制相移键控中,通常用相位0 和π来分别表示“0”或“1”。
2PSK 已调信号的时域表达式为s2psk(t)=s(t)cosωct, 2PSK移相键控中的基带信号与频移键控和幅度键控是有区别的,频移键控和幅度键控为单极性非归零矩形脉冲序列,移相键控为为双极性数字基带信号,就模拟调制法而言,与产生2ASK 信号的方法比较,只是对s(t)要求不同,因此2PSK 信号可以看作是双极性基带信号作用下的DSB 调幅信号。
在二进制数字调制中，当正弦载波的相位随二进制数字基带信号离散变化时，则产生二进制移相键控(2PSK)信号。
通常用已调信号载波的 0°和 180°分别表示二进制数字基带信号的 1 和 0。
二进制移相键控信号的时域表达式为e2PSK(t)=[nna g(t-nT s)]cosw c t其中, an与2ASK和2FSK时的不同，在2PSK调制中，an应选择双极性。
1, 发送概率为Pan=-1, 发送概率为1-P若g(t)是脉宽为Ts, 高度为1的矩形脉冲时，则有cosωct, 发送概率为Pe2PSK(t)=-cosωct, 发送概率为1-P由上式(6.2-28)可看出，当发送二进制符号1时，已调信号e2PSK(t)取0°相位，发送二进制符号0时，e2PSK(t)取180°相位。
若用φn表示第n个符号的绝对相位，则有0°, 发送 1 符号φn=180°, 发送 0 符号由于在2PSK信号的载波恢复过程中存在着180°的相位模糊，所以2PSK信号的相干解调存在随机的“倒π”现象，从而使得2PSK 方式在实际中很少采用。
为了解决2PSK 信号解调过程的反向工作问题， 提出了二进制差分相位键控(2DPSK)，这里不再详述。
通信原理DPSKFSKASKPSK检波及波形实现的MATLAB程序
通信原理DPSKFSKASKPSK检波及波形实现的MATLAB程序十、附录程序程序一2ASK调制解调程序及注释clear all close all i=10;%10个码元j=5000; t=linspace(0,5,j);%0-5之间产生5000个点行矢量，即分成5000份fc=10;%载波频率fm=i/5;%码元速率%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%产生基带信号x=(rand(1,i))%rand函数产生在0-1之间随机数，共1-10个figure(2) plot(x) a=round(x);%随机序列，round取最接近小数的整数figure(3) stem(a)%火柴梗状图st=t; for n=1:10 if a(n) for m=j/i*(n-1)+1:j/i*n st(m)=0; end else for m=j/i*(n-1)+1:j/i*n st(m)=1; end end end figure(1); subplot(421); plot(t,st); axis([0,5,-1,2]); title('基带信号st'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%载波s1=cos(2*pi*fc*t); subplot(422); plot(s1); title('载波信号s1'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%调制e_2ask=st.*s1; subplot(423); plot(t,e_2ask); title('已调信号'); noise =rand(1,j); e_2ask=e_2ask+noise;%加入噪声subplot(424);plot(t,e_2ask); title('加入噪声的信号'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解调at=e_2ask.*cos(2*pi*fc*t); at=at-mean(at);%因为是单极性波形，还有直流分量，应去掉subplot(425); plot(t,at); title('与载波相乘后信号'); [f,af] = T2F(t,at);%通过低通滤波器[t,at] = lpf(f,af,2*fm); subplot(426); plot(t,at);title('相干解调后波形'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽样判决for m=0:i-1; if at(1,m*500+250)+0.5 for j=m*500+1:(m+1)*500; at(1,j)=0; end else for j=m*500+1:(m+1)*500; at(1,j)=1; end end end subplot(427); plot(t,at); axis([0,5,-1,2]); title('抽样判决后波形')程序二2FSK调制解调程序及注释clear all close all i=10;%基带信号码元数j=5000; a=round(rand(1,i));%产生随机序列t=linspace(0,5,j); f1=10;%载波1频率f2=5;%载波2频率fm=i/5;%基带信号频率%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%产生基带信号st1=t; for n=1:10 if a(n) for m=j/i*(n-1)+1:j/i*n st1(m)=0; end else for m=j/i*(n-1)+1:j/i*n st1(m)=1; end end end st2=t; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%基带信号求反for n=1:j; if st1(n) st2(n)=0; else st2(n)=1; end end; figure(1);subplot(411); plot(t,st1); title('基带信号st1'); axis([0,5,-1,2]); subplot(412); plot(t,st2); title('基带信号反码st2'); axis([0,5,-1,2]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%载波信号s1=cos(2*pi*f1*t) s2=cos(2*pi*f2*t) subplot(413),plot(s1); title('载波信号s1'); subplot(414),plot(s2); title('载波信号s2'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%调制F1=st1.*s1;%加入载波 1 F2=st2.*s2;%加入载波 2 figure(2); subplot(411); plot(t,F1); title('F1=s1*st1'); subplot(412); plot(t,F2); title('F2=s2*st2'); e_fsk=F1+F2; subplot(413); plot(t,e_fsk); title('2FSK信号')%键控法产生的信号在相邻码元之间相位不一定连续nosie=rand(1,j); fsk=e_fsk+nosie; subplot(414); plot(t,fsk); title('加噪声后信号') %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解调st1=fsk.*s1;%与载波1相乘[f,sf1] = T2F(t,st1);%通过低通滤波器[t,st1] = lpf(f,sf1,2*fm); figure(3); subplot(311); plot(t,st1); title('与s1相乘后波形'); st2=fsk.*s2;%与载波2相乘[f,sf2] = T2F(t,st2);%通过低通滤波器[t,st2] = lpf(f,sf2,2*fm); subplot(312); plot(t,st2); title('与s2相乘后波形'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽样判决for m=0:i-1; if st1(1,m*500+250)st2(1,m*500+250); for j=m*500+1:(m+1)*500; at(1,j)=0; end else forj=m*500+1:(m+1)*500; at(1,j)=1; end end end; subplot(313); plot(t,at); axis([0,5,-1,2]); title('抽样判决后波形')程序三2PSK调制解调程序及注释clear all close all i=10; j=5000; fc=4;%载波频率fm=i/5;%码元速率B=2*fm; t=linspace(0,5,j); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%产生基带信号a=round(rand(1,i));%随机序列,基带信号figure(3); stem(a); st1=t; for n=1:10 if a(n) for m=j/i*(n-1)+1:j/i*n st1(m)=0; end else for m=j/i*(n-1)+1:j/i*n st1(m)=1; end end end figure(1); subplot(411); plot(t,st1); title('基带信号st1'); axis([0,5,-1,2]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%基带信号求反%由于PSK中的是双极性信号，因此对上面所求单极性信号取反来与之一起构成双极性码st2=t; for k=1:j; if st1(k) st2(k)=0; else st2(k)=1; end end; subplot(412); plot(t,st2); title('基带信号反码st2'); axis([0,5,-1,2]); st3=st1-st2; subplot(413); plot(t,st3); title('双极性基带信号st3'); axis([0,5,-2,2]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%载波信号s1=sin(2*pi*fc*t); subplot(414); plot(s1); title('载波信号s1'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%调制e_psk=st3.*s1; figure(2); subplot(511); plot(t,e_psk);title('e_2psk'); noise=rand(1,j); psk=e_psk+noise;%加入噪声subplot(512); plot(t,psk); title('加噪后波形'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解调psk=psk.*s1;%与载波相乘subplot(513);plot(t,psk); title('与载波s1相乘后波形'); [f,af] = T2F(t,psk);%%%%%%%%%%%通过低通滤波器[t,psk] = lpf(f,af,B); subplot(514); plot(t,psk); title('低通滤波后波形'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽样判决for m=0:i-1; if psk(1,m*500+250) for j=m*500+1:(m+1)*500; psk(1,j)=0; end else for j=m*500+1:(m+1)*500; psk(1,j)=1; end end end subplot(515); plot(t,psk); axis([0,5,-1,2]); title('抽样判决后波形')程序四2DPSK调制解调程序及注释clear all close all i=10; j=5000; fc=4;%载波频率fm=i/5;%码元速率B=2*fm; t=linspace(0,5,j); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%产生基带信号a=round(rand(1,i)); figure(4); stem(a); st1=t; for n=1:10 if a(n) for m=j/i*(n-1)+1:j/i*n st1(m)=0; end else for m=j/i*(n-1)+1:j/i*n st1(m)=1; end end end figure(1); subplot(321); plot(t,st1); title('绝对码'); axis([0,5,-1,2]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%差分变换b=zeros(1,i);%%%%%%%%全零矩阵b(1)=a(1); for n=2:10 if a(n)if b(n-1)=1 b(n)=0; elseb(n)=1; end else b(n)=b(n-1); end end st1=t; for n=1:10 if b(n) for m=j/i*(n-1)+1:j/i*n st1(m)=0; end else for m=j/i*(n-1)+1:j/i*n st1(m)=1; end end end subplot(323); plot(t,st1); title('相对码st1'); axis([0,5,-1,2]); st2=t; for k=1:j; if st1(k) st2(k)=0; else st2(k)=1; end end; subplot(324); plot(t,st2); title('相对码反码st2'); axis([0,5,-1,2]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%载波信号s1=sin(2*pi*fc*t); subplot(325); plot(s1); title('载波信号s1'); s2=sin(2*pi*fc*t+pi); subplot(326); plot(s2); title('载波信号s2'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%信号调制d1=st1.*s1; d2=st2.*s2; figure(2); subplot(411); plot(t,d1); title('st1*s1'); subplot(412); plot(t,d2); title('st2*s2'); e_dpsk=d1+d2; subplot(413); plot(t,e_dpsk); title('调制后波形'); noise=rand(1,j); dpsk=e_dpsk+noise;%%%%%%%%%%加入噪声subplot(414); plot(t,dpsk); title('加噪声后信号');%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解调dpsk=dpsk.*s1;%%%%与载波s1相乘figure(3); subplot(411); plot(t,dpsk); title('与载波相乘后波形'); [f,af]=T2F(t,dpsk);%%%%通过低通滤波器[t,dpsk]=lpf(f,af,B); subplot(412); plot(t,dpsk); title('低通滤波后波形'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽样判决st=zeros(1,i);%%全零矩阵for m=0:i-1; if dpsk(1,m*500+250)st(m+1)=0; for j=m*500+1:(m+1)*500; dpsk(1,j)=0; end else for j=m*500+1:(m+1)*500; st(m+1)=1; dpsk(1,j)=1; end end end subplot(413); plot(t,dpsk); axis([0,5,-1,2]); title('抽样判决后波形') %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%码反变换dt=zeros(1,i);%%全零矩阵dt(1)=st(1); for n=2:10; if (st(n)-st(n-1))=0(st(n)-st(n-1)) dt(n)=0; else dt(n)=1; end end st=t; for n=1:10 if dt(n) for m=j/i*(n-1)+1:j/i*n st(m)=0; end else for m=j/i*(n-1)+1:j/i*n st(m)=1; end end end subplot(414); plot(t,st); axis([0,5,-1,2]); title('码反变换后波形');。
DPSK调制解调的Matlab程序
%DPSK调制解调程序清单%仿真DPSK信号的产生及相干解调过程%将成形滤波器系数量化为10比特后，写入tra_lpf.coe文件中%将接收滤波器系数量化为10比特后，写入rec_lpf.coe文件中%将DPSK已调数据量化为8比特后，写入Dpsk.txt文件中ps=1*10^6; %码速率为1MHza=0.8; %成形滤波器系数为0.8B=(1+a)*ps; %中频信号处理带宽Fs=8*10^6; %采样速率为8MHzfc=2*10^6; %载波频率为2MHzN=20000; %仿真数据的长度t=0:1/Fs:(N*Fs/ps-1)/Fs;%产生长度为N,频率为fs的时间序列s=randint(N,1,2); %产生随机数据作为原始数据,%并将绝对码变换为相对码ds=ones(1,N);for i=2:Nif s(i)==1ds(i)=-ds(i-1);elseds(i)=ds(i-1);endend%对相对码数据以Fs频率采样Ads=upsample(ds,Fs/ps);%设计平方升余弦滤波器n_T=[-2 2];rate=Fs/ps;T=1;Shape_b = rcosfir(a,n_T,rate,T);%figure(4);freqz(Shape_b)%对采样后的数据进行升余弦滤波;rcos_Ads=filter(Shape_b,1,Ads);%产生载频信号f0=sin(2*pi*fc*t);%产生DPSK已调信号dpsk=rcos_Ads.*f0;%与相干载波相乘，实现相干解调demod_mult=dpsk.*f0;%设计接收端低通滤波器fc=[ps 3.1*10^6]; %过渡带mag=[1 0]; %窗函数的理想滤波器幅度—dev=[0.01 0.01]; %纹波[n,wn,beta,ftype]=kaiserord(fc,mag,dev,Fs) %获取凯塞窗参数fpm=[0 fc(1)*2/Fs fc(2)*2/Fs 1]; %firpm函数的频段向量magpm=[1 1 0 0]; %firpm函数的幅值向量rec_lpf=firpm(n,fpm,magpm); %firpm函数返回的最优滤波器系数%对乘法运算后的数据进行低通滤波，输出解调后的基带信号demod_lpf=filter(rec_lpf,1,demod_mult);figure(1)%绘制成形滤波后信号频谱、DPSK信号频谱、DPSK信号时域波形m_rcos_Ads=20*log10(abs(fft(rcos_Ads,1024)));m_rcos_Ads=m_rcos_Ads-max(m_rcos_Ads);m_dpsk=20*log10(abs(fft(dpsk,1024)));m_dpsk=m_dpsk-max(m_dpsk);%设置幅频响应的横坐标单位为MHzx_f=[0:(Fs/length(m_dpsk)):Fs/2];x_f=x_f/10^6;%只显示正频率部分的幅频响应mrcos_Ads=m_rcos_Ads(1:length(x_f));mdpsk=m_dpsk(1:length(x_f));%设置时域波表的横坐标单位为usLen=100;%设置时域波形显示的点数x_t=1:Len;%产生长度为Len的时间序列x_t=x_t/Fs*10^6;%显示所需的频谱及时域波形subplot(311); plot(x_f,mrcos_Ads);legend('成形滤波后信号频谱');xlabel('频率(MHz)');ylabel('幅度(dB)');grid on;subplot(312); plot(x_f,mdpsk);legend('DPSK已调信号频谱');xlabel('频率(MHz)');ylabel('幅度(dB)');grid on;subplot(313); plot(x_t,dpsk(101:Len+100));grid on;legend('DPSK时域信号波形');xlabel('时间(us)');ylabel('幅度(V)');grid on;figure(2)%对相对码重复Fs/ps倍采样，便于绘图比较s_Ads=rectpulse(ds,Fs/ps);Len=500;%设置时域波形显示的点数x_t=1:Len;%产生长度为Len的时间序列x_t=x_t/Fs*10^6;%绘制DPSK解调前后时域波形delay=18;%为便于观察，对解调后的基带波形超前显示delay个点。
通信原理课程设计基于MATLAB的2PSK和2FSK调制仿真
通信原理A课程设计报告题目：基于MATLAB的2PSK和2FSK调制仿真院系：自动化与信息工程学院专业：通信工程班级：学号：姓名：指导教师：职称：讲师2012年12月24日－2012年12月28日一、设计任务编写2PSK和2FSK调制程序，任意给定一组二进制数，计算经过这两种调制方式的输出信号。
程序书写要规范，加必要的注释；经过程序运行的调制信号波形要与理论计算出的波形一致。
分步实施：1 ）熟悉2PSK和2FSK调制原理；2 ）编写2PSK和2FSK调制程序；3 ）画出原信号和调制信号的波形图。
课程设计的最后成果是提交一份实验报告，内容包括：1)2PSK和2FSK调制原理；对给定信号画出理论调制波形；2)程序设计思想，画出流程图；3)源程序代码（需打印）；4)测试结果（需打印）和理论计算结果对比是否一致；5)小结。
六、参考文献【1】冯象初,甘小冰. 数值泛函与小波理论西安：西安电子科技大学出版社,2003.5【2】樊昌信,曹丽娜. 通信原理（第六版）北京：国防工业出版社, 2010.6【3】罗建军，扬琦.精讲多练MATLAB（第2版）西安：西安交通大学出版社，2009.7附录：源程序代码clear allclose alli=10; %基带信号码元数j=5000;a=round(rand(1,i)); %产生随机序列t=linspace(0,5,j);f1=4; %2FSK载波1频率 2PSK载波频率f2=8; %2FSK载波2频率fm=i/5; %基带信号频率%%%%%%%%%%产生基带信号st1=t;for n=1:10if a(n)<1;for m=j/i*(n-1)+1:j/i*nst1(m)=0;endelsefor m=j/i*(n-1)+1:j/i*nst1(m)=1;endend如有你有帮助，请购买下载，谢谢！endfigure(1);subplot(311);plot(t,st1);title('基带信号st1');axis([0,5,-1,2]);%%%%%%%%%%基带信号求反st2=t;for n=1:j;if st1(n)>=1;st2(n)=0;elsest2(n)=1;endend;%%%%%%%%%%构成双极性码st3=st1-st2;%%%%%%%%%%载波信号s1=sin(2*pi*f1*t)s2=sin(2*pi*f2*t)%subplot(312),plot(s1);%title('载波信号s1');%subplot(313),plot(s2);%title('载波信号s2');%%%%%%%%%%%调制%figure(2);F1=st1.*s1; %加入载波1 (2FSK)F2=st2.*s2; %加入载波2 (2FSK)e_fsk=F1+F2;subplot(312);plot(t,e_fsk);title('2FSK调制信号');e_psk=st3.*s1; %加入载波 (2PSK)subplot(313);plot(t,e_psk);title('2PSK调制信号');如有你有帮助，请购买下载，谢谢！四、程序运行结果及分析00.51 1.52 2.53 3.54 4.55-112基带信号st100.51 1.52 2.53 3.54 4.55-112FSK 调制信号00.51 1.52 2.53 3.54 4.55-1012PSK 调制信号。
ask、psk、fsk的调制与解调原理
调制和解调是现代通信系统中至关重要的过程，它们可以实现信息的传输和接收。
在数字通信中，有三种常见的调制和解调技术，分别是ask、psk和fsk。
本文将详细讨论这三种调制和解调技术的原理和应用。
一、ASK调制与解调原理1. ASK调制ASK（Amplitude Shift Keying）调制是一种将数字信号转换为模拟信号的调制技术。
在ASK调制中，数字信号被用来控制载波的振幅，当输入信号为1时，振幅为A；当输入信号为0时，振幅为0。
ASK 调制一般用于光纤通信和无线电通信系统。
2. ASK解调ASK解调是将接收到的模拟信号转换为数字信号的过程。
它通常是通过比较接收到的信号的振幅与阈值来实现的。
当信号的振幅高于阈值时，输出为1；当信号的振幅低于阈值时，输出为0。
ASK解调在数字通信系统中有着广泛的应用。
二、PSK调制与解调原理1. PSK调制PSK（Phase Shift Keying）调制是一种将数字信号转换为模拟信号的调制技术。
在PSK调制中，不同的数字信号会使载波的相位发生变化。
常见的PSK调制方式有BPSK（Binary Phase Shift Keying）和QPSK（Quadrature Phase Shift Keying）。
PSK调制在数字通信系统中具有较高的频谱效率和抗噪声性能。
2. PSK解调PSK解调是将接收到的模拟信号转换为数字信号的过程。
它通常是通过比较接收到的信号的相位与已知的相位来实现的。
PSK解调需要根据已知的相位来判断传输的是哪个数字信号。
PSK调制技术在数字通信系统中被广泛应用，特别是在高速数据传输中。
三、FSK调制与解调原理1. FSK调制FSK（Frequency Shift Keying）调制是一种将数字信号转换为模拟信号的调制技术。
在FSK调制中，不同的数字信号对应着不同的载波频率。
当输入信号为1时，载波频率为f1；当输入信号为0时，载波频率为f2。
FSK调制常用于调制通联方式线路和调制调制解调器。
通信原理课程设计(基于MATLAB的PSK,DPSK仿真)
通讯原理课程设计报告题目鉴于Matlab的2PSK,2DPSK仿真学院电子信息工程学院专业学生姓名学号年级指导教师职称讲师2013年12月20日设计报成功绩（依据优、良、中、及格、不及格评定）指导教师考语：指导教师（署名）年月日说明：指导教师评分后，设计报告交院实验室保留。
鉴于 Matlab 的 2PSK,2DPSK仿真专业：学号：学生：指导老师：纲要：现代通讯系统要求通讯距离远、通讯容量大、传输质量好，作为其重点技术之一的调制技术向来是研究的一个重要方向。
本设计主要表达了数字信号的调制方式，介绍了2PSK数字调制方式的基来源理，功率谱密度，并运用MATLAB软件对数字调制方式2PSK进行了编程仿真切现，在MATLAB 平台上成立2PSK和 2DPSK调制技术的仿真模型。
进一步学习了MATLAB编程软件，将 MATLAB与通讯系统中数字调制知识联系起来，为此后在通讯领域学习和研究打下了基础在计算机上，运用MATLAB 软件来实现对数字信号调制技术的仿真。
重点词：数字调制与解调；MA TLAB ； 2PSK； 2DPSK ；目录第 1 章绪论 (1)1.1 调制方式 (1)1.2 设计要求 (1)设计内容 (1)设计仪器 (1)第 2 章 2PSK,2DPSK原理 (2)2.1 2PSK 原理 (2)2PSK 基来源理 . (2)2PSK 调制原理 . (2)2PSK 解调原理 . (3)2.2 2DPSK 原理 (4)2DPSK 基来源理 . (4)2DPSK 调制原理 . (5)2DPSK 解调原理 . (6)第 3 章实验过程 (8)3.1 2PSK 仿真部分 (8)2PSK 仿真图 . (8)2PSK 模块的参数设置： . (8)3.2 2DPSK 仿真部分 (9)2DPSK 仿真图 . (9)2DPSK 模块的参数设置： . (10)第 4 章仿真结果 (15)4.1 2PSK 仿真结果 (15)4.2 2DPSK 仿真结果 (15)总结 . (16)参照文件 . (17)道谢 . (18)第1章绪论1.1调制方式数字通讯系统 ,按调制方式能够分为基带传输和带通传输。
psk和fsk调制与解调matlab
clc;clear;close all;%%psk调制与解调M=4;N=20;fc=100;fs=2000;Rb=40;dphi=2*pi/M;phi0=pi/M;%%初相位Ts=fs/Rb;an=randint(1,N,M);T=length(an)/Rb;t=0:1/fs:T-1/fs;for i=1:length(an)for k=1:Tsst((i-1)*Ts+k)=an(i);endendfigure,subplot(211),plot(st);axis([1 length(st) -0.1 3.1]),title('原始信号st');sig=cos(2*pi*fc*t+phi0+st.*dphi);subplot(212),plot(sig);axis([1 length(sig) -1.1 1.1]),title('PSK调制信号sig'); %%解调sig=hilbert(sig);sig = sig.*exp(-j*2*pi*fc*t);%%频谱搬移phi=angle(-sig)-pi/M;ax=max(phi);in=min(phi);sg=(ax-in)/(M-1);for i=1:Nphi1=phi((i-1)*Ts+Ts/2);if phi1>in-0.5*sg& phi1<in+0.5*sgbsg(i)=0;elseif phi1<in+1.5*sg& phi1>in+0.5*sgbsg(i)=1;elseif phi1>in+1.5*sg& phi1<in+2.5*sgbsg(i)=2;elsebsg(i)=3;endendfigure,stem(bsg);%%%%%Fsk调制与解调clc;clear;close all;M=4;N=20;% A=2;fc=100;fs=2000;Rb=40;f=40;Ts=fs/Rb;an=randint(1,N,M);T=length(an)/Rb;t=0:1/fs:T-1/fs;for i=1:length(an)for j=1:Tsst((i-1)*Ts+j)=an(i);endendfigure,subplot(211),plot(st);axis([1 length(st) -0.1 3.1]),title('原始信号st');sig=cos(2*pi*(fc+st*f).*t);subplot(212),plot(sig);axis([1 length(sig) -1.1 1.1]),title('FSK调制信号sig'); %%%解调sig=hilbert(sig);phi=unwrap(angle(sig));%%求相位fsig = diff(phi);%%求频率。
《通信原理实验》ASK、PSK、BFSK等实验报告
《通信原理实验》ASK、PSK、BFSK等实验报告《通信原理》实验报告⼀、实验⽬的1、掌握⽤键控法产⽣ASK、FSK信号的⽅法。
2、掌握ASK、FSK⾮相⼲解调的原理。
3、掌握BFSK调制和解调的基本原理。
4、掌握BFSK数据传输过程，熟悉典型电路。
5、了解数字基带波形时域形成的原理和⽅法，掌握滚降系数的概念。
6、熟悉BPSK调制载波包络的变化。
7、掌握BFSK载波恢复特点与位定时恢复的基本⽅法。
⼆、实验器材1、主控&信号源模块，9号、13号模块各⼀块2、双踪⽰波器⼀台3、连接线若⼲三、实验原理1、ASK调制及解调实验原理框图2、FSK调制及解调实验原理框图3、BPSK调制及解调实验原理框图四、实验步骤实验项⽬⼀ASK调制１、分别观测调制输⼊和调制输出信号：以9号模块TH1为触发，⽤⽰波器同时观测9号模块TH1和模块TH4，验证ASK调制原理。
调制输⼊信号和调制输出信号：由图可知，当输⼊为“1”时，输出为正弦信号；输⼊为“0”时，输出信号为0。
注：CH1（上⾯的波形）为调制输⼊信号，CH2（下⾯的波形）为调制输出信号。
调制输⼊信号频谱：调制输出信号频谱：2、将PN序列输出频率改为64KHz，观察载波个数是否发⽣变化。
调制输⼊信号和调制输出信号：将图与题1中的图作⽐较，可以发现，PN序列的输出频率改为64KHz时，载波的个数没有发⽣变化。
可以得出，ASK调制时，PN序列输出频率的改变，不会对载波产⽣影响。
注：CH1（上⾯的波形）为调制输⼊信号，CH2（下⾯的波形）为调制输出信号。
调制输⼊信号频谱：调制输出信号频谱：实验项⽬⼆ASK解调1、对⽐观测调制信号输⼊以及解调输出：以9号模块TH1为触发，⽤⽰波器同时观测9号模块TH1和TH6，调节W1直⾄⼆者波形相同；再观测TP4（整流输出）、TP5（LPF-ASK）两个中间过程测试点，验证ASK解调原理。
解调信号输⼊和解调输出：整流输出和LPF-ASK：注：CH1（上⾯的波形）为调制输⼊信号，CH2（下⾯的波形）为调制输出信号；CH1（上⾯的波形）为整流输出，CH2（下⾯的波形）为LPF-ASK从调制输⼊信号和输出信号的波形对⽐来看，两个的波形⼀致，但是存在这相位差。
dpsk信号matlab实现
clear all；close all;fs=4000000; %设定系统的抽样频率k=20000；％设定数字基带信号的频率fc=200000；％设定正弦载波频率t=0:1/fs：4000/fs；%仿真时间范围p=21；s=randint（1,p,2); ％设定需要产生的码元个数m=s（ceil（k＊t+0.01）)；%将基带生成时域信号figure（1) subplot（311）plot(t,m）；axis（[0 10e-4 —0。
2 1。
2］）;grid on；title（’数字基带信号'）;b=randint(1，p，2);％将生成的基带转换为差分码for i=1：pif (i==1)if (s（i)==0)b(i）=0;elseb(i）=1；endelseif （s（i)==b(i-1））b(i）=0；elseb（i）=1；endendn=b（ceil(k*t+0。
01）); %将差分码生成时域信号subplot（312)plot(t,n）；axis([0 10e-4 -0.2 1。
2］）;grid on;title(’差分码'）x=(n-0.5).*2car=sin(2＊pi*fc＊t)；％定义载波dpsk=x。
*car; ％2dpsk信号的载波调制subplot（313);plot（t，car);axis（［0 10e—4 —1.2 1.2］)；title('正弦载波’）;figure（2) subplot（311)；plot（t，dpsk);axis(［0 10e-4 -1.2 1。
2]);title('2DPSK信号'）;grid on；vn=0。
05；noise=vn.*(randn(size（t））)；％产生噪音subplot（312）；plot(t，noise）;grid on;title（’噪音信号')；axis（[0 10e—4 —0。
基于MATLAB实现的PSK和DPSK信号系统仿真设计
文档从互联网中收集，已重新修正排版，word格式支持编辑，如有帮助欢迎下载支持。
目录摘要 ................................................................................................................... 错误!未定义书签。
Abstract............................................................................................................. 错误!未定义书签。
1.设计目的与要求 ............................................................................................ 错误!未定义书签。
2.方案的选择 .................................................................................................... 错误!未定义书签。
2.1调制部分............................................................................................. 错误!未定义书签。
2.2解调部分............................................................................................. 错误!未定义书签。
3.单元电路原理和设计 .................................................................................... 错误!未定义书签。
通信原理实验教程MATLAB
通信原理实验教程MATLAB通信原理是一个非常重要的学科，它涉及到人类社会中所有的信息传递和交流。
在通信原理实验中，MATLAB是一个广泛应用的软件工具，它可以帮助我们更好地理解和分析各种信号处理、调制和调解技术。
本文将介绍通信原理实验教程MATLAB，包括实验内容、实验步骤和实验效果。
一、实验内容通信原理实验教程MATLAB主要包括以下内容：1.信号处理实验：涉及对不同类型的信号进行采样、量化、编码、解码和滤波等处理。
2.模拟调制实验：涉及常见的调制方式，如AM、FM、PM等，通过信号模拟生成带载波的调制信号。
3.数字调制实验：涉及数字调制方式，如ASK、FSK、PSK 等，通过二进制数字信号生成带载波的数字调制信号。
4.调解实验：涉及不同的调解技术，如干扰消除、正交调解等，可以帮助我们更好地理解信号的编码和解码过程。
二、实验步骤通信原理实验教程MATLAB的实验步骤如下：1.信号处理实验：利用MATLAB实现信号处理算法，包括采样、量化、编码、解码和滤波等过程。
通过图形界面展示处理后的信号波形和频谱，来验证算法的正确性。
2.模拟调制实验：利用MATLAB生成正弦波载波和调制信号，利用Mod函数进行AM、FM和PM模拟调制操作。
通过图形界面展示调制信号的波形和频谱，以及载波和调制信号的相位、频率和幅值信息。
3.数字调制实验：生成二进制数字信号，利用MATLAB实现ASK、FSK、PSK等数字调制算法，利用图形界面展示调制信号的波形和频谱，以及载波和数字信号的相位和频率信息。
4.调解实验：利用MATLAB实现干扰消除和正交调解等调解技术，通过图形界面展示编码和解码过程的波形和频谱信息，并比较不同技术之间的性能差异。
三、实验效果通信原理实验教程MATLAB的实验效果如下：1.信号处理实验：通过MATLAB实现信号处理算法，可以快速准确地分析和优化不同类型的信号，从而保证通信系统的稳定性和可靠性。
2.模拟调制实验：通过MATLAB模拟AM、FM和PM模拟调制操作，可以深入了解不同调制方式的优缺点以及应用场景，从而更好地选择调制方式。
通信原理 matlab代码
通信原理matlab代码以下为一篇关于通信原理的1500-2000字文章，主要讨论了使用MATLAB编写通信原理相关代码的步骤。
通信原理是现代通信系统的基础，通过使用各种信号处理技术和模型，使信息能够在发送和接收设备之间进行传输和解码。
在通信原理的研究中，MATLAB是一种常用的工具，它提供了一系列函数和工具箱，可用于模拟和分析各种通信系统。
在开始使用MATLAB编写通信原理相关代码之前，首先需要安装MATLAB软件。
安装完成后，可以打开MATLAB环境，通过编写和执行代码实现通信原理的模拟和分析。
第一步是导入需要用到的工具箱和相关函数。
MATLAB提供了多个通信系统工具箱，包括通信系统工具箱、数字信号处理工具箱和无线通信工具箱等。
使用`import`关键字可以导入特定的工具箱和函数。
matlabimport commsys.*;import dsp.*;import wireless.*;第二步是设置通信系统的参数。
通信系统的参数包括发送和接收设备的特性、信道特性和数据传输速率等。
通过定义这些参数，可以在整个通信系统中进行一致的模拟和分析。
matlabsampling_frequency = 10e6; 采样频率carrier_frequency = 2.4e9; 载波频率data_rate = 1e6; 数据传输速率snr = 10; 信噪比第三步是生成发送信号。
发送信号的生成主要是根据具体的通信系统规范和调制方案来确定。
比如，要生成一个使用QPSK调制方案的信号，可以使用`comm.QPSKModulator`工具箱生成调制器，并通过调用`step`函数传入待调制的数据。
matlabmodulator = comm.QPSKModulator();data = randi([0 3], 1000, 1); 生成随机的待调制数据modulated_signal = step(modulator, data);第四步是加入信道特性和噪声。
PSK原理及MATLAB仿真
基基 信信
Uc(t)
t
再码1 0 1 10 0
未调 载滤
S(t)
t
绝相 调鉴
t
图3绝对移相的波形示意图
在同步解调的PSK系统中，由于收端载波恢复 存在相位含糊的问题，即恢复的载波可能与未调载 波同相，也可能反相，以至使解调后的信码出现 “0”、“1”倒置，发送为“1”码，解调后得到 “0”码；发送为“0”码，解调后得到“1”码。 这是我们所不希望的，为了克服这种现象，人们提 出了相对移相方式。 相对移相的调制规律是：每一个码元的载波相 位不是以固定的未调载波相位作基准的，而是以相 邻的前一个码元的载波相位来确定其相位的取值。 例如，当某一码元取“1”时，它的载波相位与前 一码元的载波同相；码元取“0”时，它的载波相 位与前一码元的载波反相。相对移相的波形如下图 所示。
= A cos( wc t + Ψ n)
其中， Ψ n表示第n个符号的绝对相位：
Ψ
n
=
{
0发送‘1’ 时 π发送‘0’ 时
• 2PSK信号波形如图5所示
图5
解调
2PSK系统的解调部分框图如下图所示
Ud1
鉴鉴鉴
低低
再再再Um2 PSK 入压压 振振鉴Ud
环环 滤滤鉴
Ud
模模 鉴相鉴
差差 译再
Um1 Ud2
2PSK调制
2PSK系统的调制部分框图如下图所示
M序序
发再鉴 差差 编再 调鉴
PSK
1MHZ 振振
外载滤 放放鉴
异 步
同 步
10 差 频
2差频
外
内
10MHZ 晶振
图2 2PSK调制部分框图
1. M序列发生器 实际的数字基带信号是随机的，为了实验 和测试方便，一般都是用M序列发生器产生 一个伪随机序列来充当数字基带信号源。按 照本原多项式f(x)=X5+X3+1组成的五级线性 移位寄存器，就可得到31位码长的M序列。 码元定时与载波的关系可以是同步的，以 便清晰观察码元变化时对应调制载波的相应 变化；也可以是异步的，因为实际的系统都 是异步的，码元速率约为1Mbt/s。
matlab设计DPSK信号调制与解调
2DPSK调制与解调系统的仿真1、 2DPSK基本原理2DPSK信号原理2DPSK方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。
现假设用Φ表示本码元初相与前一码元初相之差，并规定：Φ＝0表示0码，Φ＝π表示1码。
则数字信息序列与2DPSK信号的码元相位关系可举例表示如2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的，在接收端只能采用相干解调，它的时域波形图如图所示。
图 2DPSK信号在这种绝对移相方式中，发送端是采用某一个相位作为基准，所以在系统接收端也必须采用相同的基准相位。
如果基准相位发生变化，则在接收端回复的信号将与发送的数字信息完全相反。
所以在实际过程中一般不采用绝对移相方式，而采用相对移相方式。
定义为本码元初相与前一码元初相之差，假设：数字信息“0”；→数字信息“1”。
→则数字信息序列与2DPSK信号的码元相位关系可举例表示如下：数字信息： 1 0 1 1 0 1 1 1 0 1DPSK信号相位：或：0 ? 2DPSK信号的调制原理一般来说，2DPSK信号有两种调试方法，即模拟调制法和键控法。
2DPSK 信号的的模拟调制法框图如图1.2.1所示，其中码变换的过程为将输入的单极性不归零码转换为双极性不归零码。
．s(t)(t)相码变载图1.2.1 模拟调制法2DPSK信号的的键控调制法框图如图1.2.2所示，其中码变换的过程为将输入的基带信号差分，即变为它的相对码。
选相开关作用为当输入为数字信息“0”时接相位0，当输入数字信息为“1”时接pi。
图1.2.2 键控法调制原理图2DPSK信号的解调原理2DPSK信号最常用的解调方法有两种，一种是极性比较和码变换法，另一种是差分相干解调法。
2DPSK信号解调的极性比较法它的原理是2DPSK信号先经过带通滤波器，去除调制信号频带以外的在信道中混入的噪声，再与本地载波相乘，去掉调制信号中的载波成分，再经过低通滤波器去除高频成分，得到包含基带信号的低频信号，将其送入抽样判决器中进行抽样判决的到基带信号的差分码，再经过逆差分器，就得到了基带信号。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
十、附录程序程序一2ASK调制解调程序及注释clear allclose alli=10;%10个码元j=5000;t=linspace(0,5,j);%0-5之间产生5000个点行矢量，即分成5000份fc=10;%载波频率fm=i/5;%码元速率%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%产生基带信号x=(rand(1,i))%rand函数产生在0-1之间随机数，共1-10个figure(2)plot(x)a=round(x);%随机序列，round取最接近小数的整数figure(3)stem(a)%火柴梗状图st=t;for n=1:10if a(n)<1;for m=j/i*(n-1)+1:j/i*nst(m)=0;endelsefor m=j/i*(n-1)+1:j/i*nst(m)=1;endendendfigure(1);subplot(421);plot(t,st);axis([0,5,-1,2]);title('基带信号st'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%载波s1=cos(2*pi*fc*t);subplot(422);plot(s1);title('载波信号s1'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%调制e_2ask=st.*s1;subplot(423);plot(t,e_2ask);title('已调信号');noise =rand(1,j);e_2ask=e_2ask+noise;%加入噪声subplot(424);plot(t,e_2ask);title('加入噪声的信号'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解调at=e_2ask.*cos(2*pi*fc*t);at=at-mean(at);%因为是单极性波形，还有直流分量，应去掉subplot(425);plot(t,at);title('与载波相乘后信号');[f,af] = T2F(t,at);%通过低通滤波器[t,at] = lpf(f,af,2*fm);subplot(426);plot(t,at);title('相干解调后波形'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽样判决for m=0:i-1;if at(1,m*500+250)+0.5<0.5;for j=m*500+1:(m+1)*500;at(1,j)=0;endelsefor j=m*500+1:(m+1)*500;at(1,j)=1;endendendsubplot(427);plot(t,at);axis([0,5,-1,2]);title('抽样判决后波形')程序二2FSK调制解调程序及注释clear allclose alli=10;%基带信号码元数j=5000;a=round(rand(1,i));%产生随机序列t=linspace(0,5,j);f1=10;%载波1频率f2=5;%载波2频率fm=i/5;%基带信号频率%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%产生基带信号st1=t;for n=1:10if a(n)<1;for m=j/i*(n-1)+1:j/i*nst1(m)=0;endelsefor m=j/i*(n-1)+1:j/i*nst1(m)=1;endendendst2=t; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%基带信号求反for n=1:j;if st1(n)>=1;st2(n)=0;elsest2(n)=1;endend;figure(1);subplot(411);plot(t,st1);title('基带信号st1');axis([0,5,-1,2]);subplot(412);plot(t,st2);title('基带信号反码st2');axis([0,5,-1,2]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%载波信号s1=cos(2*pi*f1*t)s2=cos(2*pi*f2*t)subplot(413),plot(s1);title('载波信号s1');subplot(414),plot(s2);title('载波信号s2'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%调制F1=st1.*s1;%加入载波1F2=st2.*s2;%加入载波2figure(2);subplot(411);plot(t,F1);title('F1=s1*st1');subplot(412);plot(t,F2);title('F2=s2*st2');e_fsk=F1+F2;subplot(413);plot(t,e_fsk);title('2FSK信号')%键控法产生的信号在相邻码元之间相位不一定连续nosie=rand(1,j);fsk=e_fsk+nosie;subplot(414);plot(t,fsk);title('加噪声后信号') %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解调st1=fsk.*s1;%与载波1相乘[f,sf1] = T2F(t,st1);%通过低通滤波器[t,st1] = lpf(f,sf1,2*fm);figure(3);subplot(311);plot(t,st1);title('与s1相乘后波形');st2=fsk.*s2;%与载波2相乘[f,sf2] = T2F(t,st2);%通过低通滤波器[t,st2] = lpf(f,sf2,2*fm);subplot(312);plot(t,st2);title('与s2相乘后波形'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽样判决for m=0:i-1;if st1(1,m*500+250)<st2(1,m*500+250);for j=m*500+1:(m+1)*500;at(1,j)=0;endelsefor j=m*500+1:(m+1)*500;at(1,j)=1;endendend;subplot(313);plot(t,at);axis([0,5,-1,2]);title('抽样判决后波形')程序三2PSK调制解调程序及注释clear allclose alli=10;j=5000;fc=4;%载波频率fm=i/5;%码元速率B=2*fm;t=linspace(0,5,j); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%产生基带信号a=round(rand(1,i));%随机序列,基带信号figure(3);stem(a);st1=t;for n=1:10if a(n)<1;for m=j/i*(n-1)+1:j/i*nst1(m)=0;endelsefor m=j/i*(n-1)+1:j/i*nst1(m)=1;endendendfigure(1);subplot(411);plot(t,st1);title('基带信号st1');axis([0,5,-1,2]);%%%%%%%%%%%&&&&&&%%%%%%%%%%%%%%%%%基带信号求反%由于PSK中的是双极性信号，因此对上面所求单极性信号取反来与之一起构成双极性码st2=t;for k=1:j;if st1(k)>=1;st2(k)=0;elsest2(k)=1;endend;subplot(412);plot(t,st2);title('基带信号反码st2');axis([0,5,-1,2]);st3=st1-st2;subplot(413);plot(t,st3);title('双极性基带信号st3');axis([0,5,-2,2]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%载波信号s1=sin(2*pi*fc*t);subplot(414);plot(s1);title('载波信号s1'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%调制e_psk=st3.*s1;figure(2);subplot(511);plot(t,e_psk);title('e_2psk');noise=rand(1,j);psk=e_psk+noise;%加入噪声subplot(512);plot(t,psk);title('加噪后波形'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解调psk=psk.*s1;%与载波相乘subplot(513);plot(t,psk);title('与载波s1相乘后波形');[f,af] = T2F(t,psk);%%%%%%%%%%%通过低通滤波器[t,psk] = lpf(f,af,B);subplot(514);plot(t,psk);title('低通滤波后波形'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽样判决for m=0:i-1;if psk(1,m*500+250)<0;for j=m*500+1:(m+1)*500;psk(1,j)=0;endelsefor j=m*500+1:(m+1)*500;psk(1,j)=1;endendendsubplot(515);plot(t,psk);axis([0,5,-1,2]);title('抽样判决后波形')程序四 2DPSK调制解调程序及注释clear allclose alli=10;j=5000;fc=4;%载波频率fm=i/5;%码元速率B=2*fm;t=linspace(0,5,j); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%产生基带信号a=round(rand(1,i));figure(4);stem(a);st1=t;for n=1:10if a(n)<1;for m=j/i*(n-1)+1:j/i*nst1(m)=0;endelsefor m=j/i*(n-1)+1:j/i*nst1(m)=1;endendendfigure(1);subplot(321);plot(t,st1);title('绝对码');axis([0,5,-1,2]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%差分变换b=zeros(1,i);%%%%%%%%全零矩阵b(1)=a(1);for n=2:10if a(n)>=1;if b(n-1)>=1b(n)=0;elseb(n)=1;endelseb(n)=b(n-1);endendst1=t;for n=1:10if b(n)<1;for m=j/i*(n-1)+1:j/i*nst1(m)=0;endelsefor m=j/i*(n-1)+1:j/i*nst1(m)=1;endendendsubplot(323);plot(t,st1);title('相对码st1');axis([0,5,-1,2]);st2=t;for k=1:j;if st1(k)>=1;st2(k)=0;elsest2(k)=1;endend;subplot(324);plot(t,st2);title('相对码反码st2');axis([0,5,-1,2]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%载波信号s1=sin(2*pi*fc*t);subplot(325);plot(s1);title('载波信号s1');s2=sin(2*pi*fc*t+pi);subplot(326);plot(s2);title('载波信号s2'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%信号调制d1=st1.*s1;d2=st2.*s2;figure(2);subplot(411);plot(t,d1);title('st1*s1');subplot(412);plot(t,d2);title('st2*s2');e_dpsk=d1+d2;subplot(413);plot(t,e_dpsk);title('调制后波形');noise=rand(1,j);dpsk=e_dpsk+noise;%%%%%%%%%%加入噪声subplot(414);plot(t,dpsk);title('加噪声后信号');%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解调dpsk=dpsk.*s1;%%%%与载波s1相乘figure(3);subplot(411);plot(t,dpsk);title('与载波相乘后波形');[f,af]=T2F(t,dpsk);%%%%通过低通滤波器[t,dpsk]=lpf(f,af,B);subplot(412);plot(t,dpsk);title('低通滤波后波形'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽样判决st=zeros(1,i);%%全零矩阵for m=0:i-1;if dpsk(1,m*500+250)<0;st(m+1)=0;for j=m*500+1:(m+1)*500;dpsk(1,j)=0;endelsefor j=m*500+1:(m+1)*500;st(m+1)=1;dpsk(1,j)=1;endendendsubplot(413);plot(t,dpsk);axis([0,5,-1,2]);title('抽样判决后波形') %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%码反变换dt=zeros(1,i);%%全零矩阵dt(1)=st(1);for n=2:10;if (st(n)-st(n-1))<=0&&(st(n)-st(n-1))>-1;dt(n)=0;elsedt(n)=1;endendst=t;for n=1:10if dt(n)<1;for m=j/i*(n-1)+1:j/i*nst(m)=0;endelsefor m=j/i*(n-1)+1:j/i*nst(m)=1;endendendsubplot(414);plot(t,st);axis([0,5,-1,2]);title('码反变换后波形');。