信号预处理--提取趋势项

ken1983414

A Detrending Method in Data Processing
We all know that spectral estimation request that the sigal is at least weakly stationary. So, detrending is just trying to remove the slow nonstationary trends from the data that need processing.
This is a detrending method using smoothness priors aproach(SPA).
Ref. An advanced detrending method with application to HRV analysis. IEEE trans. on Biomedical Engineering
An example for detrending. the sinusoidal signal that have lower frequency is used as the trend.

[ 本帖最后由 ken1983414 于 2009-8-7 14:09 编辑 ]

luoye919

LZ你想说什么呢？(⊙_⊙)？

ChaChing

楼主是分享其程序!?

ken1983414

想分享这个做数据预处理，提取趋势项的方法，觉得效果不错。

aini2009

那楼主为何不把程序贴出来？

李元

大侠你好，可不可以麻烦你几秒钟给我说一下工程上对于振动信号的预处理方法（步骤）呢？谢谢:loveliness:

ken1983414

主要是 ，分享信号预处理--提取趋势项的方法
发现这个方法效果不错

pengchen

LZ发的帖子 在论坛其他地方看见过啊
clear;
fs=1000;    %fs为采样频率;
t=1/fs:1/fs:1;
x1=2*sin(100*pi*t);
subplot(411);plot(x1);
title('原始信号')
x2=3*sin(pi*t);
subplot(412);plot(x2);
title('趋势项')
x=x1+x2;
subplot(413);plot(x);
title('原始信号＋趋势项')

m=2;
a=polyfit(t,x,m);     %计算多项式待定系数向量a
y=x-polyval(a,t);     %用x减去多项式系数a生成的趋势项
title('original signal');
grid on;
subplot(2,1,2);
subplot(414);plot(t,y);
title('去除趋势项后得到的信号')
这就是源程序啊

hitalfa

找的辛苦

wickedperson

好好好好好

rogen

回复 6 # 李元 的帖子

一般的信号处理的书上就有介绍..........