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