基于切比雪夫滤波器设计及matlab仿真

基于切比雪夫滤波器设计及matlab仿真内容摘要：

r polynomial coefficients % a = denominator polynomial coefficients % N = order of the elliptic filter % Rp = passband ropple in dB。 Rp 0 % Omegac = cutoff frequency in radians/sec % [z， p， k] = cheblap（ N， Rp）； % a = real(poly(p))。 % aNn = a(N+1)。 % p = p*Omegac。 % aNu = a(N+1)。 % k = k*aNu/aNn。 % b0 = k。 % B = real(poly(z))。 % b = k*B； 滤波器设计方程 已知p， s ，pR和 sA ，有三个参数就可以确定一个切比雪夫 I 型滤波器，根据相对线性标尺中的 21 与 22 式我们可得到： 1010 1pR 和 2020sAA 根据上面 以前讨论的性质有： c = p 和psr  阶 N 给出为 2 21Ag  210210l o g 1l o g 1ggNrr     4 滤波器 设计实现 滤波器设计 过程 设计一个低通的切比雪夫 I 型低通滤波器满足 ： 通带截止频率： p = ；通带波纹： pR =1dB 阻带 截止频率： *110 1 = ；阻带波纹： sA =16dB 我们首先计算必要的参数： *110 1 =， A= 162020 = c = p = , r =  = 2 21Ag  =,N=4 现在能确定 ()aHs，  = 1 + 211  = a = ( 1N N  ) = b = ( 1N N  ) = ()aHs有四个极点： =（ a c ） cos 28(b c ) sin 28 =  = （ a c ） 3cos 28 (b c ) 3sin 28 =  因此， 30()()a kkKHssp = 220 . 8 9 1 2 5 0 . 1 1 0 3 0 . 3 8 9 5( 0 . 1 7 5 4 0 . 3 8 9 5 ) ( 0 . 4 2 3 4 0 . 1 1 0 3 )s s s s    分子是要使得有： 21( 0) 1aHj   = MATLAB 的实现 已知滤波器的技术要求，利用 Uchblap 函数，我们提供了一个称为 afdchb1的函数用于设计一个模拟切比雪夫 I 型低通滤波器。 这个函数给出如后面的目录所示。 为了展示模拟滤波器的频域图，我们提供的一个称为 freqs_m的函数，它是有 MATLAB 提供的函数 freqs 的修正形式。 这个函数以绝对值和相对值 dB 尺度计算幅度响应以及相位响应，在 freqs_m 函数中响应是一 直要计算到最大频率max。 freqs_m函数见附录。 这个模拟滤波器的脉冲响应 ha(t)是用 MATLAB 的 impulde 函数计算出来的。 通过 MATLAB。