相位调制 matlab

% PM调制 figure('Name','信号调制过程中波形及其频谱','NumberTitle','off') a0=2;f0=10;fc=50;fs=1000;snr=5; t= linspace(-20,20,60001); pm1=cos(2*pi*f0*t); %信息信号 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(3,2,1); plot(t(19801:20200),pm1(19801:20200)); title('信息信号波形');

% PM调制
figure（‘Name‘‘信号调制过程中波形及其频谱‘‘Numbertitle‘‘off‘）
a0=2;f0=10;fc=50;fs=1000;snr=5;
t= linspace（-202060001）;
pm1=cos（2*pi*f0*t）; %信息信号
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（‘信息信号波形‘）;
subplot（322）; plot（ffftshift（abs（PM1）））; title（‘信息信号频谱‘）;
subplot（323）; plot（t（19801:20200）t1（19801:20200））; title（‘载波信号‘）;
subplot（324）; plot（ffftshift（abs（T1）））; title（‘载波信号频谱‘）;
subplot（325）; plot（t（19801:20200）s_pm（19801:20200））; title（‘已调信号‘）;
subplot（326）; plot（ffftshift（abs（S_PM）））; title（‘已调信号频谱‘）;
 
%产生噪声
figure（‘Name‘‘添加噪声及带通滤波过程波形及其频谱‘‘Numbertitle‘‘off‘）;
y=awgn（s_pmsnr）;
fsamp = 1000;         %采样频率为 1MHz 
fcuts = [95 105];   %通带截止频率为 95Hz，阻带截止频率为 105Hz 
mags = [1 0]; 
devs = [0.01 0.05];    %  通带波动 1%，阻带波动 5% 
[nWnbetaftype] =