【正文】 ..........15 FIR設(shè)計(jì)步驟和過程....................................................................................................154 設(shè)計(jì)程序的調(diào)試和運(yùn)行結(jié)果.............................................................................................16...........................................................16...........................................................18...........................................................20 hamming濾波器程序的調(diào)試和運(yùn)行結(jié)果...................................................................235 總結(jié).......................................................................................................................................25參考文獻(xiàn)...................................................................................................................................27附錄…………………………………………………………………………………...………28附錄 A..............................................................................................................................28附錄 B..............................................................................................................................33第 43 頁 共 56 頁1 緒論 課題背景隨著軟硬件技術(shù)的發(fā)展，儀器的智能化與虛擬化已成為未來實(shí)驗(yàn)室及研究機(jī)構(gòu)的發(fā)展方向。 語音信號(hào)的數(shù)字濾波處理(九) CHANGSHA UNIVERSITY OF SCIENCE amp。虛擬儀器技術(shù)的優(yōu)勢(shì)在于可由用戶定義自己的專用儀器系統(tǒng)，且功能靈活，很容易構(gòu)建，所以應(yīng)用面極為廣泛。數(shù)字濾波器(Digital Filter, DF)在信號(hào)的過濾、檢測(cè)和參數(shù)估計(jì)等方面起著重要的作用。 課程設(shè)計(jì)預(yù)習(xí) 卷積運(yùn)算的演示(1) 線性卷積序列x1(n)=[2 0 1 2 5 7 0 5 0 1 1 2]，序列x2(n)=[2 0 1 2 5 7 0 5]。)ylabel(39。)xn2=[2 0 1 2 5 7 0 5]。)yn(1)=sum((xm1.*xm2)39。)for N=1:17xm3=[zeros(1,N) fliplr(xn2) zeros(1,M)]。yn(N+1)=sum((xm1.*xm3)39。)endxm3=[zeros(1,18) fliplr(xn2)]figure()stem(xm3)。)。(2) 循環(huán)卷積序列x1(n)=[2 0 0 7 5 7 1 7 0 1 1 9]，序列x2(n)=[2 0 0 7 4 7 0 3]，N=12。)ylabel(39。xm2=[xn2 zeros(1,length(xn1)length(xn2))]。yn(1)=sum((xn1.*xm2)39。)for N=1:11xm1=[fliplr(xn1(1:N)) fliplr(xn1(N+1:12))]。)。stem(n,xm1)xlabel(39。figure()stem(n,yn)xlabel(39。(3) 聲音文件線性卷積序列x1(n)=[2 0 1 2 5 7 0 5],讀取一段聲音數(shù)據(jù)（Windows XP 菜單命令），當(dāng)循環(huán)卷積長(zhǎng)度大于或等于兩序列長(zhǎng)度之和時(shí)，循環(huán)卷積等于線性卷積。 zeros(1,7)]。)endxm3=[zeros(1,686) fliplr(xn2)]。n=N39。X=fft(xn2,1024)。 n=12。 subplot(5,1,1) plot(x1)。 subplot(5,1,5) plot(yn)； 采樣定理演示1② 要求輸入采樣頻率fs（根據(jù)程序處理需要指定范圍）后，在時(shí)域演示信號(hào)波形、采樣脈沖及采樣后信號(hào)；在頻域演示不同采樣頻率下對(duì)應(yīng)信號(hào)的頻譜。 subplot(5,1,1) plot(x1)。 subplot(5,1,5) plot(yn)；. 采樣定理演示2 課程設(shè)計(jì)原理 頻譜分析原理時(shí)域分析只能反映信號(hào)的幅值隨時(shí)間的變化情況，除單頻率分量的簡(jiǎn)單波形外，很難明確提示信號(hào)的頻率組成和各頻率分量大小，而頻譜分析能很好的解決此問題。以頻率f為橫坐標(biāo)，為縱坐標(biāo)，可以得到幅值譜；以頻率f為橫坐標(biāo)，為縱坐標(biāo)，可以得到相位譜；以頻率f為橫坐標(biāo)，為縱坐標(biāo)，可以得到實(shí)頻譜；以頻率f為橫坐標(biāo)，為縱坐標(biāo)，可以得到虛頻譜。脈沖響應(yīng)不變法是從S平面映射到z平面，這種映射不是簡(jiǎn)單的代數(shù)映射，而是S平面的每一條寬為2π/T的橫帶重復(fù)地映射到整個(gè)z平面。在MATLAB中，實(shí)現(xiàn)布拉克曼窗的函數(shù)為blackman，調(diào)用格式如下： w=blackman (N)w=blackman (N,39。 瓣雙線性Z變換滿足以上4個(gè)條件的映射關(guān)系，其變換公式為 FIR課程設(shè)計(jì)步驟和過程(1)根據(jù)對(duì)阻帶衰減及過渡帶的指標(biāo)要求，選擇窗函數(shù)的類型，并估計(jì)窗口長(zhǎng)度N；(2)構(gòu)造希望逼近的頻率響應(yīng)函數(shù)Hd(ejw)，即Hd(ejw)=Hdg(w)ejw（N1）/2；(3)計(jì)算hd(n)；(4)加窗得到設(shè)計(jì)結(jié)果：h(n)=hd(n)w(n)。當(dāng)無相移時(shí)，N為奇數(shù)時(shí)，所設(shè)計(jì)的FIR數(shù)字帶通濾波器為Ⅰ型濾波器；當(dāng)N為偶數(shù)時(shí)，為Ⅱ型濾波器。 原語音信號(hào)時(shí)域圖和頻譜圖。其設(shè)計(jì)程序如下：wpz=2*4500/fs。)。 濾波后語音信號(hào)的時(shí)域圖和頻譜圖 切比雪夫高通程序的調(diào)試和運(yùn)行結(jié)果 。x=[rot90(x1,3),rot90(x1,3)]。wsz=2*500/fs。[b,a]=cheby2(N,rs,wc,39。其濾波程序如下：Yk=Z.*Hk。x1=*sin(2*pi*100*t)+*sin(2*pi*220*t)+ *sin(2*pi*6000*t)。 rp=20。[N,wc]=cheb2ord(wp,ws,rp,rs,39。 帶通濾波器的幅頻圖。通過老師的指導(dǎo)，和自己的思考，總結(jié)如下：?jiǎn)栴}1：Subscript indices must either be real positive integers or logicals.解決辦法：在訪問矩陣（包括向量、二維矩陣、多維數(shù)組，下同）的過程中，下標(biāo)索引要么從 0 開始，要么出現(xiàn)了負(fù)數(shù)。解決方法：檢查打開文件目錄，一般打開文件需要在MATBLE下的work內(nèi)。(7) 數(shù)組元素求和：sum直方圖：hist(10) 矩陣分析：左右翻轉(zhuǎn)：fliplr線性變化數(shù)組：linspace，聚合矩陣：cat/horzcat/vertcat[8]二項(xiàng)分布：binorndG:\Windows 39。y39。xlabel(39。)t=0::。)。subplot(2,1,2)。title(39。)。)。rp=2。[bz,az]=bilinear(b,a,fs)。)。f=22050*(0:16384)/32768。)。n39。 (2) 切比雪夫高通：[y,fs,nbits]=wavread(39。ylabel(39。plot(f,abs(Y(1:16385)))。原語音信號(hào)頻譜圖39。n39。f=22050*(0:16384)/32768。)。)。s39。figure(3)Hk=freqz(bz,az,16384,fs)。title(39。subplot(2,1,2)。title(39。)。G:\Windows 39。y39。xlabel(39。)t=0::。)。subplot(2,1,2)。title(39。wpz=[2*700/fs,2*4500/fs]。)。f39。Hk=[H fliplr(H)]。f39。subplot(2,1,1)。)。plot(x1)。title(39。title(39。)。hng=fir1(Ng1,wcg,39。)。Yg=fft(yg)。)。)。wpd=pi*。subplot(2,2,1)。title(39。濾波后信號(hào)頻譜39。subplot(2,2,3)。加噪聲聲音濾波后波形39。Mdai=ceil(*pi/Btdai)。)。Ydai=fft(ydai)。)。)。盡我所知，除文中特別加以標(biāo)注和致謝的地方外，不包含其他人或組織已經(jīng)發(fā)表或公布過的研究成果，也不包含我為獲得 及其它教育機(jī)構(gòu)的學(xué)位或?qū)W歷而使用過的材料。本人授權(quán)　　 　 　大學(xué)可以將本學(xué)位論文的全部或部分內(nèi)容編入有關(guān)數(shù)據(jù)庫進(jìn)行檢索，可以采用影印、縮印或掃描等復(fù)制手段保存和匯編本學(xué)位論文。作者簽名： 日期： 年 月 日導(dǎo)師簽名： 日期： 年 月 日注 意 事 項(xiàng)（論文）的內(nèi)容包括：1）封面（按教務(wù)處制定的標(biāo)準(zhǔn)封面格式制作）2）原創(chuàng)性聲明3）中文摘要（300字左右）、關(guān)鍵詞4）外文摘要、關(guān)鍵詞 5）目次頁（附件不統(tǒng)一編入）6）論文主體部分：引言（或緒論）、正文、結(jié)論7）參考文獻(xiàn)8）致謝9）附錄（對(duì)論文支持必要時(shí)）：理工類設(shè)計(jì)（論文）正文字?jǐn)?shù)不少于1萬字（不包括圖紙、程序清單等）。作 者 簽 名：　　 　 　　日　 期：　　 　 　　 指導(dǎo)教師簽名：　　 　 　　 日　　期：　　 　 　　 使用授權(quán)說明本人完全了解 大學(xué)關(guān)于收集、保存、使用畢業(yè)設(shè)計(jì)（論文）的規(guī)定，即：按照學(xué)校要求提交畢業(yè)設(shè)計(jì)（論文）的印刷本和電子版本；學(xué)校有權(quán)保存畢業(yè)設(shè)計(jì)（論文）的印刷本和電子版，并提供目錄檢索與閱覽服務(wù)；學(xué)?？梢圆捎糜坝?、縮印、數(shù)字化或其它復(fù)制手段保存論文；在不以贏利為目的前提下，學(xué)?？梢怨颊撐牡牟糠只蛉?jī)?nèi)容。plot(ydai)。Amplitude39。plot(abs(Ydai))。Normaliaed Magnitude(dB)39。hndai=fir1(Mdai,wcdai,