完整版)通信原理课程设计——DSB调制解调系统设计与仿真通信原理

1、HEJ业1 liMltD LMVtJ祁HT通信原理课程设计设计题目：DSB调制解调系统设计与仿真通信原理班 级:学生姓名：学生学号:指导老师:引言31、课程设计目的32、课程设计要求3一、DSB调制解调模型的建立 41、DSB信号的模型42、DSB信号调制过程分析 43、高斯白噪声信道特性分析 64、DSB解调过程分析95、DSB调制解调系统抗噪声性能分析 10二、仿真过程13三、心得体会15四、参考文献15引言本课程设计用于实现 DSB信号的调制解调过程。信号的调制与解调在通信系统中具有重要的作 用。调制过程是一个频谱搬移的过程，它是将低频信号的频谱搬移到载频位置。解调是调制的逆过 程，即是

2、将已调制的信号还原成原始基带信号的过程。信号的接收端就是通过解调来还原已调制信 号从而读取发送端发送的信息。因此信号的解调对系统的传输有效性和传输可靠性有着很大的影响。 调制与解调方式往往决定了一个通信系统的性能。双边带DSB信号的解调采用相干解调法，这种方式被广泛应用在载波通信和短波无线电话通信中。1、课程设计目的本课程设计是实现 DSB的调制解调。在此次课程设计中，我们将通过多方搜集资料与分析，来理解DSB调制解调的具体过程和它在MATLAB的实现方法。预期通过这个阶段的研习，更清晰地认识DSB的调制解调原理，同时加深对MATLAB这款通信仿真软件操作的熟练度，并在使用中去感受MATLAB

3、勺应用方式与特色。利用自主的设计过程来锻炼自己独立思考，分析和解决问题的能力，为 我们今后的自主学习研究提供具有实用性的经验。2、课程设计要求(1) 熟悉MATLAB M文件的使用方法，掌握 DSB信号的调制解调原理，以此为基础用M文件 编程实现DSB信号的调制解调。(2) 绘制出SSB信号调制解调前后在时域和频域中的波形，观察两者在解调前后的变化，通过 对分析结果来加强对 DSB信号调制解调原理的理解。(3) 对信号分别叠加大小不同的噪声后再进行解调，绘制出解调前后信号的时域和频域波形， 比较未叠加噪声时和分别叠加大小噪声时解调信号的波形有何区别，由所得结果来分析噪声对信号 解调造成的影响。

