DSB信号的仿真分析

1、武汉理工大学MATLAB课程设计报告M ATLAB课程设计报告题目基于MATLAB的DSB调制与解调分析专业班级：通信1104班学生姓名:指导教师:武汉理工大学MATLAB课程设计报告武汉理工大学MATLAB课程设计报告MATLAB课程设计任务书学生姓名：专业班级：指导教师：工作单位：题目：基于MATLA啲DSB调制与解调分析设计内容和要求DSB言号的仿真分析调制信号：分别为300Hz正弦信号和矩形信号；载波频率：30kHz; 解调：同步解调；要求：画出以下三种情况下调制信号、已调信号、解调信号的波形、频 谱以及解调器输入输出信噪比的关系曲线；1）调制信号幅度=0.8 X载波幅度；2）调制信号

2、幅度=载波幅度；3）调制信号幅度=1.5 X载波幅度；时间安排2013年12月25日：复习DSB勺原理，初步构想设计的流程。2013年12月26日至28日：程序编写及调试。2013年12月29日：写报告。指导教师签名:武汉理工大学MATLAB课程设计报告目录目录1摘要1Abstract 21. DSB调制与解调原理 31.1DSB调制原理 31.2DSB解调原理与抗噪性能 52. DSB调制解调分析的 MATLAB实现72.1正弦波调制72.1.1调制信号幅度=0.8 X载波幅度 102.1.2调制信号幅度=载波幅度122.1.3调制信号幅度=1.5*载波幅度 142.2矩形波调制152.2.

3、1调制信号幅度=0.8 X载波幅度 182.2.2调制信号幅度=载波幅度202.2.3调制信号幅度=1.5*载波幅度213. 模拟仿真结果分析 234. 小结与体会 24参考文献25武汉理工大学MATLAB课程设计报告摘要调制在通信系统中有十分重要的作用。通过调制，不仅可以进行频谱搬移，把调制信号的频谱搬移到所希望的位置上，从而将调制信号转换成适合于传播的已调信号，而且它对系统的传输有效性和传输的可靠性有着很大的影响，调制方式往往决 定了一个通信系统的性能。MATLAB软件广泛用于数字信号分析，系统识别，时序分析与建模，神经网络、动态仿真等方面有着广泛的应用。本课题利用MATLAB软件对DSB

4、调制解调系统进行模拟仿真，分别利用 300HZ正弦波和矩形波，对30KHZ正弦波进行调 制，观察调制信号、已调信号和解调信号的波形和频谱分布， 并在解调时引入高斯白噪声， 对解调前后信号进行信噪比的对比分析，估计DSB调制解调系统的性能。1武汉理工大学MATLAB课程设计报告AbstractModulati on in com muni cati on systems have an importa nt role. Through the modulati on, not only can move the spectrum, the modulated signal spectrum mo

5、ve to the desired position, which will convert into a modulated signal suitable for transmission of modulated signals, and that its transmission system, the effectiveness and reliability of transmission has a great impact, the modulati on method is ofte n decided on a com muni cati on system perform

6、a nee. MATLAB software is widely used in digital signal analysis, system identification, time series analysis and modeling, neural networks, dynamic simulation have a wide range of applications. This topic using MATLAB software DSB modulation and demodulation system simulation, use, respectively, 30

7、0HZ sine wave and rectangular wave, sine wave modulation of the 30KHZ observed modulated sig nal modulated sig nal and demodulate the sig nal waveform and spectrum distributi on, and in the soluti on white Gaussia n no ise in troduced whe n adjusted for demodulat ing the signal-noise ratio before an

8、d after the comparative analysis, it is estimated DSB modulation and demodulatio n performa nee of the system.#武汉理工大学MATLAB课程设计报告1. DSB调制与解调原理1.1DSB调制原理DSB调制属于幅度调制。幅度调制是用调制信号去控制咼频载波的振幅，使其按调制 信号的规律而变化的过程。设正弦型载波c(t)=Acos( ct),式中：A为载波幅度，A；为载波角频率。根据调制定义，幅度调制信号(已调信号)一般可表示为：Sm (t)=Am(t)cosC c t)(公式1-1 ),其

