MATLAB: FIR Parks-McClellan

fir parks-mcclellan

Hi,

I have a set of FIR Parks-McClellan coefficient :

Numerator: FR_SV = [3.91719822748777E-02, 0.103629842929331, 0.171922134825388, 0.221881476438683, 0.221881476438683, 0.171922134825388, 0.103629842929331, 3.91719822748777E-02];

Denominator: [1.000];

Can i use spectrogram function to plot and see the filter response? I have browse through the help file and could not find Parks-McClellan window. Anyone can help?

Regards, Newbie

Best Answer

Parks-McClellan is an algorithm for filter design. It does not use window. Spectrogram is normally used to view the signal since it explores both time and frequency domain. You have a set of fixed filter coefficients so there is no time dependence. You can certainly use spectrogram, but it won't give you additional information than freqz does.

Related Solutions

MATLAB: Blackman FIR bandpass filter design

Hi olga,

There is a subtle difference between the specification of an IIR filter and an FIR filter. Because you are designing an FIR filter, you need to use 'N,Fc1,Fc2' instead of 'N,F3dB1,F3dB2'. Here is an example

h = fdesign.bandpass
h.Specification = 'N,Fc1,Fc2'
Hd = window(h,'window',@blackman)
fvtool(Hd)

HTH,

Honglei

MATLAB: Convert Time Domain Signal Data into Frequency Domain, How to handle the imaginary terms

You want to plot the magnitude and phase separately for the complex-valued data. Think of the polar form of a complex number. A few other things: you want to create a frequency vector to use in plotting the frequency domain data (DFT). You may or may not want to center 0 frequency in your Fourier transform, I do this below. Because the mean of your time data is so large, you are going to get a large 0 frequency magnitude in your Fourier transform.

     value = [42007935 111212895 184546560 238219725 238219725 184546560 111212895   42007935];
     time = 0:7;
     stem(time,value,'markerfacecolor',[0 0 1])
     title('Time Data'); xlabel('Time'); ylabel('Amplitude');
     dftvalue = fft(value);
     freq = -pi:(2*pi)/length(value):pi-(2*pi)/length(value);
     figure;
     stem(freq,abs(fftshift(dftvalue)),'markerfacecolor',[0 0 1])
     title('Magnitude of DFT'); xlabel('Frequency'); ylabel('Magnitude');
     figure;
     stem(freq,angle(fftshift(dftvalue)),'markerfacecolor',[0 0 1])
     title('Phase of DFT');
     xlabel('Frequency'); ylabel('Magnitude');
     % or stem(freq,unwrap(angle(dftvalue)),'markerfacecolor',[0 0 1])

Best Answer

Related Solutions

MATLAB: Blackman FIR bandpass filter design

MATLAB: Convert Time Domain Signal Data into Frequency Domain, How to handle the imaginary terms

Related Question