基于matlab的AM和DSB信号的调制与解调精品资料.doc

《基于matlab的AM和DSB信号的调制与解调精品资料.doc》由会员分享，可在线阅读，更多相关《基于matlab的AM和DSB信号的调制与解调精品资料.doc（29页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、塔里木大学信息工程学院课程设计目录1前言12正文12.1设计的目的12.2设计的内容12.3设计的原理12.3.1 AM信号的调制原理12.3.2 AM信号的解调原理22.3.3 DSB信号的调制原理32.3.4 DSB信号解调原理32.4AM信号与DSB信号的仿真52.4.1 AM信号的仿真52.4.2 DSB信号的仿真62.5结论83致谢84参考文献85附录91前言调制就是使一个信号（如光、高频电磁振荡等）的某些参数（如振幅、频率等）按照另一个欲传输的信号（如声音、图像等）的特点变化的过程。用所要传播的语言或音乐信号去改变高频振荡的幅度，使高频振荡的幅度随语言或音乐信号的变化而变化，这个控

2、制过程就称为调制。其中语言或音乐信号叫做调制信号，调制后的载波就载有调制信号所包含的信息，称为已调波。解调是调制的逆过程，它的作用是从已调波信号中取出原来的调制信号。对于幅度调制来说，解调是从它的幅度变化提取调制信号的过程。对于频率调制来说，解调是从它的频率变化提取调制信号的过程。频率解调要比幅度解调复杂，用普通检波电路是无法解调出调制信号的，必须采用频率检波方式，如各类鉴频器电路。关于鉴频器电路可参阅有关资料，这里不再细述。本课题利用MATLAB软件对AM信号DSB信号调制解调系统进行模拟仿真，分别对正弦波进行调制，观察调制信号、已调信号和解调信号的波形和频谱分布，并在解调时引入高斯白噪声，

3、对解调前后信号进行信噪比的对比分析，估计AM信号和DSB信号的调制解调系统的性能。2正文2.1设计的目的1.熟悉掌握AM与DSB信号的调制与解调原理。2.利用Matlab软件进行振幅调制和解调程序设计，输出显示调制信号、以调信号以及解调信号波形，并输出显示三种信号频谱图。3.对产生波形进行分析，并通过参数的改变，观察波形变化，分析实验现象。2.2设计的内容1.设计AM与DSB信号实现的Matlab程序，输出调制信号、载波信号以及已调信号波形以及频谱图，并改变参数观察信号变化情况，进行实验分析。2.设计AM与DSB信号解调实现的Matlab程序，输出并观察解调信号波形，分析实验现象。2.3设计的

4、原理2.3.1 AM信号的调制原理 AM是指调制信号去控制高频载波的幅度，使其随调制信号呈线性变化的过程为基带信号，它可以是确知信号，也可以是随机信号，但通常认为它的平均值为0.载波为：上式中，为载波振幅，为载波角频率为载波的初始相位。2.3.2 AM信号的解调原理解调是将位于载波的信号频谱再搬回来，并且不失真的恢复出原始基带信号。解调的方式有两种：相干解调与非相干解调。相干解调是利用乘法器，输入一路与载频相干（同频同相）的参考信号与载频相乘；适用于各种线性调制系统，非相干解调一般适用幅度调制（AM）信号。本文中用到的AM信号解调是相干解调图中为已调信号，n(t)为加性高斯白噪声。 和n(t)

5、首先经过一带通滤波器，滤出有用信号，滤除带外的噪声。经过带通滤波器后到达解调器输入端的信号为，噪声为高斯窄带噪声，显然解调器输入端的噪声带宽与已调信号的带宽是相同的。最后经解调器解调输出的有用信号为，噪声为。AM调制系统的相干解调模型如下图所示。图2.1 AM调制系统的相干解调模型AM信号的时域表达式为 通过分析可得AM信号的平均功率为又已知输入功率, 其中B表示已调信号的带宽。由此可得AM信号在解调器的输入信噪比为AM信号经相干解调器的输出信号为因此解调后输出信号功率为在上图中输入噪声通过带通滤波器之后，变成窄带噪声，经乘法器相乘后的输出噪声为经LPF后， 因此解调器的输出噪声功率为可得AM