9、中，m(t)为基带调制信号。设调制信号 m(t)的频谱为M&；)，则由公式1-1不难A得到已调信号(t)的频谱()： SmC )=- M(+c)+M(一c)。由以上表示式可见，在波形上，幅度已调信号随基带信号的规律呈正比地变化；在频谱结构上，它的频谱完全是基带信号频谱在频域内的简单搬移。标准振幅就是常规双边带调制，简称调幅(AM )。假设调制信号m(t)的平均值为0， 将其叠加一个直流偏量后与载波相乘，即可形成调幅信号。其时域表达式为：Sam (t)=Ao+m(t)cos( ct)式中：为外加的直流分量；m(t )可以是确知信号，也可以是随机信号。若为确知信号，则AM信号的频谱为1Sa

10、m ()二二 Ao( 'c) ( ' - c)石 側(c) M ( JAM信号的频谱由载频分量、上边带、下边带三部分组成。AM信号的总功率包括载波功率和边带功率两部分。只有边带功率才与调制信号有关，也就是说，载波分量并不携 带信息。因此，AM信号的功率利用率比较低。AM调制典型波形和频谱如图1-1所示：图1-1 AM调制典型波形和频谱如果在AM调制模型中将直流去掉，即可得到一种高调制效率的调制方式一抑制载波双 边带信号(DSB SC)，简称双边带信号。其时域表达式为 SDSB(t)=m(t)cos(ct) 式中，假设的平均值为0。DSB的频谱与AM的谱相近，只是没有了在处的函数

11、，即1SdsbC - )=2 M( c) M ( - c)-0 O)(旳八仪co轴O) 缁©sdsb(3AA-CQO)Gc3武汉理工大学MATLAB课程设计报告#武汉理工大学MATLAB课程设计报告图1-2 DSB调制典型波形和频谱#武汉理工大学MATLAB课程设计报告与AM信号比较，因为不存在载波分量，DSB信号的调制效率是100 ，即全部效率都 用于信息传输1.2DSB解调原理与抗噪性能解调是调制的逆过程，其作用是从接收的已调信号中恢复原基带信号(即调制信号) 解调的方法可分为两类：相干解调和非相干解调(包络检波)。相干解调，也称同步检波，为了无失真地恢复原基带信号，接收端必须提

12、供一个与接 收的已调载波严格同步(同频同相)的本地载波(称为相干载波)，它与接受的已调信号相乘后，经低通滤波器取出低频分量，即可得到原始的基带调制信号。包络检波器就是直接从已调波的幅度中提取原调制信号，通常由半波或全波整流器和 低通滤波器组成。由于DSB信号的包络不再与调制信号的变化规律一致，因而不能采用简单的包络检波 来恢复调制信号。DSB信号解调时需采用相干解调。DSB相干解调性能分析模型如图1-3所示：图1-3 DSB相干解调性能分析模型设解调器输入信号为Sm (t)= m(t)cos匕t)，与相干载波cos(t)相乘后，得2 1 1 1 m(t)cos ct m(t) m(t)cos

13、2 ct,经低通滤波器后，输出信号为：m°(t) m(t)。1 因此，解调器输出端的有用信号功率为So =m2(t) = -m2(t)4解调DSB信号时，接收机中的带通滤波器的中心频率与调制频率相同，因此解调器输入端的窄带噪声nt二n c t cos - 'ct -ns t sin t ,它与相干载波cos( ct)相乘后，1 1 _ _得 ni t cos ct =2氏 t ' 2 II nc t cos 2 J n$ t sin 2 ct1经低通滤波器后，解调器最终输出噪声为 心=扌亿t1 ii故输出噪声功率为 N。二n/(t)二M二n°B444式中，B

14、=2,为DSB的带通滤波器的带宽，no为噪声单边功率谱密度。1 :解调器输入信号平均功率为S =sm2(t) =m(t)cosi*ct 2 =?m2(t)可得解调器的输入信噪比-更，解调器的输出信噪比Ni n°BSom2(t)Non°B因此制度增益为Gdsb =2，也就是说,DSB信号的解调器使信噪比改善一倍5武汉理工大学MATLAB课程设计报告2. DSB调制解调分析的MATLA实现信号DSB调制采用MATLAB函数modulate实现，其函数格式为：Y = MODULATE(X,fc,fs,METHOD,OPT)X为基带调制信号，fc为载波频率，fs为抽样频率，METH

