% 幅度调制 AM 信号输出和功率谱 程序清单
dt=0.001; %时间采样频谱
fmax=1; %信源最高频谱
fc=10; %载波中心频率
T=5; %信号时长
N=T/dt;
t=[0:N-1]*dt;
mt=sqrt(2)*cos(2*pi*fmax*t); %信源
A=2;
s_am=(A+mt).*cos(2*pi*fc*t);
[f,Xf]=FFT_SHIFT(t,s_am); %调制信号频谱
PSD=(abs(Xf).^2)/T; %调制信号功率谱密度
figure(1)
subplot(211);
plot(t,s_am);hold
on; %画出AM信号波形
plot(t,A+mt,'r--'); %表示AM包络
title('AM调制信号及其包络');
xlabel('t');
subplot(212); %画出功率谱图形
plot(f,PSD);
axis([-2*fc 2*fc 0 1.5*max(PSD)]);
title('AM信号功率谱');
xlabel('f');
% ----------------------------------------