2PSK系统设计课程设计报告材料

1、华南理工大学通信原理课程设计报告题目：2PSK系统仿真专业：班级：某某:学号:日期：20XX年XX月、实验需要材料MATLAB软件二、实验要求完成规定系统的MATLAB程以与simulink的仿真，根本内容包括：输入信号，系统中各个关键模块的输出情况。并调整仿真的参数得到不同的仿真结果。三、设计原理2PSK汉语全称：二进制相移键控。2PSK是相移键控的最简单的一种形式，它用两个初相相隔为180的载波来传递二进制信息。所以也被称为BPSKSimulink简介：Simulink是Mathworks公司推出的基于Matlab平台的著名仿真环境Simulin作为一种专业和功能强大且操作简单的仿真工具，

2、目前已被越来越多的工程技术人员所青睐，它搭建积木式的建模仿真方式既简单又直观，而且已经在各个领域得到了广泛的应用。数字信号的传输方式分为基带传输和带通传输，在实际应用中，大多数信道具有带通特性而不能直接传输基带信号。为了使数字信号在带通信道中传输，必须使用数字基带信号对载波进展调制，以使信号与信道的特性相匹配。这种用数字基带信号控制载波，把数字基带信号变换为数字带通信号的过程称为数字调制。数字调制技术的两种方法： 利用模拟调制的方法去实现数字式调制，即把数字调制看成是模拟调制的一个特例，把数字基带信号当做模拟信号的特殊情况处理。利用数字信号的离散取值特点通过开关键控载波，从而实现数字调制。这种

3、方法通常称为键控法，比如对载波的相位进展键控，便可获得相移键控PSK1根本的调制方式。图1相应的信号波形的示例101调制原理:在二进制数字调制中，当正弦载波的相位随二进制数字基带信号离散变化时，如此产生二进制移相键控(2PSK)信号。2PSK信号调制有两种方法，即模拟调制法和键控法。通常用已调信号载波的0。和180。分别表示二进制数字基带信号的1和0,模拟调制法用两个反相的载波信号进展调制。2PSK以载波的相位变化作为参考基准的，当基带信号为0时相位相对于初始相位为0°,当基带信号为1时相对于初始相位为180°。键控法，是用载波的相位来携带二进制信息的调制方式。通常用0。和

