（通信课程设计）DSB调制与解调电路的MATLAB实现及调制性能分析

1、dsb调制与解调电路的matlab实现及调制性能分析 题 目： dsb调制与解调电路的matlab实现及调制性能分析专业班级: 学生姓名： 指导教师： dsb调制与解调电路的matlab实现及调制性能分析学生姓名： 专业班级： 指导教师：_ 工作单位： _题目: dsb调制与解调电路的matlab实现及调制性能分析 设计内容和要求dsb信号的仿真分析调制信号：分别为300hz正弦信号和矩形信号；载波频率：30khz；解调：同步解调；要求：画出以下三种情况下调制信号、已调信号、解调信号的波形、频 谱以及解调器输入输出信噪比的关系曲线；1）调制信号幅度=0.8载波幅度；2）调制信号幅度=载波幅度；

2、3）调制信号幅度=1.5载波幅度； 初始条件 1）matlab软件 2）数字信号处理基础知识指导教师签名： 目录摘要4abstract51dsb调制与解调原理61.1dsb调制原理61.2dsb解调原理与抗噪性能82dsb调制解调分析的matlab实现102.1正弦波调制102.1.1调制信号幅度=0.8载波幅度142.1.2调制信号幅度=载波幅度162.1.3调制信号幅度=1.5*载波幅度182.2矩形波调制202.2.1调制信号幅度=0.8载波幅度232.2.2调制信号幅度=载波幅度252.2.3调制信号幅度=1.5*载波幅度273模拟仿真结果分析304小结与体会315参考文献32摘要 调

3、制在通信系统中有十分重要的作用。通过调制，不仅可以进行频谱搬移，把调制信号的频谱搬移到所希望的位置上，从而将调制信号转换成适合于传播的已调信号，而且它对系统的传输有效性和传输的可靠性有着很大的影响，调制方式往往决定了一个通信系统的性能。matlab软件广泛用于数字信号分析，系统识别，时序分析与建模，神经网络、动态仿真等方面有着广泛的应用。本课题利用matlab软件对dsb调制解调系统进行模拟仿真，分别利用300hz正弦波和矩形波，对30khz正弦波进行调制，观察调制信号、已调信号和解调信号的波形和频谱分布，并在解调时引入高斯白噪声，对解调前后信号进行信噪比的对比分析，估计dsb调制解调系统的性

4、能。abstractmodulation in communication systems have an important role. through the modulation, not only can move the spectrum, the modulated signal spectrum move to the desired position, which will convert into a modulated signal suitable for transmission of modulated signals, and that its transmissi

5、on system, the effectiveness and reliability of transmission has a great impact, the modulation method is often decided on a communication system performance. matlab software is widely used in digital signal analysis, system identification, time series analysis and modeling, neural networks, dynamic

6、 simulation have a wide range of applications. this topic using matlab software dsb modulation and demodulation system simulation, use, respectively, 300hz sine wave and rectangular wave, sine wave modulation of the 30khz observed modulated signal modulated signal and demodulate the signal waveform

7、and spectrum distribution, and in the solution white gaussian noise introduced when adjusted for demodulating the signal-noise ratio before and after the comparative analysis, it is estimated dsb modulation and demodulation performance of the system.1dsb调制与解调原理1.1dsb调制原理dsb调制属于幅度调制。幅度调制是用调制信号去控制高频载波

8、的振幅，使其按调制信号的规律而变化的过程。设正弦型载波c(t)=acos(ct),式中：a为载波幅度，c为载波角频率。根据调制定义，幅度调制信号（已调信号）一般可表示为：sm(t)=am(t)cos(ct)（公式1-1）,其中，m(t)为基带调制信号。设调制信号m(t)的频谱为m()，则由公式1-1不难得到已调信号sm(t)的频谱sm()：sm()=a2m(+c)+m(-c)。由以上表示式可见，在波形上，幅度已调信号随基带信号的规律呈正比地变化；在频谱结构上，它的频谱完全是基带信号频谱在频域内的简单搬移。标准振幅就是常规双边带调制，简称调幅（am）。假设调制信号m(t)的平均值为0，将其叠加一

