基于matlab的am调制与仿真,基于MATLAB的AM调制及解调系统仿真分解

图7 滤波前的解调信号波形及频谱

图8滤波后的解调信号波形及频谱

3.6 AM波的功率

3.6.1 仿真程序

function [] =Gonglv()

m=-1:0.01:1; %调制度扫描范围 Ucm=5; %载波信号幅值 RL=1000; %负载电阻

13

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); %新建一个图形窗口8 plot(m,e); %做出m与e的关系曲线 xalbel(‘调制度m’)； % 设定横纵坐标显示

yxabel(‘双边带总功率与平均总功率之比’)； grid on

3.6.2仿真波形

图9 双边带功率与总功率之比与调制度的关系曲线

3.7 调制度m对AM调制的影响

3.7.1 仿真程序

function [ ] =m_yingxiang()

t=-1:0.00001:1; %t扫描范围-1到1

14

U1=4; %载波信号幅度 U2=2; %调制信号幅度 f1=3000; %载波信号频率 f2=3; %调制信号频率

m1=0; %调制度为0 m2=0.4; %调制度为0.3 m3=0.7; %调制度为0.6 m4=1; %调制度为1 m5=1.3; %调制度为1.3 m6=3; %调制度为3

u1=U1*(1+m1.*cos(2*pi*f2*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*f1*t); %调制度为3的调制信号

figure; %新建一个图形窗口

subplot(3,2,1) %将图形窗口1分为6，把第一个作为当前图形窗口 plot(t,u1); %绘制调制度为0时的波形 xlabel('t');ylabel('u1'); title('m=0时AM调制信号'); axis([0,1,-10,10]); grid on;

subplot(3,2,2) %将第2个子图形窗口作为当前图形窗口 plot(t,u2); %绘制调制度为0.4时的波形 xlabel('t');ylabel('u2');

15

title('m=0.4时AM调制信号'); grid on;

subplot(3,2,3) %将第3个子图形窗口作为当前图形窗口 plot(t,u3); %绘制调制度为0.7时的波形 xlabel('t');ylabel('u3'); title('m=0.7时AM调制信号'); grid on;

subplot(3,2,4) plot(t,u4); xlabel('t');ylabel('u4'); title('m=1时AM调制信号'); grid on;

subplot(3,2,5) plot(t,u5); xlabel('t');ylabel('u5');

title('m=1.3时AM调制信号'); grid on;

subplot(3,2,6) plot(t,u6); xlabel('t');ylabel('u6'); title('m=3时AM调制信号'); grid on;

3.7.2仿真波形

%将第4个子图形窗口作为当前图形窗口 %绘制调制度为1时的波形 %将第5个子图形窗口作为当前图形窗口 %绘制调制度为1.3时的波形 %将第6个子图形窗口作为当前图形窗口 %绘制调制度为3时的波形 16