matlab快速傅里叶变换

Question

t=-10:0.01:10;
A0=100; %载波信号振幅
A1=0.05; %调制信号振幅
f=50; %载波信号频率
Id=17;
fc=1000+100*Id; %调制信号频率
mes=A1*sin(f*2*pi*t)./(f*2*pi*t);%调制信号
uam=A0*(mes).*cos(2*pi*fc*t);%AM已调信号
subplot(221);
plot(t,mes);
xlabel('t'),title('调制信号');
axis([-0.1,0.1,-0.1,0.1])

subplot(222);
plot(t,uam);
grid on;
title('AM已调信号波形');
axis([-0.01,0.01,-10,10]) ;

subplot(223);
Y1=fft(mes); % 对调制信号进行傅里叶变换
plot(abs(Y1));
title('调制信号频谱');
axis([-500,500,0,0.001]);

subplot(224);
Uam=fft(uam); % 对AM已调信号进行傅里叶变换
plot(abs(Uam));
title('AM已调信号频谱');
axis([-20000,20000,0,0.01])
为什么得不到想要的信号的傅里叶变换呢。。

Answer 1

几个问题，

（1）t设置的采样间隔有点大，改的再小一写。即将

t=-10:0.01:10;

改为

t=-10:0.001/pi:10;

（2）matlab中的fft之后的范围是（0-2π），要使用fftshift转换到（-π -- π），即将

Y1=fft(mes);

改为

Y1=fft(mes);         
Y1 = fftshift(Y1);

（3）频域的范围咋能按照时域的写呢，你得换算成频域的范围。

可以将

axis([-500,500,0,0.001]);

去掉，换算好之后再画。

追问

得出的图还是不对啊，和理论计算差太多了。。

追答

哪一块不对，跟理论差别大？