前图形窗口
plot(t,u3); %绘制已调信号波形
xlabel('t');ylabel('u3'); %横坐标为t，纵坐标显示为u3 title('已调信号波形'); %标题为‘已调信号波形’ axis([0,1,-15,15]); %设置显示范围
Y3=fft(u3); %对u3进行傅里叶变换
subplot(2,1,2); %将第二个子图形窗口作为当前图形窗口 plot(abs(Y3)); %绘制Y3的图形
xlabel('w');ylabel('Y3'); %横坐标为t，纵坐标显示为u3 title('已调信号频谱'); %标题为‘已调信号频谱’ axis([5900,6100,0,600000]); %设置显示范围
3.3.2仿真波形
图4 已调信号波形及频谱
3.4 AM波解调（包络检波法）
3.4.1 仿真程序
function [] =Jietiao1()
9
%经过AM调制产生已调信号u3
env=abs(hilbert(u3)); %找出已调信号的包络 u4=18*(env-U1)*m; %去掉直流分量并重新缩放 figure(4); %新建一个图形窗口4
subplot(2,1,1); %将图形窗口一分为二，并把第一个作为当前图
形窗口
plot(t,u4); %绘制解调波形
xlabel('t');ylabel('u4'); %横坐标为t，纵坐标显示为u4
title('AM已调信号的包络检波波形') %标题为‘AM已调信号的包络检波波形’ axis([0,1,-15,15]); %设置显示范围
Y4=fft(u4); %对u4进行傅里叶变换
subplot(2,1,2); %将第二个子图形窗口作为当前图形窗口 plot(abs(Y4)); xlabel('w');ylabel('Y4');
title('AM已调信号的包络检波频谱'); %标题为‘AM已调信号的包络检波频谱’ axis([0,250,0,2000000]); %设置显示范围
3.4.2仿真波形
图5 AM已调信号的包络检波波形及频谱
10
3.5 AM波解调（同步乘积型检波法）
3.5.1 仿真程序
function [ ] =Jietiao2()
%经过AM调制产生已调信号u3
u5=u3.*cos(w1*t); %滤波前的解调信号 figure(5); %新建图形窗口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]); 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)); figure(6) [h,f]=freqz(b,1,512);
%将图形窗口一分为二 %绘制滤波前的解调信号波形 %横坐标为t，纵坐标显示为u5 %设置标题 %对u5进行傅里叶变换
%将第二个子图形窗口作为当前图形窗口 %绘制Y5的波形 %巴特沃斯低通滤波器 %待滤波信号频率 %采样频率 %定义过度带宽 %定义窗函数的长度 %定义滤波器的阶数 %使用fir1函数设计滤波器
%输入的参数分别是滤波器的阶数和截止频率 %新建图形窗口6 %滤波器的幅频特性图
%[H,W]=freqz(B,A,N)当N是一个整数时函数返回N点的频率向量和幅频响应向量
11
plot(f*fs/(2*pi),20*log10(abs(h))) %参数分别是频率与幅值 xlabel('频率/赫兹');ylabel('增益/分贝'); title('滤波器的增益响应'); grid
u6=filter(b,1,u5); %滤波后的解调信号 figure(7); %新建图形窗口7
subplot(2,1,1); %将图形窗口一分为二，并把第一个作为当前
图形窗口
plot(t,u6); %绘制滤波后的解调信号波形 xlabel('t');ylabel('u6'); title('滤波后的解调信号波形');
Y6=fft(u6); %对u6进行傅里叶变换
subplot(2,1,2); %将第二个子图形窗口作为当前图形窗口 plot(abs(Y6)); %绘制Y6的波形 xlabel('w');ylabel('Y6'); title('滤波后的解调信号频谱'); axis([0,250,0,600000]);
3.5.2仿真波形
图6 滤波器的增益响应
12