一阶Butterworth带通滤波

预言

谢谢，很有意义

呆妹小霸王

我想问一下，这里边的阶数该如何选择呢？我不太理解它的阶数是什么意思，使用不同阶数计算出来的都不太一样。

伽蓝鸟

呆妹小霸王 发表于 2020-8-6 10:01
我想问一下，这里边的阶数该如何选择呢？我不太理解它的阶数是什么意思，使用不同阶数计算出来的都不太一样 ...

Murakami (1979, Mon Wea Rev)用的是1阶，Wei et al. (2019, GRL)用的是6阶，我最近用的也是6阶。阶数高的话，滤波效果会好一点；但是也不能太高了。你可以画响应函数看看。

呆妹小霸王

伽蓝鸟 发表于 2020-8-6 10:40
Murakami (1979, Mon Wea Rev)用的是1阶，Wei et al. (2019, GRL)用的是6阶，我最近用的也是6阶。阶数高 ...

什么是响应函数，我搜了下还是不知道是啥

伽蓝鸟

呆妹小霸王 发表于 2020-8-6 19:07
什么是响应函数，我搜了下还是不知道是啥

加我微信或者QQ吧，气象家园MATLAB交流群里能找到我

伽蓝鸟

经常有人问我Butterworth用几阶比较好，额，Murakami 1979 MWR用的是1阶，Wei et al. 2019 GRL用的是6阶。我之前用的是1阶，最近用的都是6阶了。滤波器的效果可以通过响应函数（滤波后的振幅/滤波前的振幅）来看，这分别是1阶、6阶、7年滑动的结果。

一诺倾情

感谢楼主！

蝴蝶振翅

楼主新年快乐！我想请问一下，如果是逐月数据滤波时如何计算采样频率呢。例如逐月数据滤2~8年的周期，用butter函数该如何写？谢谢楼主！

伽蓝鸟

蝴蝶振翅 发表于 2021-2-15 23:33
楼主新年快乐！我想请问一下，如果是逐月数据滤波时如何计算采样频率呢。例如逐月数据滤2~8年的周期，用but ...

[b,a]=butter(6,[1/96,1/24]/0.5,'bandpass');
y=filtfilt(b,a,y);
抱歉才看到消息，我觉得是这样的？

蝴蝶振翅

伽蓝鸟 发表于 2021-2-17 11:12
=butter(6,[1/96,1/24]/0.5,'bandpass');
y=filtfilt(b,a,y);
抱歉才看到消息，我觉得是这样的？

十分感谢您！