用matlab设计陷波器的程序中，NLen是什么意思

Question

画陷波器的幅频特性时，已知Fs=4000Hz,f0=500Hz,陷波器的传递函数形式为H(Z)=(1+a/z+1/z^2)/(1+ab/z+(b/z)^2),怎么用matlab编写程序，画出幅频特性曲线。