ex0905.m

来自「信号与系统分析及MATLAB实现一书中的所有源代码」· M 代码 · 共 38 行

%例9.5调制信号的波形及其功率谱；
Fs=1000;
Fc=400;
N=1000;
n=0:N-2;
t=n/Fs;
x=sin(2*pi*50*t);
subplot(221)
plot(t,x);
xlabel('t(s)');
ylabel('x');
title('被调信号');
axis([0 0.1 -1 1])
Nfft=1024;
window=hamming(512);
noverlap=256;
dflag='none';
[Pxx,f]=psd(x,Nfft,Fs,window,noverlap,dflag);
subplot(222)
plot(f,Pxx)
xlabel('频率(Hz)');
ylabel('功率谱(X)');
title('被调信号的功率谱')
grid
y=modulate(x,Fc,Fs,'am');
subplot(223)
plot(t,y)
xlabel('t(s)');
ylabel('y');
axis([0 0.1 -1 1])
title('已调信号')
[Pxx,f]=psd(y,1024,Fs,window,noverlap,dflag);
subplot(224)
plot(f,Pxx)
xlabel('频率(Hz)');
ylabel('功率谱(Y)');
title('已调信号的功率谱');
grid