離散信号系统分析.docVIP

离散信号与离散系统分析 一、实验目的 1、掌握双线性变换法及脉冲响应不变法设计IIR数字滤波器的具体设计方法及原理，熟悉用双线性变换法及脉冲响应不变法设计低通高通带通IIR数字滤波器的编程。 2、观察双线性变换及脉冲响应不变法设计的滤波器的频域特性，了解双线性变换法及脉冲响应不变法的特点。 3、熟悉巴特沃思滤波器、切比雪夫滤波器和椭圆滤波器的频率特性。 二、实验原理 [N, Wn] = buttord(Wp, Ws, Rp, Rs, s): 求巴特沃思滤波器的阶数N和3dB边界角频率Wn。Rp和Rs分别为通带波动δ和阻带衰减At，单位均为dB；s表示模拟滤波器设计，如无此参数，则表示数字滤波器设计；Wp和Ws分别表示通带边界角频率Ωc和阻带边界角频率Ωr，其值为标量(低通和高通)或双元素向量(带通和带阻)，见下表。 Filter Type Wp Ws Passband Stopband Lowpass WpWs (0, Wp) (Ws, ∞) Highpass WpWs (Wp, ∞) (0, Ws) Bandpass 0Ws(1)Wp(1)Wp(2)Ws(2) (Wp(1), Wp(2)) (0, Ws(1)) (Ws(2), ∞) Bandstop 0Wp(1)Ws(1)Ws(2)Wp(2) (0, Wp(1)) (Wp(2), ∞) (Ws(1), Ws(2)) [b, a] = butter(N, Wn, ftype, s): 根据N和Wn求巴特沃思滤波器的系统函数H(s)，b和a分别为H(s)分子和分母多项式的系数。‘ftype’指定滤波器类型，值为‘low’、‘high’或‘stop’， ‘low’为缺省值，表示设计低通或带通滤波器(取决于Wn是标量还是双元素向量)，‘high’或‘stop’分别表示设计高通和带阻滤波器。低通和高通滤波器的阶数为N，带通和带阻滤波器的阶数为2N。 [z, p, k] = butter(N, Wn, ftype, s): 根据N和Wn求巴特沃思滤波器的零、极点和增益。z、p和k分别为H(s)的零点向量、极点向量和增益。系统函数H(s)与b和a或z、p和k的关系为： [N, Wp] = cheb1ord(Wp, Ws, Rp, Rs, s): 求切比雪夫I型滤波器的阶数N。输入参数的含义与巴特沃思滤波器相同，输出参数Wp等于输入参数Wp。 [b, a] = cheby1(N, Rp, Wp, ftype, s): 求切比雪夫I型滤波器的系统函数H(s)，b和a分别为H(s)分子和分母多项式的系数。 [z, p, k] = cheby1(N, Rp, Wp, ftype, s): 求切比雪夫I型滤波器的零、极点和增益。 [N, Wp] = ellipord(Wp, Ws, Rp, Rs,s): 求椭圆滤波器的阶数N。 [b, a] = ellip(N, Rp, Rs, Wp, ftype, s): 求椭圆滤波器的系统函数H(s)，b和a分别为H(s)分子和分母多项式的系数。 [z, p, k] = ellip(N, Rp, Rs, Wp, ftype, s): 求椭圆滤波器的零、极点和增益。 [bz, az] = impinvar(b, a, fs): 通过脉冲响应不变法求数字滤波器的系统函数H(z)。b和a分别为模拟滤波器系统函数H(s)的分子和分母多项式系数；fs为采样频率；bz和az分别为H(z)的分子和分母多项式系数。 [bz, az] = bilinear(b, a, fs): 通过双线性变换法求数字滤波器的系统函数H(z)。输入、输出参数的含义与impinvar函数相同。系统函数H(z)与bz和az的关系为： 三、实验内容 (1)，，，,;设计一切比雪夫高通滤波器，观察其通带损耗和阻带衰减是否满足要求 （2），，，,;分别用脉冲响应不变法及双线性变换法设计一巴特沃思数字低通滤波器，观察所设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。比较这两种方法的优缺点。 （3）利用双线性变换法分别设计满足下列指标的巴特沃思型、切比雪夫型和椭圆型数字低通滤波器，并作图验证设计结果：，，，,。 （4）分别用脉冲响应不变法及双线性变换法设计一巴特沃思数字带通滤波器，已知，其等效的模拟滤波器指标为，，，，，。 （5）利用双线性变换法设计满足下列指标的切比雪夫型数字带阻滤波器，并作图验证设计当时，，当以及时，；采样频率 四、实验结果与分析 第一题： bz =[0.0262 -0.1047 0.1570 -0.1047 0.0262] az =[1.0000 1.5289 1.6537 0.94