长度为6的滑动平均滤波器的增益响应
clf;
M = 6;%设定滤波器的长度
num = ones(1,M)/M;
[g,w] = gain(num,1);%求增益函数
plot(w/pi,g);grid
axis([0 1 -50 0.5])
xlabel('\omega /\pi');ylabel('Gain in dB');
title(['M = ', num2str(M)])
gain函数:
function [g,w] = gain(num,den)
w = 0:pi/255:pi;
h = freqz(num,den,w);%求频率响应
g = 20*log10(abs(h));%求频率响应h的增益响应