基于MATLAB的AM调制及解调系统仿真Word下载.docx

1、式中，ka为比例系数，一般由调制电路确定，故又称为调制灵敏度。由此可得调幅信号的表达式 （2-5） 上面的分析是在单一正弦信号作为调制信号的情况下进行的，而一般传送的信号并非为单一频率的信号，例如是一连续频谱信号f（t）,这时，可用下式来描述调幅波: （2-6）式中，f（t）是均值为零的归一化调制信号，|f（t）|max=1。若将调制信号分解为： （2-7）则调幅波表示式为： （2-8）2.1.2 调幅波的频谱 在单一频率的正弦信号的调制情况下，调幅波如前所描述。将其用三角公式展开，可得: （2-9）可见,单一频率信号调制的调幅波包含三个频率分量, 由三个高频正弦波叠加而成, 调制信号的幅度及

2、频率信息只含在边频分量中。图1 单音调制时已调波的频谱（a）调制信号频谱;（b）载波信号频谱;（c）AM信号频谱2.1.3 调幅波的功率AM调幅信号： （2-10）设上式代表的调幅信号传输至负载电阻RL上，那么调幅波各频率分量在RL上消耗的功率分别为：（1）RL上消耗的载波功率： （2-11）（2） 上、下边频分量所消耗的平均功率： （2-12）（3）在调制信号的一个周期内，调幅信号的平均总功率： （2-13）由此可得边频功率、载波功率与平均总功率之间的关系比为： （2-14） （2-15）2.2调幅信号的解调从高频已调信号中恢复出调制信号的过程称为解调，又称为检波。对于振幅调制信号，解调就是

3、从它的幅度变化上提取调制信号的过程。解调是调制的逆过程。振幅波解调方法可分为包络检波和同步检波两大类。包络检波是指解调器输出电压与输入已调波的包络成正比的检波方法。由于AM信号的包络与调制信号成线性关系，因此包络检波只适用于AM波。同步检波可分为乘积型同步检波和叠加型。对于乘积型同步检波，首先使调幅波与本地恢复信号相乘，最后经低通滤波器输出。当恢复载波与发射载波同频同相时，同步检波器可以无失真地将调制信号恢复出来。 叠加型同步检波是将DSB或SSB信号插入恢复载波，使之成为或近似为AM信号，再利用包络检波器将调制信号恢复出来。2.3巴特沃斯滤波器巴特沃斯滤波器的特点是通频带内的频率响应曲线最大

4、限度平坦，没有起伏，而在阻频带则逐渐下降为零。在振幅的对数对角频率的波特图上,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无穷大。一阶巴特沃斯滤波器的衰减率为每倍频6分贝，每十倍频20分贝。二阶巴特沃斯滤波器的衰减率为每倍频12分贝、三阶巴特沃斯滤波器的衰减率为每倍频18分贝、如此类推。巴特沃斯滤波器的振幅对角频率单调下降，并且也是唯一的无论阶数，振幅对角频率曲线都保持同样的形状的滤波器。只不过滤波器阶数越高，在阻频带振幅衰减速度越快。其他滤波器高阶的振幅对角频率图和低级数的振幅对角频率有不同的形状。巴特沃斯低通滤波器可用如下振幅的平方对频率的公式表示： （2-16）其中，n=

5、滤波器的阶数，=截止频率=振幅下降为-3dB时的频率，=通频带边缘频率3. MATLAB仿真3.1 载波信号3.1.1 仿真程序function = Zaiboxinhao（）U1=5; %载波幅值为5f1=3000; %载波频率为3000t=-1:0.00001:1; % t扫描范围为-1到1w1=2*pi*f1; %载波信号角频率u1=U1*cos（w1*t）; %载波信号表达式figure（1）; %新建一个图形窗口1subplot（2,1,1）; %将图形窗口一分为二，并把第一个作为当 前图形窗口 plot（t,u1）; %绘制载波信号波形xlabel（t）;ylabel（u1 %横坐

6、标为t，纵坐标显示为u1title（载波信号波形 %标题为载波信号波形axis（0,0.01,-10,10）; %设置显示范围Y1=fft（u1）; %对u1进行傅里叶变换subplot（2,1,2）; %将第二个子图形窗口作为当前图形窗口plot（abs（Y1）; %绘制Y1的图形wY1 %横坐标为w，纵坐标显示为Y1载波信号频谱 %标题为载波信号频谱axis（5800,6200,0,600000）; %设置显示范围3.1.2仿真波形图2 载波信号波形及频谱3.2调制信号3.2.1 仿真程序function =Tiaozhixinhao（）U2=3;f2=3; %t扫描范围为-1到1w2=2

