如何在 MATLAB 中制作多陷波滤波器？

Question

我正在设计一个过滤器来消除一段音频中的噪音。我发现噪声频率位于频谱中的 745 和 1965 Hz，但我不知道如何制作多陷波滤波器来去除这两个特定频率

这是我的代码。我只能删除音频中的一个频率。有什么方法可以在 MATLAB 中对两个滤波器进行卷积吗？

%Reading first sample file
[x1,fs1] = audioread('sample.wav');

%Creating the time span for the file
t1=(0:length(x1)-1)/fs1;

%Creating the frequency span for the file
k1 = 0:length(x1)-1;
f1=k1*fs1/length(x1);
wg=[744.5*2/fs1 745.5*2/fs1 ];

%Creating  filter
[b1,a1] = butter(2,wg,'stop'); 

%Performing filtering on file
x1f = filter(b1,a1,x1);

Answer 1

正如您已经提到的，您必须对两个过滤器进行卷积以将它们组合起来。这可以使用 conv 函数来完成。

% Design first filter
wg1 = [744.5*2/fs1, 745.5*2/fs1];
[b1,a1] = butter(2,wg1,'stop');

% Design second filter
wg2 = [1964.5*2/fs1, 1965.5*2/fs1];
[b2,a2] = butter(2,wg2,'stop');

% Convolve filters
a0 = conv(a1, a2);
b0 = conv(b1, b2);

% Plot filter
freqz(b0, a0, 4096, fs1);