基于Matlab的AM振幅调制与解调仿真

1、基于Matlab的AM振幅调制与解调仿真摘要：本次高频电子电路大作业的设计，我组所选的题目为振幅调制电路（am及解调。在本课程设计报告中，首先说明了进行此次课程设计的目的、内容及要求；阐明了标准振幅调制与解调的基本原理以及操作方法，同时也对滤波电路的原理加以说明。接着叙述了利用Matlab软件对振幅调制、解调以及滤波器等所设计编写的程序，并附上了调试后输出的载波信号、调制信号、AM已调信号及滤波前后的解调信号等的波形图和频谱图，另外还附上了滤波器的增益响应和双边带总功率与平均总功率之比。报告的最后，是个人对本次大作业结果的分析、过程反思以及总结。关键词：振幅调制解调AMMatlab仿真Abst

2、ract:InThehigh-frequencyelectroniccircuitdesigningjob,ourgroupselectedthetopicasamplitudemodulationcircuit(AM)anddemodulation.Inthiscoursedesignreport,firstexplainsthepurpose,contentandrequirementsofthecurriculumdesign;clarifythebasicprinciplesandmethodsofoperationstandardamplitudemodulationanddemod

3、ulation,andalsotoillustratetheprinciplesofthefiltercircuit.ThendescribestheuseofMatlabandotheramplitudemodulation,demodulationandfilterdesignprogramwritten,alongwiththecarriersignaldebuggingoutputmodulationsignal,AMmodulatedanddemodulatedsignalwaveformsignalbeforeandafterfiltering,etc.mapandspectrum

4、,alsoattachedatotalpowerandaveragepowerratioofthetotalgainresponseandbilateralbandfilter.Attheendofthereportistoanalyzetheindividualresultsofthislargeoperation,processreflectionandsummary.Keywords:amplitudemodulation,demodulation,Matlabsimulation基于Matlab的AM振幅调制与解调仿真引言：无线通信系统中，信号通过一定的传输介质在发射机和接受机之间进行

5、传送时，信号的原始形式一般不适合传输。因此，为了将信号从发射端传输到接收端，必须转换它们的形式。而把信息加载到信息载体上以便传输的处理过程便是调制。通常称代表信息的信号为调制信号，称信息载体为载波信号，称调制后的频带信号为已调波信号。调制的种类很多，分类方法各不同，按调制信号的形式，可分为模拟调制和数字调制；模拟调制中的正弦波调制又分为振幅调制、频率调制和相位调制。其中的振幅调制是由调制信号去控制载波的振幅，使之按调制信号的规律变化，严格地讲，是使高频振荡的振幅与调制信号成线性关系，其他参数（频率和相位）不变。这是使高频振荡的振幅有消息的调制方式。振幅调制分为三种方式：普通的调幅方式（AM、抑

6、制载波的双边带调制（DSB及抑制载波的单边带调制（SSB方式。解调是与调制相反的过程，即从接收到的已调波信号中恢复原调制信息的过程。与振幅调制、频率调制和相位调制相对应，有振幅解调、频率解调和相位解调，并分别称为检波、鉴频和鉴相。本次课程设计所使用软件Matlab是1984年由美国Mathworks公司推向市场。它是一种科学计算软件，专门以矩阵的形式处理数据。Matlab将高性能的数值计算和可视化集成在一起，并提供了大量的内置函数，从而被广泛的应用于科学计算、控制系统和信息处理等领域的分析、仿真和设计工作。Matlab软件包括五大通用功能，数值计算功能（Nemeric）、符号运算功能（Symb

7、olic）数据可视化功能（Graphic）数字图形文字统一处理功能（Notebook）和建模仿真可视化功能（Simulink）。该软件有三大特点，一是功能强大；二是界面友善、语言自然；三是开放性强。Matlab的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故Matlab来解算问题要比用C,Fortran等语言完相同的事情简捷得多，并且mathwork也吸收了像Maple等软件的优点，使Matlab成为一个强大的数学软件。在新的版本中也加入了对C,Fortran,C+,Java的支持。可以直接调用，用户也可以将自己编写的实用程序导入到MATLA函数库中方便自己以后调用，此

