对AM信号FFT的matlab仿真

%已知基带信号为1hz，载波为64hz，调制系数ma=0.3，采样频率1024hz，FFT变换区间N为2048

clear;     
fs=1024;     
f=1;     %1hz基带信号
fc=64;   %64hz载波
t=0:1/fs:2-0.5/fs;    %采样时间为0到2秒
x=(1+0.3*cos(2*pi*f*t)).*cos(2*pi*fc*t);  %调幅表达式
figure(1);   %在第一个窗口中画图
plot(t,x);     
grid;       %画方格线
%画出采样时间内的AM信号

y=fft(x);
N=length(x);
figure(2);
stem((0:N-1)*fs/N,abs(y),‘.‘);   %画出以模拟频率为横坐标，以FFT后的幅度为纵坐标
%画出AM信号做傅里叶变化后的幅度和频率图


%接下来过滤掉一半的AM信号。
for n=1:N
    if  x(n)<0
        x(n)=0;     %对应电路即是通过二极管来整流滤波
    end
end


figure(3);
y=fft(x);
stem((0:N-1)*fs/N,abs(y),‘.‘);
grid;
%画出截掉一半信号的傅里叶变换的幅度与频率图