15、OD为调制方式选择，DSB 调制时为'am', OPT在DSB调制时可不选，fs需满足fs > 2*fc + BW，BW为调制信号带 宽。DSB信号解调采用MATLAB函数demod实现，其函数使用格式为：X = DEMOD(Y ,fc,fs,METHOD,OPT)丫为DSB已调信号，fc为载波频率，fs为抽样频率，METHOD为解调方式选择，DSB 解调时为'am', OPT在DSB调制时可不选。观察信号频谱需对信号进行傅里叶变换，采用 MATLAB函数fft实现，其函数常使用格 式为：Y=FFT(X,N)，X为时域函数，N为傅里叶变换点数选择，一般取值

16、。频域变换后， 对频域函数取模，格式：丫仁ABS(Y)，再进行频率转换，转换方法：f=(0:length(Y)-1) ' f*ength(Y)分析解调器的抗噪性能时，在输入端加入高斯白噪声，采用 MATLAB函数awgn实现， 其函数使用格式为：丫 =AWGN(X,S_N)，加高斯白噪声于X中，S_N为信噪比，单位为dB， 其值在假设X的功率为OdBM的情况下确定。信号的信噪比为信号中有用的信号功率与噪声功率的比值，根据信号功率定义，采用 MATLAB函数var实现，其函数常使用格式为：Y =VAR(X)，返回向量的方差，则信噪比 为：S_N=VAR(X1)/VAR(X2)。绘制曲线采

