clc;
close?all;
clear?all;
%%?20181120
%?by呂冬治
%?信號的FM調制與解調
%%?信號參數
fm?=?3;????
fc?=?10;????
fs?=?200;????
Am?=?1;???????
A?=?2;????????
N?=?1e3;??????
K?=?N-1;
n?=?0:K;
f?=?n*fs/N;???
t?=?0:1/fs:K/fs;????????
wm?=?2*pi*fm;?
wc?=?2*pi*fc;?
thm?=?0;???????
thc?=?0;
B?=?2*fm;
KFM?=?10;
beatFM?=?KFM*Am/wc;
ct?=?A*cos（wc*t）;%???
fft_ct=?abs（fft（ct））;
figure（1）
subplot（211）;
plot（tct）;
xlabel（‘時間/s‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘頻率為10Hz的載波信號byLDM‘）;
subplot（212）;
plot（f（1:N/2）fft_ct（1:N/2））;
xlabel（‘頻率/Hz‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘載波信號信號頻譜‘）;

%%?調制FM
sig_fm?=?A*cos（wc*t+beatFM*sin（wm*t））;
fft_sig_fm?=?abs（fft（sig_fm））;
figure（2）
subplot（211）;
plot（tsig_fm）;
xlabel（‘時間/s‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘調制信號FM‘）;
subplot（212）;
plot（f（1:N/2）fft_sig_fm（1:N/2））;
xlabel（‘頻率/Hz‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘調制信號頻譜‘）;
%%?加噪聲
sig_noise_fm?=?awgn（sig_fm40）;?????%?添加信噪比為6dB的高斯白噪聲
fft_sig_noise_fm?=?abs（fft（sig_noise_fm））;?????????%?對AM解調信號進行傅里葉變換??
sig_noise_fm2?=?awgn（sig_fm20）;?????%?添加信噪比為15dB的高斯白噪聲
fft_sig_noise_fm2?=?abs（fft（sig_noise_fm））;?????????%?對AM解調信號進行傅里葉變換??
figure（3）
subplot（211）;
plot（tsig_noise_fm）;
xlabel（‘時間/s‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘添加信噪比為40dB噪聲后的FM信號byLDM‘）;
subplot（212）;
plot（f（1:N/2）fft_sig_noise_fm（1:N/2））;
xlabel（‘頻率/Hz‘）;
ylabel（‘幅值‘）;
grid;
title（‘添加信噪比為40dB噪聲后的FM信號頻譜‘）;

figure（4）
subplot（211）;
plot（tsig_noise_fm2）;
xlabel（‘時間/s‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘添加信噪比為20dB噪聲后的FM信號byLDM‘）;
subplot（212）;
plot（f（1:N/2）fft_sig_noise_fm2（1:N/2））;
xlabel（‘頻率/Hz‘）;
ylabel（‘幅值‘）;
grid;
title（‘添加信噪比為20dB噪聲后的FM信號頻譜‘）;
%%?解調
%?無噪聲
diff_sig_fm?=?diff（sig_fm）;?????????%?調頻信號經過一個微分器
d_m?=?abs（hilbert（diff_sig_fm））;????%?用希爾伯特變換能夠實現包絡檢波
fft_dm?=?abs（fft（d_m））;???????????%?對AM解調信號進行傅里葉變換??
%?3?dB高斯白噪聲
diff_sig_fm1?=?diff（sig_noise_fm）;?????????%?調頻信號經過一個微分器
d_m2?=?abs（hilbert（diff_sig_fm1））;????%?用希爾伯特變換能夠實現包絡檢波
fft_dm2?=?abs（fft（d_m2））;???????????%?對AM解調信號進行傅里葉變換??
%?20?dB高斯白噪聲
diff_sig_fm2?=?diff（sig_noise_fm2）;?????????%?調頻信號經過一個微分器
d_m3?=?abs（hilbert（diff_sig_fm2））;????%?用希爾伯特變換能夠實現包絡檢波
fft_dm3?=?abs（fft（d_m3））;???????????%?對AM解調信號進行傅里葉變換??
%?diff_sig_fm2?=?diff（sig_noise_fm2）;?????????%?調頻信號經過一個微分器
%?d_m3?=?abs（hilbert（diff_sig_fm2））;????%?用希爾伯特變換能夠實現包絡檢波
%?fft_dm3?=?abs（fft（d_m3））;???????????%?對AM解調信號進行傅里葉變換??
%?diff_sig_fm2?=?diff（sig_noise_fm2）;?????????%?調頻信號經過一個微分器
%?d_m3?=?abs（hilbert（diff_sig_fm2））;????%?用希爾伯特變換能夠實現包絡檢波
%?fft_dm3?=?abs（fft（d_m3））;???????????%?對AM解調信號進行傅里葉變換??
%?diff_sig_fm2?=?diff（sig_noise_fm2）;?????????%?調頻信號經過一個微分器
%?d_m3?=?abs（hilbert（diff_sig_fm2））;????%?用希爾伯特變換能夠實現包絡檢波
%?fft_dm3?=?abs（fft（d_m3））;???????????%?對AM解調信號進行傅里葉變換??
%?diff_sig_fm2?=?diff（sig_noise_fm2）;?????????%?調頻信號經過一個微分器
%?d_m3?=?abs（hilbert（diff_sig_fm2））;????%?用希爾伯特變換能夠實現包絡檢波
%?fft_dm3?=?abs（fft（d_m3））;???????????%?對AM解調信號進行傅里葉變換??
%?diff_sig_fm2?=?diff（sig_noise_fm2）;?????????%?調頻信號經過一個微分器
%?d

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????5429??2018-11-25?16:34??FM.m

-----------?---------??----------?-----??----

?????????????????5429????????????????????1