为我讲解一下代码

好的，我可以对这个代码段逐行进行解释：

wc=2*1000*tan(2*pi*400/(2*1000)); 
wt=2*1000*tan(2*pi*317/(2*1000));

这两行代码定义了两个截止频率，`wc`是滤波器的通带截止频率，`wt`是滤波器的阻带截止频率。其中`2*1000`代表采样率，`tan`是正切函数，`2*pi*400/(2*1000)`代表通带的中心频率，`2*pi*317/(2*1000)`代表阻带的中心频率。

[N,wn]=cheblord(wc,wt,0.5,19,'s');

这行代码使用Chebyshev Type I滤波器设计方法计算出了所需的滤波器参数，其中`N`是滤波器的阶数，`wn`是滤波器的截止频率。

[B,A]=cheby1(N,0.5,wn,'high','s');

这行代码使用Chebyshev Type I滤波器设计方法生成了滤波器的分子项和分母项系数，`B`是分子项系数，`A`是分母项系数。

[num,den]=bilinear(B,A,1000);

这行代码使用双线性变换将连续时间域中的滤波器转换为离散时间域中的数字滤波器，`num`是数字滤波器的分子项系数，`den`是数字滤波器的分母项系数。

[h,w]=freqz(num,den);

这行代码计算数字滤波器的幅频响应，`h`是滤波器的幅度响应，`w`是对应的频率。

f=w/pi*500; 
plot(f,20*log10(abs(h))); 
axis([0,500,- 80,10]); 
grid; 
xlabel('频率/Hz') 
ylabel('幅度/dB')

这几行代码将计算出的数字滤波器的幅频响应进行绘制，其中`f`是频率，`20*log10(abs(h))`是滤波器的幅度响应转换为分贝单位后的值。`axis`函数设置图像的坐标轴范围，`grid`函数添加网格线，`xlabel`和`ylabel`分别用于添加X轴和Y轴的标签。最终绘制出的图像显示了数字滤波器的幅频响应。