基于matlab的fir濾波器設(shè)計(jì)與仿真-文庫(kù)吧
2025-06-08 15:53 本頁(yè)面
【正文】 即觀察到改變后的仿真結(jié)果。本文將介紹在Simulink環(huán)境下應(yīng)用數(shù)字濾波器設(shè)計(jì)模塊進(jìn)行信號(hào)處理（濾波）的可視化設(shè)計(jì)技術(shù)。由于調(diào)用了Matlab高度集成化和可視化的計(jì)算機(jī)仿真環(huán)境，用戶可以不用編寫(xiě)任何程序語(yǔ)句，僅僅通過(guò)鼠標(biāo)拖放和點(diǎn)擊操作，就可完成復(fù)雜的數(shù)字濾波器的設(shè)計(jì)工作。下面通過(guò)一個(gè)應(yīng)用示例介紹這種基于可視化的信號(hào)濾波技術(shù)，用數(shù)字濾波器設(shè)計(jì)模塊（Digital Filter Design Block）設(shè)計(jì)、分析和實(shí)現(xiàn)濾波器。三、實(shí)驗(yàn)結(jié)果與分析制作數(shù)據(jù)源：我用電腦錄制了同學(xué)的一段語(yǔ)音，內(nèi)容是：現(xiàn)代語(yǔ)音處理技術(shù)。該語(yǔ)音大概有4，5秒之長(zhǎng)。然后用wavread函數(shù)將語(yǔ)音文件讀入到matlab中，其實(shí)現(xiàn)的語(yǔ)句如下：[a,fs]=wavread(39。39。)。%讀入語(yǔ)音a=a(:,1)。從中我們知道fs語(yǔ)音的采樣頻率。信號(hào)中混入隨機(jī)噪聲：noise_mu = 0。 noise_var =。 a0= randn(size(a)) .* sqrt(noise_var) + noise_mu。a1=a+a0。事先取一個(gè)參考值為0，，噪聲強(qiáng)度不宜過(guò)大或過(guò)小，適中即可。然后通過(guò) randn(size(a)) .* sqrt(noise_var) + noise_mu。語(yǔ)句產(chǎn)生與原始語(yǔ)音長(zhǎng)度相同的隨機(jī)噪聲，最后把這個(gè)噪聲添加到原始語(yǔ)音中，得到a1是信號(hào)中混入隨機(jī)噪聲的的語(yǔ)音。設(shè)計(jì)一個(gè)FIR濾波器：、FIR濾波器設(shè)計(jì)的原理：①先對(duì)原始語(yǔ)音進(jìn)行FFT變換，從FFT圖中得到passband和stopband從中我們可以知道了FIR濾波器的截止頻率。確定我們選擇是用低通，高通還是帶通來(lái)處理加噪語(yǔ)音。②然后我們通過(guò)：fedge=input(39。Band edges in Hz=39。)。mval=input(39。Desired magnitude values in each band=39。)。dev=input(39。Desired ripple in each band=39。)。FT=input(39。Sampling frequency in Hz=39。)。輸入從FFT上得到的數(shù)據(jù)，passband和stopband，低通或高通，帶通，通帶紋波和阻帶紋波，以及從wavread函數(shù)中得到的采樣頻率fs。通過(guò)：[N,fpts,mag,wt]=remezord(fedge,mval,dev,FT)。b=remez(N,fpts,mag,wt)。產(chǎn)生濾波器的階數(shù)N和濾波器的系數(shù)b。③然后把加噪的語(yǔ)音a1與b進(jìn)行卷積：a11=conv(b,a1) 得到的是濾波后的語(yǔ)音信號(hào)。、FIR濾波器設(shè)計(jì)的程序：[a
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1