8、外许多的Matlab爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。Matlab的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用Matlab函数集）扩展了Matlab环境，以解决这些应用领域内特定类型的问题。基于Matlab的AM振幅调制与解调仿真、作业题目1 .类型：设计性课程作业。2 .题目：基于Matlab的振幅调制（AM与解调。二、作业要求1 .掌握振幅调制和解调的原理。2 .学会Matlab仿真软件在振幅调制和解调中的应用。3 .掌握参数设置方法和性能分析方法。4 .通过实验中波

9、形的变换，学会分析实验现象。5 .建立较为完善的信号链路模型，能够较好的描述电路的工作过程;6 .正确分析输入输出信号的特征，关键步骤有相关图形输出。三、设计原理3.1振幅调制的一般模型振幅是用调制信号去控制高频正弦载波的幅度，使其按调制信号的规律变化的过程。幅度调制器的一般模型如图3.1所示。图3.1振幅调制的一般模型图中，脚。）为调制信号，为已调信号，为为滤波器的冲激响应，则已调信号的时域和频域一般表达式分别为：5.=地）cos/水航外（3-1）s架T露）+肠3-小）向口3-2）式中，”3）为调制信号切的频谱，见。岭/为载波角频率。由以上表达式可见，对于幅度调制信号，在波形上它的幅度随基带

10、信号规律而变化；在频谱结构上它的频谱完全是基带信号频谱在频域内的简单搬移。由于这种搬移是线性的，因此幅度调制通常又称为线性调制，相应地，幅度调制系统也称为线性调制系统。在图3.1的一般模型中，适当选择滤波器的特性丹（此，便可得到各种幅度调制信号，例如：常规双边带调幅（AM、抑制载波双边带调幅（DSB-SC、单边带调制（SSB和残留边带调制（VSB信号等。基于Matlab的AM振幅调制与解调仿真3.2 标准振幅调制(AM3.2.1 AM信号的表达式、频谱及带宽在图3.1中，若假设滤波器为全通网络(H=1),调制信号陶叠加直流4后再与载波相乘，则输出的信号就是常规双边带调幅(AM信号。AM调制器模

11、型如图3.2.1所示。小图3.2.1AM调制器模型wtOTA+祖Qc口汽AM信号的时域和频域表示式分别为：»一”；(3-3)加3)=飞回岸+M)+通-%)4卅3+/3-，。八2(3-4)式中，4为外加的直流分量；脚可以是确知信号也可以是随机信号，但通常认为其平均值为0,即丽-0。AM信号的典型波形和频谱分别如图3-3(a)、(b)所示，图中假定调制信号谣的上限频率为皿&。显然，调制信号脚©的带宽为勾=%。(a)基于Matlab的AM振幅调制与解调仿真由图3-3（a）可见，AM信号波形的包络与输入基带信号烟S成正比，故用包络检波的方法很容易恢复原始调制信号。但为了保证

12、包络检波时不发生失真，必须满足4'M否则将出现过调幅现象而带来失真。AM信号的频谱是由载频分量和上、下两个边带组成（通常称频谱中画斜线的部分为上边带，不画斜线的部分为下边带）。上边带的频谱与原调制信号的频谱结构相同，下边带是上边带的镜像。显然，无论是上边带还是下边带，都含有原调制信号的完整信息。故AM信号是带有载波的双边带信号，它的带宽为基带信号带宽的两倍，即：Ejna=2为=2为（3-5）式中，日僵为为调制信号脚S的带宽，i为调制信号的最高频率。3.2.2 AM信号的功率分配及调制效率AM信号在1口电阻上的平均功率应等于的均方值。当脚为确知信号时,$3的均方值即为其平方的时间平均，即