4、(4) 在老师的指导下，独立完成课程设计的全部内容，并按要求编写课程设计论文，文中能正 确阐述和分析设计和实验结果。UMItU UHUSflTDSB调制解调模型的建立1、DSB信号的模型在AM信号中，载波分量并不携带信息，信息完全由边带传送。如果将载波抑制，只需在将直流Ao去掉，即可输出抑制载波双边带信号，简称双边带信号(DSB调制器模型如图1所示。DSB。DSB调制器模型其中，设正弦载波为C(t)A cos( cto)式中，A为载波幅度；c为载波角频率;0为初始相位(假定0为0)。P低通滤波器调制过程是一个频谱搬移的过程，它是将低频信号的频谱搬移到载频位置。而解调是将位于载 频的信号频谱再搬

5、回来，并且不失真地恢复出原始基带信号。双边带解调通常采用相干解调的方式，它使用一个同步解调器，即由相乘器和低通滤波器组成。2所示:在解调过程中，输入信号和噪声可以分别单独解调。相干解调的原理框图如图COS(0訂+仍图2相干解调器的数学模型信号传输信道为高斯白噪声信道，其功率为2、DSB信号调制过程分析假定调制信号 m(t)的平均值为0,与载波相乘，即可形成DSB信号，其时域表达式为sdsb m(t)cos ct式中，m(t)的平均值为0。DSB的频谱为1SDsb( )M (c) M (c)2DSB信号的包络不再与调制信号的变化规律一致，因而不能采用简单的包络检波来恢复调制信号，需180的突变。

6、除了不再含有载频分量离散谱外，DSB信号的频谱与AM信号的频谱完全相同，仍由上下对称的两个边带组成。所以 DSB信号的带宽与 AM信号的带宽相同，也为基带信号带宽的两倍，即BdsbBam2 fH式中，fH为调制信号的最高频率。调制信号产生的代码及波形为elf;%清除窗口中的图形ts=0.01;%定义变量区间步长t0=2;%定义变量区间终止值t=-t0+0.0001:ts:t0;%定义变量区间fe=10;%给出相干载波的频率A=1;%疋义输入信号幅度fa=1;%定义调制信号频率mt=A*cos(2*pi*fa.*t);%输入调制信号表达式ct=cos(2*pi*fc.*t);%输入调制信号表达式

7、psn t=mt.*cos(2*pi*fc.*t);%输出调制信号表达式subplot(3,1,1);%划分画图区间plot(t,mt,g);%画出输入信号波形title（输入信号波形）;xlabel(Variable t);ylabel(Variable mt);subplot(3,1,2);plot(t,ct,b);%画出输入信号波形title（输入载波波形）;xlabel(Variable t);ylabel(Variable ct);subplot(3,1,3);plot(1:le ngth(ps nt),ps nt,丫)；%length用于长度匹配title（已调信号波形）;%画出已

8、调信号波形xlabel(Variable t);ylabel(Variable psn t);运行结果：13tc ATD F a V10.50-0.5-1-1.5-1输入信号波形-0.500.5Variable t输入载波波形11.52t m FbFav-210.50-0.5-1-2-1.5210-1-2050-1-0.500.5Variable t已调信号波形100150200Variable t25011.52300350400tns ptorborav图3调制信号、载波、已调信号波形3、高斯白噪声信道特性分析在实际信号传输过程中，通信系统不可避免的会遇到噪声，例如自然界中的各种电磁波噪声

9、和 设备本身产生的热噪声、散粒噪声等，它们很难被预测。而且大部分噪声为随机的高斯白噪声，所 以在设计时引入噪声，才能够真正模拟实际中信号传输所遇到的问题，进而思考怎样才能在接受端 更好地恢复基带信号。信道加性噪声主要取决于起伏噪声，而起伏噪声又可视为高斯白噪声，因此 我在此环节将对双边带信号添加高斯白噪声来观察噪声对解调的影响情况。为了具体而全面地了解噪声的影响问题，我们将分别引入大噪声（信噪比为20dB）与小噪声（信噪比为2dB）作用于双边带信号，再分别对它们进行解调，观察解调后的信号受到了怎样的影响。在此过程中，我用函数 randn来添加噪声，此函数功能为向信号中添加噪声功率为其方差的高

10、斯白噪声。正弦波通过加性高斯白噪声信道后的信号为r(t) Acos( ct)n(t)故其有用信号功率为A2噪声功率为22Httl UMilbD UMVtfiSlT信噪比S N满足公式B 1Ologi0(SN)则可得到公式AB2?1O10我们可以通过这个公式方便的设置高斯白噪声的方差。为了便于比较，我们显示了双边带信号加入两种噪声后的时频波形图。实现代码和波形如图4elf;ts=0.01;t0=2;t=-tO+O.OOO1:ts:tO;fe=10;A=1;fa=1;mt=A*cos(2*pi*fa.*t);xzb=2;sn r=10.A(xzb/10);h,l=size(mt);fan geha

11、=A*A./(2*s nr);n it=sqrt(fa ngcha).*ra ndn( h,l);psmt=mt.*eos(2*pi*fe.*t);psn t=psmt+ nit;xzb=20;sn r1=10.A(xzb/10);h,l=size(mt);fan geha1=A*A./(2*s nr1);n it1=sqrt(fa ngcha1).*ra ndn( h,l);psn t1=psmt+ nit1;subplot(2,2,1);plot(t, nit,g);title(小信噪比高斯白躁声);xlabel(Variable t);ylabel(Variable nit);%清除窗口

12、中的图形%定义变量区间步长%定义变量区间终止值%定义变量区间%给出相干载波的频率%定义输入信号幅度%定义调制信号频率%输入调制信号表达式%输入小信躁比(dB)%求调制信号的维数%由信躁比求方差%产生小信噪比高斯白躁声%输出调制信号表达式%输出叠加小信噪比已调信号波形%输入大信躁比(dB)%求调制信号的维数%由信躁比求方差%产生大信噪比高斯白躁声%输出已调信号波形%划分画图区间%画出输入信号波形2大小不同信噪比的解调波形，如图6:plot(t,ps nt,b);title(叠加小信噪比已调信号波形);xlabel(Variable t);%length用于长度匹配%画出输入信号与噪声叠加波形yl

13、abel(Variable psn t); subplot(2,2,3);plot(t ,n it1,r);title(大信噪比高斯白躁声);xlabel(Variable t);ylabel(Variable nit);subplot(2,2,4);%画出输出信号波形plot(t,ps nt1,k);title(叠加大信噪比已调信号波形);xlabel(Variable t);ylabel(Variable23tn torbaLravIILI1 rI IUUjlUUili I Jr? r,i. j Ja-2-1.5-1-0.500.5Variable t11.52小信噪比高斯白躁声-2-0

14、03 1111EC-2-1.5-1-0.500.511.52Variable t大信噪比咼斯白躁声0.03 ,.,2 1o oo oo-31.5-2-1.5-1-0.500.511.52Variable t叠加小信噪比已调信号波形2 10-2tnsp XUSFav-2-1.5-1-0.500.511.52Variable t叠加大信噪比已调信号波形msp toTDOrav-1.5psmt);图4不同信噪比的噪声及含噪声的已调波形可以清晰地看出，加大噪声后，解调信号的波形杂乱无章，起伏远大于加小噪声时的波形。造成此现象的原因是当信噪比较小时，噪声的功率在解调信号中所占比重较大，所以会造成杂 波较

15、多的情况；而信噪比很大时，噪声的功率在解调信号中所占比重就很小了，噪声部分造成的杂 乱波形相对就不是很明显，甚至可以忽略。综上所述，叠加噪声会造成解调信号的失真，信噪比越小，失真程度越大。所以当信噪比低于4、DSB解调过程分析实际应用中，应该尽量减少噪声的产生。所谓相干解调是为了从接收的已调信号中，不失真地恢复原调制信号，要求本地载波和接收信 号的载波保证同频同相。相干解调的一般数学模型如图所示。一低通滤波器cos (仏 f+g)图5 DSB相干解调模型设图四的输入为 DSB信号Sm(t) SDSB(t) m(t)COS( ct0)乘法器输出为(t)SDSB(t)m(t)cos( cto)co

16、s( ct)1m(t)cos(o) cos(2 cto )2通过低通滤波器后1mo(t) m(t)cos( o )2当o常数时，解调输出信号为1mo(t)m(t)FCitikJ Uhvltu WIVUHsrx大信噪比解调信号波形50-50200250300350400450500550Variable t小信噪比解调信号波形50600PKbavlav-50200250300350400450500550600Variable t图6不同信噪比解调波形5、DSB调制解调系统抗噪声性能分析由于加性噪声只对已调信号的接收产生影响，因而调制系统的抗噪声性能主要用解调器的抗噪声性能来衡量。为了对不同调制

17、方式下各种解调器性能进行度量，通常采用信噪比增益G (又称调制制度增益)来表示解调器的抗噪声性能。有加性噪声时解调器的数学模型如图7所示。艸)图7有加性噪声时解调器的数学模型图7中Sm(t)为已调信号,Sm(t)、噪声为高斯窄带n(t)为加性高斯白噪声。Sm(t)和n(t)首先经过带通滤波器，滤出有用信号，滤除带外的噪声。经过带通滤波器后到达解调器输入端的信号为噪声ni(t)，显然解调器输入端的噪声带宽与已调信号的带宽是相同的。最后经解调器解调输出的有用信号为mo(t)，噪声为no(t)。Sm(t) m(t)COS ct与相干载波COS ct相乘后，得m(t)cos211ctm(t) m(t)

18、cos 2 ct22经低通滤波器后，输出信号为mo(t)-m(t)2因此，解调器输出端的有用信号功率为So212m(t)-m (t)4解调DSB信号时，接收机中的带通滤波器的中心频率o与调制载频c相同，因此解调器输出端的窄带噪声ni(t)可表示为ni(t)n c(t)cos ctns(t )sin或它与相干载波相乘后，得n i(t)cos ct n c(t)cos 11nc(t)nc(t)cos 222经低通滤波器后，解调器最终的输出噪声为ctctno(t)1严)ns(t )sinctns(t)si n 2 或故输出噪声功率为No n o2(t)、c2(t)4这里，B 2fH，为DSB信号的带

