声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2517|回复: 7

[FFT] FFT变换之前,非平稳数据怎么去趋势项,急啊!

[复制链接]
发表于 2012-2-26 21:52 | 显示全部楼层 |阅读模式

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

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

x
大家好:
最近在做一些数据处理工作,就是把时域的采样数据转化为频域数据,可是这些数据本身就是非平稳数据。
那么怎么去掉非平稳数据的趋势项呢?比如,设数据的周期为T,一个周期包含N个数据,共采样M个周期,总的数据量为N*M个,采样足够多,帖子里有人提到用最小二乘法或者高通滤波来求趋势项,可是怎么具体怎么做呢?希望有人能帮忙解答,最好给个例子。
回复
分享到:

使用道具 举报

 楼主| 发表于 2012-2-26 22:05 | 显示全部楼层
在论坛里搜索了一下,用MATLAB程序可以实现趋势项的处理
http://forum.vibunion.com/viewth ... 6amp%3Btypeid%3D192
问题是:我现在用的是FORTRAN,要自己编程实现趋势项处理,FORTRAN里好像没有相关的吧?所以想搞清楚去除趋势项的基本原理,希望大家赶快给予帮助!
发表于 2012-5-23 13:56 | 显示全部楼层
不知道现在还在不在弄,先要分清楚趋势主要有什么趋势项,然后再针对性地去除,一般的话就是去掉直流偏置,然后做下高通滤波去掉低于采样率的低频成分,线性趋势的话可以用线性拟合的方式去除,然后有些非线性的,我自己是用EMD算法去的,效果还行,就是计算量大一点。
发表于 2012-6-2 10:19 | 显示全部楼层
刚接触EMD感觉复杂
 楼主| 发表于 2012-6-25 16:14 | 显示全部楼层
回复 3 # 渠爷 的帖子

我的问题可能要用高通滤波去掉,这个工作已经做过,尝试了几种高通滤波方法,非线性问题这个倒是没有注意到,EMD算法还没有接触过,我试着看看。谢谢
发表于 2012-7-19 10:16 | 显示全部楼层
建议你看看这个,里面有个omega算法,其实就是fft逆变换,一开始对时域信号进行滤波和去DC之后fft再对频域进行积分。blog.prosig.com
发表于 2012-7-19 11:01 | 显示全部楼层
编辑滤波器,这个要按照你自己的要求了。
 楼主| 发表于 2012-8-15 14:43 | 显示全部楼层
回复 6 # 渠爷 的帖子

好的,谢谢。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-11 05:13 , Processed in 0.125696 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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