基于MATLAB的2ASK2FSK和2PSK的调制仿真Word文档下载推荐.docx

1、2FSK信号的调制方法主要有两种。第一种是用二进制基带矩形脉冲信号去调制一个调频器，使其能够输出两个不同频率的码元。第二种方法是用一个受基带脉冲控制的开关电路去选择两个独立频率源的振荡作为输出。 2FSK信号的接收也分为相关和非相关接收两类。相关接收根据已调信号由两个载波f1、f2调制而成，则先用两个分别对f1、f2带通的滤波器对已调信号进行滤波，然后再分别将滤波后的信号与相应的载波f1、f2相乘进行相干解调，再分别低通滤波、用抽样信号进行抽样判决器即可。原理图如下：非相关接收经过调制后的2FSK数字信号通过两个频率不同的带通滤波器f1、f2滤出不需要的信号，然后再将这两种经过滤波的信号分别通

2、过包络检波器检波，最后将两种信号同时输入到抽样判决器同时外加抽样脉冲，最后解调出来的信号就是调制前的输入信号。其原理图如下图所示：图5（b）非相干方式33、二进制相移键控（2PSK）相移键控是利用载波的相位变化来传递数字信息，而振幅和频率保持不变。在2PSK中通常用初始相位0和 分别表示二进制“1”和“0”。因此，2PSK信号的时域表达式 其中， 表示第n个符号的绝对相位： 0 发送“0”时 发送“1”时因此，式子可以改写为概率为P概率为1-P典型波形如图所示。由于表示信号的两种码元的波形相同，极性相反，故2PSK信号一般可以表述为一个双极性全占空矩形脉冲序列与一个正弦载波的相乘，即其中这里，

3、g（t）是脉宽为Ts的单个矩形脉冲，而 的统计特性为1 概率为P-1 概率为1-P2PSK信号的调制原理框图如图所示。2ASK信号的产生方法比较知识对s（t）的要求不同，在2ASK中s（t）是单极性的，而在2PSK中s（t）是双极性的基带信号。2PSK信号的解调通常采用相干解调法，解调器原理框图如图所示。四仿真源程序和代码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）%ra

4、nd函数产生在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）1; for m=j/i*（n-1）+1:j/i*n st（m）=0; end else st（m）=1;endfigure（1）;subplot（421）;plot（t,st）;axis（0,5,-1,2）;title（基带信号st）;%载波s1=cos（2*pi*fc*t）;subplot（422）;plot（s1）;载波信号s1%调制e_2ask=st.*s1;s

5、ubplot（423）;plot（t,e_2ask）;已调信号noise =rand（1,j）;e_2ask=e_2ask+noise;%加入噪声subplot（424）;加入噪声的信号%相干解调at=e_2ask.*cos（2*pi*fc*t）;at=at-mean（at）;%因为是单极性波形，还有直流分量，应去掉subplot（425）;plot（t,at）;与载波相乘后信号f,af = T2F（t,at）;%通过低通滤波器t,at = lpf（f,af,2*fm）;subplot（426）;相干解调后波形%抽样判决for m=0:i-1; if at（1,m*500+250）+0.5=1

6、; st2（n）=0; st2（n）=1;end;subplot（411）;plot（t,st1）;基带信号st1subplot（412）;plot（t,st2）;基带信号反码st2%载波信号s1=cos（2*pi*f1*t）s2=cos（2*pi*f2*t）subplot（413）,plot（s1）;subplot（414）,plot（s2）;载波信号s2F1=st1.*s1;%加入载波1F2=st2.*s2;%加入载波2figure（2）;plot（t,F1）;F1=s1*st1plot（t,F2）;F2=s2*st2e_fsk=F1+F2;subplot（413）;plot（t,e_fs

7、k）;2FSK信号）%键控法产生的信号在相邻码元之间相位不一定连续nosie=rand（1,j）;fsk=e_fsk+nosie;subplot（414）;plot（t,fsk）;加噪声后信号st1=fsk.*s1;%与载波1相乘f,sf1 = T2F（t,st1）;t,st1 = lpf（f,sf1,2*fm）;figure（3）;subplot（311）;与s1相乘后波形st2=fsk.*s2;%与载波2相乘f,sf2 = T2F（t,st2）;t,st2 = lpf（f,sf2,2*fm）;subplot（312）;与s2相乘后波形 if st1（1,m*500+250） st2（k）=

