4-1AF滤波器设计.ppt

例：设计一个满足下列指标BW型模拟低通滤波器 p 1 . 0 = w p ， p 4 . 0 = w s ， dB A p 1 ? ， dB A s 10 ? 取N=2，将N=2带入满足通带的方程 解： (1) 计算N和?c 例：设计一个满足下列指标BW型模拟低通滤波器 p 1 . 0 = w p ， p 4 . 0 = w s ， dB A p 1 ? ， dB A s 10 ? 解： (2) 确定滤波器的极点 (3) 确定模拟低通滤波器的系统函数HL (s) 验证：Ap=1dB 满足指标；As= 18.3 dB 存在裕量 [num,den]=butter(N,Wc,s) 确定阶数为N、3db截频为Wc(radian/s)的Butterworth 滤波器分子和分母多项式。 s 表示模拟域。 [N,Wc]=buttord(Wp,Ws,Ap,As,s) 确定模拟Butterworth filter的阶数N和3db截频Wc。Wc是由阻带参数确定的。 MATLAB设计模拟Butterworth滤波器 [z,p,k]=buttap(N) 确定N阶归一化的Butterworth filter 的零点、极点和增益(gain) Wp=2*pi*1000;Ws=2*pi*5000;Ap=1;As=40; [N,Wc]=buttord(Wp,Ws,Ap,As,s); % filter order； fprintf(Order of the filter=%.0f\n,N) [num,den] = butter(N,Wc,s); fprintf(num=%.4e\n,num); fprintf(den=%.4e\n,den); %Compute Ap and As of designed filter omega=[Wp Ws]; h = freqs(num,den,omega); fprintf(Ap= %.4f\n,-20*log10(abs(h(1)))); fprintf(As= %.4f\n,-20*log10(abs(h(2)))); % Compute and plot the frequency response omega=[0:200:12000*pi];h= freqs(num,den,omega); gain = 20*log10(abs(h)); plot (omega/(2*pi),gain); xlabel(Frequency in Hz); ylabel(Gain in dB) 实现例子的MATLAB程序 Order of the filter=4 Numerator polynomial 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 9.7414e+015 Denominator polynomial 1.0000e+000 2.5961e+004 3.3698e+008 2.5623e+012 9.7414e+015 Ap= 0.1098，As= 40.0000 运行结果 4.1.3 Chebyshev(CB,Chebyshev I型)低通滤波器 e ：通带波纹 wc：通带截频 N：阶数，由阻带指标确定 选择等波纹特性的逼近函数 Chebyshev多项式 Chebyshev I 性质： w ) ( w j H 1 c w N=2 N=3 N=7 1)当 ， 2)当 ， 单调下降。N越大，下降越快。 Chebyshev I 性质： 3) ，e 控制了通带衰耗 归一化Cheybshev filter的极点 N为偶时 N为奇时 CB AF 设计步骤 (2)通带指标确定e (3)由阻带指标确定N(N越大阻带衰耗就越大) (1)通带截频确定wc=wp CB AF 设计步骤 (4)极点确定：取左半平面极点 k=1，2，~ ，2N 例：设计一个满足下列指标的模拟CB I型低通滤波器 推导过程 (1)求模拟滤波器的?c、e 取N=3 (2)求模拟滤波器的阶数N 推导过程(续) (3)确定极点，求系统函数 MATLAB设计模拟type I Chebyshev filter [z,p,k]=cheb1ap(N,Ap); 确定N阶归一化的Chebyshev filter 的零点、极点和增益(gain) [num,den]=cheby1(N,Ap,Wc,s) 确定阶数为N，通带截频为Wc(radian/s)的Chebyshev filter。s 表示模拟域 [N,Wc]=cheb1ord(Wp,Ws,Ap,As,s) 确定模拟Chebyshev filter的