我想在此代码中找到H(z)并打印出来。但我不知道我该怎么办?
n=-11:11;
lp=(1/3)*sinc((1/3)*n);
wh=rectwin(23);
b=lp.*wh';
figure(1);
k=0:22;
stem(k,b);
title('N=23 FIR Filter impulse Response');
xlabel('Time');
ylabel('Mag');
figure(2);
[h,w]=freqz(b,1,1024);
plot(w/pi,20*log10(abs(h)));
grid;
title('LPF FIR -rectwin N=23 frequency response');
axis([0 1 -100 10]);
答案 0 :(得分:1)
对于具有系数b(1)
,b(2)
,...的FIR滤波器,传递函数是k从1到len(b)的b(k)*z^(-k+1)
之和。请参阅https://ccrma.stanford.edu/~jos/fp/FIR_Transfer_Function.html(z的幂+1
的差异是因为Matlab数组从1开始编制索引,但链接中的向量b
从b_0
开始。)