matlab滤波

夏夜

这个问题有段时间了，来求助大家：

目的：想实现20-90day的带通滤波，附件中是一组涡动能数据（>0）.

工具：matlab中的filtfilt滤波

问题：滤波的结果出现了负值。原始数据都是大于0的，滤波结果的负值怎么解释？

谢谢各位。

程序如下：
load ekei.mat
fh=90; %90day
hl=20;  %
n=5;    %阶数
dt=1; %采样间隔1day
wn_bandpass=[1/fh/(1/(2*dt)) 1/fl/(1/(2*dt))]; %转为频率
[filter_bandpassa,filter_bandpassb]=butter(n,wn_bandpass); % bandpass
fff_bandpass=filtfilt(filter_bandpassa,filter_bandpassb,ekei); %滤波结果，时间序列


plot(ekei)
hold on
plot(fff_bandpass,'r')





另，
大家可以尝试下附件中的lanczosfilter，可以做低通或者高通滤波，带通可以用低通减去高通得到。可以与已有结果比较，有问题欢迎继续讨论。

夏夜

数据太长 已删除

lysx

这是一堆杀啊

rs06

如果滤波结果反映的是周期很小的波动情况时，很有可能出现负值，我之前对空气污染序列滤波时，短期分量就有负值。

夏夜

lysx 发表于 2014-10-23 21:29
这是一堆杀啊

抱歉，mat数据上传后需要金钱下载，复制后这么一堆，又不知怎么删掉。

夏夜

rs06 发表于 2014-10-24 08:54
如果滤波结果反映的是周期很小的波动情况时，很有可能出现负值，我之前对空气污染序列滤波时，短期分量就有 ...

嗯 这种负值该怎么解释呢

小淼在梦游№

我的原始数据是降雨量滤波后也会有负值，因为权重系数有负的，而且滤波后反应的是周期变化，应该有正负位相

夏夜

小淼在梦游№ 发表于 2014-10-24 14:57
我的原始数据是降雨量滤波后也会有负值，因为权重系数有负的，而且滤波后反应的是周期变化，应该有正负位相

滤波结果大小不是表示信号振幅吗

桥土豆

还是没搞清楚

zhaoyesi

感谢楼主分享，正在学习中