通信仿真实验报告（高频）(2)

1）设计AM振幅调制与解调制仿真电路，要求调制信号的幅度A和频率F可变；载波信号的幅度A和频率F可变，调制度ma可变；

2）绘制调制信号u11、载波信号uc1和已调波信号uam1的时域波形图和频谱图（要求谱线清晰）；

3）要求调制信号为三个以上正弦波信号的合成，幅度和频率均可变，绘制调制信号u12、载波信号uc1和已调波信号uam2的时域波形图和频谱图（要求谱线清晰）； 4）用同步检波对已调波信号uam1进行解调制，在同一幅图中绘制原调制信号u11和解调后的信号y11，并绘制这两个信号的频谱图；

5）用同步检波对已调波信号uam2进行解调制，在同一幅图中绘制原调制信号u12和解调后的信号y12，并绘制这两个信号的频谱图。

1.单调制信号的AM波

function [u11,uc1,uam1] =dantiaozhixinhaoAMbo (a1,a2,f1,f2,ma) ? 调制信号的幅度; ? 调制信号的频率; ￠ 载波信号的幅度; ò 载波信号的频率; %ma 调制度

fs=20*max(f1,f2);

t=0:1/fs:2/min(f1,f2);

u11=a1*cos(2*pi*f1*t); %设调制信号频率 uc1=a2*cos(2*pi*f2*t); %设载波频率 uam1=(1+ma*u11).*uc1; %AM振幅调制 X=fftshift(fft(u11)); Y=fftshift(fft(uc1)); Z=fftshift(fft(uam1));

F=linspace(-fs/2,fs/2,length(t)); subplot(3,1,1);plot(t,u11);grid on subplot(3,1,2);plot(t,uc1);grid on subplot(3,1,3);plot(t,uam1);grid on figure;

subplot(3,1,1);plot(F,abs(X));grid on subplot(3,1,2);plot(F,abs(Y));grid on subplot(3,1,3);plot(F,abs(Z));grid on s=uam1.*uc1;

Rp=3;%信号衰减幅度 Rs=60;%信号衰减幅度 Wp=40/500;%通带截止频率 Ws=150/500;%阻带截止频率

[n,Wn]=ellipord(Wp,Ws,Rp,Rs);%阶数n

[b,a]=ellip(n,Rp,Rs,Wn);%传递函数分子分母b，a y11=50*filter(b,a,s); A=fftshift(fft(u11)); B=fftshift(fft(y11));

F=linspace(-fs/2,fs/2,length(t)); figure;

subplot(4,1,1);plot(t,u11);grid on subplot(4,1,2);plot(t,y11);grid on

subplot(4,1,3);plot(F(1800:2250),abs(A(1800:2250)));grid on subplot(4,1,4);plot(F(1800:2250),abs(B(1800:2250)));grid on 在命令窗口出入[u11,uc1,uam1] =dantiaozhixinhaoAMbo (1,3,1000,100000,0.3) 得到：

单个调制信号的调制波形

频谱图

检波波形

2.三个输入的AM波

Function [u11,uc1,uam1]

=duogeshurudeAMbo(a11,a12,a13,a2,f11,f12,f13,f2,ma) ?1 调制信号1的幅度; ?2 调制信号2的幅度; ?3 调制信号3的幅度; ?1 调制信号1的频率; ?2 调制信号2的频率; ?3 调制信号3的频率; ￠ 载波信号的幅度; ò 载波信号的频率; %ma 调制度 fs=20*f2;

t=0:1/fs:2/f11;

u11=a11*cos(2*pi*f11*t)+ a12*cos(2*pi*f12*t)+ a13*cos(2*pi*f13*t); %设调制信号频率 uc1=a2*cos(2*pi*f2*t); %设载波频率 uam1=(1+ma*u11).*uc1; %AM振幅调制 X=fftshift(fft(u11)); Y=fftshift(fft(uc1)); Z=fftshift(fft(uam1));

F=linspace(-fs/2,fs/2,length(t)); subplot(3,1,1);plot(t,u11);grid on subplot(3,1,2);plot(t,uc1);grid on subplot(3,1,3);plot(t,uam1);grid on figure;

subplot(3,1,1);plot(F,abs(X));grid on subplot(3,1,2);plot(F,abs(Y));grid on subplot(3,1,3);plot(F,abs(Z));grid on s=uam1.*uc1;

Rp=3;%信号衰减幅度 Rs=60;%信号衰减幅度 Wp=40/500;%通带截止频率 Ws=150/500;%阻带截止频率

[n,Wn]=ellipord(Wp,Ws,Rp,Rs);%阶数n

[b,a]=ellip(n,Rp,Rs,Wn);%传递函数分子分母b，a y11=5*filter(b,a,s); A=fftshift(fft(u11)); B=fftshift(fft(y11));

F=linspace(-fs/2,fs/2,length(t)); figure;

subplot(4,1,1);plot(t,u11);grid on subplot(4,1,2);plot(t,y11);grid on subplot(4,1,3);plot(F,abs(A));grid on subplot(4,1,4);plot(F,abs(B));grid on

在命令窗口出入：[u11,uc1,uam1] =duogeshurudeAMbo(1, 2, 3,1,300,500,1000,10000,0.3) 得到：

三个调制信号的调制波形

频谱图