模拟滤波器设计.ppt

归一低通到带通的变换为取+号时，无论w怎样变化，w均为正取-号时，无论w怎样变化，w均为负讨论取+号的情况：(2)-(1)得：(2)*(1)得：*第三章IIR数字滤波器的设计IIR滤波器设计的基本思想：1)将数字滤波器的设计为模拟滤波器的设计。2)设计满足技术指标的模拟滤波器。3)将模拟滤波器转换为数字滤波器。本章的主要内容1)模拟低通滤波器设计2)模拟域频率变换3)冲激响应不变法4)双线性变换法模拟滤波器的技术要求pw:通带截止频率wst:阻带截止频率dp:通带波动ds:阻带波动通带衰减(db)(passbandAttenuation)阻带衰减(db)(stopbandAttenuation)3.2模拟低通滤波器的设计|H(jw)|10通带过渡带阻带pwstwsdpd-1wA(w)=20log10|H(jw)|dB滤波器的衰减函数wc10N=1N=3N=50.7073.2.1ButterworthLowpassfilterN:滤波器阶数性质：2)幅度响应单调下降(monotonicallydecreasing)3)|H(jw)|2在w=0点1到2N-1阶导数零。称为最大平坦性。(maximallyflatmagnitudefilter)1)|H(j0)|=1,|H(j￥)|=0,-20log10|H(jwc)|?3dbwc:3db截频,当wc=1时，称其为归一化的BWF在w=0点做Taylorseries展开N增大通带衰减减小，阻带衰减加大。归一化的Butterworth滤波器(BWF)任意的BWF和归一化BWF的关系归一化的Butterworth滤波器的极点：条件：h(t)是实的H(jw)=H*(-jw)极点：共有2N个极点，为了保证系统的稳定，选左半平面的N个极点。为左半平面的N个极点当N为偶数时例：N=2,=p/4;k=1例：N=4,=p/8,3p/8;k=1,2当N为奇数时例：N=1N=3例：设计一个满足下列指标BW型模拟滤波器p1.0=wp，p4.0=wst，dbAp1=，dbAs10=取N=2，将N=2带入通带满足的方程通带满足指标，阻带超过指标验证：Ap=0.9999db;As=18.2795dbMATLAB设计模拟Butterworthfilter[z,p,k]=buttap(N)确定N阶归一化的Butterworthfilter的零点、极点和增益(gain)[num,den]=zp2tf(z,p,k)将零、极点和增益表示的H(s)转换为多项式系数表示的H(s)。[num,den]=butter(N,Wc,s)确定阶数为N，3-db截频为Wc(radian/s)的Butterworthfilter分子和分母多项式。s表示模拟域。[N,Wc]=buttord(Wp,Ws,Ap,As,s)确定模拟Butterworthfilter的阶数N和3-db截频Wc。Wc是由阻带参数确定的。例p3_1：设计满足下列条件的模拟Butterworth低通滤波器Wp=2*pi*1000;Ws=2*pi*5000;Ap=1;As=40;[N,Wc]=buttord(Wp,Ws,Ap,As,s);fprintf(Orderofthefilter=%.0f\n,N)[num,den]=butter(N,Wc,s);disp(Numeratorpolynomial);fprintf(%.4e\n,num);disp(Denominatorpolynomial);fprintf(%.4e\n,den);omega=[WpWs];h=freqs(num,den,omega);fprintf(Ap=%.4f\n,-20*log10(abs(h(1))));fprintf(As=%.4f\n,-20*log10(abs(h(2))));omega=[0:200:12000*pi];h=freqs(num,den,omega);gain=20*log10(abs(h));p