%?PM調(diào)制
figure（‘Name‘‘信號(hào)調(diào)制過(guò)程中波形及其頻譜‘‘Numbertitle‘‘off‘）
a0=2;f0=10;fc=50;fs=1000;snr=5;
t=?linspace（-202060001）;
pm1=cos（2*pi*f0*t）;?%信息信號(hào)
t1=cos（2*pi*fc*t）;?%載波
s_pm=cos（2*pi*fc*t+1*pm1）;
PM1=fft（pm1）;?T1=fft（t1）;?S_PM=fft（s_pm）;
f=（0:60000）*fs/60001-fs/2;
subplot（321）;?plot（t（19801:20200）pm1（19801:20200））;?title（‘信息信號(hào)波形‘）;
subplot（322）;?plot（ffftshift（abs（PM1）））;?title（‘信息信號(hào)頻譜‘）;
subplot（323）;?plot（t（19801:20200）t1（19801:20200））;?title（‘載波信號(hào)‘）;
subplot（324）;?plot（ffftshift（abs（T1）））;?title（‘載波信號(hào)頻譜‘）;
subplot（325）;?plot（t（19801:20200）s_pm（19801:20200））;?title（‘已調(diào)信號(hào)‘）;
subplot（326）;?plot（ffftshift（abs（S_PM）））;?title（‘已調(diào)信號(hào)頻譜‘）;
?
%產(chǎn)生噪聲
figure（‘Name‘‘添加噪聲及帶通濾波過(guò)程波形及其頻譜‘‘Numbertitle‘‘off‘）;
y=awgn（s_pmsnr）;
fsamp?=?1000;?????????%采樣頻率為?1MHz?
fcuts?=?[95?105];???%通帶截止頻率為?95Hz，阻帶截止頻率為?105Hz?
mags?=?[1?0];?
devs?=?[0.01?0.05];????%??通帶波動(dòng)?1%，阻帶波動(dòng)?5%?
[nWnbetaftype]?=?