%****FM調(diào)制余弦函數(shù)*******
%調(diào)制信號頻率：10Hz；調(diào)制信號幅度：1V；載波信號頻率：100Hz；載波信號幅度：1V；
%調(diào)頻比例常數(shù)kf：10；采樣頻率fs:?500Hz。
dt=0.001;??%設(shè)定時間步長
t=0:dt:1.5;?%產(chǎn)生時間向量
am=1;%設(shè)定調(diào)制信號幅度
fm=10;%設(shè)定調(diào)制信號頻率
mt=am*cos（2*pi*fm*t）;%生成調(diào)制信號
fc=100;%設(shè)定載波頻率
Vc=cos（2*pi*fc*t）;%生成載波
kf=1;%設(shè)定調(diào)頻指數(shù)
sfm=cos（2*pi*fc*t+am*kf*sin（2*pi*fm*t））;%調(diào)制，產(chǎn)生已調(diào)信號
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變換，求絕對值得到瞬時幅度（包絡(luò)檢波）
?zero=（max（diff_Vt2）+min（diff_Vt2））/2;
?diff_Vt3=diff_Vt2-zero;
figure（1）
subplot（221）;plot（tmt）;%繪制調(diào)制信號的時域圖
title（‘調(diào)制信號的時域圖‘）;
?subplot（223）;plot（tV