DSB调制与解调的MATLAB实现及

1、DSB信号的调制与解调信号的调制与解调 摘要摘要本次课程设计用于实现模拟调制信号经DSB调制后的解调过程。信号的调制与解调在通信系统中具有重要的作用。调制过程实际上是一个频谱搬移的过程，即是将低频信号的频谱（调制信号）搬移到载频位置（载波）。而解调是调制的逆过程，即是将已调制信号还原成原始基带信号的过程。信号的接收端就是通过解调来还原已调制信号从而读取发送端发送的信息。因此信号的解调对系统的传输有效性和传输可靠性有着很大的影响。调制与解调方式往往能够决定一个通信系统的性能。双边带DSB调制信号的解调采用相干解调法（即是将已调信号与相同载波频率相乘），这种方式被广泛应用在载波通信和短波无线电话通

2、信中。但是由于在信道传输过程中必将引入高斯白噪声，虽然经过带通滤波器后会使其转化为窄带噪声，但它依然会对解调信号造成影响，即使其失真，而这种失真是不可避免的。DSB调制与解调原理调制与解调原理DSB调制原理调制原理tccos)(tm)(tsDSB图3-7 DSB调制器模型DSB调制器模型ttmtscDSBcos)()( )()(21)(ccDSBMMS图形cos0tOttOm (t)sDSB(t)OtOccM ()OHHSDSB()Occ载波反相点2HDSB解调原理图解调原理图带通滤波器sm(t)sm(t)n(t)ni(t)mo(t)no(t)低通滤波器cosctDSB是AM调制的一种，AM信

3、号通过信道后自然会叠加有燥声，经过接收天线进入带通滤波器。 BPF的作用有两个，一是让AM信号直接通过，二是滤出带外噪声。AM信号通过BFP后与本地载波相乘后。进入LPF，LPF的截止频率设定为一个定值，它不允许频率大于截止频率的成分通过，因此LPF的输出仅为与要的信号。原理简述DSB解调原理 基于基于MATLAB仿真结果仿真结果 无高斯白噪声时 DSB信号调制过程 ts=0.01;%定义变量区间步长t0=2;%定义变量区间终止值t=-t0:ts:t0;%定义变量区间取值情况fc=10;%给出相干载波的频率A=1;%定义调制信号幅度fa=1;%定义调制信号频率mt=A*cos(2*pi*fa.

4、*t);%输入调制信号表达式ct=cos(2*pi*fc.*t);%输入载波信号表达式psnt=mt.*ct;%输出调制信号表达式subplot(5,1,1);%划分画图区间plot(t,mt,g);%画出调制信号波形title(输入信号波形);xlabel(Variable t);ylabel(Variable mt);subplot(5,1,2);plot(t,ct,b);%画出载波信号波形title(输入载波波形);DSB信号调制过程 xlabel(Variable t);ylabel(Variable ct);subplot(5,1,3);plot(1:length(psnt),psn

5、t,r);%length用于长度匹配title(已调信号波形);%画出已调信号波形xlabel(Variable t);ylabel(Variable psnt);mt1=fftshift(fft(mt);mt2=abs(mt1.2);df=1/(2*t0);ff=length(mt1);f=-ff/2*df:df:ff/2*df-df;subplot(5,1,4);plot(f,mt2);%调制信号功率谱密度psnt1=fftshift(fft(psnt);psnt2=abs(psnt1.2);subplot(5,1,5);plot(f,psnt2);%已调信号功率谱密度高斯白噪声信道特性分

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

7、来添加噪声，此函数功能为向信号中添加噪声功率为其方差的高斯白噪声。正弦波通过加性高斯白噪声信道后的信号为： 故其有用信号功率为： 噪声功率为：信噪比 满足公式: ( )cos()( )cr tAtn t22AS 2NSN1 01 0 lo g()SBN高斯白噪声信道特性分析则可得到公式: 2210210BA我们可以通过这个公式方便的设置高斯白噪声的方差。高斯白噪声信道特性：高斯白噪声信道特性：分别在大信噪比和小信噪比两种情况下，已调信号的波形仿真结果如图 -2-1012-2-1012小 信 噪 比 高 斯 白 躁 声Variable tVariable nit-2-1012-4-202叠 加

8、小 信 噪 比 已 调 信 号 波 形Variable tVariable psnt-2-1012-0.2-0.100.10.2大 信 噪 比 高 斯 白 躁 声Variable tVariable nit-2-1012-2-1012叠 加 大 信 噪 比 已 调 信 号 波 形Variable tVariable psmt高斯白噪声信道特性高斯白噪声信道特性 clf;%清除窗口中的图形ts=0.01;%定义变量区间步长t0=2;%定义变量区间终止值t=-t0+0.0001:ts:t0;%定义变量区间fc=10;%给出相干载波的频率A=1;%定义输入信号幅度fa=1;%定义调制信号频率mt=A

9、*cos(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);%

10、由信噪比求方差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(Variable t);ylabel(Variable psnt);subplot(2,2,

11、3);plot(t,nit1,r);%length用于长度匹配title(大信噪比高斯白躁声);%画出输入信号与噪声叠加波形xlabel(Variable t);ylabel(Variable nit);subplot(2,2,4);plot(t,psnt1,k);title(叠加大信噪比已调信号波形);%画出输出信号波形xlabel(Variable t);ylabel(Variable psmt);高斯白噪声信道特性：高斯白噪声信道特性： 由于在解调过程中，当噪声相对于有用信号较大时会会对有用信号造成失真，从而使解调出来的信号不能完全反应调制信号的信息，而信噪比较小时造成的失真比较小，从而

12、反应调制信息也比较准确。调制解调仿真过程：调制解调仿真过程：-2-1012-1-0.500.51输 入 信 号 波 形Variable tVariable mt-2-1012-0.2-0.100.10.2输 入 噪 声 波 形Variable tVariable nit0100200300400-2-1012输 入 信 号 与 噪 声 叠 加 波 形Variable tVariable snit-2-1012-1-0.500.51输 出 信 号 波 形Variable tVariable psmt-2-1012-0.2-0.100.10.2输 出 噪 声 波 形Variable tVariab

13、le pnit-2-1012-2-1012输 出 信 号 与 输 出 噪 声 叠 加 波 形Variable tVariable psnt0204060050100150低 通 滤 波 器 频 域 波 形Variable wVariable htw150200250-20-1001020低 通 滤 波 器 时 域 波 形Variable tVariable psnt200300400500600-50050输 出 信 号 与 输 出 噪 声 叠 加 波 形Variable tVariable jt调制解调仿真过程 clf;%清除窗口中的图形ts=0.01;%定义变量区间步长t0=2;%定义变量

14、区间终止值t=-t0+0.0001: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=sqrt(fangcha).*randn(h,l); %产生高斯白噪声snit=mt+nit; %调制信号与噪声叠加psmt=mt.*cos(2*pi*fc.*t); %输出调制信号表达式pnit=nit

15、.*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); %解调信号的时域表达式subplot(3,3,1); %划分画图区间plot(t,mt,g); %画出输入信号波形title(输入信号波形);xlabel(Variable t);ylabe

16、l(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),snit,r); %length用于长度匹配title(输入信号与噪声叠加波形); %画出输入信号与噪声叠加波形xlabel(Variable t);ylabel(Variable snit);subplot(3,3,4);plot(t,psmt,k);title(输出信号波形); %画出输出信号波形xlabel(

17、Variable t);ylabel(Variable psmt);subplot(3,3,5);plot(t,pnit,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);title(低通滤波器频

18、域波形); %画出低通滤波器频域波形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);title(输出信号与输出噪声叠加波形);%画出输出信号与输出噪声叠加波形xlabel(Variable t);ylabel(Variable jt);axis(200 600 -50 50);调制解调仿真过程：调制解调仿真过程： 调制信号经过调制解调过程之后引入白噪声，分别经过带通滤波器和低通滤波