模拟角度调制系统仿真,实现FM与PM调制

fc = 150; %载波频率

fs = 550; %采样频率

t = (0 :0.001:0.18); %时间区域

x = sin(3*pi*50*t); %调制信号

y = modulate(x,fc,fs,'FM');%对调制信号FM调制

z = modulate(x,fc,fs,'PM');%对调制信号PM调制

y1 = y + awgn(y,10,0);%y = awgn(x,SNR) 在信号x中加入高斯白噪声。信噪比SNR以dB为单位。x的强度假定为0dBW。如果x是复数，就加入复噪声。

%y = awgn(x,SNR,SIGPOWER) 如果SIGPOWER是数值，则其代表以dBW为单位的信号强度；如果SIGPOWER为'measured'，则函数将在加入噪声之前测定信号强度。

z1 = z + awgn(z,10,0);%已调信号

plot(t,x);

title('调制信号波形图');

xlabel('t');

ylabel('x');

plot(t,y);

xlabel('t(s)');

ylabel('y');

title('FM波形图');

plot(t,z);

xlabel('t(s)');

ylabel('y');

title('PM波形图');

a=fft(x,1024); %对x进行傅利叶变换

f1=(0:length(a)-1)*fs/length(a) -fs/2;

%figure

b=fft(y1,1024); %对y1进行傅利叶变换--plot(f2,abs(b));

f2=(0:length(b)-1)*fs/length(b) -fs/2;

%figure

c=fft(z1,1024); %对y2进行傅利叶变换

f3=(0:length(c)-1)*fs/length(c) -fs/2;

plot(f1,abs(a));

xlabel('Frequence(Hz)');

ylabel('powerSpectrum(x)');

title('调制信号频谱图');

plot(f2,abs(b));

xlabel('Frequence(Hz)'); ylabel('Power Spectrum(y)');

title('FM频谱图');

plot(f3,abs(c));

xlabel('Frequence(Hz)'); ylabel('Power Spectrum(y)');

title('PM频谱图');