am.m

信号与系统时域、频域分析及MATLAB软件的应用及仿真m程序

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('f(Hz)');
ylabel('功率谱(X)');
title('被调信号的功率谱')
grid
y=modulate(x,Fc,Fs,'am');
subplot(233)
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('f(Hz)');
ylabel('功率谱(Y)');
title('已调信号的功率谱');
grid