[Matlab]模拟相位调制寄解调

clear
close all
clc
Fs=100;
t=[0:2*Fs+1]'/Fs;
x=sin(2*pi*t)+sin(4*pi*t);
Fc=10;
phasedev=1;  %调制度
y=pmmod(x,Fc,Fs,phasedev);
%y=awgn(y,10,'measured',103);
z=pmdemod(y,Fc,Fs,phasedev);
figure;
plot(t,x,'k-',t,z,'g-');
figure
plot(t,y);
hold on;
y2=cos(2*pi*Fc*t+x*phasedev);
plot(t,y2,'r');
z2=pmdemod(y2,Fc,Fs,phasedev);
figure;

plot(t,x,'k-',t,z2,'g-');