裁缝Butter+filtfilt滤波函数

一夜荒原

参考论坛里的之前交流的用MATLAB设计滤波的帖子，用butter和filtfilt函数裁缝一个普适性比较强的（个人认为）滤波函数，最高可用于3维矩阵的滤波。

这个函数的优点在于：（1）对于三维的变量场，不需要写循环，运行速度会很快；（2）矩阵中含有NaN的格点可以自动排除掉，省去判断if ~isnan(y)的过程。

需要注意的是：

如果你的数据是二维的，必须把时间维放到第一个维度。因为filtfilt这个函数只支持第一维度的滤波（不像filter函数可以指定维度）。

自己和写循环的方式对比过，结果无差异，应该没什么问题吧。

函数使用过程如果有什么问题，或函数本身有什么本质上的错误可随时联系我！谢谢！

474620961

请教一下，n（阶数）应该怎么选；
如果是用bandpass，参数应该如何设置？

一夜荒原

474620961 发表于 2023-3-5 09:46
请教一下，n（阶数）应该怎么选；
如果是用bandpass，参数应该如何设置？

阶数选取我也不太清楚原则；
bandpass的Hz可以设置为[1/n 1/m]就好