神经网络第二章.ppt

bilinear.m：双线性变换，由模拟滤波器01得到数字滤波器。02[Bz,Az]=bilinear(B,A,Fs)03式中B,A分别是G（s）的分子、分母多项式04的系数向量，Bz,Az分别是H(z)的分子、分05母多项式的系数向量，Fs是抽样频率。06butter.m用来直接设计Butterworth数字滤波器，实际上它把buttord.m,buttap.m,lp2lp.m,bilinear.m等文件都包含了进去，从而使设计过程更简捷。格式(1)～(3)用来设计数字滤波器，B，A分别是H(z)的分子、分母多项式的系数向量，Wn是通带截止频率，范围在0～1之间。若Wn是标量，(1)用来设计低通数字滤波器，若Wn是1×2的向量，则(1)用来设计数字带通滤波器；(2)用来设计数字高通滤波器；(3)用来设计数字带阻滤波器，显然，这时的Wn是1×2的向量；格式(4)用来设计模拟滤波器。[B,A]=butter(N,Wn);(2)[B,A]=butter(N,Wn,’high’);[B,A]=butter(N,Wn,’stop’);(4)[B,A]=butter(N,Wn,’s’)例6.7.1(例6.5.1)clearall;fp=100;fs=300;Fs=1000;rp=3;rs=20;wp=2*pi*fp/Fs;ws=2*pi*fs/Fs;Fs=Fs/Fs;%letFs=1%Firstlytofinishfrequencyprewarping;wap=tan(wp/2);was=tan(ws/2);%[n,wn]=buttord(wap,was,rp,rs,s)%Note:s![z,p,k]=buttap(n);%[bp,ap]=zp2tf(z,p,k)%[bs,as]=lp2lp(bp,ap,wap)%%Note:s=(2/Ts)(z-1)/(z+1);Ts=1,thatis2Fs=1,Fs=0.5;[bz,az]=bilinear(bs,as,Fs/2)%[h,w]=freqz(bz,az,256,Fs*1000);plot(w,abs(h));gridon;设计IIRLPDF，例6.7.1(例6.5.1)clearall;wp=.2*pi;ws=.6*pi;Fs=1000;rp=3;rs=20;%%Firstlytofinishfrequencyprewarping;wap=2*Fs*tan(wp/2);was=2*Fs*tan(ws/2);[n,wn]=buttord(wap,was,rp,rs,s);%Note:s![z,p,k]=buttap(n);[bp,ap]=zp2tf(z,p,k);[bs,as]=lp2lp(bp,ap,wap)w1=[0:499]*2*pi;h1=freqs(bs,as,w1);[bz,az]=bilinear(bs,as,Fs)%Note:z=(2/ts)(z-1)/(z+1);[h2,w2]=freqz(bz,az,500,Fs);plot(w1/2/pi,abs(h1),w2,abs(h2),k);gridon;设计IIRLPDF，例6.7.1(例6.5.1)clearall;wp=.2*pi;ws=.6*pi;Fs=1000;rp=3;rs=20;[n,wn]=buttord(wp/pi,ws/pi,rp,rs);[bz,az]=butter(n,wp/pi)[bz1,az1]=butter(n,wn)[h,w]=freqz(bz,az,128,Fs);[h1,w1]=freqz(bz1,az1,128,Fs);plot(w,abs(h),w1,abs(h1),g.);gridon;设计IIRLPDF，非线性关系设计的AF并不是按给定的技术指标，但当再由变回后，保证了DF的技术要求。又称为频率的预变形（Freq.Warping)。例如：抽样频率设计出7.6数字高通,带通及带阻滤波器的设计得到模拟低通的技术要求得到模拟高通转移函数得到模拟高通的技术要求数字高通滤波器设计步骤给出数字高通的技术要求最后得到数字高通转移函数对带通（BP）、带阻（BS）数字滤波器的设计，只需改变图中Step2和Step4：带阻带通要求：按上述转换办法，可以求出：例6.6.2：设计一IIRBPDF