clc;
close?all;
clear?all;
%%?20181120
%??by呂冬治
%?信號的PM調(diào)制與解調(diào)
%?此程序給出了信號的PM調(diào)制的與解調(diào)的仿真過程，并分析了在不同信噪比條件下解調(diào)的性能
%%?信號參數(shù)
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?=?5;
beatFM?=?KFM*Am;
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（‘載波信號信號頻譜‘）;
%%?調(diào)制PM
sig_PM?=?cos（wc*t+beatFM*cos（2*pi*wm*t））;%PM已調(diào)制信號
fft_sig_PM?=?abs（fft（sig_PM））;
figure（2）
subplot（211）;
plot（tsig_PM）;
xlabel（‘時間/s‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘調(diào)制信號tPM‘）;
subplot（212）;
plot（f（1:N/2）fft_sig_PM（1:N/2））;
xlabel（‘頻率/Hz‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘調(diào)制信號頻譜‘）;
%%?加噪聲
sig_noise_PM1?=?awgn（sig_PM3）;?????%?添加信噪比為6dB的高斯白噪聲
fft_sig_noise_PM1?=?abs（fft（sig_noise_PM1））;?????????%?對AM解調(diào)信號進行傅里葉變換??
sig_noise_PM2?=?awgn（sig_PM20）;?????%?添加信噪比為15dB的高斯白噪聲
fft_sig_noise_PM2?=?abs（fft（sig_noise_PM2））;?????????%?對AM解調(diào)信號進行傅里葉變換??
figure（3）
subplot（211）;
plot（tsig_noise_PM1）;
xlabel（‘時間/s‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘添加信噪比為40dB噪聲后的FM信號byLDM‘）;
subplot（212）;
plot（f（1:N/2）fft_sig_noise_PM1（1:N/2））;
xlabel（‘頻率/Hz‘）;
ylabel（‘幅值‘）;
grid;
title（‘添加信噪比為40dB噪聲后的FM信號頻譜‘）;

figure（4）
subplot（211）;
plot（tsig_noise_PM2）;
xlabel（‘時間/s‘）;
ylabel（‘幅值‘）;
grid?on;?
title（‘添加信噪比為20dB噪聲后的FM信號byLDM‘）;
subplot（212）;
plot（f（1:N/2）fft_sig_noise_PM2（1:N/2））;
xlabel（‘頻率/Hz‘）;
ylabel（‘幅值‘）;
grid;
title（‘添加信噪比為20dB噪聲后的FM信號頻譜‘）;
%?%%?解調(diào)
%?%?無噪聲
%?diff_sig_PM?=?diff（sig_PM）;?????????%?調(diào)頻信號經(jīng)過一個微分器
%?d_m?=?abs（hilbert（diff_sig_PM））;????%?用希爾伯特變換能夠?qū)崿F(xiàn)包絡(luò)檢波
%?fft_dm?=?abs（fft（d_m））;???????????%?對AM解調(diào)信號進行傅里葉變換??
%?%?3?dB高斯白噪聲
%?diff_sig_PM1?=?diff（sig_noise_PM1）;?????????%?調(diào)頻信號經(jīng)過一個微分器
%?d_m2?=?abs（hilbert（diff_sig_PM1））;????%?用希爾伯特變換能夠?qū)崿F(xiàn)包絡(luò)檢波
%?fft_dm2?=?abs（fft（d_m2））;???????????%?對AM解調(diào)信號進行傅里葉變換??
%?%?20?dB高斯白噪聲
%?diff_sig_fm2?=?diff（sig_noise_fm2）;?????????%?調(diào)頻信號經(jīng)過一個微分器
%?d_m3?=?abs（hilbert（diff_sig_fm2））;????%?用希爾伯特變換能夠?qū)崿F(xiàn)包絡(luò)檢波
%?fft_dm3?=?abs（fft（d_m3））;???????????%?對AM解調(diào)信號進行傅里葉變換??
%?
%?figure（5）
%?subplot（211）;
%?plot（t（1:end-1）d_m）;
%?xlabel（‘時間/s‘）;
%?ylabel（‘幅值‘）;
%?grid?on;?
%?title（‘無噪聲條件下解調(diào)信號波形FMbyLDM‘）;
%?subplot（212）;
%?plot（f（1:N/2）fft_dm（1:N/2））;
%?xlabel（‘頻率/Hz‘）;
%?ylabel（‘幅值‘）;
%?axis（[0?100?0?100]）;
%?grid;
%?title（‘FM解調(diào)信號頻譜byCL‘）;
%?
%?figure（6）
%?subplot（211）;
%?plot（t（1:end-1）d_m2）;
%?xlabel（‘時間/s‘）;
%?ylabel（‘幅值‘）;
%?grid?on;?
%?title（‘信噪比為40dB條件下解調(diào)信號波形FMbyLDM‘）;
%?subplot（212）;
%?plot（f（1:N/2）fft_dm2（1:N/2））;
%?xlabel（‘頻率/Hz‘）;
%?ylabel（‘幅值‘）;
%?axis（[0?100?0?100]）;
%?grid;
%?title（‘FM解調(diào)信號頻譜byCL‘）;
%?
%?figure（7）
%?subplot（211）;
%?plot（t（1:end-1）d_m3）;
%?xlabel（‘時間/s‘）;
%?ylabel（‘幅值‘）;
%

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

?????文件???????4591??2018-11-25?16:36??pm_.m

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

?????????????????4591????????????????????1