%****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