【正文】
? [b,a]=ellip(n,Rp,Rs,Wn) ? [b,a]=ellip(n,Rp,Rs,Wn,’ftype’) ? [b,a]=ellip(n,Rp,Rs,Wn,’s’) ? [b,a]=ellip(n,Rp,Rs,Wn,’ftype’,’s’) ? [z,p,k]=ellip(…) ? [A,B,C,D]=ellip(…) ? 參數(shù)含義: n-濾波器階 數(shù) Rp-通帶波紋系數(shù) (dB) Rs-阻帶波紋系數(shù) (dB) Wn-歸一化截止頻率 (若為二元向量則給出 2n階數(shù)字帶通濾波器) ftype=high-高通數(shù)字濾波器 ftype=stop,Wn為二元向量-帶阻數(shù)字濾波器 ? 返回變量 b,a-濾波器分子分母多項式系數(shù) z,p,k-濾波器零極點及增益 A,B,C,D-濾波器的狀態(tài)空間型描述參數(shù) ellipord: 切比雪夫 I型濾波器的階數(shù)計算函數(shù) ? [n,Wn]=ellipord(Wp,Ws,Rp,Rs) ? [n,Wn]=ellipord(Wp,Ws,Rp,Rs,’s’) ? 參數(shù)含義: Wp-通帶轉(zhuǎn)折頻率 Ws-阻帶轉(zhuǎn)折頻率 Rp-通帶波紋系數(shù) Rs-阻帶衰減系數(shù) 說明 ? MATLAB設(shè)計函數(shù)的規(guī)范化頻率設(shè)置不是通常的歸一化概念! ? 當(dāng) 設(shè)計的濾波器轉(zhuǎn)折頻率為 Wp(rad)和Ws(rad)時, 它們 必須被 π除 。(*fs/2π=*fNyq/π) ? 例如:濾波器技術(shù)指標(biāo)要求當(dāng)抽樣頻率為5kHz時,系統(tǒng)截止頻率為 1kHz時,在MATLAB設(shè)計函數(shù)中 必須 以 : ωc=2π(1000/5000)=2π()= 典型 IIR濾波器的比較 類型 特點 |H(jΩ)|(Ω=1) butter Ω=0和 Ω=∞處提供最好的 Taylor級數(shù)逼近,整體平滑 單調(diào) cheby1 理想與實際頻響絕對誤差最小,阻帶響應(yīng)最平坦 10Rp/20 cheby2 優(yōu)點:理想 與實際頻響絕對誤差最小,通帶響應(yīng)無波紋最 平坦 缺點:阻帶逼近零比 cheby1慢,且階次 n為偶數(shù)時不逼近零 10Rs/20 ellip 滿足相同要求所需階次最低,相同 指標(biāo)( n,Rp,Rs)過渡帶 最窄 10Rp/20 yulewalk: 基于最小均方意義的遞歸濾波器設(shè)計函數(shù) ? [b,a]=yulewalk(n,f,m) ? designs recursive IIR digital filters using a leastsquares fit to a specified frequency response. ? f is a vector of frequency points, specified in the range between 0 and 1, where 1 corresponds to half the sample frequency (the Nyquist frequency). The first point of f must be 0 and the last point 1, with all intermediate points in increasing order. Duplicate frequency points are allowed, corresponding to steps in the frequency response. ? m is a vector containing the desired magnitude response at the points specified in f.(給定頻率點的理想幅值響應(yīng) ) ? f and m must be