< 我可以提供代码给你，参考如下：

请参考一下：
函数besself用于设计贝塞尔模拟滤波器，该函数不适于数字滤波器设计。函数调用格式为
   [b,a]=besself(n,wn)      [b,a]=besself(n,wn,'ftype')    [z,p,k]=besself(.....)    [A,B,C,D]=besself(...)
其中,n为滤波器除数，wn为滤波器截止频率（-3dB）
函数impinvar是基于脉冲不变法实现模拟滤波器至数字滤波器的变换，调用格式为
   [bz,ba]=impinvar(b,a,Fs)
     [bz,ba]=impinvar(b,a)
其中b,a分别为模拟滤波器的分子和分母多项式系数向量；Fs为采样频率，HZ，缺省情况下Fs=1HZ；bz,az
分别为数字滤波器分子和分母多项式系数向量。
函数bilinear采用双线性变换法实现模拟s域至数字z域的映射，直接用于模拟滤波器变换为数字滤波器。
   其中零极点增益形式为
            [zd,Pd,Kd]=bilinear(z,p,k,Fs)
                        [zd,Pd,Kd]=bilinear(z,p,k,Fs,Fp)
式中，z,p分别为模拟滤波器零，极点列向量；k为模拟滤波器增益；Fs为采样频率，Hz.Fp为预畸变频率，Hz，函数选择项。zd,Pd,Kd为函数滤波器的零极点增益。