模拟角度调制系统仿真,实现FM与PM调制
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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频谱图');。