基于matlab的fir濾波器設(shè)計(jì)與仿真-免費(fèi)閱讀
【正文】 開始用MATLAB來(lái)編程實(shí)現(xiàn)對(duì)我來(lái)說(shuō)還是有點(diǎn)難度的，但是通過(guò)查看相關(guān)資料以及認(rèn)真的閱讀老師提供的資料，我慢慢的對(duì)實(shí)現(xiàn)FIR濾波器設(shè)計(jì)有了初步的程序輪廓，然后我花了將近一天的時(shí)間來(lái)編程，中間出現(xiàn)過(guò)很多困難，都是自己對(duì)函數(shù)的使用不熟悉，但在老師的幫助下我都逐個(gè)解決了，并最后實(shí)現(xiàn)FIR濾波器功能，對(duì)加噪語(yǔ)音進(jìn)行濾波，輸出的語(yǔ)音與原始語(yǔ)音在聽覺效果上分辨不出有何差別，很好的實(shí)現(xiàn)了濾波效果。用Simulink設(shè)計(jì)FIR濾波器模型并進(jìn)行仿真，所得結(jié)果進(jìn)行分析因?yàn)闉V波器的階數(shù)是N=27階，則我通過(guò)Simulink設(shè)計(jì)了一個(gè)27階的濾波器，其圖形如下，下圖是simulink總的設(shè)計(jì)模型，其中我們可以看到示波器可以顯示三路的波形，第一路是原始語(yǔ)音信號(hào)的波形，第二路是加噪語(yǔ)音信號(hào)的波形，第三路是加噪語(yǔ)音信號(hào)進(jìn)過(guò)濾波后的語(yǔ)音信號(hào)輸出，通過(guò)三路波形的同時(shí)顯示我們可以很好的比較原始語(yǔ)音信號(hào)與濾波后的語(yǔ)音信號(hào)之間的差別，從而可以判斷濾波器的濾波效果是否合格。)。title(39。a139。plot(t,a1)。a39。plot(t,a)。)。grid on。幅頻特性39。b=remez(N,fpts,mag,wt)。)。abs(x)39。title(39。、FIR濾波器設(shè)計(jì)的程序：[a,fs]=wavread(39。FT=input(39。Band edges in Hz=39。 a0= randn(size(a)) .* sqrt(noise_var) + noise_mu。該語(yǔ)音大概有4，5秒之長(zhǎng)。有限長(zhǎng)FIR數(shù)字濾波器的設(shè)計(jì)方法主要是窗函數(shù)設(shè)計(jì)法。然后通過(guò) randn(size(a)) .* sqrt(noise_var) + noise_mu。語(yǔ)句產(chǎn)生與原始語(yǔ)音長(zhǎng)度相同的隨機(jī)噪聲，最后把這個(gè)噪聲添加到原始語(yǔ)音中，得到加噪語(yǔ)音信號(hào)。常用的窗函數(shù)有以下幾種, 矩形窗, 三角窗,漢寧窗, 海明窗, 布拉克曼窗, 凱澤窗(β =7. 865) ,各種窗函數(shù)基本參數(shù)比較如下表1 。然后用wavread函數(shù)將語(yǔ)音文件讀入到matlab中，其實(shí)現(xiàn)的語(yǔ)句如下：[a,fs]=wavread(39。a1=a+a0。)。Sampling frequency in Hz=39。39。原始信號(hào)a的fft圖39。)。dev=input(39。disp(39。)。%畫出相頻特性subplot(2,1,2)pha=angle(h)。ylabel(39。title(39。)。title(39。)。經(jīng)過(guò)濾波器后的信號(hào)a1139。grid on。上圖的Filter濾波器是下圖所示的子系統(tǒng)封裝得到的，該子系統(tǒng)設(shè)計(jì)了階數(shù)為27階的FIR濾波器，其中的濾波器系數(shù)是
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1