傅里叶变换中的一个问题

dingdingysu · 发表于 2012-3-4 19:32

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

本帖最后由 dingdingysu 于 2012-3-4 20:07 编辑

我要求信号的傅里叶变换，程序如下。但是求出的幅频图是对称的。而且在200和800处有幅值，请问真是怎么回事？

N=1000;
fs=2000;
t=0:1/fs:(N-1)/fs;
x=sin(400*pi*t)
y=fft(x,N);%进行fft变换
mag=abs(y);%求幅值
f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换
figure(1);
subplot(232);
plot(f,mag);%做频谱图

321forever · 发表于 2012-3-4 20:46

本帖最后由 321forever 于 2012-3-4 20:51 编辑

在200的地方是实际的频率处，而800是对称的，这个是对离散点的傅立叶变化产生的，可以只取一半的频率进行研究就可以了，具体可以参考wiki aliasing

这幅图来源于http://en.wikipedia.org/wiki/Aliasing
从图中就可以看出，离散点可以表示两种频率的sin函数，所以fft后出来的是在200 和800处出现了对称的峰值

Ericwjc · 发表于 2012-3-5 14:28

0~N/2点为正频率，N/2~N-1点为负频率，它们关于中点对称。所以实际中只画出0~N/2点点即可，将其幅值乘以2（正负频率相应点的幅值相等）。

dingdingysu · 发表于 2012-3-6 16:42

回复 2 # 321forever 的帖子

O(∩_∩)O谢谢

dingdingysu · 发表于 2012-3-6 16:42

回复 3 # Ericwjc 的帖子

我得好好学习啊...

账号		自动登录	找回密码
密码			我要加入

[综合讨论] 傅里叶变换中的一个问题

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

评分

评分