如何在这个matlab代码中找到H（z）

我想在这段代码中找到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]);

---

对于系数

b（1）

，

b（2）

，…，传递函数是k从1到len（b）的

b（k）*z^（-k+1）

之和。请参阅（z的幂与

+1

的区别是，Matlab数组的索引从1开始，但是链接中的向量

b

从

b_0

开始）

请注意，Matlab是区分大小写的。