【正文】
Rs=。dev=[Rs Rp Rs]。 [N,fo,ao,w] = remezord(f,a,dev)。 h = remez(N,fo,ao,w)。 w=linspace(0,pi,1000)。 mag=freqz(h,[1],w)。 plot(w/pi,20*log10(abs(mag)))。 xlabel(39。Normalized frequency39。)。 ylabel(39。Gain, dB39。)。 用 Remez算法實(shí)現(xiàn)的 MATLAB程序: 0 1 ?40 ?20 0 Normalized frequency G a i n , d B M=25時(shí)( remezord 確定的濾波器階數(shù))濾波器的增益響應(yīng)。 上例的設(shè)計(jì)結(jié)果 M=25 ?Ws1=,Wp1=,Wp2=, Ws1=, ?p=,?s=) As??20log10(?s)=40dB 阻帶衰減 不滿足指標(biāo) 上例的設(shè)計(jì)結(jié)果 M=30 (Ws1=, Wp1=, Wp2=, Ws1=, ?p=, ?s=) 0 1 ?40 ?20 0 Normalized frequency G a i n , d B M=30,由 remez得出的濾波器的增益響應(yīng) As??20log10(?s)=40dB 阻帶衰減 滿足指標(biāo)