寻求振动数据快速带通滤波方案：bandpass过慢、fftfilt失效

我需要对振动数据（变量为accH）进行带通滤波：

• 采样率48kHz，需提取40个连续频段（如20-600Hz、600-1200Hz…）
• 用bandpass()速度太慢，无法满足需求
• 用fftfilt()时滤波后波形振幅为0，但频谱对应频段有信号

现有代码 #

accH=randn(48000*5,1);% 用randn模拟振动数据

%% 慢的bandpass滤波
LB=0:600:24000-600;LB(1)=20;
UB=600:600:24000;UB(end)=24000-20;
A1=zeros(length(accH),length(LB));
 for k=1:length(UB)
 A1(:,k)=bandpass(accH,[LB(k) UB(k)],Fs);
end

%% 无效的fftfilt滤波
A=zeros(length(accH),length(LB));
for j=1:length(UB)
 bpfilt = designfilt('bandpassfir', ...
 'FilterOrder',20,'CutoffFrequency1',LB(k), ...
 'CutoffFrequency2',UB(k),'SampleRate',Fs);
 A(:,k)=fftfilt(bpfilt,accH);
end

问题现象 #

对比两种方法的结果：6000-6600Hz频段中，bandpass()输出的波形正常，fftfilt()输出的波形振幅为0，明显异常。需要无需陡峭过渡带的快速滤波方案。

解决方案 #

1. 修复fftfilt的直接错误

你的fftfilt代码存在变量名误用：循环变量是j，但调用了LB(k)、UB(k)，这是导致输出异常的核心原因。修正后fftfilt可正常工作，且速度远快于bandpass()。

修正后的fftfilt代码：

A=zeros(length(accH),length(LB));
Fs = 48000; % 补充定义采样率
for j=1:length(UB)
 bpfilt = designfilt('bandpassfir', ...
 'FilterOrder',20,'CutoffFrequency1',LB(j), ...
 'CutoffFrequency2',UB(j),'SampleRate',Fs);
 A(:,j)=fftfilt(bpfilt,accH);
end

2. 频域一次性滤波（更快方案）

由于你的频段连续且覆盖0~24kHz（除首尾边缘），可以直接通过FFT将信号转到频域，给每个频段生成掩码后逆FFT得到结果，速度远优于循环FIR滤波。

代码示例：

Fs = 48000;
N = length(accH);
f = Fs/2 * linspace(0,1,N/2+1);

% 对原始信号做一次FFT
accH_fft = fft(accH);

A2 = zeros(N, length(LB));
for k=1:length(LB)
    % 生成频域掩码
    mask = zeros(size(accH_fft));
    % 定位目标频段的频率索引
    idx_low = find(f >= LB(k), 1, 'first');
    idx_high = find(f <= UB(k), 1, 'last');
    % 给正频率区间赋值
    mask(idx_low:idx_high) = 1;
    % 给负频率区间（对称部分）赋值
    if idx_high < N/2+1
        mask(N - idx_high + 2 : N - idx_low + 2) = 1;
    end
    % 逆FFT得到滤波后信号，取实部消除数值误差
    A2(:,k) = real(ifft(accH_fft .* mask));
end

该方案优势：

• 仅做一次FFT，循环仅处理掩码和逆FFT，计算效率大幅提升
• 无需陡峭过渡带时，直接用矩形掩码即可，完全匹配需求

3. 额外优化：减少重复计算

若频段为均匀划分（仅首尾特殊），可预先生成所有频段的索引范围，避免循环中重复查找频率索引，进一步加快速度。

内容的提问来源于stack exchange，提问作者Sebastian Carta