17、用MATLAB函数plot实现，其函数常使用格式：PLOT(X，丫)，X 为横轴变量，丫为纵轴变量，坐标范围限定AXIS(x1 x2 y1 y2)，轴线说明XLABEL(')和 YLABEL('。2.1正弦波调制用频率300HZ正弦波调制频率30KHZ的正弦波，采用同步解调，观察调制信号、已调信号、解调信号的波形、频谱以及解调器输入输出信噪比的关系MATLAB源程序如下：fc=30000;fs=100000;N=1000;n=0:N-1;t=n/fs;x=A*si n( 2*pi*300*t);y=modulate(x,fc,fs,'am' );%抑制双边带振幅

18、调制fft1=fft(x,N);%傅里叶变换mag 仁 abs(fft1);%取模f1=(0:le ngth(fft1)-1)'*fs/le ngth(fft1);fft2=fft(y,N);mag2=abs(fft2);f2=(0:le ngth(fft2)-1)'*fs/le ngth(fft2);figure(1);subplot(2,2,1); plot(t,x);xlabel('调制信号波形)subplot(2,2,2); plot(f1,mag1); axis(0 600 0 600);xlabel('调制信号频谱)subplot(2,2,3); p

19、lot(t,y);xlabel('已调信号波形)subplot(2,2,4); plot(f2,mag2); axis(28000 32000 0 400);xlabel('已调信号频谱)%yn=awg n(y,4);%加入高斯白噪声znn=demod(y,fc,fs,'am' );%无噪声已调信号解调zn=demod(yn,fc,fs,'am' );%加噪声已调信号解调fft3=fft(z nn ,N);mag3=abs(fft3);f3=(0:le ngth(fft3)-1)'*fs/le ngth(fft3);figure(2);s

20、ubplot(3,1,1);plot(t,z n);xlabel('加噪声解调信号波形)subplot(3,1,2);plot(t,z nn);xlabel('无噪声解调信号波形)subplot(3,1,3); plot(f3,mag3); axis(O 500 0 500);xlabel('解调信号频谱')%yn仁 awg n(y,8);yn 2=awg n(y,12);yn 3=awg n(y,16);yn 4=awg n(y,20);zn 1=demod( yn 1,fc,fs,'am');zn 2=demod( yn 2,fc,fs,&#

21、39;am');zn 3=demod( yn 3,fc,fs,'am');zn4=demod(y n4,fc,fs,'am');dyi=yn_y;%高斯白噪声s_ni=var(y)/var(dyi);%输入信噪比dyo=z n-znn;%解调后噪声s_no=var(z nn)/var(dyo);%俞出信噪比dyi1= yn 1-y;s_n i仁var(y)/var(dyi1);dyo1=z n1-z nn;s_n o仁var(z nn)/var(dyol);dyi2=y n2-y;s_ni 2=var(y)/var(dyi2);dyo2=z n2-z n

22、n;s_no 2=var(z nn)/var(dyo2);dyi3=y n3-y;s_ni 3=var(y)/var(dyi3);dyo3=z n3-z nn;s_no 3=var(z nn)/var(dyo3);dyi4=yn 4-y;s_n i4=var(y)/var(dyi4);dyo4=z n4-znn;s_no4=var(z nn)/var(dyo4);in=s_ ni,s_ ni1,s_ ni 2,s_ni3,s_ ni4;out=s_ no,s_n o1,s_ no 2,s_ no 3,s_ no 4;figure(3);plot(i n, out,'*')hol

23、d onplot(i n, out)xlabel('输入信噪比');ylabel( '输出信噪比')2.1.1调制信号幅度=0.8 X载波幅度调用程序，程序中A=0.8o调制信号、已调信号的波形、频谱如图2-1所示:图2-1调制信号、已调信号的波形、频谱图解调信号的波形、频谱如图2-2所示：图2-2解调信号的波形、频谱图输入输出信噪比关系曲线如图2-3所示：13武汉理工大学MATLAB课程设计报告#武汉理工大学MATLAB课程设计报告图2-3输入输出信噪比关系曲线2.1.2调制信号幅度=载波幅度调用函数，函数中A=1。调制信号、已调信号的波形、频谱如图 2-4所

24、示:图2-4调制信号、已调信号的波形、频谱图15武汉理工大学MATLAB课程设计报告解调信号的波形、频谱如图2-5所示:图2-5解调信号的波形、频谱图#武汉理工大学MATLAB课程设计报告#武汉理工大学MATLAB课程设计报告输入输出信噪比关系曲线如图2-6所示:Q Figure 3105File Edit VieM Insert Tools Desktop window Help D曰k|鬓叫幼SI宾"J剋 DI5 0 5 0 5 0 54 4 3 3 2 2 151015202530辙入信噪比#武汉理工大学MATLAB课程设计报告#武汉理工大学MATLAB课程设计报告图2-6输入

25、输出信噪比关系曲线#武汉理工大学MATLAB课程设计报告17武汉理工大学MATLAB课程设计报告2.1.3调制信号幅度=15*载波幅度调用程序，程序中A=1.5。调制信号、已调信号的波形、频谱如图2-7所示:图2-7调制信号、已调信号的波形、频谱图解调信号的波形、频谱如图2-8所示：图2-8解调信号的波形、频谱图#武汉理工大学MATLAB课程设计报告输入输出信噪比关系曲线如图2-9所示:Figure 319武汉理工大学MATLAB课程设计报告View Insert Tools Desktop Window Helpa口匡i#武汉理工大学MATLAB课程设计报告#武汉理工大学MATLAB课程设计

26、报告输入廳比00804020o20#武汉理工大学MATLAB课程设计报告#武汉理工大学MATLAB课程设计报告图2-9输入输出信噪比关系曲线2.2矩形波调制用频率300HZ矩形波调制频率30KHZ的正弦波，采用同步解调，观察调制信号、已 调信号、解调信号的波形、频谱以及解调器输入输出信噪比的关系。MATLAB源程序如下：fc=30000;%载波频率fs=100000;%抽样频率N=10000;n=0:N-1;t=n/fs;x=A*square(2*pi*300*t,50);y=modulate(x,fc,fs,'am' );%抑制双边带振幅调制fft1=fft(x,N);%傅里

27、叶变换%取模mag1=abs(fft1);f1=(0:le ngth(fft1)-1)'*fs/le ngth(fft1);fft2=fft(y,N);mag2=abs(fft2);f2=(0:le ngth(fft2)-1)'*fs/le ngth(fft2);figure(1);subplot(2,2,1); plot(t,x); axis(0 0.1 -2 2);xlabel('调制信号波形)subplot(2,2,2); plot(f1,mag1); axis(0 5000 0 8000);xlabel('调制信号频谱)subplot(2,2,3); p

28、lot(t,y); axis(0 0.01 -2 2);xlabel('已调信号波形)subplot(2,2,4); plot(f2,mag2); axis(0 50000 0 8000);xlabel( '已调信号频谱)%yn=awg n(y,4);%加入高斯白噪声znn=demod(y,fc,fs,'am' );%无噪声已调信号解调zn=demod(yn,fc,fs,'am' );%加噪声已调信号解调fft3=fft(z nn ,N);mag3=abs(fft3);f3=(0:le ngth(fft3)-1)'*fs/le ngth(

29、fft3);figure(2);subplot(3,1,1); plot(t,z n);xlabel('加噪声解调信号波形)subplot(3,1,2); plot(t,z nn);xlabel('无噪声解调信号波形)subplot(3,1,3); plot(f3,mag3); axis(0 5000 0 4000);xlabel('解调信号频谱)%yn 仁awg n(y,8);yn 2=awg n(y,12);yn 3=awg n(y,16);'am');'am');'am');'am');%高斯白噪声%

30、俞入信噪比%解调后噪声%俞出信噪比yn 4=awg n(y,20);zn 仁demod( yn 1,fc,fs,zn 2=demod( yn 2,fc,fs, zn 3=demod( yn 3,fc,fs, zn4=demod(y n4,fc,fs, dyi=yn_y;s_ni=var(y)/var(dyi); dyo=z n-znn;s_no=var(z nn)/var(dyo); dyi1= yn 1-y;s_n i仁var(y)/var(dyi1); dyo1=z n1-z nn;s_n o仁var(z nn)/var(dyo1);dyi2=y n2-y;s_n i2=var(y)/va

31、r(dyi2);dyo2=z n2-z nn;s_no 2=var(z nn)/var(dyo2);dyi3=y n3-y;s_ni 3=var(y)/var(dyi3);dyo3=z n3-z nn;s_no 3=var(z nn)/var(dyo3);dyi4=yn 4-y;s_n i4=var(y)/var(dyi4);dyo4=z n4-znn;s_no4=var(z nn)/var(dyo4);in=s_ ni,s_ ni1,s_ ni 2,s_ni3,s_ ni4;out=s_ no,s_n o1,s_ no 2,s_ no 3,s_ no 4; figure(3);plot(i

32、n, out,'*')hold on输出信噪比')plot(i n,o ut)xlabel('输入信噪比');ylabel(2.2.1调制信号幅度=0.8 X载波幅度调用程序，程序中A=0.8。调制信号、已调信号的波形、频谱如图2-10所示:23武汉理工大学MATLAB课程设计报告#武汉理工大学MATLAB课程设计报告图2-10调制信号、已调信号的波形、频谱图#武汉理工大学MATLAB课程设计报告#武汉理工大学MATLAB课程设计报告解调信号的波形、频谱如图2-11所示：#武汉理工大学MATLAB课程设计报告图2-12输入输出信噪比关系曲线#武汉理工大学

33、MATLAB课程设计报告图2-12输入输出信噪比关系曲线#武汉理工大学MATLAB课程设计报告图2-12输入输出信噪比关系曲线#武汉理工大学MATLAB课程设计报告图2-11解调信号的波形、频谱图输入输出信噪比关系曲线如图2-12所示：Q Figure 3豆图2-12输入输出信噪比关系曲线#武汉理工大学MATLAB课程设计报告图2-12输入输出信噪比关系曲线#武汉理工大学MATLAB课程设计报告File Edit view Insert Tools Desktop Window Help2 j j k v x o ®'a n s 口图2-12输入输出信噪比关系曲线#武汉理工大

34、学MATLAB课程设计报告图2-12输入输出信噪比关系曲线#武汉理工大学MATLAB课程设计报告5101520输入桶比3035504 32OJ1图2-12输入输出信噪比关系曲线25武汉理工大学MATLAB课程设计报告#武汉理工大学MATLAB课程设计报告2.2.2调制信号幅度=载波幅度调用程序，程序中A=1。调制信号、已调信号的波形、频谱如图 2-13所示:图2-13调制信号、已调信号的波形、频谱图解调信号的波形、频谱如图2-14所示：E Figure 2图2-14解调信号的波形、频谱图27武汉理工大学MATLAB课程设计报告#武汉理工大学MATLAB课程设计报告输入输出信噪比关系曲线如图2-15所示:图2-15输入输出信噪比关系曲线2.2.3调制信号幅度=15*载波幅度调用程序，程序中A=1.5。调制信号、已调信号的波形、频谱如图2-16所示:图2-16调制信号、已调信号的波形、频谱图29武汉理工大学