7、*pi*f2;u2=1.2*U2*cos（w2*t）+U2*cos（2*w2*t）+1.8*U2*cos（3*w2*t）; %调制信号表达式figure（2）; %新建一个图形窗口2 %将图形窗口一分为二，并把第一个作为当前图形窗口 plot（t,u2）;u2 %横坐标为t，纵坐标显示为u2调制信号波形 %标题为调制信号波形axis（0,1,-15,15）;Y2=fft（u2）; %对u2进行傅里叶变换plot（abs（Y2）; %绘制Y2的图形Y2 %横坐标为w，纵坐标显示为u1调制信号频谱 %标题为载波信号频谱axis（0,250,0,2000000）;3.2.2仿真波形图3调制信号波形及

8、频谱3.3 AM调制3.3.1 仿真程序function =Tiaozhi（） %t扫描范围-1到1 %载波信号幅度 %载波信号频率m=0.1; %调制度为0.1 %调制信号u3=U1*（1+m*u2）.*cos（w1）*t）; %AM已调信号figure（3）; %新建一个图形窗口3 %将图形窗口一分为二，并把第一个作为当plot（t,u3）; %绘制已调信号波形u3 %横坐标为t，纵坐标显示为u3已调信号波形 %标题为已调信号波形Y3=fft（u3）; %对u3进行傅里叶变换 %将第二个子图形窗口作为当前图形窗口plot（abs（Y3）; %绘制Y3的图形Y3已调信号频谱 %标题为已调信号

9、频谱axis（5900,6100,0,600000）;3.3.2仿真波形图4 已调信号波形及频谱3.4 AM波解调（包络检波法）3.4.1 仿真程序 function =Jietiao1（）%经过AM调制产生已调信号u3env=abs（hilbert（u3）; %找出已调信号的包络 u4=18*（env-U1）*m; %去掉直流分量并重新缩放figure（4）; %新建一个图形窗口4 %将图形窗口一分为二，并把第一个作为当前图形窗口 plot（t,u4）; %绘制解调波形u4 %横坐标为t，纵坐标显示为u4AM已调信号的包络检波波形） %标题为AM已调信号的包络检波波形Y4=fft（u4）;

10、%对u4进行傅里叶变换plot（abs（Y4）;Y4AM已调信号的包络检波频谱 %标题为AM已调信号的包络检波频谱3.4.2仿真波形图5 AM已调信号的包络检波波形及频谱3.5 AM波解调（同步乘积型检波法）3.5.1 仿真程序function =Jietiao2（）u5=u3.*cos（w1*t）; %滤波前的解调信号 figure（5）; %新建图形窗口5 %将图形窗口一分为二 plot（t,u5）; %绘制滤波前的解调信号波形u5 %横坐标为t，纵坐标显示为u5滤波前的解调信号波形 %设置标题Y5=fft（u5）; %对u5进行傅里叶变换plot（abs（Y5）; %绘制Y5的波形 Y5

11、滤波前的解调信号频谱axis（187900,188100,0,600000）; %巴特沃斯低通滤波器f1=100;f2=200; %待滤波信号频率fs=2000; %采样频率m=（0.3*f1）/（fs/2）; %定义过度带宽M=round（8/m）; %定义窗函数的长度N=M-1; %定义滤波器的阶数b=fir1（N,0.5*f2/（fs/2）; %使用fir1函数设计滤波器 %输入的参数分别是滤波器的阶数和截止频率figure（6） %新建图形窗口6h,f=freqz（b,1,512）; %滤波器的幅频特性图 %H,W=freqz（B,A,N）当N是一个整数时函数返回N点的频率向量和幅频响

12、应向量plot（f*fs/（2*pi）,20*log10（abs（h） %参数分别是频率与幅值频率/赫兹增益/分贝滤波器的增益响应gridu6=filter（b,1,u5）; %滤波后的解调信号 figure（7）; %新建图形窗口7 %将图形窗口一分为二，并把第一个作为当前 图形窗口plot（t,u6）; %绘制滤波后的解调信号波形u6滤波后的解调信号波形Y6=fft（u6）; %对u6进行傅里叶变换plot（abs（Y6）; %绘制Y6的波形 Y6滤波后的解调信号频谱axis（0,250,0,600000）;3.5.2仿真波形图6 滤波器的增益响应图7 滤波前的解调信号波形及频谱图8滤波后

13、的解调信号波形及频谱3.6 AM波的功率3.6.1 仿真程序function =Gonglv（）m=-1:0.01: %调制度扫描范围Ucm=5; %载波信号幅值 RL=1000; %负载电阻Pc=1/2*Ucm*Ucm/RL; %负载上消耗的载波功率 Pu=（1/2*m*Ucm）.*（1/2*m*Ucm）/（2*RL）; %上边频分量所消耗的平均功率Pl=Pu; %下边频分量所消耗的平均功率PAM=Pc+Pu+Pl; %在调制信号的一个周期内，调幅信号的平均总功率e=（Pu+Pl）./PAM; %双边带总功率与平均总功率之比 figure（8）; %新建一个图形窗口8plot（m,e）; %

