基于MATLAB的2ASK、2FSK和2PSK的调制仿真

1实验报告（一）一、实验名称：基于MATLAB的2ASK、2FSK和2PSK的调制仿真二、实验目的：（1）熟悉2ASK、2FSK和2PSK的调制原理。（2）学会运用Matlab编写2ASK、2FSK和2PSK调制程序。（3）会画出原信号和调制信号的波形图。（4）掌握数字通信的2ASK、2FSK和2PSK的调制方式。三、实验原理分析3.1二进制振幅键控（2ASK）振幅键控是利用载波的幅度变化来传递数字信息，而其频率和初始相位保持不变。在2ASK中，载波的幅度只有两种变化状态，分别对应二进制信息“0”或“1”。二进制振幅键控的表达式为：s(t)=A(t)cos(w0+θ)0＜t≤T式中，w0=2πf0为载波的角频率；A(t)是随基带调制信号变化的时变振幅，即A(t)=0A典型波形如图所示：2ASK信号的产生方法通常有两种：相乘法和开关法，相应的调制器如图2。图2（a）就是一般的模拟幅度调制的方法，用乘法器实现；图2（b）是一种数字键控法，其中的开关电路受s(t)控制。在接收端，2ASK有两种基本的解调方法：非相干解调（包络检波法）和相干解调（同步检测法），相应的接收系统方框图如图：当发送“1”时当发送“0”时23.2、二进制频移键控（2FSK）二进制频移键控信号码元的“1”和“0”分别用两个不同频率的正弦波形来传送，而其振幅和初始相位不变。故其表达式为：)(st时0发送“),cos(”时1发送“),cos21(nntAtA图42FSK信号时间波形由图可见，2FSK信号的波形（a）可以分解为波形（b）和波形（c）,也就是说，一个2FSK信号可以看成是两个不同载频的2ASK信号的叠加。2FSK信号的调制方法主要有两种。第一种是用二进制基带矩形脉冲信号去调制一个调频器，使其能够输出两个不同频率的码元。第二种方法是用一个受基带脉冲控制的开关电路3去选择两个独立频率源的振荡作为输出。2FSK信号的接收也分为相关和非相关接收两类。相关接收根据已调信号由两个载波f1、f2调制而成，则先用两个分别对f1、f2带通的滤波器对已调信号进行滤波，然后再分别将滤波后的信号与相应的载波f1、f2相乘进行相干解调，再分别低通滤波、用抽样信号进行抽样判决器即可。原理图如下：非相关接收经过调制后的2FSK数字信号通过两个频率不同的带通滤波器f1、f2滤出不需要的信号，然后再将这两种经过滤波的信号分别通过包络检波器检波，最后将两种信号同时输入到抽样判决器同时外加抽样脉冲，最后解调出来的信号就是调制前的输入信号。其原理图如下图所示：3．3、二进制相移键控（2PSK）相移键控是利用载波的相位变化来传递数字信息，而振幅和频率保持不变。在2PSK中通常用初始相位0和分别表示二进制“1”和“0”。因此，2PSK信号的时域表达式其中，表示第n个符号的绝对相位：图5（b）非相干方式)(2teFSK)cos(nctAn40发送“0”时发送“1”时因此，式子可以改写为概率为P概率为1-P典型波形如图所示。由于表示信号的两种码元的波形相同，极性相反，故2PSK信号一般可以表述为一个双极性全占空矩形脉冲序列与一个正弦载波的相乘，即其中这里，g(t)是脉宽为Ts的单个矩形脉冲，而的统计特性为1概率为P-1概率为1-P2PSK信号的调制原理框图如图所示。2ASK信号的产生方法比较知识对s(t)的要求不同，在2ASK中s(t)是单极性的，而在2PSK中s(t)是双极性的基带信号。2PSK信号的解调通常采用相干解调法，解调器原理框图如图所示。n)(2teFSKtAccostAccos)(2teFSKttsccos)(nnnTtga)(nanas(t)码型变换双极性不归零乘法器e2PSK(t)cosct(a)cosct0°开关电路e2PSK(t)180°移相s(t)(b)带通滤波器e2PSK(t)a相乘器c低通滤波器dbe抽样判决器输出cosct定时脉冲)(ts5四．仿真源程序和代码2ASK调制解调程序及注释clearallclosealli=10;%10个码元j=5000;t=linspace(0,5,j);%0-5之间产生5000个点行矢量，即分成5000fc=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;forn=1:10ifa(n)1;form=j/i*(n-1)+1:j/i*nst(m)=0;endelseform=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('相干解调后波形');%抽样判决form=0:i-1;ifat(1,m*500+250)+0.50.5;forj=m*500+1:(m+1)*500;at(1,j)=0;endelseforj=m*500+1:(m+1)*500;at(1,j)=1;endendendsubplot(427);plot(t,at);axis([0,5,-1,2]);title('抽样判决后波形')2FSK调制解调程序及注释clearallclosealli=10;%基带信号码元数j=5000;a=round(rand(1,i));%产生随机序列t=linspace(0,5,j);f1=10;%载波1频率f2=5;%载波2频率fm=i/5;%基带信号频率6%产生基带信号st1=t;forn=1:10ifa(n)1;form=j/i*(n-1)+1:j/i*nst1(m)=0;endelseform=j/i*(n-1)+1:j/i*nst1(m)=1;endendendst2=t;%基带信号求反forn=1:j;ifst1(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相乘后波形');%抽样判决form=0:i-1;ifst1(1,m*500+250)st2(1,m*500+250);forj=m*500+1:(m+1)*500;at(1,j)=0;endelseforj=m*500+1:(m+1)*500;at(1,j)=1;endendend;subplot(313);plot(t,at);axis([0,5,-1,2]);title('抽样判决后波形')72PSK调制解调程序及注释clearallclosealli=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;forn=1:10ifa(n)1;form=j/i*(n-1)+1:j/i*nst1(m)=0;endelseform=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;fork=1:j;ifst1(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('低通滤波后波形');%抽样判决form=0:i-1;ifpsk(1,m*500+250)0;forj=m*500+1:(m+1)*500;psk(1,j)=0;endelseforj=m*500+1:(m+1)*500;psk(1,j)=1;endendend8subplot(515);plot(t,psk);axi