DSB调制与解调

1、DSB调制与解调1 课程设计目的本课程设计是实现DSB 的调制解调。在此次课程设计中，我将通过多方搜集资料与分析，来理解DSB调制解调的具体过程和它在MATLAB中的实现方法。预期通过这个阶段的研习，更清晰地认识DSB 的调制解调原理，同时加深对MATLAB这款通信仿真软件操作的熟练度，并在使用中去感受MATLAB的应用方式与特色。利用自主的设计过程来锻炼自己独立思考，分析和解决问题的能力，为我今后的自主学习研究提供具有实用性的经验。2 课程设计要求（1）熟悉 MATLAB中 M 文件的使用方法，掌握DSB信号的调制解调原理，以此为基础用 M 文件编程实现 DSB信号的调制解调。（2）绘制出

2、SSB信号调制解调前后在时域和频域中的波形，观察两者在解调前后的变化，通过对分析结果来加强对DSB信号调制解调原理的理解。（3）对信号分别叠加大小不同的噪声后再进行解调，绘制出解调前后信号的时域和频域波形， 比较未叠加噪声时和分别叠加大小噪声时解调信号的波形有何区别，由所得结果来分析噪声对信号解调造成的影响。（4）在老师的指导下，独立完成课程设计的全部内容，并按要求编写课程设计论文，文中能正确阐述和分析设计和实验结果。3 相关知识在 AM 信号中，载波分量并不携带信息，信息完全由边带传送。如果将载波抑制，只需在将直流 A0 去掉，即可输出抑制载波双边带信号，简称双边带信号（ DSB）。 DSB

3、调制器模型如图 1 所示。图 1DSB调制器模型其中，设正弦载波为c(t) A cos( ct0)式中， A 为载波幅度；c 为载波角频率；0 为初始相位（假定0为 0）。调制过程是一个频谱搬移的过程，它是将低频信号的频谱搬移到载频位置。而解调是将位于载频的信号频谱再搬回来，并且不失真地恢复出原始基带信号。双边带解调通常采用相干解调的方式，它使用一个同步解调器， 即由相乘器和低通滤波器组成。 在解调过程中， 输入信号和噪声可以分别单独解调。相干解调的原理框图如图2 所示：图 2相干解调器的数学模型信号传输信道为高斯白噪声信道，其功率为2 。4 课程设计分析DSB信号调制过程分析假定调制信号 m

4、(t ) 的平均值为 0，与载波相乘， 即可形成 DSB信号，其时域表达式为sDSBm(t )cosct式中，m(t )的平均值为0。DSB的频谱为sDSB( )1M (c)M (c )2DSB信号的包络不再与调制信号的变化规律一致，因而不能采用简单的包络检波来恢复调制信号，需采用相干解调 (同步检波 )。另外，在调制信号m(t ) 的过零点处，高频载波相位有180°的突变。除了不再含有载频分量离散谱外，DSB信号的频谱与 AM 信号的频谱完全相同，仍由上下对称的两个边带组成。 所以 DSB信号的带宽与 AM 信号的带宽相同，也为基带信号带宽的两倍，即BB2 fHDSBAM式中， f

5、H 为调制信号的最高频率。调制信号产生的代码及波形为clf;%清除窗口中的图形ts=;%定义变量区间步长t0=2;%定义变量区间终止值t=-t0+:ts:t0;%定义变量区间fc=10;%给出相干载波的频率A=1;%定义输入信号幅度fa=1;%定义调制信号频率mt=A*cos(2*pi*fa.*t);%输入调制信号表达式ct=cos(2*pi*fc.*t);%输入调制信号表达式psnt=mt.*cos(2*pi*fc.*t);%输出调制信号表达式subplot(3,1,1);%划分画图区间plot(t,mt,'g');%画出输入信号波形title(' 输入信号波形

6、9;);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:length(psnt),psnt,'r');%length 用于长度匹配title(' 已调信号波形');%画出已调信号波形xlabe

7、l('Variable t');ylabel('Variable psnt');运行结果：输入信号波形1t 0.5melb0airaV-0.5-1-1.5-1-0.500.511.52-2Variable t1输入载波波形t0.5celb0airaV-0.5-1-1.5-1-0.500.511.52-2Variable t2已调信号波形t1nspel0bairaV -1 -2050100150200250300350400Variable t图 3调制信号、载波、已调信号波形高斯白噪声信道特性分析在实际信号传输过程中， 通信系统不可避免的会遇到噪声， 例如自然界

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

