基于matlab的fir濾波器設(shè)計與仿真-全文預覽
【正文】 %讀入語音a=a(:,1)。grid on。)。經(jīng)過濾波器后的信號a1139。m=0:ts:M。)。t39。title(39。 a0= randn(size(a)) .* sqrt(noise_var) + noise_mu。)。t39。title(39。ta=(length(a)1)/fs。ylabel(39。)。%畫出相頻特性subplot(2,1,2)pha=angle(h)。ylabel(39。)。%畫出幅頻特性figure(2)subplot(2,1,1)。disp(39。Sampling frequency in Hz=39。dev=input(39。)。)。f39。原始信號a的fft圖39。f=0:fs/FS:(FS1)*fs/FS。39。b=remez(N,fpts,mag,wt)。Sampling frequency in Hz=39。dev=input(39。)。設(shè)計一個FIR濾波器：、FIR濾波器設(shè)計的原理：①先對原始語音進行FFT變換，從FFT圖中得到passband和stopband從中我們可以知道了FIR濾波器的截止頻率。a1=a+a0。從中我們知道fs語音的采樣頻率。然后用wavread函數(shù)將語音文件讀入到matlab中，其實現(xiàn)的語句如下：[a,fs]=wavread(39。由于調(diào)用了Matlab高度集成化和可視化的計算機仿真環(huán)境，用戶可以不用編寫任何程序語句，僅僅通過鼠標拖放和點擊操作，就可完成復雜的數(shù)字濾波器的設(shè)計工作。常用的窗函數(shù)有以下幾種, 矩形窗, 三角窗,漢寧窗, 海明窗, 布拉克曼窗, 凱澤窗(β =7. 865) ,各種窗函數(shù)基本參數(shù)比較如下表1 。FIR 濾波器由于單位沖激響應是有限長的,可以用快速傅立葉變換( FFT) 算法來實現(xiàn)過濾信號,從而大大提高運算效率。語句產(chǎn)生與原始語音長度相同的隨機噪聲，最后把這個噪聲添加到原始語音中，得到加噪語音信號。學號：054100603 姓名：陳文珠 班級：電信051基于Matlab的FIR濾波器設(shè)計與仿真實驗地點：曹光彪信息樓219 實驗日期：7月13日至 7月18日一、實驗目的掌握基本的MATLAB編程方法；理解FIR濾波器的設(shè)計原理；學會用MATLAB來編程實現(xiàn)FIR濾波器；掌握基本的simulink交互式仿真，并對FIR濾波器模型并進行仿真；學會對所得的結(jié)果進行分析。然后通過 randn(size(a)) .* sqrt(noise_var) + noise_mu。只要經(jīng)過一定的延時,任何非因果有限長序列都能變成因果的有限長序列, 因而總能用因果系統(tǒng)來實現(xiàn)
法律信息相關(guān)推薦
鄂ICP備17016276號-1