14、做出m与e的关系曲线xalbel（调制度m）； % 设定横纵坐标显示yxabel（双边带总功率与平均总功率之比）； grid on 3.6.2仿真波形图9 双边带功率与总功率之比与调制度的关系曲线3.7 调制度m对AM调制的影响3.7.1 仿真程序function =m_yingxiang（） %t扫描范围-1到1U1=4; %载波信号幅度U2=2; %调制信号幅度 %调制信号频率m1=0; %调制度为0m2=0.4; %调制度为0.3m3=0.7; %调制度为0.6m4=1; %调制度为1m5=1.3; %调制度为1.3m6=3; %调制度为3u1=U1*（1+m1.*cos（2*pi*f2

15、*t）.*cos（2*pi*f1*t）; %调制度为0的调制信号u2=U1*（1+m2.*cos（2*pi*f2*t）.*cos（2*pi*f1*t）; %调制度为0.4的调制信号u3=U1*（1+m3.*cos（2*pi*f2*t）.*cos（2*pi*f1*t）; %调制度为0.7的调制信号u4=U1*（1+m4.*cos（2*pi*f2*t）.*cos（2*pi*f1*t）; %调制度为1的调制信号u5=U1*（1+m5.*cos（2*pi*f2*t）.*cos（2*pi*f1*t）; %调制度为1.3的调制信号u6=U1*（1+m6.*cos（2*pi*f2*t）.*cos（2*pi*

16、f1*t）; %调制度为3的调制信号figure; %新建一个图形窗口subplot（3,2,1） %将图形窗口1分为6，把第一个作为当前图形窗口 %绘制调制度为0时的波形m=0时AM调制信号axis（0,1,-10,10）;grid on;subplot（3,2,2） %将第2个子图形窗口作为当前图形窗口 %绘制调制度为0.4时的波形m=0.4时AM调制信号subplot（3,2,3） %将第3个子图形窗口作为当前图形窗口 %绘制调制度为0.7时的波形m=0.7时AM调制信号subplot（3,2,4） %将第4个子图形窗口作为当前图形窗口 %绘制调制度为1时的波形m=1时AM调制信号sub

17、plot（3,2,5） %将第5个子图形窗口作为当前图形窗口 %绘制调制度为1.3时的波形m=1.3时AM调制信号subplot（3,2,6） %将第6个子图形窗口作为当前图形窗口 %绘制调制度为3时的波形m=3时AM调制信号3.7.2仿真波形图10 调制系数对调制过程的影响3.8仿真结果分析（1）由调制解调仿真结果可知，原调制信号的频率分量在0-50之间，调制过后，全在6000左右，解调过后，又恢复到0-50之间。AM调制解调的过程实际上是完成信号频谱的线性搬移。调制时，将调制信号的频谱搬移到载波信号的频段上；解调时，把低频的调制信号从高频的已调波上恢复出来。（2）调幅波幅度随调制信号的大小

18、线性变化。（3）调制系数m对AM波的调制有着直接的影响。当m=0时，无调幅作用；当m的范围为0-1时，随着调制度的增大，信号调制的深度越大；当m=1时，调制波形达到临界状态；当m大于1时，已调波的包络形状与调制信号不一样，发生了过量调幅，产生了严重的包络失真，且m越大，失真越厉害。（4）巴特沃斯低通滤波器的上限截止频率为100左右。（5）在利用包络检波解调时，直接利用了Matlab里自带的hilbert（）函数，没有产生各类失真情况；在利用乘积型同步检波解调时，与本地恢复波相乘之后，通过了巴特沃斯低通滤波器，也很好地还原出了调制信号。（6）由功率比与调制系数的曲线可得：m在0-1范围内，随调制

19、系数的减小，双边带功率与调幅波总功率的比值越小，当100%调制（即m=1）时双边带功率只有调幅波总功率的1/3，由于有用信息只携带在边频带内，而载波本身并不携带信息，但它的功率却占用了整个调幅波功率的大部分，因此AM调幅波的功率浪费大，效率低。4.总结通过本次大作业，让我再一次体会到了Matlab软件强大的功能。在整个过程中，通过查书和上网查阅资料，我温习了Matlab编程的基本语法格式及一些基本函数的用法，学会了建立M文件和实现函数编辑与调用，以及知道怎么去建立完整的信号仿真模型。同时，我对AM调制与解调的过程有了一个更深的认识。信号的调制与解调实际上就是完成信号频谱的线性搬移，以便于信号的传送。同时我了解了