9、个直流偏量a0后与载波相乘，即可形成调幅信号。其时域表达式为: sam（t）=(a0+m(t)cos(ct)式中：a0为外加的直流分量；m(t)可以是确知信号，也可以是随机信号。若为确知信号，则am信号的频谱为sm()=a0(+c)+(-c)+ 12m(+c)+m(-c)am信号的频谱由载频分量、上边带、下边带三部分组成。am信号的总功率包括载波功率和边带功率两部分。只有边带功率才与调制信号有关，也就是说，载波分量并不携带信息。因此，am信号的功率利用率比较低。am调制典型波形和频谱如图1-1所示：图1-1 am调制典型波形和频谱如果在am调制模型中将直流a0去掉，即可得到一种高调制效率的调制

10、方式抑制载波双边带信号(dsbsc)，简称双边带信号。 其时域表达式为sdsb(t)=m(t)cos(ct)式中，假设的平均值为0。dsb的频谱与am的谱相近，只是没有了在处的函数，即sm()=12m(+c)+m(-c)其典型波形和频谱如图1-2所示：图1-2 dsb调制典型波形和频谱与am信号比较，因为不存在载波分量，dsb信号的调制效率是100，即全部效率都用于信息传输。1.2dsb解调原理与抗噪性能解调是调制的逆过程，其作用是从接收的已调信号中恢复原基带信号（即调制信号）。解调的方法可分为两类：相干解调和非相干解调（包络检波）。相干解调，也称同步检波，为了无失真地恢复原基带信号，接收端必

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

12、ct),经低通滤波器后，输出信号为: m0(t)= 12m(t)。因此，解调器输出端的有用信号功率为s0=m02(t)=14m2(t)解调dsb信号时，接收机中的带通滤波器的中心频率c与调制频率c相同，因此解调器输入端的窄带噪声ni(t)= nc(t)cos(ct)- ns(t)sin(ct),它与相干载波cos(ct)相乘后，得 ni(t)cos(ct)= 12nc(t)+ 12nc(t)cos(2ct)- ns(t)sin(2ct)经低通滤波器后，解调器最终输出噪声为n0(t)= 12nc(t)故输出噪声功率为 n0=n02(t)=14nc2(t)=14ni=14n0b式中，b=2fh,为

13、dsb的带通滤波器的带宽，n0为噪声单边功率谱密度。解调器输入信号平均功率为si=12m2(t)可得解调器的输入信噪比 sini=12m2(t)n0b，解调器的输出信噪比s0n0=m2(t)n0b因此制度增益为gdsb=s0/n0si/ni=2，也就是说，dsb信号的解调器使信噪比改善一倍。2dsb调制解调分析的matlab实现信号dsb调制采用matlab函数modulate实现，其函数格式为：y = modulate(x,fc,fs,method,opt)x为基带调制信号，fc为载波频率，fs为抽样频率，method为调制方式选择，dsb调制时为am，opt在dsb调制时可不选，fs需满足

14、fs 2*fc + bw，bw为调制信号带宽。dsb信号解调采用matlab函数demod实现，其函数使用格式为：x = demod(y,fc,fs,method,opt)y为dsb已调信号，fc为载波频率，fs为抽样频率，method为解调方式选择，dsb解调时为am，opt在dsb调制时可不选。观察信号频谱需对信号进行傅里叶变换，采用matlab函数fft实现，其函数常使用格式为：y=fft(x,n)，x为时域函数，n为傅里叶变换点数选择，一般取值2n。频域变换后，对频域函数取模，格式：y1=abs(y)，再进行频率转换，转换方法：f=(0:length(y)-1)*fs/length(y

15、)分析解调器的抗噪性能时，在输入端加入高斯白噪声，采用matlab函数awgn实现，其函数使用格式为：y =awgn(x,snr)，加高斯白噪声于x中，snr为信噪比，单位为db，其值在假设x的功率为0dbm的情况下确定。信号的信噪比为信号中有用的信号功率与噪声功率的比值，根据信号功率定义，采用matlab函数var实现，其函数常使用格式为：y =var(x)，返回向量的方差，则信噪比为:snr=var(x1)/var(x2)。 绘制曲线采用matlab函数plot实现，其函数常使用格式：plot（x，y），x为横轴变量，y为纵轴变量，坐标范围限定axis(x1 x2 y1 y2)，轴线说明x