8、0; st2（k）=1;st3=st1-st2;plot（t,st3）;双极性基带信号st3axis（0,5,-2,2）;s1=sin（2*pi*fc*t）;e_psk=st3.*s1;subplot（511）;plot（t,e_psk）;e_2psknoise=rand（1,j）;psk=e_psk+noise;subplot（512）;plot（t,psk）;加噪后波形psk=psk.*s1;%与载波相乘subplot（513）;与载波s1相乘后波形f,af = T2F（t,psk）;%通过低通滤波器t,psk = lpf（f,af,B）;subplot（514）;低通滤波后波形 if p

9、sk（1,m*500+250）0; psk（1,j）=0; psk（1,j）=1;subplot（515）;用到的傅立叶T2F函数%利用FFT计算信号的频谱并与信号的真实频谱的抽样比较。%脚本文件T2F.m定义了函数T2F，计算信号的傅立叶变换。function f,sf= T2F（t,st）dt = t（2）-t（1）;T=t（end）;df = 1/T;N = length（st）;f=-N/2*df:df:N/2*df-df; sf = fft（st）;sf = T/N*fftshift（s）；用到的低通滤波器函数function t,st=lpf（f,sf,B）df = f（2）-f（

10、1）;T = 1/df;hf = zeros（1,length（f）;%全零矩阵bf = -floor（ B/df ）: floor（ B/df ） + floor（ length（f）/2 ）;hf（bf）=1;yf=hf.*sf;t,st=F2T（f,yf）;st = real（st）;用到的反傅立叶函数%脚本文件F2T.m定义了函数F2T，计算信号的反傅立叶变换。function t,st=F2T（f,sf）%This function calculate the time signal using ifft function for the input%signals spectrum

11、 df = f（2）-f（1）;Fmx = （ f（end）-f（1） +df）;dt = 1/Fmx;N = length（sf）;T = dt*N;%t=-T/2:dt:T/2-dt;t = 0:T-dt;sff = fftshift（sf）;st = Fmx*ifft（sff）;五、仿真结果及分析5.1 2ASK分析通过编写M文件程序，产生随机信号，编程后。程序中注有需注意语句及解释。运行程序，实现2ASK的调制与解调过程。本次设计采用模拟调制法（相乘器法）和相干解调法。仿真后调制过程及解调过程的图形分别如图1和图2所示。图 1. 2ASK调制过程仿真图图2. 2ASK解调过程仿真图由图

12、可以看出，产生的数字随机信号为“0100100001”，经载波调制后信号为“1”的时间内有正弦波形，信号为“0”的时间内无波形。经过加随机噪声，相干解调后，恢复出原始信号，与基带信号一致，因此达到本次设计目的。5.2 2FSK分析本次设计中采用键控法调制法和相干解调法。仿真后调制过程及解调过程的图形分别如图3、图4和图5所示。图 3 2FSK键控法调制过程仿真图图 4 2FSK键控法调制过程仿真图图 5 2FSK相干解调过程仿真图由图可以看出，产生的随机信号为“1011001001”，经过反相产生反码，并分别与两个载波相乘，经过加入随机噪声后波形如图2-15所示。在解调时，分别与对应的载波相乘。解调出基带信号，可以看出实现了本次设计目的。5.3.2PSK仿真结果及分析运行程序，实现2PSK的调制与解调过程。本次设计采用模拟调制法和相干解调法。仿真后调制过程及解调过程的图形分别如图6和图7所示。图 6 2PSK模拟调制方法过程仿真图图 7 2PSK相干解调过程仿真图由图可以看出，产生的随机信号为“0011001100”，经过反相产生反码，并将原码跟反码一起合成双极性码，与载波相乘后加入随机噪声。在解调时，与对应的载波相乘经过低通滤波、抽样判决后，解调出基带信号与原基带信号一致，可以看出实现了本次设计目的。