武汉理工大学《信号分析处理》课程设计说明书

图4.7解调信号频谱

图4.8加噪声解调信号频谱

4.3解调信号频谱分析

调制前信号频谱与解调后信号频谱比较：

18

武汉理工大学《信号分析处理》课程设计说明书

图4.9频谱比

调制后的信号频谱与调制器前的相似

图4.10调制信号与解调后的信号

由上图通过比较解调出来的信号与调制前的信号，可知二者信号大致相同，但是解调出来的有相位的偏移，但是不影响信号的质量。

19

武汉理工大学《信号分析处理》课程设计说明书

5 MATLAB代码仿真

5.1程序流程图

初始化程序FM调制程序添加高斯白噪声程序FM解调程序绘制频谱图程序显示程序 图5.1程序流程图

5.2 程序代码

%FM调制解调系统.m

%频率调制与解调的Matlab演示源程序 %电信1006班 石勇

%*****************初始化****************** echo off close all clear all clc

%****************FM调制*******************

dt=0.001; %设定时间步长 t=0:dt:1.5;

%产生时间向量

20

武汉理工大学《信号分析处理》课程设计说明书

am=5; %设定调制信号幅度 fm=5; %设定调制信号频率 mt=am*cos(2*pi*fm*t); %生成调制信号 fc=50; %设定载波频率 ct=cos(2*pi*fc*t); %生成载波 kf=10; %设定调频指数 int_mt(1)=0; for i=1:length(t)-1

int_mt(i+1)=int_mt(i)+mt(i)*dt; %求信号m(t)的积分 end %调制，产生已调信号 sfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt); %调制信号 %*************添加高斯白噪声**************

sn1=10; %设定信躁比(小信噪比) sn2=30; %设定信躁比(大信噪比) sn=0; %设定信躁比(无信噪比) db=am^2/(2*(10^(sn/10))); %计算对应的高斯白躁声的方差 n=sqrt(db)*randn(size(t)); %生成高斯白躁声

nsfm=n+sfm; %生成含高斯白躁声的已调信号（信号通

%过信道传输）

%****************FM解调*******************

for i=1:length(t)-1 %接受信号通过微分器处理 diff_nsfm(i)=(nsfm(i+1)-nsfm(i))./dt; end

diff_nsfmn = abs(hilbert(diff_nsfm)); %hilbert变换，求绝对值得到瞬时幅度（包络检波） zero=(max(diff_nsfmn)-min(diff_nsfmn))/2; diff_nsfmn1=diff_nsfmn-zero;

%**************时域到频域转换**************

ts=0.001; %抽样间隔

21