4、180。来分别代表0和1O其时域表达式为：e2PSKang(tnTDcosct其中，2PSK的调制中an必须为双极性码。两种方法原理图分别如图1-1和图1-2所示。<6®不归菱COSfU!图1-1模拟调制法原理图解调原理2PSK信号的解调方法是相干解调法，如图2-1所示。由于PSK信号本身就是利用相位传递信息的，所以在接收端必须利用信号的相位信息来解调信号。由于2PSK的幅度是恒定的，必须进展相干解调。经过带通滤波的信号在相乘器中与本地载波相乘，然后用低通滤波器滤除高频分量，在进展抽样判决。判决器是按极性来判决的。即正抽样值判为1,负抽样值判为0。0PSK(t)图2-12PSK

5、相干解调法原理图0相位，如图2-2所假设相干载波的基准相位与2PSK信号的基准一致通常默认为示。但是由于2PSK信号的载波回复过程中存在着180°的相位模糊，即恢复的本地载波与所需相干载波可能一样，也可能相反，这种相位关系的不确定性将会造成解调出的数字基带信号与发送的基带信号正好相反，即“1"变成"0",“0"变成"1",判决器输出数字信号全部出错。这种现象称为2PSK方式的“倒兀"现象或“反相工作"。但在本次仿真中是直接给其同频同相的载波信号，所以不存在此问题。abcde图2-2相干解调中各点波形图相关

6、公式：2PSK信号在一个码元的持续时间Ts内可以表示为u1T(t)发送“1"时Sr(t)=uoT(t)=-u1T(t)发送“0"时其中Acos3ct0<t<Tsu1T(t)=0其他设发送端发出的信号如上式所示，如此接收端带通滤波器输出波形y(t)为a+nc(t)cos3ct-ns(t)sinwct,发送"T'时y(t)=-a+nc(t)cos3ct-ns(t)sinwct，发送"0”时y(t)经过想干解调相乘一低通后，送入抽样判决器的输入波形为a+nc(t),发送“1"时x(t)=-a+nc(t)，发送“0"时由最

7、优判决门限分析可知，在发送“1"和“0"概率相等时，即P(1)=P(0)时，最优门限b*=0.此时，发“1"而错判为“0"的概率为P(0/1)=P(x三0)=/0-fi(x)dx=1/2erfc(yfT)式中：r=a2/2b2n同理，发“0而错判为“1"的概率为P(1/0)=P(x>0)=/°。f0(x)dx=1/2erfc(y/Y)故2PSK信号相干解调系统的总误码率为Pe=P(1)P(0/1)+P(0)P(0/1)=1/2erfc()在大信噪比(r>>1)的条件下，上式可近似为四、源程序与相应实验结果1、相干调解

8、法i=10;j=5000;fc=4;%载波频率fm=i/5;%码元速率B=2*fm;t=linspace(0,5,j);a=round(rand(1,i);%随机序列，基带信号figure(1);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(2);subplot(411);plot(t,st1);title('基带信号st1');axis(0,5,-1,2);%由于PS仲的是双极性信号，

9、因此对上面所求单极性信号取反来与之一起构成双极性码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

10、=st3.*s1;figure(3);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相乘后波形');lpf=fir1(101,2/1000,52/1000);%低通滤波器设置st=filter(lpf,1,psk);%经过低通滤波器后的信号subp

11、lot(514);plot(t,st);title('低通滤波后波形');form=0:i-1;ifst(1,m*500+250)<0;forj=m*500+1:(m+1)*500;st(1,j)=0;endelseforj=m*500+1:(m+1)*500;st(1,j)=1;endendendsubplot(515);plot(t,st);axis(0,5,-1,2);title('抽样判决后波形')2、模拟调制法clearallclosealli=10;j=5000;fc=4;%载波频率fm=i/5;%码元速率B=2*fm;t=linspace(0

12、,5,j);a=round(rand(1,i);figure(4);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(5);subplot(321);plot(t,st1);title('绝对码');axis(0,5,-1,2);b=zeros(1,i);%全零矩阵b(1)=a(1);forn=2:10ifa(n)>=1;ifb(n-1)>=1b(n)=0;elseb(n)=1;

13、endelseb(n)=b(n-1);endendst1=t;forn=1:10ifb(n)<1;form=j/i*(n-1)+1:j/i*nst1(m)=0;endelseform=j/i*(n-1)+1:j/i*nst1(m)=1;endendendsubplot(323);plot(t,st1);title('相对码st1');axis(0,5,-1,2);st2=t;fork=1:j;ifst1(k)>=1;st2(k)=0;elsest2(k)=1;endend;subplot(324);plot(t,st2);title('相对码反码st2

14、9;);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(4

15、13);plot(t,e_dpsk);title('调制后波形');noise=rand(1,j);参加噪声dpsk=e_dpsk+noise;subplot(414);plot(t,dpsk);title('加噪声后信号');与载波s1相乘dpsk=dpsk.*s1;figure(3);subplot(411);plot(t,dpsk);title('与载波相乘后波形');低通滤波器设置lpf=fir1(101,2/1000,52/1000);%k=filter(lpf,1,dpsk);%经过低通滤波器后的信号subplot(412);plot

16、(t,k);title('低通滤波后波形');st=zeros(1,i);%全零矩阵form=0:i-1;ifk(1,m*500+250)<0;st(m+1)=0;forj=m*500+1:(m+1)*500;k(1,j)=0;endelseforj=m*500+1:(m+1)*500;st(m+1)=1;k(1,j)=1;endendendsubplot(413);plot(t,k);axis(0,5,-1,2);title('抽样判决后波形')dt=zeros(1,i);%全零矩阵dt(1)=st(1);forn=2:10;if(st(n)-st(n-

17、1)<=0&&(st(n)-st(n-1)>-1;dt(n)=0;elsedt(n)=1;endendst=t;forn=1:10ifdt(n)<1;form=j/i*(n-1)+1:j/i*nst(m)=0;endelseform=j/i*(n-1)+1:j/i*nst(m)=1;endendendsubplot(414);plot(t,st);axis(0,5,-1,2);title('码反变换后波形');st=real(st);运行示例：方法一效果图Figure?XFileEditViewnsertTookDesktop理indo*Hel

18、pD总敛SC®基带信号mD50010001SOO2000250030003500400045005000方法二效果图：ffj|Figure4XFileEditView£hs«rtTookDesktopWindowHelp窖巨LI昌|O®|曰罔1BO绝对码2345相对码all相对玛反玛st22I112345%123戴波信号$1载波信号K'firllliimm-IWml'WWillknllJFigure5XFileEditViewMsertTookDesktopVWndowHelpD喝LI写饪鼠忌©®_LSB"。200040006000"200040006000参考文献1李晓峰.通信原理：清华大学五、个人体会与感想这次课程设计对我的自身能力有了进一步了解。第一点，这进一步端正了我的学习态度，学会了实事求是、严谨的作风，提高了动手能力。也要对自己严格要求，不能够一知半解，要力求明明白白。浮躁的性格对丁搞设计来说或者学习是致命的，一定要静下心来，踏实的做事。第二点，我觉得动手之前，头脑里必须活楚应该怎么做，这一点是很重要的，所谓三思而后行。第三点，在这两周的时间里，深深感觉到我对所学的知识的了解还是很局限、很肤浅的，在以后的时问里还要继续加深对这门课程的理解，加强知