系数向量；Fs为采样频率；bd为
数字滤波器的分子多项式系数；ad为数字滤波器的分母多项式系数。
②双线性变换法：bdadbili
earbaFs参数含义与impi
var一致。
2、ChebyshevⅠ型数字滤波器设计
ChebyshevⅠ型滤波器为通带纹波控制器：在通带呈现纹波特性，在阻带单调衰减。其格式
如下：
NOmegaccheb1ordOmegapOmegasPpAs‘s’
z0p0k0cheb1apNRp
参数含义与bueeordbuttap中参数一致。
3、ChebyshevⅡ型滤波器设计
ChebyshevⅡ型滤波器为阻带纹波滤波控制器：在阻带呈现纹波特性，在通带单调衰减。其
格式如下：
NOmegaccheb2ordOmegapOmegasRpAs‘s’
z0p0k0cheb2apNAs
其参数含义同上。
4、椭圆滤波器设计
椭圆滤波器在通带和阻带都呈现纹波特性。
NOmegacellipordOmegapOmegasRpAs‘s’
z0p0k0ellipapNRpAs
其参数含义同上。
三，实验内容与方法
【例341】利用脉冲响应不变法，用巴特沃斯原型设计一个数字低通滤波器，其中滤波器
的技术指标如下：wp04、Rp05dB、ws06、As50dB，滤波器采样频率为
1kHz。MATLAB程序如下：wp04piws06piRp05As50Fs1000OmegaPwpFsOmegaSwsFs
fepsqrt10Rp101Ripplesqrt11epepAtt
110As20NOmegaCbuttordOmegaPOmegaSRpAssz0p0k0buttapNpp0OmegaCzz0OmegaCkk0OmegaCNbakrealpolyzaarealpolypbdadimpi
varbaaaFsHwfreqzbdad1000wholeHH11501ww11501magabsHdb20log10magepsmaxmagphaa
gleHgrdgrpdelaybdadwsubplot221plotwpimagtitle幅度响应xlabelylabelHaxis01011setgcaXTickModema
ualXTick004061grido
setgcaYTickmodema
ualYTick0Att
Ripple1subplot223plotwpidbtitle幅度（单位：dB）xlabel频率（单位：pi）ylabel分贝axis0115050setgcaXTickModema
ualXTick004061setgcaYTickmodema
ualYTick15050050grido
subplot222plotwpiphapititle相位响应xlabelylabel单位：piaxis0111setgcaXTickModema
ualXTick004061setgcaYTickmodema
ualYTick101grido
subplot224plotwpigrdtitle群延迟xlabel频率（单位：pi）ylabel样本axis01020setgcaXTickModema
ualXTick004061setgcaYTickmodema
ualYTick0220grido
setgcfcolorw程序运行结果如图所示。
f【例342】按照例341的技术指标，即wp04、Rp05dB、ws06、As50dB，
采样频率为1kHz。利用脉冲响应不变法，要求用切比雪夫Ⅰ型来实现数字低通滤波器。MATLAB程序如下：
wp04piws06piRp05As50r