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