How to Plot Absolute Frequency Response of Filter?

조회 수: 9 (최근 30일)
Matlab Student
Matlab Student 2016년 5월 15일
댓글: Star Strider 2016년 5월 16일
I am currently using the freqz() command to plot the frequency response of an FIR filter designed in fdatool and exported to the workspace as a variable.This plots the correct frequency response (as seen in FDAtool), but with normalized frequency on the x-axis.I understand that the freqz() command returns 512 points from 0 to pi.I divide the index by 512 and multiply it with FS/2, which "corrects" the x-axis but when I plot this against 20*log10(abs(filter)), the amplitude is not as seen in FDA tool.
mag=(freqz(Filter1,512));
plot(((1:512)/512)*(44100/2),20*log10(abs(mag(1:512))));

답변 (1개)

Star Strider
Star Strider 2016년 5월 15일
You can output the frequency in Hz if you give freqz the sampling frequency as an argument (at least in more recent MATLAB releases):
[mag,frq] = freqz(Filter1,2048,FS);
I’ve also noticed that fdatool and freqz give different results. I design my own filters with the individual filter functions, and analyse them with freqz. (I have nothing against using fdatool and the others, but when I learned signal processing, we had to design our own filters at every step, and I just got used to designing them that way.)
  댓글 수: 2
Matlab Student
Matlab Student 2016년 5월 16일
Thank you for your reply! I'd like to add that the correct frequency response can be obtained with
mag=(freqz(Filter1,512));
plot(((1:512)/512)*(44100/2),20*log10(abs(mag(1:512))));
Star Strider
Star Strider 2016년 5월 16일
My pleasure.
If you use freqz with no outputs, it plots the frequency response by default:
figure(1)
freqz(Filter1,2048,FS)

댓글을 달려면 로그인하십시오.

카테고리

Help CenterFile Exchange에서 Filter Design에 대해 자세히 알아보기

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by