16、label( )和ylabel( )。2.1正弦波调制用频率300hz正弦波调制频率30khz的正弦波，采用同步解调，观察调制信号、已调信号、解调信号的波形、频谱以及解调器输入输出信噪比的关系。matlab源程序如下：fs=100000; %抽样频率 fc=30000; %载波频率n=1000; %fft长度n=0:n-1;t=n/fs; %截止时间和步长x=a*sin(2*pi*300*t); %基带调制信号 y=modulate(x,fc,fs,am); %抑制双边带振幅调制yn=awgn(y,4); %加入高斯白噪声yn1=awgn(y,10);yn2=awgn(y,15);yn3=aw

17、gn(y,20);yn4=awgn(y,25);y1=demod(y,fc,fs,am); %无噪声已调信号解调yyn=demod(yn,30000,fs,am); %加噪声已调信号解调yyn1=demod(yn1,30000,fs,am);yyn2=demod(yn2,30000,fs,am);yyn3=demod(yn3,30000,fs,am);yyn4=demod(yn4,30000,fs,am);dy1=yn-y; %高斯白噪声snr1=var(y)/var(dy1); %输入信噪比dy2=yyn-y1; %解调后噪声snr2=var(y1)/var(dy2); %输出信噪比 dy1

18、1=yn1-y; snr11=var(y)/var(dy11); dy21=yyn1-y1; snr21=var(y1)/var(dy21); dy12=yn2-y; snr12=var(y)/var(dy12); dy22=yyn2-y1; snr22=var(y1)/var(dy22); dy13=yn3-y; snr13=var(y)/var(dy13); dy23=yyn3-y1; snr23=var(y1)/var(dy23); dy14=yn4-y; snr14=var(y)/var(dy14); dy24=yyn4-y1; snr24=var(y1)/var(dy24); in=

19、snr1,snr11,snr12,snr13,snr14; out=snr2,snr21,snr22,snr23,snr24;ff1=fft(x,n); %傅里叶变换mag1=abs(ff1); %取模f1=(0:length(ff1)-1)*fs/length(ff1); %频率转换ff2=fft(y,n);mag2=abs(ff2);f2=(0:length(ff2)-1)*fs/length(ff2);ff3=fft(y1,n);mag3=abs(ff3);f3=(0:length(ff3)-1)*fs/length(ff3);figure(1);subplot(221) %绘制曲线pl

20、ot(t,x)xlabel(调制信号波形)subplot(222)plot(f1,mag1)axis(0 1000 0 1000)xlabel(调制信号频谱)subplot(223)plot(t,y)xlabel(已调信号波形)subplot(224)plot(f2,mag2)axis(0 40000 0 500)xlabel(已调信号频谱)figure(2);subplot(311)plot(t,yyn)xlabel(加噪声解调信号波形)subplot(313)plot(f3,mag3)axis(0 1000 0 600)xlabel(解调信号频谱)subplot(312)plot(t,y1

21、)xlabel(无噪声解调信号波形)figure(3);plot(in,out,*)hold onplot(in,out)xlabel(输入信噪比)ylabel(输出信噪比)2.1.1调制信号幅度=0.8载波幅度调用程序，程序中a=0.8。调制信号、已调信号的波形、频谱如图2-1所示：图2-1 调制信号、已调信号的波形、频谱图解调信号的波形、频谱如图2-2所示：图2-2解调信号的波形、频谱图输入输出信噪比关系曲线如图2-3所示：图2-3 输入输出信噪比关系曲线2.1.2调制信号幅度=载波幅度调用函数，函数中a=1。调制信号、已调信号的波形、频谱如图2-4所示：图2-4调制信号、已调信号的波形、