9、中添加噪声功率为其方差的高斯白噪声。正弦波通过加性高斯白噪声信道后的信号为r (t )A cos( ct)n(t )故其有用信号功率为S噪声功率为A222N信噪比 S N 满足公式B10log 10(S N )则可得到公式2 A2B2 ?1010我们可以通过这个公式方便的设置高斯白噪声的方差。为了便于比较， 我显示了双边带信号加入两种噪声后的时频波形图。实现代码和波形如图 4：clf;%清除窗口中的图形ts=;%定义变量区间步长t0=2;%定义变量区间终止值t=-t0+:ts:t0;%定义变量区间fc=10;%给出相干载波的频率A=1;%定义输入信号幅度fa=1;%定义调制信号频率mt=A*c

10、os(2*pi*fa.*t);%输入调制信号表达式xzb=2;%输入小信躁比 (dB)snr=10.(xzb/10);h,l=size(mt);%求调制信号的维数fangcha=A*A./(2*snr);%由信躁比求方差nit=sqrt(fangcha).*randn(h,l);%产生小信噪比高斯白躁声psmt=mt.*cos(2*pi*fc.*t);%输出调制信号表达式psnt=psmt+nit;%输出叠加小信噪比已调信号波形xzb=20;%输入大信躁比(dB)snr1=10.(xzb/10);h,l=size(mt);%求调制信号的维数fangcha1=A*A./(2*snr1);%由信躁

