matlab - 巴特沃斯带通滤波器的命令是什么

Question

我有信号数据要使用频率为 10-45s、45-150s 和 150-600s 的带通滤波器进行过滤。如何在matlab中创建带通滤波器命令来过滤该频率？

Answer 1

假设您有信号处理工具箱（对于butter命令 - 如果butter不可用，请参阅下面的代码），并且如果您知道采样频率fs，您可以使用低频和高频（以 Hz 为单位）-3dB 点制作一个简单的递归巴特沃斯滤波器下面的代码。较高order的 s 将为您提供更好的异频抑制，但代价是更长的脉冲响应和更多的计算费用。

[b,a] = butter(order, [lowFreq hiFreq]/(fs/2), 'bandpass');
 y = filter(b,a,x)

输出信号是通过使用该命令生成的系数y对输入信号进行滤波而获得的。xbutter

如果您想同时获得所有 3 个波段的输出，您可以将 3 个单独的滤波器操作的输出相加，每个波段一个。但是，如果您想保留波段之间的相位关系，您可能希望在执行此操作时使用该filtfilt命令。filter

如果您没有信号处理工具箱，您可以使用下面的代码创建二阶带通巴特沃斯系数，其中dt = 1/fs和是低截止频率fl和fu高截止频率。

function [ b, a ] = butterTwoBp( dt, fl, fu ) 
q=pi*dt*(fu-fl);
r=pi*dt*(fu+fl);
N = (tan(q)^2) + sqrt(2)*tan(q) + 1;
M = (tan(q)^2) / N; %M after N because it depends on N
O = -cos(r) * (2*sqrt(2)*tan(q) + 4) / ((cos(q))*N);
P = (-2*(tan(q)^2) + ((  (2*cos(r))   /  (cos(q))   )^2) + 2 )  /   N;
Q = cos(r)*(2*sqrt(2)*tan(q) - 4)/(cos(q)*N);
R = (   (tan(q)^2) - sqrt(2)*tan(q) + 1   )  /  N;

b=[M 0 -2*M 0 M];
a=[1 O P Q R];