22、频谱图解调信号的波形、频谱如图2-5所示：图2-5解调信号的波形、频谱图输入输出信噪比关系曲线如图2-6所示：图2-6 输入输出信噪比关系曲线2.1.3调制信号幅度=1.5*载波幅度调用程序，程序中a=1.5。调制信号、已调信号的波形、频谱如图2-7所示：图2-7调制信号、已调信号的波形、频谱图解调信号的波形、频谱如图2-8所示：图2-8解调信号的波形、频谱图输入输出信噪比关系曲线如图2-9所示：图2-9输入输出信噪比关系曲线2.2矩形波调制用频率300hz矩形波调制频率30khz的正弦波，采用同步解调，观察调制信号、已调信号、解调信号的波形、频谱以及解调器输入输出信噪比的关系。matlab源

23、程序如下：clear;f0=300; w0=2*pi*f0; %基带调制信号频率fs=100000; %抽样频率n=10000; %fft长度n=0:n-1;t=n/fs; %截止时间和步长m=a *square(w0*t,50); %基带调制信号y1=fft(m,n); %进行fft变换mag1=abs(y1); %求幅值f1=(0:length(y1)-1)*fs/length(y1); %进行对应的频率转换y=modulate(m,30000,fs,am); %信号抑制载波双边带幅度调制yn=awgn(y,5); %加高斯白噪声于y中yn1=awgn(y,10);yn2=awgn(y,1

24、5);yn3=awgn(y,20);yn4=awgn(y,25);dy1=yn-y; %高斯白噪声snr1=var(y)/var(dy1); %输入信噪比yyn=demod(yn,30000,fs,am); %加噪声已调信号解调yyn1=demod(yn1,30000,fs,am);yyn2=demod(yn2,30000,fs,am);yyn3=demod(yn3,30000,fs,am);yyn4=demod(yn4,30000,fs,am);yy=demod(y,30000,fs,am); %无噪声已调信号dy2=yyn-yy; %解调后输出噪声snr2=var(yy)/var(dy2)

25、; %输出信噪比dy11=yn1-y; snr11=var(y)/var(dy11); dy21=yyn1-yy; snr21=var(yy)/var(dy21); dy12=yn2-y; snr12=var(y)/var(dy12); dy22=yyn2-yy; snr22=var(yy)/var(dy22); dy13=yn3-y; snr13=var(y)/var(dy13); dy23=yyn3-yy; snr23=var(yy)/var(dy23); dy14=yn4-y; snr14=var(y)/var(dy14); dy24=yyn4-yy; snr24=var(yy)/var

26、(dy24); %输出信噪比in=snr1,snr11,snr12,snr13,snr14; out=snr2,snr21,snr22,snr23,snr24; %输入输出信噪比关系y2=fft(y,n); %进行fft变换mag2=abs(y2); %求幅值f2=(0:length(y2)-1)*fs/length(y2); %进行对应的频率转换yy2=fft(yy,n); %进行fft变换mag3=abs(yy2); %求幅值f3=(0:length(yy2)-1)*fs/length(yy2); %进行对应的频率转换figure(1); %绘制曲线subplot(221) plot(t,

27、m)axis(0 0.1 -2 2)xlabel(调制信号波形)subplot(222)plot(f1,mag1)axis(0 5000 0 10000)xlabel(调制信号频谱)subplot(223)plot(t,y)axis(0 0.004 -2 2)xlabel(已调信号波形)subplot(224)plot(f2,mag2)axis(0 50000 0 5000)xlabel(已调信号频谱)figure(2);subplot(311)plot(t,yy)axis(0 0.05 -2 2)xlabel(无噪声解调信号波形)subplot(312)plot(t,yyn)axis(0 0

28、.05 -2 2)xlabel(有噪声解调信号波形)subplot(313)plot(f3,mag3)axis(0 5000 0 5000)xlabel(解调信号频谱)figure(3);plot(in,out,*)hold onplot(in,out)xlabel(输入信噪比)ylabel(输出信噪比)2.2.1调制信号幅度=0.8载波幅度调用程序，程序中a=0.8。调制信号、已调信号的波形、频谱如图2-10所示：图2-10 调制信号、已调信号的波形、频谱图解调信号的波形、频谱如图2-11所示：图2-11 解调信号的波形、频谱图输入输出信噪比关系曲线如图2-12所示：图2-12 输入输出信噪比关系曲线2.2.2调制信号幅度=载波幅度调用程序，程序中a=1。调制信号、已调信号的波形、频谱如图2-13所示：图2-13 调制信号、已