13、：为+跋co"/=封COS2G*+（J）COSaG/+24两。）。0GJ因为调制信号不含直流分量，即丽=口，且cqFg/=1/2,所以：(3-6)式中，尺二用/2为载波功率；尺为边带功率，它是调制信号功率外=标的一半。由此可见，常规双边带调幅信号的平均功率包括载波功率和边带功率两部分。只有边带功率分量与调制信号有关，载波功率分量不携带信息，定义调制效耳Pw率：(3-7)3.3 AM信号的解调调制过程的逆过程叫做解调。振幅波解调方法可分为包络检波和同步检波两大类。包络检波是指解调器输出电压与输入已调波的包络成正比的检波方法。由基于Matlab的AM振幅调制与解调仿真于AM信号的包络与调

14、制信号成线性关系，因此包络检波只适用于AM波。AM信号的解调是把接收到的已调信号还原为调制信号透。图3.3.1包络检波原理图下图是二极管峰值包络检波器的原理电路。它是由输入回路、二极管VD和RC低通滤波器组成。图3.3.2二极管峰值包络检波器（a）原理电路；（b）二极管导通；（c）二极管截止当RC两足条件：玛包络检波器的输出与输入信号的包络十分相近，即：叫任”4+双。（3-9）根据对电路的分析，可以知道：检波过程就是信号源通过二极管给电容充电与电容对电阻R放电的交替重复过程；由于RC时常数远大于输入电压载波周期，放电慢，使得二极管负极永远处于正的较高的电位（因为输出电压接近于高频正基于Matl

15、ab的AM振幅调制与解调仿真弦波的峰值，即UgUm）二极管电流iD包含平均分量（此种情况为直流分量）Iav及高频分量。Uo(t”0t图3.3.3输入为AM信号时检波器的输出波形图包络检波法特点是：解调效率高；解调电路简单，特别是接收端不需要与发送端同频同相位的载波信号，大大降低实现难度。故几乎所有的调幅（AM式接收机都采用这种电路。3.4 巴特沃斯低通滤波器巴特沃斯滤波器的特点是同频带内的频率响应曲线最为平坦，没有起伏，而在组频带则逐渐下降为零。在振幅的对数对角频率的波特图上，从某一边界见频率开始，振幅随着角频率的增加而逐渐减少，趋向于负无穷大。一阶巴特沃斯滤波器的衰减率为每倍频20分贝，二阶

16、巴特沃斯滤波器的衰减率为每倍频12分贝，三阶的衰减率为每分贝18分贝，如此类推，巴特沃斯滤波器的振幅对角频率单调下降，并且滤波器的结束越高，在组频带振幅衰减速度越快，其他滤波器高阶的振幅对角频率图和低阶数的振幅对角频率有不同的形状。1H(s)H(-s)=21(S)NJ2)-c上述函数的特点是等距离分布在半径为C的圆上因此，极点用下式表示为基于Matlab的AM振幅调制与解调仿真jnj(2k1)Sk=jce»eNk:0,1,2,N-1Ha(s)的表小式:H小"Ha(s)-NjI1(s-Sk)k：0为了使设计公式和图表统一，将频率归一化。巴特沃斯滤波器采用3dB截1sSk-)c

