声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 4892|回复: 24

[FFT] 在采样点少的情况下傅里叶变换后的曲线怎么不是sinc函数

[复制链接]
发表于 2007-11-12 20:01 | 显示全部楼层 |阅读模式

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

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

x
傅里叶变换的频谱应该是圆滑的曲线,即在采样点少的情况下应该是sinc函数的样子,但我在采样点少的情况下得到的频谱图是这样的,
还望赐教什么原因
fft.jpg
回复
分享到:

使用道具 举报

发表于 2007-11-12 22:29 | 显示全部楼层
什么意思,把程序帖出来看看。

跟采样点多少无关,要细化才能看到sinc的形状

[ 本帖最后由 yangzj 于 2007-11-12 22:55 编辑 ]
 楼主| 发表于 2007-11-13 10:44 | 显示全部楼层

程序如下

interferogramoerror1.m

1.52 KB, 下载次数: 65

 楼主| 发表于 2007-11-13 10:46 | 显示全部楼层
整个程序是采样之后进行量化处理,然后进行傅里叶变换的,改变采样点可以看到细化程度不同,但欠采样的时候,为什么是这样的形状呢,跟矩形截断的理论不一致呢??
发表于 2007-11-13 10:57 | 显示全部楼层
单频成分的信号(或者是间隔较远的多频信号)在信号频率附近做频域细化(不是用FFT,用离散时间傅立叶变换DTFT)才能看到sinc的形状。
程序里X11这样的信号不适合用傅立叶变换
 楼主| 发表于 2007-11-14 09:45 | 显示全部楼层

yangzj:频域与空域信号是互为傅里叶变换关系的,怎么不可以用FFT呀?

这个是信号的传统的理解吧,频域与空域互为傅里叶变换关系,对空域信号采样之后进行量化,然后进行傅里叶变换,就得到频域信息,可为什么不适合用FFT呢,我看到这方面的理论都是用基2--FFT的呀,为什么呢?
发表于 2007-11-14 10:46 | 显示全部楼层

回复 #6 xiaokongliuliu 的帖子

我昨天弄错了,看到你的信号只是在末端有个冲击.
刚才看了下,信号基本上平稳,可以用傅立叶变换.
基2是由于FFT算法的要求
发表于 2007-11-14 14:47 | 显示全部楼层
FFT是计算连续傅立叶变换的离散而又快速的方式。相当于离散积分。

当时欠采样的时候,相当于采样步长太宽,使得数值积分的误差太大而导致画出的图形面目可憎。
 楼主| 发表于 2007-11-14 15:02 | 显示全部楼层
但是矩形函数的截断效应具体表现在哪里呢?
还有信号的切趾不是为了减小截断效应带来的频谱混叠么,可是我作出来的变换后的频谱图没有体现出来矩形函数的变换sinc函数的影响呀,那这样以来怎么采用切趾函数修正呢?
还望仔细指教!!
 楼主| 发表于 2007-11-14 15:16 | 显示全部楼层
还原后的谱线怎么没有”旁瓣‘?
发表于 2007-11-14 16:54 | 显示全部楼层
我不完全明白楼主的意思,用楼主的程序我得的图和楼主1层给出的不一样(见笫1图)。如果把谱图取对数,则能很明显地看到旁瓣(见笫2图),说明是由sinc函数产生的效果。
xk2a.jpg
xk2b.jpg
 楼主| 发表于 2007-11-14 18:32 | 显示全部楼层

回复 #11 songzy41 的帖子

谢谢你的回复,你的结果就是我想要的结果,可我怎么没作出来,就是把我的程序中的采样点数增加之后,我的频谱图跟你的也不一样的呢,
 楼主| 发表于 2007-11-14 18:33 | 显示全部楼层
不知道哪里是不对的呢

[ 本帖最后由 xiaokongliuliu 于 2007-11-14 18:37 编辑 ]

2.fig

7.9 KB, 下载次数: 11

 楼主| 发表于 2007-11-14 18:35 | 显示全部楼层
刚才图的格式错了,这次更正过来,我作出来的谱图是这样的,就是没有sinc函数的影响,不知道在哪里出错了,还望楼上的指出来,非常感谢!!
2.jpg
发表于 2007-11-14 19:12 | 显示全部楼层

回复 #14 xiaokongliuliu 的帖子

改下纵座标,songzy41用的是对数座标,突出了小值.
如果你想要看的是sinc函数的形状的话,我觉得这还不是你要的结果.
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-17 12:47 , Processed in 0.070972 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表