How to apply an IIR-Filter correctly on a PCM-Signal?

Question

Hannes Broschk 2022-7-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1753805-how-to-apply-an-iir-filter-correctly-on-a-pcm-signal

回答： Jan 2022-7-5

I don't know much about filters yet and am just trying things out. I wanted to apply a "pink filter" to my measurement signal (maximum-length-sequence). I got the parameters for the IIR from this page: ccrma.stanford.edu. When I plot the filter, everything looks correct to me. However, when I apply the filter to the measurement signal, it is so quiet that you almost can't hear it anymore. In the original, the measurement signal is in full scale and jumps between 1 and -1. How can this be, what am I doing wrong?

Here's my code:

    In=readtable("idealMLS32768.txt"); %In is an Array filled with 32768 points
Error using readtable
Unable to find or open 'idealMLS32768.txt'. Check the path and filename or file permissions.
    Dirac = [1;zeros(32768-1,1)]; %Dirac could be used to demonstrate the filter instead of the Variable In
    %Fs=48000.0;
    
    
    b = [0.049922035 -0.095993537 0.050612699 -0.004408786];
    a = [1 -2.494956002   2.017265875  -0.522189400];
    iir = dsp.IIRFilter('Numerator',b,'Denominator',a);
    
    Out=iir(In.Var1); % here i want to apply the filter to the pcm-audio input.
    plot(Out)
    %fvtool(iir,'Fs',Fs);
     Y=fft(Out);
     P = sqrt(real(Y).^2+imag(Y).^2);
     P=20*(log(P)/log(10));
     plot(f,P(1:n/2+1)) 
     set(gca, 'XScale', 'log')
     title('FFT')
     xlabel('Frequency (f)')
     ylabel('|P(f)|')