高頻課程設(shè)計(jì)報(bào)告--基于matlab的帶通濾波器設(shè)計(jì)(編輯修改稿)
2025-04-19 12:51 本頁面
【文章內(nèi)容簡介】 最小衰減。凱瑟給出的估算β和濾波器階數(shù)N的公式如下： 式中，=是數(shù)字濾波器過渡帶寬度。應(yīng)當(dāng)注意，因?yàn)樯鲜綖殡A數(shù)估算，所以必須對設(shè)計(jì)結(jié)果進(jìn)行檢驗(yàn)。 ： 根據(jù)題目要求。因此，我設(shè)置的阻帶和通帶截止頻率分別為wsl=,wpl=,wph=,wsh=,Rs=60dB,故過渡帶寬為DB=wplwsl=,（）。由奈奎斯特抽樣定理得，fs=2*320=640(Hz)，這里為了得到更好的抽樣效果，選取fs=1000Hz。 matlab程序： %kaiser濾波器%wsl= *pi。 %臨界頻率采用模擬角頻率表示 wsh= *pi。 %臨界頻率采用模擬角頻率表示wpl= *pi。 %臨界頻率采用模擬角頻率表示 wph=*pi。 %臨界頻率采用模擬角頻率表示rs=60。 %阻帶最小衰減DB=wplwsl。 %過渡帶寬B=*()。 %計(jì)算凱瑟窗控制參數(shù)M=ceil((rs8)/)。%計(jì)算所需h(n)長度M，ceil(x)取大于等于x的最小整數(shù)wc=[(wpl+wsl)/2/pi,(wph+wsh)/2/pi]。%計(jì)算理想帶通濾波器通帶截止頻率（關(guān)于π歸一化）h=fir1(M1,wc,kaiser(M,B))。 %調(diào)用fir1計(jì)算帶通FIRDF的h(n)Fs=1000。 %抽樣頻率[H,w1]=freqz(h,1,M,Fs)。 %求濾波器幅度響應(yīng)，設(shè)置最大幅度為1plot(w1,abs(H))。 %畫出濾波器幅度響應(yīng)grid。xlabel(39。頻率/Hz39。)。ylabel(39。幅值39。)。title(39。濾波器幅頻響應(yīng)39。)。 %設(shè)置圖像窗口標(biāo)題figure(2)。 %創(chuàng)建圖像窗口（2）freqz(h)。 %畫圖，濾波器幅度響應(yīng)(db)和相位響應(yīng) %信號(hào)的處理%t=0::。 %設(shè)置t變量范圍，和步長n=2000。 %抽樣點(diǎn)數(shù)Fs=1000。 %抽樣頻率y=cos(2*pi*65*t)+*cos(2*pi*150*t)+2*sin(2*pi*230*t)+*sin(2*pi*320*t)。%輸入信號(hào)y1=fft(y)。 %輸入信號(hào)的傅里葉變換y2=fftshift(y1)。 %調(diào)整頻譜圖f=(0:1999)*Fs/nFs/2。 %計(jì)算頻率ffigure(3)。 %創(chuàng)建圖像窗口（3）plot(f,abs(y2)/Fs)。 %畫出輸入信號(hào)頻譜圖title(39。輸入信號(hào)頻譜圖39。)。 %設(shè)置圖像窗口標(biāo)題 %輸出的信號(hào)amp。G=fftfilt(h,y)。 %輸入信號(hào)y通過濾波器G1=fft(G)。 %濾波后，輸出信號(hào)傅里葉變換G2=fftshift(G1)。 %輸出信號(hào)傅里葉變換重新排布，使數(shù)據(jù)與頻率對應(yīng)figure(4)。 %創(chuàng)建圖像窗口（4）plot(f,abs(G2)/Fs)。 %畫出輸出信號(hào)頻譜圖title(39。輸出信號(hào)頻譜圖39。)。 %設(shè)置圖像窗口標(biāo)題
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號(hào)-1