MATLAB进行FM调制与解调

%****FM调制余弦函数*******
%调制信号频率：10Hz；调制信号幅度：1V；载波信号频率：100Hz；载波信号幅度：1V；
%调频比例常数kf：10；采样频率fs: 500Hz。
dt=0.001;  %设定时间步长
t=0:dt:1.5; %产生时间向量
am=1;%设定调制信号幅度
fm=10;%设定调制信号频率
mt=am*cos（2*pi*fm*t）;%生成调制信号
fc=100;%设定载波频率
Vc=cos（2*pi*fc*t）;%生成载波
kf=1;%设定调频指数
sfm=cos（2*pi*fc*t+am*kf*sin（2*pi*fm*t））;%调制，产生已调信号
Vt1=sfm;
for i=1:length（t）-1 %接受信号通过微分器处理
 diff_Vt1（i）=（Vt1（i+1）-Vt1（i））./dt;
 end
 diff_Vt2=abs（hilbert（diff_Vt1））;%hilbert变换，求绝对值得到瞬时幅度（包络检波）
 zero=（max（diff_Vt2）+min（diff_Vt2））/2;
 diff_Vt3=diff_Vt2-zero;
figure（1）
subplot（221）;plot（tmt）;%绘制调制信号的时域图
title（‘调制信号的时域图‘）;
 subplot（223）;plot（tV