Designing an Arbitrary FIR Filter

Question

Krishna Ganesan 2020-2-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/505945-designing-an-arbitrary-fir-filter

评论： Star Strider 2020-2-19

I need to design an Arbitrary FIR Filter whose response looks like this,

Frequency = [100 125 160 200 250 315 400 500 630 800 1000 1250 1600 2000 2500 3150 4000 5000 6300 8000 10000]

Required_Filter_Response = [55,38 55,74 52,13 52,64 51,25 49,04 46,50 44,99 42,95 40,73 39,47 37,53 35,36 34,33 32,33 30,36 28,01 25,98 23,19 20,42 18,65]

My professor suggested me using designfilt and fvtool for realizing it. I am not new to MATLAB, but to filters and digital signal processing, I am.

Any suggestions on how to achieve it would be greatly appreciated. Thanks!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2020-2-17

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/505945-designing-an-arbitrary-fir-filter#answer_415982

在 MATLAB Online 中打开

This is probably as close as you can get:

Frequency = [100  125  160  200  250  315  400  500  630  800  1000  1250  1600  2000  2500  3150  4000  5000  6300  8000  10000];
Required_Filter_Response = [55.38  55.74  52.13  52.64  51.25  49.04  46.50  44.99  42.95  40.73  39.47  37.53  35.36  34.33  32.33  30.36  28.01  25.98  23.19  20.42 18.65];
Fs = 20000;                                                                 % Create Sampling Frequency (Not Supplied, Must Be > 20000 Here)
Fn = Fs/2;
Fnew = linspace(min(Frequency), max(Frequency), fix(numel(Frequency)/2)*2);
Rnew = interp1(Frequency, Required_Filter_Response, Fnew);
b = firls(100, Fnew/Fn, Rnew/max(Rnew));
figure
freqz(b, 1, 2^14, Fs)
[h,f] = freqz(b, 1, 2^14, Fs);
figure
plot(Fnew, 20*log10(Rnew/max(Rnew)))
hold on
plot(f, 20*log10(abs(h)))
hold off

Note that the filter and the desired response are essentially the same in the second plot.

6 个评论
显示 4更早的评论隐藏 4更早的评论

Star Strider 2020-2-17

编辑：Star Strider 2020-2-17

在 MATLAB Online 中打开

My pleasure.

It is always good to know the sampling frequency annd the intended application.

That problem is then relatively straightforward to solve. It helps to estimate the parameters of mathematical model of the frrequency response, in order to see what the correct extrapolation values would be, out to the Nyquist frequency. This appears to be a reasonable model:

fcn = @(b,f) b(1).*f.^b(2) + b(3);
B = fminsearch(@(b) norm(Required_Filter_Response - fcn(b,Frequency)), [25; -0.15; 10])

producing

B =
          345.886358298793
       -0.0285538321269912
         -245.322882425467

So evaluating the function at the Nyquist frequency (22050 Hz) and adding those to the original vectors:

Frequency = [100  125  160  200  250  315  400  500  630  800  1000  1250  1600  2000  2500  3150  4000  5000  6300  8000  10000 22050];
Required_Filter_Response = [55.38  55.74  52.13  52.64  51.25  49.04  46.50  44.99  42.95  40.73  39.47  37.53  35.36  34.33  32.33  30.36  28.01  25.98  23.19  20.42 18.65 14.57];
Fs = 44100;                                                                 % Sampling Frequency
Fn = Fs/2;
Fnew = linspace(min(Frequency), max(Frequency), fix(numel(Frequency)/2)*2);
Rnew = interp1(Frequency, Required_Filter_Response, Fnew);
b = firls(60, Fnew/Fn, Rnew/max(Rnew));
figure
freqz(b, 1, 2^14, Fs)
[h,f] = freqz(b, 1, 2^14, Fs);
figure
plot(Fnew, 20*log10(Rnew/max(Rnew)))
hold on
plot(f, 20*log10(abs(h)))
hold off

So with the added frequency point and a slightly different (and more efficient) filter order, the new transfer function appears to almost exactly match ‘Required_Filter_Response’. It appears to be a linear-phase filter, however I would still recommend using it with filtfilt rather than filter.

EDIT — (17 Feb 2020 at 17:25)

Added Bode plot figure —

and code creating it —

[h,f] = freqz(b, 1, 2^14, Fs);
figure
subplot(2,1,1)
plot(Fnew, 20*log10(Rnew/max(Rnew)),'.r')
hold on
plot(f, 20*log10(abs(h)),'-b')
hold off
xlabel('Frequency (Hz)')
ylabel('Magnitude (dB)')
xlim([min(f) max(f)])
legend('Desired Filter Response','Actual Filter Response')
grid
subplot(2,1,2)
plot(f, rad2deg(unwrap(angle(h))))
xlabel('Frequency (Hz)')
ylabel('Phase (degrees)')
xlim([min(f) max(f)])
grid

This is simply an augmented freqz plot. It is not part of the filter design code.

Krishna Ganesan 2020-2-18

编辑：Krishna Ganesan 2020-2-18

在 MATLAB Online 中打开

Dear Star Strider,

Thank you very much for providing the code for the sampling frequency that I wanted. Still, after running and studying the code, I have the following doubts.

--------------------------------------------------------------------------------------------------------------------------------------------------

Doubt 1 - The filter response that I need is from 55 dB at 100 Hz to 20 dB at 10,000 Hz. I understand that the frequency is normalized using,

plot(Fnew, 20*log10(Rnew/max(Rnew)))

and I understand that normalized frequencies are always adjusted from 0 to 1. But, when I run the code, the filter response is from 0 dB at 0 Hz to -12 dB at 20,000 Hz as presented below,

What my doubt is, should some sort of scaling be done? Or any gain be added to boost the filter response to make it to the scaling (from 55 dB at 100 Hz to 20 dB at 10,000 Hz) that I want?

--------------------------------------------------------------------------------------------------------------------------------------------------

Doubt 2 - I have an audio file with two channels on which the filter has to be applied. I read in the audio file using,

[y,Fs]    = audioread('Input.wav'); %Audio measurement recorded using Binaural headset microphones
Left_Ear  = y(:,1);
Right_Ear = y(:,2);

My doubt here is, how do I apply the filter on the audio file?

I understand that the filter is stored in the variable 'b'

b = firls(60, Fnew/Fn, Rnew/max(Rnew));

Would this be a correct approach to apply the filter on the audio file? As my professor wanted me to apply the filter in time-domain.

Filtered_Left_Ear  = conv(Left_Ear,b);
Filtered_Right_Ear = conv(Right_Ear,b);

I understand that the filter has to be applied to each channels separately and then they have to be appended together to save the filtered audio file.

Filtered_Audio_File = [Filtered_Left_Ear, Filtered_Right_Ear];
audiowrite('Output.wav', Filtered_Audio_File, Fs);

--------------------------------------------------------------------------------------------------------------------------------------------------

Kindly help me solve the issue

Krishna Ganesan 2020-2-19

Thank you very much for the clarification!

Star Strider 2020-2-19

My pleasure.

If my Answer helped you solve your problem, please Accept it!

请先登录，再进行评论。

Designing an Arbitrary FIR Filter

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

6 个评论
显示 4更早的评论隐藏 4更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

Designing an Arbitrary FIR Filter

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

6 个评论 显示 4更早的评论隐藏 4更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

WeChat

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

6 个评论
显示 4更早的评论隐藏 4更早的评论