基于切比雪夫?yàn)V波器設(shè)計(jì)及matlab仿真(文件)
【正文】 MATLAB 提供了 [Z， P， K]=cheblap（ N， Rp）函數(shù)來設(shè)計(jì) N 階 通帶波紋為 Rp 的歸一化切比雪夫 I 型模擬原型濾波器 ， 并得到數(shù)組 Z 和 P（即零點(diǎn)和極點(diǎn)）以及增益 值 K。 新的增益 K 利用 式子： ( 0)aHj=1， N為奇數(shù)或21( 0) 1aHj ?? ?， N 為偶數(shù)確定這可通過將原來的 K 倍乘以非歸一化分母多項(xiàng)式對歸一化分母多項(xiàng)式的比在 s=0 的求值來完成。 Rp 0 % Omegac = cutoff frequency in radians/sec % [z， p， k] = cheblap（ N， Rp）； % a = real(poly(p))。 % k = k*aNu/aNn。這個函數(shù)給出如后面的目錄所示。 這個模擬濾波器的脈沖響應(yīng) ha(t)是用 MATLAB 的 impulde 函數(shù)計(jì)算出來的。通常我們用濾波器是來降低噪聲，估計(jì)信號的，不過濾波器也可以用來處理一些其它的問題，如圖象壓縮等。這三種頻域?yàn)V波器設(shè)計(jì)的原理是根據(jù)噪聲與有用信號的頻率分布是基本分開的，通過確定一個頻率開關(guān)值，將噪聲的頻率去掉，再進(jìn)行傅立葉反變換就得到有用信號?；具_(dá)到了了解 信號處理課程設(shè)計(jì) 的任務(wù)，明確了 濾波器設(shè)計(jì) 的基本原則，掌握了 濾波器設(shè)計(jì)的基本方法與任務(wù)。Ws=*pi。Attn=10^(As/20)。 %Plots subplot(2,2,1)。) xlabel(39。|H|39。,39。,[0,])。manual39。grid subplot(2,2,2)。) xlabel(39。分貝數(shù) 39。,39。,[0,])。manual39。grid set(gca,39。,39。,39。]) subplot(2,2,3)。) xlabel(39?；《?39。,39。,[0,])。manual39。grid subplot(2,2,4)。) xlabel(39。ha(t)39。 OmegaC = Wp。 Fprintf(‘ \n***切比雪夫 I型濾波器的階次 =%\n’ ,N) [b,a] = cheblap_o(N,Rp,Rp,Rs) 。 g = sqrt(A*A1)/ep。axis([0,max(t),min(ha),max(ha)]) Function[b,a]=afd_cheb1(Wp,Ws,Rp,Rs) %切比雪夫 I型模擬低通濾波器的設(shè)計(jì) %b為 ()aHs分子多項(xiàng)式的系數(shù) %a為 ()aHs分母多項(xiàng)式的系數(shù) %Wp 為以弧度 /秒為單位的通帶邊緣頻率； Wp0 %Ws 為以弧度 /秒為單位的阻帶邊緣頻率； WsWp0 %Rp 為通帶中的振幅波動的 dB數(shù) %As 為阻帶衰減的 dB數(shù) %if Wp=0 Error (‘ 通帶必須大于 0’ ) end if Ws=Wp Error (‘ 阻帶邊緣必須大于通帶邊緣 ’ ) end if (Rp=0)|(As0) Error (‘ 通帶 波紋或阻 帶衰減必須大于 0’ ) end %ep = sqrt(10^(Rp/10)1)。)。title(39。YTick39。YTickmode39。,39。axis([0,1,1]) set(gca,39。)。title(39。,39。,[30,39。,39。YTick39。YTickmode39。,39。axis([0,30,5]) set(gca,39。)。title(39。YTick39。YTickmode39。,39。axis([0,0,]) set(gca,39。)。title(39。
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號-1