6、信号经过解调器后的输出信噪比为由上面分析的解调器的输入、输出信噪比可得AM信号的信噪比增益为2.3.3 DSB信号的调制原理DSB调制属于幅度调制。幅度调制是用调制信号去控制高频载波的振幅，使其按调制信号的规律而变化的过程。设正弦型载波,式中：为载波幅度, 为载波角频率。根据调制定义，幅度调制信号（已调信号）一般可表示为： 其中，为基带调制信号。设调制信号的频谱为，则得到已调信号的频谱 由以上表示式可见，在波形上，幅度已调信号随基带信号的规律呈正比地变化；在频谱结构上，它的频谱完全是基带信号频谱在频域内的简单搬移。2.3.4 DSB信号解调原理 解调是调制的逆过程，其作用是从接收的已调信号中恢

7、复原基带信号（即调制信号）。解调的方法可分为两类：相干解调和非相干解调（包络检波）。相干解调，也称同步检波，为了无失真地恢复原基带信号，接收端必须提供一个与接收的已调载波严格同步（同频同相）的本地载波（称为相干载波），它与接受的已调信号相乘后，经低通滤波器取出低频分量，即可得到原始的基带调制信号。包络检波器就是直接从已调波的幅度中提取原调制信号，通常由半波或全波整流器和低通滤波器组成。由于DSB信号的包络不再与调制信号的变化规律一致，因而不能采用简单的包络检波来恢复调制信号。DSB信号解调时需采用相干解调。DSB相干解调性能分析模型如图1.3所示：图2.2 DSB相干解调性能分析模型 设解调器

8、输入信号为，与相干载波相乘后，得,经低通滤波器后，输出信号为: 。因此，解调器输出端的有用信号功率为 解调DSB信号时，接收机中的带通滤波器的中心频率与调制频率c相同，因此解调器输入端的窄带噪声,它与相干载波相乘后，得: 经低通滤波器后，解调器最终输出噪声为 故输出噪声功率为 式中，,为DSB的带通滤波器的带宽，为噪声单边功率谱密度。解调器输入信号平均功率为 可得解调器的输入信噪比 解调器的输出信噪比 因此制度增益为，也就是说，DSB信号的解调器使信噪比改善一倍。2.4AM信号与DSB信号的仿真2.4.1 AM信号的仿真 图2.3 AM信号的调制信号与已调信号图 2.4 滤波前后的频谱图 2.

9、5 调制信号与解调信号2.4.2 DSB信号的仿真图 2.6 DSB信号的调制信号与已调信号图 2.7 滤波前后的频谱图 2.8 调制信号与解调信号2.5结论 在设计AM与DSB信号的调制与解调时选取的载波频率为25Hz，调制频率为5Hz。 通过用MATLAB对AM频谱的分析，AM信号由载频分量，上边带，下边带三部分组成。上边带的频谱结构与原调制信号的频谱结构相同，下边带是上边带的镜像，因此，AM信号时有载波分量的双边带信号，其带宽是基带信号带宽的两倍。通过MATLAB对DSB调制和解调系统的模拟仿真，观察各波形和频谱，在波形上，已调信号的幅度随基带信号的规律呈正比地变化；在频谱结构上，它的频

10、谱完全是基带信号频谱在频域内的简单搬移。通过在已调信号中加入高斯白噪声，通过解调器解调，根据对输入输出信噪比关系曲线绘制观察，在理想情况下，输出信噪比为输入信噪比的二倍，即DSB信号的解调器使信噪比改善一倍；不同的调制信号对系统性能有一定的影响。在频谱图上可以看出AM信号与DSB信号的的相同点是AM信号与DSB信号都是双边带的。不同点就是AM信号的调制信号与已调信号的频谱上都有载频分量，而DSB信号的调制信号频谱与已调信号的频谱图上没有载频分量。3致谢两周的课程设计已经接近尾声，虽然时间很短，但还是收获很多。我做的是AM与DSB信号的调制与解调，在这次设计中，不仅对课本上的知识有了更深层次的了