19、通滤波器的带宽。解调器输入信号平均功率为Ss；(t) m(t)cos ct21 2 尹(t)NinoB同时可得解调器的输出信噪比1 2777So卩(t)-Ni4m2(t)因此制度增益为NonoBGdsbSNi由此可见，DSB调制系统的制度增益为 2。也就是说DSB信号的解调器使信噪比改善了一倍。%清除窗口中的图形%定义变量区间步长%定义变量区间终止值%定义变量区间%给出相干载波的频率%定义输入信号幅度%定义调制信号频率%输入调制信号表达式%输入信噪比(dB)%求调制信号的维数%由信躁比求方差%产生高斯白噪声%调制信号与噪声叠加%输出调制信号表达式%输出噪声表达式%输出已调信号波形%调制信号乘以

20、相干载波%低通滤波器的时域表达式%低通滤波器的频域表达式%解调信号的时域表达式%划分画图区间%画出输入信号波形elf;ts=0.01;t0=2;t=-tO+O.OOO1:ts:tO;fe=10;A=1;fa=1;mt=A*cos(2*pi*fa.*t);xzb=20;sn r=104(xzb/10);h,l=size(mt);fan geha=A*A./(2*s nr);n it=sqrt(fa ngcha).*ra ndn( h,l);sn it=mt+nit;psmt=mt.*eos(2*pi*fe.*t);pn it= ni t.*eos(2*pi*fe.*t);psn t=psmt+p

