FM的MATLAB调制与解调仿真程序，可以直接运行

FM的MATLAB调制与解调仿真程序，可以直接运行，需要的下载

clc; clear;
t0=-0.1; t1=0.1; %调制信号的时域范围 
fs=1000; ts=1/fs; %采样率
t=t0:ts:t1;
m=sinc（100*t）; %调制信号fm=50;
fc=300; %载频
x=cos（2*pi*fc*t）; %载波信号
kf=100; %偏差常数
u=GetFMSignal（mt0t1tsfckf）; %已调制信号
%画调制信号和已调制信号的波形
subplot（321）;plot（tm）;title（‘基带信号‘）;
subplot（322）;plot（tu）;title（‘已调制信号‘）;

df=0.25; %FFT时的频率分辨率
f=[0:df:df*（length（t）-1）]-fs/2; %与t对应的频率向量
[Mm_newdf_new_m]=GetFFT（mt0t1tsdf）;
[Uu_newdf_new_u]=GetFFT（ut0t1tsdf）;
%FFT变换后，对两个频率向量重新赋值
f1=[0:df_new_m:df_new_m*（length（m_new）-1）]-fs/2;
f2=[0:df_new_u:df_new_u*（length（u_new）-1）]-fs/2;
subplot（323）;plot（f1abs（fftshift（M）））;xlabel（‘基带信号频谱‘）;
subplot（324）;plot（f2abs（fftshift（U）））;xlabel（‘已调制信号频谱‘）;

phase=GetDEM（ut0t1tsfc）; %解调，得到相角
phase_new=unwrap（phase）; %对相角平滑处理
m_dem=1/（2*pi*kf）*diff（phase_new）;
m_dem_new=[m_dem0]; %求微分后少了一个点，补0
subplot（325）;plot（tm）;xlabel（‘原基带信号‘）;
subplot（326）;plot（tm_dem_new）;xlabel（‘解调后基带信号‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1077  2008-12-07 14:47  FM\FM.m

     文件        126  2008-12-07 14:48  FM\GetDEM.m

     文件        347  2008-12-07 14:49  FM\GetFFT.m

     文件        201  2008-12-07 09:50  FM\GetFMSignal.m

     文件        135  2008-12-07 14:13  FM\GetQuaSignal.m

     文件        143  2008-12-07 09:49  FM\m_INT.m

     文件     137592  2018-11-17 18:08  fm.docx

     目录          0  2008-12-12 18:30  FM

----------- ---------  ---------- -----  ----

               139621                    8