17、止频率Cc归一化，归一化后的系统函数为:Ga(-二NT",cIT(kz0令p="+j7u=sCc&=Q/Cc,人称为归一化频率，p称为归一化复变量,这样巴特沃斯滤波器的归一化低通原型系统函数为：G二1aN1i【(p-Pk)k=0式中，Pk=s/Gc,为归一化极点，用下式表示Pk=e._.J2k1、j.()22Nk=0,1,2,N-1四、MATLAB4.1 载波信号4.1.1 仿真程序function=Zaiboxinhao()U1=5;f1=3000;t=-1:0.00001:1;%载波幅值为5%载波频率为3000%t扫描范围为-1到1%载波信号角频率%载波信号表达

18、式%新建一个图形窗口1%将图形窗口一分为二，并把第一个作为当前图形窗口%绘制载波信号波形%横坐标为t,纵坐标显示为u1%标题为载波信号波形%设置显示范围%对u1进行傅里叶变换%将第二个子图形窗口作为当前图形窗口%绘制Y1的图形%横坐标为w,纵坐标显示为Y1%标题为载波信号频谱%设置显示范围基于Matlab的AM振幅调制与解调仿真w1=2*pi*f1;u1=U1*cos(w1*t);figure(1);subplot(2,1,1);plot(t,u1);xlabel('t');ylabel('u1');title('载波信号波形');axis(0,

19、0.01,-10,10);Y1=fft(u1);subplot(2,1,2);plot(abs(Y1);xlabel('w');ylabel('Y1');title('载波信号频谱');axis(5800,6200,0,600000);4.1.2仿真波形基于Matlab的AM振幅调制与解调仿真图4.1.2载波信号波形及频谱figure(2);subplot(2,1,1);plot(t,u2);xlabel('t');ylabel('u2');title('调制信号波形');4.2调制信号4.2.1仿

20、真程序function=Tiaozhixinhao()U2=3;f2=3;t=-1:0.00001:1;%t扫描范围为-1到1w2=2*pi*f2;u2=1.2*U2*cos(w2*t)+U2*cos(2*w2*t)+1.8*U2*cos(3*w2*t);%调制信号表达式%新建一个图形窗口2%将图形窗口一分为二，并把第一个作为当前图形窗口%绘制载波信号波形%横坐标为t,纵坐标显示为u2%标题为调制信号波形%对u2进行傅里叶变换%将第二个子图形窗口作为当前图形窗口%绘制Y2的图形%横坐标为w,纵坐标显示为u1%标题为载波信号频谱%设置显示范围基于Matlab的AM振幅调制与解调仿真axis(0,

21、1,-15,15);%设置显示范围Y2=fft(u2);subplot(2,1,2);plot(abs(Y2);xlabel('w');ylabel('Y2');title('调制信号频谱');axis(0,250,0,2000000);4.2.2仿真波形图4.2.2调制信号波形及频谱4.3 AM调制4.3.1 仿真程序%t扫描范围-1到1function=Tiaozhi()t=-1:0.00001:1;%调制信号%AM已调信号%新建一个图形窗口3%将图形窗口一分为二，并把第一个作为当前图形窗口%绘制已调信号波形%横坐标为t,纵坐标显示为u3%标

22、题为巴调信号波形%设置显示范围%对u3进行傅里叶变换%将第二个子图形窗口作为当前图形窗口%绘制Y3的图形%横坐标为t,纵坐标显示为u3%标题为巴调信号频谱%设置显示范围基于Matlab的AM振幅调制与解调仿真U1=5;%载波信号幅度U2=3;f1=3000;%载波信号频率f2=3;m=0.1;%调制度为0.1w1=2*pi*f1;%载波信号角频率w2=2*pi*u2=1.2*U2*cos(w2*t)+U2*cos(2*w2*t)+1.8*U2*cos(3*w2*t);u3=U1*(1+m*u2).*cos(w1)*t);figure(3);subplot(2,1,1);plot(t,u3);x

23、label('t');ylabel('u3');title('已调信号波形');axis(0,1,-15,15);Y3=fft(u3);subplot(2,1,2);plot(abs(Y3);xlabel('w');ylabel('Y3');title('已调信号频谱');axis(5900,6100,0,600000);4.3.2 仿真波形基于Matlab的AM振幅调制与解调仿真图4.3.2已调信号波形及频谱4.4 AM波解调（包络检波法）4.4.1 仿真程序function=Jietiao1()

24、%经过AM调制产生已调信号u3env=abs(hilbert(u3);u4=18*(env-U1)*m;figure(4);subplot(2,1,1);plot(t,u4);xlabel('t');ylabel('u4');title('AM已调信号的包络检波波形')axis(0,1,-15,15);%找出已调信号的包络%去掉直流分量并重新缩放%新建一个图形窗口4%将图形窗口一分为二，并把第一个作为当前图形窗口%绘制解调波形%横坐标为t,纵坐标显示为u4%标题为AME调信号的包络检波波形，%设置显示范围丫4=fft(u4);%对u4进行傅里叶变

25、换基于Matlab的AM振幅调制与解调仿真subplot(2,1,2);%将第二个子图形窗口作为当前图形窗口plot(abs(Y4);xlabel('w');ylabel('Y4');title('AM已调信号的包络检波频谱');%标题为AM已调信号的包络检波频谱axis(0,250,0,2000000);%设置显示范围4.4.2 仿真波形图4.4.2AM已调信号的包络检波波形及频谱4.5 AM波解调(同步乘积型检波法)4.5.1 仿真程序function=Jietiao2()%经过AM调制产生已调信号u3基于Matlab的AM振幅调制与解调仿真

26、u5=u3.*cos(w1*t);figure(5);subplot(2,1,1);plot(t,u5);xlabel('t');ylabel('u5');title('滤波前的解调信号波形');Y5=fft(u5);subplot(2,1,2);plot(abs(Y5);xlabel('w');ylabel('Y5');title('滤波前的解调信号频谱');axis(187900,188100,0,600000);%滤波前的解调信号%新建图形窗口5%将图形窗口一分为二，并把第一个作为当前图形窗口

27、%绘制滤波前的解调信号波形%横坐标为t,纵坐标显示为u5%设置标题%对u5进行傅里叶变换%将第二个子图形窗口作为当前图形窗口%绘制Y5的波形f1=100;f2=200;fs=2000;m=(0.3*f1)/(fs/2);M=round(8/m);N=M-1;b=fir1(N,0.5叶2/(fs/2);figure(6)h,f=freqz(b,1,512);%低通滤波器%待滤波信号频率%采样频率%定义过度带宽%定义窗函数的长度%定义滤波器的阶数%使用fir1函数设计滤波器%输入的参数分别是滤波器的阶数和截止频率%新建图形窗口6%滤波器的幅频特性图%H,W=freqz(B,A,N)当N是一个整数时

28、函数返回N点的频率向量和幅频基于Matlab的AM振幅调制与解调仿真响应向量plot(f*fs/(2*pi),20*log10(abs(h)%参数分别是频率与幅值xlabel('频率/赫兹');ylabel('增益/分贝');title('滤波器的增益响应');gridu6=filter(b,1,u5);figure(7);subplot(2,1,1);plot(t,u6);xlabel('t');ylabel('u6');title('滤波后的解调信号波形');Y6=fft(u6);subplot

29、(2,1,2);plot(abs(Y6);xlabel('w');ylabel('Y6');title('滤波后的解调信号频谱');axis(0,250,0,600000);%滤波后的解调信号%新建图形窗口7%将图形窗口一分为二，并把第一个作为当前图形窗口%绘制滤波后的解调信号波形%对u6进行傅里叶变换%将第二个子图形窗口作为当前图形窗口%绘制Y6的波形4.5.2仿真波形基于Matlab的AM振幅调制与解调仿真滤淞前的解弱信号波形x图滤波前的解调信号波形及频谱'滤波器的噌益响应1002003004005'006007

30、0060D9001000频率,赫兹0000-0000022468024-.-T-111图滤波器的增益响应基于Matlab的AM振幅调制与解调仿真图滤波后的解调信号波形及频谱4.6 AM波的功率4.6.1 仿真程序function=Gonglv()m=0:0.01:1;Ucm=5;RL=1000;Pc=1/2*Ucm*Ucm/RL;%调制度扫描范围%载波信号幅值%负载电阻%负载上消耗的载波功率%下边频分量所消耗的平均功率%在调制信号的一个周期内，调幅信号的平均总功率%双边带总功率与平均总功率之比%新建一个图形窗口8%做出m与e的关系曲线%设定横纵坐标显示Pu=(1/2

31、*m*Ucm).*(1/2*m*Ucm)/(2*RL);%上边频分量所消耗的平均功率Pl=Pu;PAM=Pc+Pu+Pl;e=(Pu+Pl)./PAM;figure(8);plot(m,e);xalbel(调制度m')；基于Matlab的AM振幅调制与解调仿真yxabel（双边带总功率与平均总功率之比）;gridon4.6.2 仿真波形1lLh;/hX1/*/1+1'h*/1«1/9tJrFn11-.*-1才FIIp1!iIfiI000102030405OG0,7。8。,9调制度m图4.6.2双边带总功率与平均总功率之比4.7 调制度m对AM调制的影响4.7.1 仿真

32、程序%t扫描范围-1到1%载波信号幅度%调制信号幅度%载波信号频率%调制信号频率%调制度为0function=m_yingxiang()t=-1:0.00001:1;U1=4;U2=2;f1=3000;f2=3;m1=0;基于Matlab的AM振幅调制与解调仿真m2=0.4;m3=0.7;m4=1;m5=1.3;m6=3;%调制度为0.3%调制度为0.6%调制度为1%调制度为1.3%调制度为0的调制信号%调制度为0.4的调制信号%调制度为0.7的调制信号%调制度为1的调制信号%调制度为1.3的调制信号%调制度为3的调制信号%调制度为3u1=U1*(1+m1.*cos(2*pi*f2*t).*c

33、os(2*pi*f1*t);u2=U1*(1+m2.*cos(2*pi*f2*t).*cos(2*pi*f1*t);u3=U1*(1+m3.*cos(2*pi*f2*t).*cos(2*pi*f1*t);u4=U1*(1+m4.*cos(2*pi*f2*t).*cos(2*pi*f1*t);u5=U1*(1+m5.*cos(2*pi*f2*t).*cos(2*pi*f1*t);u6=U1*(1+m6.*cos(2*pi*f2*t).*cos(2*pi*f1*t);figure;%新建一个图形窗口subplot(3,2,1)%将图形窗口1分为6,并把第一个作为当前图形窗口plot(t,u1);%

34、绘制调制度为0时的波形xlabel('t');ylabel('u1');title('m=0时AM调制信号');axis(0,1,-10,10);gridon;subplot(3,2,2)%将第2个子图形窗口作为当前图形窗口plot(t,u2);%绘制调制度为0.4时的波形xlabel('t');ylabel('u2');title('m=0.4时AM调制信号');gridon;subplot(3,2,3)%将第3个子图形窗口作为当前图形窗口基于Matlab的AM振幅调制与解调仿真plot(t,u3

35、);xlabel('t');ylabel('u3');title('m=0.7时AM调制信号');gridon;subplot(3,2,4)plot(t,u4);xlabel('t');ylabel('u4');title('m=1时AM调制信号');gridon;subplot(3,2,5)plot(t,u5);xlabel('t');ylabel('u5');title('m=1.3时AM调制信号');gridon;subplot(3,2,6)plot(t,u6);xlabel('t');ylabel('u6');title('m=3时AM调制信号');gridon;%绘制调制度为0.7时的波形%将第4个子图形窗口作为当前图形窗口%绘制调制度为1时的波形%将第5个子图形窗口作为当前图形窗口%绘制调制度为1.3时的波形%将第6个子图形窗口作为当前图形窗口%绘制调制度为3时的波形4.7.2仿真波形基于Matlab的AM振幅调制与解调仿真ttm=07时AM调制信号1-0.500.51tm=1