11、比求方差nit1=sqrt(fangcha1).*randn(h,l);%产生大信噪比高斯白躁声psnt1=psmt+nit1;%输出已调信号波形subplot(2,2,1);%划分画图区间plot(t,nit,'g');%画出输入信号波形title(' 小信噪比高斯白躁声');xlabel('Variable t');ylabel('Variable nit');subplot(2,2,2);plot(t,psnt,'b');title(' 叠加小信噪比已调信号波形');xlabel('V

12、ariable t');ylabel('Variable psnt');subplot(2,2,3);plot(t,nit1,'r');%length 用于长度匹配title(' 大信噪比高斯白躁声');%画出输入信号与噪声叠加波形xlabel('Variable t');ylabel('Variable nit');subplot(2,2,4);plot(t,psnt1,'k');title(' 叠加大信噪比已调信号波形');%画出输出信号波形xlabel('Var

13、iable t');ylable(Variable psmt);2小信噪比高斯白躁声3叠加小信噪比已调信号波形1.521t1tni0.5snepleb0l0abiariarV-0.5aV-1-1-1.5-2-2-1.5-1-0.500.511.52-3-1.5-1-0.500.511.52-2-2Variable tVariable t0.03大信噪比高斯白躁声1.5叠加大信噪比已调信号波形0.021t0.01t0.5minsepleb0l0abiarairV-0.01a-0.5V-0.02-1-0.03-1.5-1-0.500.511.52-1.5-1.5-1-0.500.511.5

14、2-2-2Variable tVariable t图 4 不同信噪比的噪声及含噪声的已调波形可以清晰地看出， 加大噪声后， 解调信号的波形杂乱无章， 起伏远大于加小噪声时的波形。造成此现象的原因是当信噪比较小时，噪声的功率在解调信号中所占比重较大，所以会造成杂波较多的情况；而信噪比很大时， 噪声的功率在解调信号中所占比重就很小了，噪声部分造成的杂乱波形相对就不是很明显，甚至可以忽略。DSB解调过程分析所谓相干解调是为了从接收的已调信号中，不失真地恢复原调制信号， 要求本地载波和接收信号的载波保证同频同相。相干解调的一般数学模型如图所示。图 5DSB相干解调模型设图四的输入为 DSB信号S (t

15、 )S (t)m(t )cos( t0)mDSBc乘法器输出为(t) S (t )m(t)cos(ct0)cos(ct)DSB1 m(t)cos(0 ) cos(2ct0)2通过低通滤波器后m0(t )1 m(t )cos(0)2当0常数时，解调输出信号为m (t )1m(t)02大小不同信噪比的解调波形，如图6：大信噪比解调信号波形50tjelba0riaV-50250300350400450500550600200Variable t小信噪比解调信号波形50t1jelbai0raV-50250300350400450500550600200Variable t图 6不同信噪比解调波形DSB

16、调制解调系统抗噪声性能分析由于加性噪声只对已调信号的接收产生影响， 因而调制系统的抗噪声性能主要用解调器的抗噪声性能来衡量。 为了对不同调制方式下各种解调器性能进行度量，通常采用信噪比增益 G（又称调制制度增益）来表示解调器的抗噪声性能。有加性噪声时解调器的数学模型如图7 所示。图 7有加性噪声时解调器的数学模型图 7 中 Sm( t ) 为已调信号， n(t ) 为加性高斯白噪声。Sm(t ) 和 n(t ) 首先经过带通滤波器，滤出有用信号， 滤除带外的噪声。 经过带通滤波器后到达解调器输入端的信号为 Sm(t )、噪声为高斯窄带噪声ni (t ) ，显然解调器输入端的噪声带宽与已调信号的

17、带宽是相同的。最后经解调器解调输出的有用信号为mo (t ) ，噪声为no (t ) 。图 8 有加性噪声时解调器的数学模型设解调器输入信号为sm(t)m(t)cosct与相干载波 cos ct 相乘后，得m(t)cos21m(t )1ctm(t)cos 2 t2c2经低通滤波器后，输出信号为1mo(t )m(t)2因此，解调器输出端的有用信号功率为So mo2 (t )1m2 (t)4解调 DSB信号时，接收机中的带通滤波器的中心频率o 与调制载频c 相同，因此解调器输出端的窄带噪声n (t) 可表示为ini(t ) nc(t)cos ctns(t )sinct它与相干载波相乘后，得ni(t

18、)cosct nc(t )cosctns(t)sinct11ctns(t )sin2ct nc(t)nc(t )cos 222经低通滤波器后，解调器最终的输出噪声为1no (t )nc(t)2故输出噪声功率为No no2 (t )1nc2 (t)1noB44这里， B 2 fH ，为 DSB信号的带通滤波器的带宽。解调器输入信号平均功率为Si sm2 (t)m(t)cosct 21m2 (t)2可得解调器的输入信噪比i1 m2 (t)2SNinoB同时可得解调器的输出信噪比12So4m(t)m 2 (t )No1on BNi4因此制度增益为SoGDSBNo2SiNi由此可见， DSB调制系统的

19、制度增益为 2。也就是说 DSB信号的解调器使信噪比改善了一倍。 这是因为采用相干解调， 使输入噪声中的正交分量 ns(t) 被消除的缘故。5 仿真源程序：clf;%清除窗口中的图形ts=;%定义变量区间步长t0=2;%定义变量区间终止值t=-t0+:ts:t0;%定义变量区间fc=10;%给出相干载波的频率A=1;%定义输入信号幅度fa=1;%定义调制信号频率mt=A*cos(2*pi*fa.*t);%输入调制信号表达式xzb=20;%输入信噪比(dB)snr=10.(xzb/10);h,l=size(mt);%求调制信号的维数fangcha=A*A./(2*snr);%由信躁比求方差nit

20、=sqrt(fangcha).*randn(h,l);%产生高斯白噪声snit=mt+nit;%调制信号与噪声叠加psmt=mt.*cos(2*pi*fc.*t);%输出调制信号表达式pnit=nit.*cos(2*pi*fc.*t);%输出噪声表达式psnt=psmt+pnit;%输出已调信号波形jic=psnt.*cos(2*pi*fc.*t);%调制信号乘以相干载波ht=(2*pi*fc.*sin(2*pi*fc.*t)./(2*pi*fc.*t)./pi;%低通滤波器的时域表达式htw=abs(fft(ht);%低通滤波器的频域表达式jt=conv(ht,jic);%解调信号的时域表达

21、式subplot(3,3,1);%划分画图区间plot(t,mt,'g');%画出输入信号波形title(' 输入信号波形');xlabel('Variable t');ylabel('Variable mt');subplot(3,3,2);plot(t,nit,'b');title(' 输入噪声波形');xlabel('Variable t');ylabel('Variable nit');subplot(3,3,3);plot(1:length(snit),sn

22、it,'r');%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,pnit,&#

23、39;k');title(' 输出噪声波形 ');%画出输出噪声波形xlabel('Variable t');ylabel('Variable pnit');subplot(3,3,6);plot(t,psnt,'k');title(' 输出信号与输出噪声叠加波形');%画出输出信号与输出噪声叠加波形xlabel('Variable t');ylabel('Variable psnt');subplot(3,3,7);plot(1:length(htw),htw,'k

24、');title(' 低通滤波器频域波形 ');%画出低通滤波器频域波形xlabel('Variable w');ylabel('Variable htw');axis(0 60 0 150);subplot(3,3,8);plot(1:length(ht),ht,'k');title(' 低通滤波器时域波形 ');%画出低通滤波器时域波形xlabel('Variable t');ylabel('Variable psnt');axis(150 250 -20 25);%给出坐标轴范围subplot(3,3,9);plot(1:length(jt),jt,'k');ti