基于MATLAB的AM调制及解调系统仿真分解 下载本文

内容发布更新时间 : 2024/11/7 20:37:26星期一 下面是文章的全部内容请认真阅读。

前图形窗口

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