11、解，也对matlab的运用有了一定的掌握。本次课题设计中遇到了一些难点，与许多同学进行了讨论，除此之外，多亏了蒋老师和张老师的认真讲解，是我最终完成了这篇论文，在这里我要对老师与同学给予我的帮助表示感谢。4参考文献1樊昌信,曹丽娜.通信原理M.北京:国防工业出版社2011.82徐幸然,程玲.现代电子技术M. 南京信息工程大学20123李建新,刘乃安,刘继平.现代通信系统分析与仿真MATLAB通信工具箱.西安:西安电子科技大学出版社，20004邓华. MATLAB通信仿真及应用实例详解.北京:人民邮电出版社20035蓝洋 浅述软件无线电中的解调调制理论 南京大学 2008（27）23-256宋辉

12、. 通信信号的特征分析 自动识别与参数提取D. 南京理工大学, 2003.7胡广书. 现代信号处理M. 北京：清华大学出版社, 2004.8罗明. 数字通信信号的自动识别与参数估计研究D. 西安电子科技大学, 2005.5附录AM信号调制程序 t=0:0.01:1;A1=3; %直流偏量m=3*sin(10*pi*t); %调制信号yk1=fft(m,512); %fuliyebianhuanyw1=2*pi*abs(fftshift(yk1);fw1=-256:255/512*100;plot(fw1,yw1)subplot(221)plot(t,m)xlabel(t/s);ylabel(幅

13、度);title(调制信号)subplot(222)plot(fw1,yw1)xlabel(频率/HZ)title(调制信号的频谱)grid;c=cos(50*pi*t); %载波信号s=(A1+m).*c; %已调信号s=awgn(s,5);yk2=fft(s,512);yw2 = 2*pi*abs(fftshift(yk2);fw2=-256:255/512*100;subplot(223)plot(t,s)xlabel(t/s);ylabel(幅度);title(已调信号)axis(0,1,-10,10);subplot(224)plot(fw2,yw2)xlabel(频率/HZ)tit

14、le(已调信号的频谱)grid;AM信号解调程序t=0:0.001:1;A=3;m=3*sin(10*pi*t); c=cos(2*25*pi*t);s=(A+m).*c; %已调信号sa=awgn(s,5); %已调信号加zao声S=sa.*c;% tuoyuan滤波器的设计Rp=0.1; %通带最大衰减为0.1dB Rs=40; %阻带最小衰减为40dB Wp=5/100; %通带截止频率Ws=10/100; %阻带起始频率n,Wn = ellipord(Wp,Ws,Rp,Rs)b,a = ellip(n,Rp,Rs,Wn);Sh=filter(b,a,S); %对信号进行滤波figure

15、(1)plot(t,m,t,Sh-A)legend(原调制信号,解调后的调制信号)title(调制信号与解调信号对比)figure(2)yk1=fft(S,512);yw1 = 2*pi*abs(fftshift(yk1);fw1=-256:255/512;yk2=fft(Sh,512);yw2 = 2*pi*abs(fftshift(yk2);fw2=-256:255/512;subplot(211)plot(fw1,yw1)legend(滤波前的频谱图)subplot(212)plot(fw2,yw2)legend(滤波后的频谱图)DSB信号调制程序t=0:0.01:1;m=3*sin(1

16、0*pi*t); %调制信号yk1=fft(m,512); %fuliyebianhuanyw1=2*pi*abs(fftshift(yk1);fw1=-256:255/512*100;plot(fw1,yw1)subplot(221)plot(t,m)xlabel(t/s);ylabel(幅度);title(调制信号)subplot(222)plot(fw1,yw1)xlabel(频率/HZ)title(调制信号的频谱)grid;c=cos(50*pi*t); %载波信号s=m.*c; %已调信号s=awgn(s,5);yk2=fft(s,512);yw2 = 2*pi*abs(fftshi

17、ft(yk2);fw2=-256:255/512*100;subplot(223)plot(t,s)xlabel(t/s);ylabel(幅度);title(已调信号)axis(0,1,-10,10);subplot(224)plot(fw2,yw2)xlabel(频率/HZ)title(已调信号的频谱)grid;DSB信号解调程序t=0:0.001:1;A=3;m=3*sin(10*pi*t); c=cos(50*pi*t);s=m.*c; %已调信号sa=awgn(s,5); %已调信号加zao声S=sa.*c;% 滤波器的设计Rp=0.1; %通带最大衰减为0.1dB Rs=40; %阻

18、带最小衰减为40dB Wp=5/100; %通带截止频率Ws=10/100; %阻带起始频率n,Wn = ellipord(Wp,Ws,Rp,Rs)b,a = ellip(n,Rp,Rs,Wn);Sh=filter(b,a,S); %对信号进行滤波figure(1)plot(t,m,t,Sh-A)legend(原调制信号,解调后的调制信号)title(调制信号与解调信号对比)figure(2)yk1=fft(S,512);yw1 = 2*pi*abs(fftshift(yk1);fw1=-256:255/512;yk2=fft(Sh,512);yw2 = 2*pi*abs(fftshift(y

19、k2);fw2=-256:255/512;subplot(211)plot(fw1,yw1)legend(滤波前的频谱图)subplot(212)plot(fw2,yw2)legend(滤波后的频谱图)附录资料：MATLAB的30个方法1 内部常数pi 圆周率 exp(1)自然对数的底数ei 或j 虚数单位Inf或 inf 无穷大 2 数学运算符a+b 加法a-b减法a*b矩阵乘法a.*b数组乘法a/b矩阵右除ab矩阵左除a./b数组右除a.b数组左除ab 矩阵乘方a.b数组乘方-a负号 共轭转置.一般转置3 关系运算符=等于大于=大于或等于=不等于4 常用内部数学函数 指数函数exp(x)以

20、e为底数对数函数log(x)自然对数，即以e为底数的对数log10(x)常用对数，即以10为底数的对数log2(x)以2为底数的x的对数开方函数sqrt(x)表示x的算术平方根绝对值函数abs(x)表示实数的绝对值以及复数的模三角函数（自变量的单位为弧度）sin(x)正弦函数cos(x)余弦函数tan(x)正切函数cot(x)余切函数sec(x)正割函数csc(x)余割函数反三角函数 asin(x)反正弦函数acos(x)反余弦函数atan(x)反正切函数acot(x)反余切函数asec(x)反正割函数acsc(x)反余割函数双曲函数 sinh(x)双曲正弦函数cosh(x)双曲余弦函数tan

21、h(x)双曲正切函数coth(x)双曲余切函数sech(x)双曲正割函数csch(x)双曲余割函数反双曲函数 asinh(x)反双曲正弦函数acosh(x)反双曲余弦函数atanh(x)反双曲正切函数acoth(x)反双曲余切函数asech(x)反双曲正割函数acsch(x)反双曲余割函数求角度函数atan2(y,x)以坐标原点为顶点，x轴正半轴为始边，从原点到点（x，y）的射线为终边的角，其单位为弧度，范围为（ ， 数论函数gcd(a,b)两个整数的最大公约数lcm(a,b)两个整数的最小公倍数排列组合函数factorial(n)阶乘函数，表示n的阶乘 复数函数 real(z)实部函数ima

22、g(z)虚部函数abs(z)求复数z的模angle(z)求复数z的辐角，其范围是（ ， conj(z)求复数z的共轭复数求整函数与截尾函数ceil(x)表示大于或等于实数x的最小整数floor(x)表示小于或等于实数x的最大整数round(x)最接近x的整数最大、最小函数max(a，b，c，)求最大数min(a，b，c，)求最小数符号函数 sign(x)5 自定义函数-调用时：“返回值列=M文件名（参数列）”function 返回变量=函数名（输入变量） 注释说明语句段（此部分可有可无）函数体语句 6进行函数的复合运算compose(f,g) 返回值为f(g(y)compose(f,g,z)

23、返回值为f(g(z)compose(f,g,x,.z) 返回值为f(g(z)compose(f,g,x,y,z) 返回值为f(g(z)7 因式分解syms 表达式中包含的变量 factor(表达式) 8 代数式展开syms 表达式中包含的变量 expand(表达式)9 合并同类项syms 表达式中包含的变量 collect(表达式，指定的变量)10 进行数学式化简syms 表达式中包含的变量 simplify(表达式)11 进行变量替换syms 表达式和代换式中包含的所有变量 subs(表达式，要替换的变量或式子，代换式)12 进行数学式的转换调用Maple中数学式的转换命令，调用格式如下：m

24、aple(Maple的数学式转换命令) 即：maple(convert(表达式，form)将表达式转换成form的表示方式 maple(convert(表达式，form, x) 指定变量为x，将依赖于变量x的函数转换成form的表示方式（此指令仅对form为exp与sincos的转换式有用） 13 解方程solve(方程，变元) 注：方程的等号用普通的等号： = 14 解不等式调用maple中解不等式的命令即可，调用形式如下： maple(maple中解不等式的命令)具体说，包括以下五种：maple( solve（不等式）) maple( solve（不等式，变元） ) maple( solv

25、e（不等式，变元） ) maple( solve（不等式，变元） ) maple( solve（不等式，变元） )15 解不等式组调用maple中解不等式组的命令即可，调用形式如下： maple(maple中解不等式组的命令) 即：maple( solve（不等式组，变元组） )16 画图方法：先产生横坐标的取值和相应的纵坐标的取值，然后执行命令： plot(x,y) 方法2：fplot(f(x)，xmin,xmax) fplot(f(x)，xmin,xmax,ymin,ymax) 方法3：ezplot(f(x) ezplot(f(x) ,xmin,xmax) ezplot(f(x) ,xmi

26、n,xmax,ymin,ymax) 17 求极限（1）极限：syms x limit(f(x), x, a) （2）单侧极限：左极限：syms x limit(f(x), x, a,left)右极限：syms x limit(f(x), x, a,right) 18 求导数diff(f(x) diff(f(x),x) 或者：syms x diff(f(x) syms x diff(f(x), x) 19 求高阶导数 diff(f(x)，n) diff(f(x),x，n)或者：syms x diff(f(x)，n)syms x diff(f(x), x，n) 20 在MATLAB中没有直接求隐函

27、数导数的命令，但是我们可以根据数学中求隐函数导数的方法，在中一步一步地进行推导；也可以自己编一个求隐函数导数的小程序；不过，最简便的方法是调用Maple中求隐函数导数的命令，调用格式如下： maple(implicitdiff(f(x,y)=0,y,x) 在MATLAB中，没有直接求参数方程确定的函数的导数的命令，只能根据参数方程确定的函数的求导公式 一步一步地进行推导；或者，干脆自己编一个小程序，应用起来会更加方便。21 求不定积分 int(f(x) int (f(x),x)或者：syms x int(f(x) syms x int(f(x), x) 22 求定积分、广义积分 int(f(x

28、),a,b) int (f(x),x,a,b)或者：syms x int(f(x),a,b) syms x int(f(x), x,a,b) 23 进行换元积分的计算自身没有提供这一功能，但是可以调用Maple函数库中的changevar命令，调用方法如下：maple( with(student) ) 加载student函数库后，才能使用changevar命令maple( changevar( m(x)=p(u), Int(f(x),x) ) ) 把积分表达式中的m(x)代换成p(u)24 进行分部积分的计算自身没有提供这一功能，但是可以调用Maple函数库中的intparts命令，调用方法如

29、下： maple( with(student) ) 加载student函数库后，才能使用intparts命令maple(intparts(Int(f(x),x),u) ) 指定u，用分部积分公式 进行计算25 对数列和级数进行求和 syms n symsum(f(n), n a ,b )26 进行连乘 maple(product(f(n),n=a.b)27 展开级数syms x taylor(f(x), x, n, a )28 进行积分变换syms s t laplace( f(t), t, s ) 拉普拉斯变换 ilaplace( F(s), s, t ) 拉普拉斯变换的逆变换 syms t

30、 fourier( f(t), t, ) 傅立叶变换 ifourier( F(), , t ) 傅立叶变换的逆变换 syms n z ztrans( f(n), n, z) Z变换 iztrans( F(z), z, n ) Z变换的逆变换 在matlab中，矩形法、梯形法和辛普森法求近似积分可以用自身的命令，也可调用Maple的相应命令。调用方法如下： maple(with(student) ) maple(Maple中求定积分近似值的命令)29 解微分方程dsolve(微分方程，自变量) dsolve(微分方程，初始条件或边界条件，自变量)30 解微分方程组dsolve(微分方程组，自变量) dsolve(微分方程组，初始条件或边界条件，自变量)第 28 页 共 29 页