21、 nit;jie=ps nt.*eos(2*pi*fe.*t);ht=(2*pi*fe.*si n( 2*pi*fe.*t)./(2*pi*fe.*t)./pi; htw=abs(fft(ht);jt=eo nv(ht,jic);subplot(3,3,1);plot(t,mt,g);title(输入信号波形); xlabel(Variable t); ylabel(Variable mt); subplot(3,3,2);plot(t ,n it,b); title(输入噪声波形); xlabel(Variable t); ylabel(Variable nit);subplot(3,3,3

22、);【 IHX点拓fHttd UKitc mtirtfisnplot(1:le ngth(s nit),s nit,丫);%length用于长度匹配title(输入信号与噪声叠加波形);xlabel(Variable t); ylabel(Variable snit); subplot(3,3,4);plot(t,psmt,k);title(输出信号波形); xlabel(Variable t); ylabel(Variable psmt); subplot(3,3,5); plot(t,p ni t,k);title(输出噪声波形);%画出输入信号与噪声叠加波形%画出输出信号波形%画出输出噪

23、声波形xlabel(Variable t);ylabel(Variable pn it);subplot(3,3,6);plot(t,ps nt,k);title(输出信号与输出噪声叠加波形);xlabel(Variable t);ylabel(Variable psn t);subplot(3,3,7);plot(1:le ngth(htw),htw,k);title(低通滤波器频域波形);xlabel(Variable w);ylabel(Variable htw);axis(0 60 0 150);subplot(3,3,8);plot(1:le ngth(ht),ht,k);title

24、(低通滤波器时域波形);%画出输出信号与输出噪声叠加波形%画出低通滤波器频域波形%画出低通滤波器时域波形xlabel(Variable t);ylabel(Variable psn t);axis(150 250 -20 25);subplot(3,3,9);plot(1:le ngth(jt),jt,k);title(输出信号与输出噪声叠加波形);%给出坐标轴范围%画出输出信号与输出噪声叠加波形xlabel(Variable t);ylabel(Variable jt); axis(200 600 -50 50);仿真结果输入信号波形15 0 5o-0t m XUQPav输入噪声波形0.2输

25、入信号与噪声叠加波形21 o 1 o 0 trn XUQPav00.5-0.5Variable tVariable t输出信号波形-0.20Variable t输出噪声波形-2Variable t-201002003001 o d thns OLnpa V400Variable t输出信号与输出噪声叠加波形2-2-2-1 0 1Variable t15020WLn OLQPav0050tnsp orbarlav500 IIE0204060Variable w-20 t1150200250Variable t图9仿真结果O-50200300400500600Variable tHEJ业1 liM

26、ltD LMVtJ祁HT三、心得体会通过这次的课程设计，我们觉得最大的收获就是既了解了噪声对信号传输的影响，又回顾了 MATLAB的相关知识。在代码的编制过程中，我分了三步走。第一步，给出一个确定的噪声信号，并观察其对有用信号的影响；第二步，对信号进行调制，采用的方法是将信号乘以一个相干载波； 第三步，对调制信号进行解调，其中涉及到一个低通滤波器的设计。虽然有函数可以直接实现这个 功能，但我们都没有用，而是自己动手设计。我们的基本思路是由低通滤波器的频域特性，通过傅 里叶逆变换，得到低通滤波器的时域表达式，进而通过频域和时域的对称性，得到解调信号，通过 仿真可直观的看出解调信号和原信号的波形。其次，在这次课程设计中，我们第一次明白了零点的重要性。尤其