声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2595|回复: 9

[滤波] 短时傅里叶变换及其逆变换问题

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

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

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

x
各位老师好,我想问下为什么对一段非平稳信号短时傅里叶变换然后再重新逆变换回来,信号大体上没变,但是在开头和结尾都有一定量的衰减呢?

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2015-5-20 13:03 | 显示全部楼层
边界效应?不知该如何解释
 楼主| 发表于 2015-5-22 08:46 | 显示全部楼层
yghit08 发表于 2015-5-20 13:03
边界效应?不知该如何解释

谢谢您的回答,我想问下这种衰减会不会和短时傅里叶加窗引起泄露有关呢?
发表于 2015-5-23 14:37 | 显示全部楼层
我对傅里叶分析的算法实现不清楚,待论坛里高人回答
发表于 2015-5-23 17:09 | 显示全部楼层
逆短时傅里叶变换,恢复出来的信号,存在边界效应。主要是由于在信号的开始与末尾,被窗函数截断时,最少时只有一半的信号会被截断。这是由于离散信号造成的。也有短时傅里叶方法可以克服边界效应。如,对恢复出来的信号,引入一补偿系数,叫Padding 方法。再如,利用矩形框截断信号,不做等时截断。再如,你把短时傅里叶系数全部加起来,取实部,仅针对单个信号。
发表于 2015-5-25 08:56 | 显示全部楼层
你所在的用户组无法查看附件
 楼主| 发表于 2015-5-26 16:18 | 显示全部楼层
yugang2010 发表于 2015-5-23 17:09
逆短时傅里叶变换,恢复出来的信号,存在边界效应。主要是由于在信号的开始与末尾,被窗函数截断时,最少时 ...

谢谢您的答复~~~~~
 楼主| 发表于 2015-5-26 16:32 | 显示全部楼层
yugang2010 发表于 2015-5-23 17:09
逆短时傅里叶变换,恢复出来的信号,存在边界效应。主要是由于在信号的开始与末尾,被窗函数截断时,最少时 ...

您能详细讲讲边界效应的原因吗,有点看不懂,谢谢您!!!!
发表于 2015-5-26 16:50 | 显示全部楼层
13869165136 发表于 2015-5-26 16:32
您能详细讲讲边界效应的原因吗,有点看不懂,谢谢您!!!!

一般情况下,短时傅里叶变换,使用移动的高斯窗"g(u-t)"去截断信号"s(u)",然后进行傅里叶变换(fft),这样就得到了时间点t对应的傅里叶变换“S(t,f)”。通过移动窗函数,不停地去截断信号,则会得到此信号在各个时间点的频率分布,即所谓的时频分布,或者时频分析。

当窗函数,在截断信号的开头与结尾时,分析信号长度不够了。此时,窗函数只能截断小于窗函数长度的信号,这种情况下,所做的fft,得到的幅值会小于其他没有截断误差的信号幅值。也就是说,STFT的开头与结尾部分的幅值,小于其他时间点的幅值。

这么说,可能还是不清楚。我建议你,利用一段正弦波,一步步地跟踪你的stft代码,观测里面数值的变化,以及窗函数的变化,就能明白了。否则,理解总是无法深入进去。stft的代码,只是看起来麻烦点,实际上挺简单的。主要代码,可能都不到十行。就是一循环,利用窗函数截断信号。然后,对截断信号,进行fft。之前的代码,主要负责各种参数的初始化,如窗函数的长度,等等。

切莫眼高手低。

评分

1

查看全部评分

 楼主| 发表于 2015-6-1 17:44 | 显示全部楼层
yugang2010 发表于 2015-5-26 16:50
一般情况下,短时傅里叶变换,使用移动的高斯窗"g(u-t)"去截断信号"s(u)",然后进行傅里叶变换(fft),这 ...

非常感谢您,大神!!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-15 16:38 , Processed in 0.075793 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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