Matlab課程設(shè)計(jì)報(bào)告-基于MATLAB有噪聲語(yǔ)音信號(hào)處理

Matlab課程設(shè)計(jì)報(bào)告題目:基于MATLAB有噪聲語(yǔ)音信號(hào)處理系（院）：計(jì)算機(jī)與信息工程學(xué)院專(zhuān)業(yè)：通信工程班級(jí)：10623102指導(dǎo)教師：學(xué)年學(xué)期:2011~2012學(xué)年第2學(xué)期簡(jiǎn)介：我們通信工程專(zhuān)業(yè)在實(shí)踐中經(jīng)常碰到需要對(duì)已接收信號(hào)進(jìn)行處理的情況，而濾波器設(shè)計(jì)在數(shù)字信號(hào)處理中占有極其重要的地位。本課題基于MATLAB有噪音語(yǔ)音信號(hào)處理的設(shè)計(jì)與實(shí)現(xiàn)，綜合運(yùn)用數(shù)字信號(hào)處理的理論知識(shí)對(duì)加噪語(yǔ)音信號(hào)進(jìn)行時(shí)域、頻域分析和濾波。通過(guò)理論推導(dǎo)得出相應(yīng)結(jié)論，再利用MATLAB作為編程工具進(jìn)行計(jì)算機(jī)實(shí)現(xiàn)。在設(shè)計(jì)實(shí)現(xiàn)的過(guò)程中，我們使用雙線性變換法設(shè)計(jì)IIR數(shù)字濾波器，對(duì)模擬加噪語(yǔ)音信號(hào)進(jìn)行低通濾波、高通濾波及帶通濾波，并利用MATLAB作為輔助工具完成設(shè)計(jì)中的計(jì)算與圖形的繪制。1緒論：數(shù)字信號(hào)處理是利用計(jì)算機(jī)或?qū)Ｓ锰幚碓O(shè)備，以數(shù)值計(jì)算的方法對(duì)信號(hào)進(jìn)行采集、抽樣、變換、綜合、估值與識(shí)別等加工處理，借以達(dá)到提取信息和便于應(yīng)用的目的。數(shù)字濾波器,是數(shù)字信號(hào)處理中及其重要的一部分。本課題采用IIR濾波器對(duì)加噪聲音信號(hào)進(jìn)行處理。IIR濾波器采用遞歸型結(jié)構(gòu)，即結(jié)構(gòu)上帶有反饋環(huán)路。IIR濾波器運(yùn)算結(jié)構(gòu)通常由延時(shí)、乘以系數(shù)和相加等基本運(yùn)算組成，可以組合成直接型、正準(zhǔn)型、級(jí)聯(lián)型、并聯(lián)型四種結(jié)構(gòu)形式，都具有反饋回路。同時(shí)，IIR數(shù)字濾波器在設(shè)計(jì)上可以借助成熟的模擬濾波器的成果，如巴特沃斯、契比雪夫和橢圓濾波器等，有現(xiàn)成的設(shè)計(jì)數(shù)據(jù)或圖表可查，在設(shè)計(jì)一個(gè)IIR數(shù)字濾波器時(shí)，我們根據(jù)指標(biāo)先寫(xiě)出模擬濾波器的公式，然后通過(guò)一定的變換，將模擬濾波器的公式轉(zhuǎn)換成數(shù)字濾波器的公式。2.原始語(yǔ)音信號(hào)采集與處理2.1語(yǔ)音信號(hào)的采集由于MATLAB只識(shí)別格式為.wav的聲音文件，我們利用PC機(jī)上的聲卡和WINDOWS操作系統(tǒng)進(jìn)行數(shù)字信號(hào)的采集。啟動(dòng)錄音機(jī)進(jìn)行錄音，以文件名“Orisound”保存入原程序所屬的文件夾中?？梢钥吹?文件存儲(chǔ)器的后綴默認(rèn)為.wav,這是WINDOWS操作系統(tǒng)規(guī)定的聲音文件存的標(biāo)準(zhǔn)。程序流程圖：開(kāi)始開(kāi)始用wavread讀取語(yǔ)音信號(hào)，并進(jìn)行采樣，建立數(shù)據(jù)文件，并用plot畫(huà)出數(shù)據(jù)文件時(shí)域波形圖和頻譜圖利用MATLAB中的隨機(jī)函數(shù)產(chǎn)生噪聲加入到語(yǔ)音信號(hào)中，模仿語(yǔ)音信號(hào)被污染，并對(duì)其進(jìn)行頻譜分析?；胤耪Z(yǔ)音信號(hào)得出所設(shè)計(jì)濾波器在語(yǔ)音處理中的優(yōu)劣并設(shè)計(jì)一個(gè)信號(hào)處理系統(tǒng)界面。結(jié)束運(yùn)用數(shù)字信號(hào)處理理論設(shè)計(jì)IIR數(shù)字濾波器，并對(duì)被噪聲污染的語(yǔ)音信號(hào)進(jìn)行濾波，分析濾波后信號(hào)的時(shí)域和頻域。特征。2.2語(yǔ)音信號(hào)的時(shí)頻分析利用MATLAB中的“wavread”命令來(lái)讀入（采集）語(yǔ)音信號(hào)，將它賦值給某一向量。再對(duì)其進(jìn)行采樣，記住采樣頻率和采樣點(diǎn)數(shù)。對(duì)語(yǔ)音信號(hào)Orisound.wav進(jìn)行采樣其程序如下：[y,fs,nbits]=wavread(‘OriSound’);%把語(yǔ)音信號(hào)加載入MATLAB仿真軟件平臺(tái)中畫(huà)出語(yǔ)音信號(hào)的時(shí)域波形，再對(duì)語(yǔ)音信號(hào)進(jìn)行頻譜分析。在本次設(shè)計(jì)中，我們利用fft函數(shù)對(duì)語(yǔ)音信號(hào)進(jìn)行快速傅里葉變換，就可以得到信號(hào)的頻譜特性。程序如下：[y,fs,nbits]=wavread(‘Orisound’);%語(yǔ)音信號(hào)的采集sound(y,fs,nbits);%語(yǔ)音信號(hào)的播放n=length(y);Y=fft(y,n);%快速傅里葉變換figure;subplot(2,1,1);plot(y);title(‘原始信號(hào)波形’,’fontweight’,’bold’);axis([7800080000-11]);grid;subplot(2,1,2);plot(abs(Y));title(‘原始信號(hào)頻譜’,’fontweight’,’bold’);axis([015000004000]);grid;程序結(jié)果如下圖：2.3語(yǔ)音信號(hào)加噪與頻譜分析利用MATLAB中的隨機(jī)函數(shù)Randn（m,n）產(chǎn)生噪聲加入到語(yǔ)音信號(hào)中，模仿語(yǔ)音信號(hào)被污染，并對(duì)其頻譜分析。主要程序如下：[y,fs,nbits]=wavread('OriSound');sound(y,fs,nbits);n=length(y);Noise=0.2*randn(n,2);s=y+Noise;sound(s);figure;subplot(2,1,1);plot(s);title('加噪語(yǔ)音信號(hào)的時(shí)域波形','fontweight','bold');axis([7800080000-11]);grid;S=fft(s);subplot(2,1,2);plot(abs(S));title('加噪語(yǔ)音信號(hào)的頻域波形','fontweight','bold');axis([015000004000]);grid;程序結(jié)果如下圖：3設(shè)計(jì)數(shù)字濾波器3.1數(shù)字濾波器設(shè)計(jì)的基本思路數(shù)字濾波器的實(shí)現(xiàn)有兩個(gè)關(guān)鍵步驟：一個(gè)從數(shù)字域到模擬域間的變換，這個(gè)變換實(shí)現(xiàn)了數(shù)字濾波器技術(shù)指標(biāo)到模擬濾波器技術(shù)指標(biāo)的轉(zhuǎn)換，同樣也實(shí)現(xiàn)了模擬濾波器系統(tǒng)函數(shù)到數(shù)字濾波器系統(tǒng)函數(shù)的轉(zhuǎn)換；另一個(gè)是從模擬濾波器技術(shù)指標(biāo)到滿(mǎn)足該指標(biāo)的模擬濾波器的設(shè)計(jì)。3.2模擬濾波器概述用模擬—數(shù)字變換法設(shè)計(jì)IIR數(shù)字濾波器，首先必須設(shè)計(jì)一個(gè)模擬濾波器，它有許多不同的類(lèi)型，主要有以下兩種類(lèi)型：（1）、巴特沃思（Botterworth簡(jiǎn)寫(xiě)B(tài)W）濾波器。BW濾波器是根據(jù)幅頻特性在通帶內(nèi)具有最拼圖特性而定義的濾波器，對(duì)一個(gè)N階低通濾波器來(lái)說(shuō)，所謂最平坦特性就是模擬函數(shù)的前（2N-1）階導(dǎo)數(shù)在處都為零。BW濾波器的另一個(gè)特性是在通帶和阻帶內(nèi)的幅頻特性始終是頻率的單調(diào)下降函數(shù)，且其模擬函數(shù)隨階次N的增大而更接近于理想低通濾波器。（2）、切比雪夫（Chbyshev簡(jiǎn)寫(xiě)為CB）濾波器。CB低通濾波器的模擬函數(shù)由切比雪夫多項(xiàng)式定義，且在通帶內(nèi)的幅頻響應(yīng)是波動(dòng)的，在阻帶則單調(diào)變化。3.3設(shè)計(jì)IIR濾波器目前IIR數(shù)字濾波器設(shè)計(jì)的最通用的方法是借助于模擬濾波器的設(shè)計(jì)方法。模擬濾波器設(shè)計(jì)已經(jīng)有了一套相當(dāng)成熟的方法，它不但有完整的設(shè)計(jì)公式，而且還有較為完整的圖表供查詢(xún)，因此，充分利用這些已有的資源將會(huì)給數(shù)字濾波器的設(shè)計(jì)帶來(lái)很大方便。IIR數(shù)字濾波器的設(shè)計(jì)步驟是：（1）、按一定規(guī)則將給出的數(shù)字濾波器的技術(shù)指標(biāo)轉(zhuǎn)換為模擬低通濾波器的技術(shù)指標(biāo)；（2）、根據(jù)轉(zhuǎn)換后的技術(shù)指標(biāo)設(shè)計(jì)模擬低通濾波器；（3）、再按一定規(guī)則將G(s)轉(zhuǎn)換成H(z)。若設(shè)計(jì)的數(shù)字濾波器是低通的，那么上述設(shè)計(jì)工作可以結(jié)束，若所設(shè)計(jì)的是高通，帶通或帶阻濾波器，那么還有步驟（4）：（4）、將高通、帶通、或帶阻數(shù)字l不去的技術(shù)指標(biāo)先轉(zhuǎn)化為低通模擬濾波器的技術(shù)指標(biāo)，然后按照上述步驟（2）設(shè)計(jì)出低通，再將轉(zhuǎn)換為所需的H(z)。利用模擬濾波器設(shè)計(jì)IIR數(shù)字低通濾波器的步驟：（1）確定數(shù)字低通濾波器的技術(shù)指標(biāo)：通帶邊界頻率、通帶最大衰減，阻帶截止頻率、阻帶最小衰減。（2）將數(shù)字低通濾波器的技術(shù)指標(biāo)轉(zhuǎn)換成相應(yīng)的模擬低通濾波器的技術(shù)指標(biāo)。（3）按照模擬低通濾波器的技術(shù)指標(biāo)設(shè)計(jì)及過(guò)渡模擬低通濾波器。（4）用雙線性變換法，模擬濾波器系統(tǒng)函數(shù)轉(zhuǎn)換成數(shù)字低通濾波器系統(tǒng)函數(shù)。程序如下：Ft=8000;Fp=1000;Fs=1200;wp=2*pi*Fp/Ft;ws=2*pi*Fs/Ft;fp=2*Ft*tan(wp/2);fs=2*Fs*tan(wp/2);[n11,wn11]=buttord(wp,ws,1,50,’s’);[b11,a11]=butter(n11,wn11,’s’);[num11,den11]=bilinear(b11,a11,0.5);[h,w]=freqz(num11,den11);figure;plot(w*8000*0.5/pi,abs(h));legend(‘IIR低通濾波器’,’Location’,’NorthWest’);grid;程序結(jié)果如下圖：3.2驗(yàn)證所設(shè)計(jì)的濾波器：為了驗(yàn)證濾波器的可使用性，我們用常用的sin函數(shù)來(lái)進(jìn)行驗(yàn)證。其具體程序及運(yùn)行結(jié)果如下：t=[0:1/1023:1];s=sin(2*pi*t);N=length(s);y=s+0.5*rand(1,N);subplot(2,1,1);plot(y);title('加噪語(yǔ)音信號(hào)的時(shí)域波形','fontweight','bold');S=fft(y);subplot(2,1,2);plot(abs(S));title('加噪語(yǔ)音信號(hào)的頻域波形','fontweight','bold');Ft=8000;Fp=1000;Fs=1200;wp=2*pi*Fp/Ft;ws=2*pi*Fs/Ft;[n11,wn11]=buttord(wp,ws,1,50,'s');%求低通濾波器的階數(shù)和截止頻率[b11,a11]=butter(n11,wn11,'s');%求S域的頻率響應(yīng)的參數(shù)[num11,den11]=bilinear(b11,a11,0.5);%利用雙線性變換實(shí)現(xiàn)頻率響應(yīng)S域到Z域的變換z11=filter(num11,den11,s);sound(z11);m11=fft(z11);%求濾波后的信號(hào)figure;subplot(2,2,1);plot(abs(S),'g');title('濾波前信號(hào)的頻譜','fontweight','bold');grid;subplot(2,2,2);plot(abs(m11),'r');title('濾波后信號(hào)的頻譜','fontweight','bold');grid;subplot(2,2,3);plot(y);title('濾波前信號(hào)的波形','fontweight','bold');grid;subplot(2,2,4);plot(z11);title('濾波后的信號(hào)波形','fontweight','bold');grid;程序結(jié)果如下圖：由所得結(jié)果可知，所設(shè)計(jì)的濾波器符合要求。4濾波用設(shè)計(jì)好的IIR低通濾波器對(duì)加噪的語(yǔ)音信號(hào)進(jìn)行濾波，程序如下：[y,fs,nbits]=wavread(‘OriSound’);%IIR低通n=length(y);%求出語(yǔ)音信號(hào)的長(zhǎng)度Noise=0.2*randn(n,2);%隨機(jī)函數(shù)產(chǎn)生噪聲s=y+Noise;%語(yǔ)音信號(hào)加入噪聲S=fft(s);Ft=8000;Fp=1000;Fs=1200;wp=2*pi*Fp/Ft;ws=2*pi*Fs/Ft;[n11,wn11]=buttord(wp,ws,1,50,’s’);%求低通濾波器的階數(shù)和截止頻率[b11,a11]=butter(n11,wn11,’s’);%求S域的頻率響應(yīng)的參數(shù)[num11,den11]=bilinear(b11,a11,0.5);%利用雙線性變換實(shí)現(xiàn)頻率響應(yīng)S域到Z域的變換z11=filter(num11,den11,s);sound(z11);m11=fft(z11);%求濾波后的信號(hào)figure;subplot(2,2,1);plot(abs(S),’g’);title(‘濾波前信號(hào)的頻譜’,’fontweight’,’bold’);axis([015000004000]);grid;subplot(2,2,2);plot(abs(m11),’r’);title(‘濾波后信號(hào)的頻譜’,’fontweight’,’bold’);axis([015000004000]);grid;subplot(2,2,3);plot(s);title(‘濾波前信號(hào)的波形’,’fontweight’,’bold’);axis([95000100000-11]);grid;subplot(2,2,4);plot(z11);title(‘濾波后的信號(hào)波形’,’fontweight’,’bold’);axis([95000100000-11]);grid;程序結(jié)果如下圖：經(jīng)過(guò)以上的加噪處理后，可在Matlab中用函數(shù)sound對(duì)聲音進(jìn)行回放。其調(diào)用格式：sound(y,Fs)，sound(y)和sound(y,Fs,bits)?？梢圆煊X(jué)濾波前后的聲音有明顯的變化。5.創(chuàng)建GUI界面總結(jié)本次設(shè)計(jì)通過(guò)選擇一個(gè)語(yǔ)音信號(hào)作為分析的對(duì)象，對(duì)其進(jìn)行頻譜分析；利用MATLAB中的隨機(jī)函數(shù)產(chǎn)生噪聲加入到語(yǔ)音信號(hào)中，模仿被噪聲干擾的語(yǔ)音信號(hào)，并對(duì)其進(jìn)行頻譜分析；運(yùn)用數(shù)字信號(hào)處理理論設(shè)計(jì)IIR數(shù)字濾波器，并對(duì)被噪聲污染的語(yǔ)音信號(hào)進(jìn)行濾波，分析濾波后信號(hào)的時(shí)域和頻域特征，回放語(yǔ)音信號(hào)。最后，設(shè)計(jì)一個(gè)信號(hào)處理系統(tǒng)界面返回。我們四個(gè)人共同完成這個(gè)課程設(shè)計(jì)，分工是很重要的，首先我們確定了大致的方向，我們要做低通，帶通以及高通三種濾波器，此外還有界面的設(shè)計(jì)，程序的整合，以及課程報(bào)告的編寫(xiě)等方面的內(nèi)容，我們的具體分工如下：聲音信號(hào)的讀入以及ppt的制作：翁淦泉1062310209低通，帶通，高通濾波器的設(shè)計(jì)及相關(guān)程序的編寫(xiě)：馮竹青1062310202，賀靜文1062301203最后整合及界面的設(shè)計(jì)：王雅青1062310208課程報(bào)告：賀靜文參考文獻(xiàn)1、張志涌楊祖櫻MATLAB教程R2011a北京航空航天大學(xué)出版社2、高西全，丁玉美.數(shù)字信號(hào)處理（第三版）[M].西安：西安電子科技大學(xué)出版社，2008.83、蔣瑞艷振動(dòng)，噪聲處理系統(tǒng)[J].大連理工大學(xué)學(xué)院學(xué)報(bào)(自然科學(xué)版)20024、馬超

高世倫基于MATLAB的噪聲信號(hào)采集與分析系統(tǒng)研究[J].華中科技大學(xué)能源與動(dòng)力工程學(xué)院院報(bào)2004.55、李正周.MATLAB數(shù)字信號(hào)處理與應(yīng)用[M].北京:清華大學(xué)出版社，20086、高萍,祖靜.基于MATLAB小波去噪技術(shù)淺析[J].科技信息(學(xué)術(shù)版),2006(11)1-3.7、劉智基于MATLAB的機(jī)械噪聲分析與處理[J].華中江師范學(xué)院學(xué)報(bào)(自然科學(xué)版)2010（4）附錄附錄（I）設(shè)計(jì)IIR數(shù)字濾波器%=========================IIR低通濾波器=======================Ft=8000;Fp=1000;Fs=1200;wp=2*pi*Fp/Ft;ws=2*pi*Fs/Ft;fp=2*Ft*tan(wp/2);fs=2*Fs*tan(wp/2);[n11,wn11]=buttord(wp,ws,1,50,’s’);[b11,a11]=butter(n11,wn11,’s’);[num11,den11]=bilinear(b11,a11,0.5);[h,w]=freqz(num11,den11);figure;plot(w*8000*0.5/pi,abs(h));legend(‘IIR低通濾波器’,’Location’,’NorthWest’);grid;程序結(jié)果如下圖：%======================IIR高通濾波器========================Ft=8000;Fp=4000;Fs=3500;wp1=tan(pi*Fp/Ft);ws1=tan(pi*Fs/Ft);wp=1;ws=wp1*wp/ws1;[n13,wn13]=cheb1ord(wp,ws,1,50,’s’);[b13,a13]=cheby1(n13,1,wn13,’s’);[num,den]=lp2hp(b13,a13,wn13);[num13,den13]=bilinear(num,den,0.5);[h,w]=freqz(num13,den13);figure;plot(w*21000*0.5/pi,abs(h));legend(‘IIR高通濾波器’,’Location’,’NorthWest’);axis([01100001.5]);grid;程序結(jié)果如下圖：%======================IIR帶通濾波器==========================Fp1=1200;Fp2=3000;Fs1=1000;Fs2=3200;Ft=8000;wp1=tan(pi*Fp1/Ft);wp2=tan(pi*Fp2/Ft);ws1=tan(pi*Fs1/Ft);ws2=tan(pi*Fs2/Ft);w=wp1*wp2/ws2;bw=wp2-wp1;%有效通帶頻率wp=1;ws=(wp1*wp2-w.^2)/(bw*w);[n12,wn12]=buttord(wp,ws,1,50,’s’);[b12,a12]=butter(n12,wn12,’s’);[num2,den2]=lp2bp(b12,a12,sqrt(wp1*wp2),bw);[num12,den12]=bilinear(num2,den2,0.5);[h,w]=freqz(num12,den12);figure;plot(w*8000*0.5/pi,abs(h));axis([0450001.5]);legend(‘IIR帶通濾波器’,’Location’,’NorthWest’);grid;程序結(jié)果如下圖：附錄（II）比較濾波前后語(yǔ)音信號(hào)的波形及頻譜%======================雙線性變換法=======================%*************************低通濾波器************************[y,fs,nbits]=wavread(‘OriSound’);%IIR低通n=length(y);%求出語(yǔ)音信號(hào)的長(zhǎng)度Noise=0.2*randn(n,2);%隨機(jī)函數(shù)產(chǎn)生噪聲s=y+Noise;%語(yǔ)音信號(hào)加入噪聲S=fft(s);Ft=8000;Fp=1000;Fs=1200;wp=2*pi*Fp/Ft;ws=2*pi*Fs/Ft;[n11,wn11]=buttord(wp,ws,1,50,’s’);%求低通濾波器的階數(shù)和截止頻率[b11,a11]=butter(n11,wn11,’s’);%求S域的頻率響應(yīng)的參數(shù)[num11,den11]=bilinear(b11,a11,0.5);%利用雙線性變換實(shí)現(xiàn)頻率響應(yīng)S域到Z域的變換z11=filter(num11,den11,s);sound(z11);m11=fft(z11);%求濾波后的信號(hào)figure;subplot(2,2,1);plot(abs(S),’g’);title(‘濾波前信號(hào)的頻譜’,’fontweight’,’bold’);axis([015000004000]);grid;subplot(2,2,2);plot(abs(m11),’r’);title(‘濾波后信號(hào)的頻譜’,’fontweight’,’bold’);axis([015000004000]);grid;subplot(2,2,3);plot(s);title(‘濾波前信號(hào)的波形’,’fontweight’,’bold’);axis([95000100000-11]);grid;subplot(2,2,4);plot(z11);title(‘濾波后的信號(hào)波形’,’fontweight’,’bold’);axis([95000100000-11]);grid;程序結(jié)果如下圖：附II-1雙線性低通濾波器比較%**********************高通濾波器*****************************[y,fs,nbits]=wavread(‘OriSound’);%IIR高通n=length(y);%求出語(yǔ)音信號(hào)的長(zhǎng)度Noise=0.2*randn(n,2);%隨機(jī)函數(shù)產(chǎn)生噪聲s=y+Noise;%語(yǔ)音信號(hào)加入噪聲S=fft(s);%傅里葉變換 Fp1=1200;Fs1=1000;Ft=8000;wp1=tan(pi*Fp1/Ft);ws1=tan(pi*Fs1/Ft);wp=1;ws=wp1*wp/ws1;[n13,wn13]=cheb1ord(wp,ws,1,50,’s’);%求模擬的低通濾波器階數(shù)和截止頻率[b13,a13]=cheby1(n13,1,wn13,’s’);%求S域的頻率響應(yīng)的參數(shù)[num,den]=lp2hp(b13,a13,wn13);%將S域低通參數(shù)轉(zhuǎn)為高通的[num13,den13]=bilinear(num,den,0.5);%利用雙線性變換實(shí)現(xiàn)頻率響應(yīng)S域到Z域轉(zhuǎn)換z13=filter(num13,den13,s);sound(z13);m13=fft(z13);%求濾波后的信號(hào)figure;subplot(2,2,1);plot(abs(S),’g’);title(‘濾波前信號(hào)的頻譜’,’fontweight’,’bold’);axis([015000004000]);grid;subplot(2,2,2);plot(abs(m13),’r’);title(‘濾波后信號(hào)的頻譜’,’fontweight’,’bold’);axis([015000004000]);grid;subplot(2,2,3);plot(s);title(‘濾波前信號(hào)的波形’,’fontweight’,’bold’);axis([95000100000-11]);grid;subplot(2,2,4);plot(z13);title(‘濾波后的信號(hào)波形’,’fontweight’,’bold’);axis([95000100000-11]);grid;程序結(jié)果如下圖：附II-2雙線性高通濾波器比較%**********************帶通濾波器*****************************[y,fs,nbits]=wavread(‘OriSound’);%IIR帶通n=length(y);%求出語(yǔ)音信號(hào)的長(zhǎng)度Noise=0.2*randn(n,2);%隨機(jī)函數(shù)產(chǎn)生噪聲s=y+Noise;%語(yǔ)音信號(hào)加入噪聲S=fft(s);%傅里葉變換Ft=8000;Fp=1000;Fs=1200;wp=2*Fp/Ft;ws=2*Fs/Ft;rp=1;rs=50;p=1-10.^(-rp/20);%通帶阻帶波紋q=10.^(-rs/20);fpts=[wpws];mag=[10];dev=[pq];[n21,wn21,beta,ftype]=kaiserord(fpts,mag,dev);%由kaiserord求濾波器的階數(shù)和截止頻率b21=fir1(n21,wn21,kaiser(n21+1,beta));%由fir1設(shè)計(jì)濾波器z21=fftfilt(b21,s);sound(z21);m21=fft(z21);%求濾波后的信號(hào)figure(4);subplot(2,2,1);plot(abs(S),’g’);title(‘濾波前信號(hào)的頻譜’,’fontweight’,’bold’);axis([015000004000]);grid;subplot(2,2,2);plot(abs(m21),’r’);title(‘濾波后信號(hào)的頻譜’,’fontweight’,’bold’);axis([015000004000]);grid;subplot(2,2,3);plot(s);title(‘濾波前信號(hào)的波形’,’fontweight’,’bold’);axis([95000100000-11]);grid;subplot(2,2,4);plot(z21);title(‘濾波后的信號(hào)波形’,’fontweight’,’bold’);axis([95000100000-11]);grid;程序結(jié)果如下圖：附II-3雙線性帶通濾波器比較基于C8051F單片機(jī)直流電動(dòng)機(jī)反饋控制系統(tǒng)的設(shè)計(jì)與研究基于單片機(jī)的嵌入式Web服務(wù)器的研究MOTOROLA單片機(jī)MC68HC（8）05PV8/A內(nèi)嵌EEPROM的工藝和制程方法及對(duì)良率的影響研究基于模糊控制的電阻釬焊單片機(jī)溫度控制系統(tǒng)的研制基于MCS-51系列單片機(jī)的通用控制模塊的研究基于單片機(jī)實(shí)現(xiàn)的供暖系統(tǒng)最佳啟停自校正（STR）調(diào)節(jié)器單片機(jī)控制的二級(jí)倒立擺系統(tǒng)的研究基于增強(qiáng)型51系列單片機(jī)的TCP/IP協(xié)議棧的實(shí)現(xiàn)基于單片機(jī)的蓄電池自動(dòng)監(jiān)測(cè)系統(tǒng)基于32位嵌入式單片機(jī)系統(tǒng)的圖像采集與處理技術(shù)的研究基于單片機(jī)的作物營(yíng)養(yǎng)診斷專(zhuān)家系統(tǒng)的研究基于單片機(jī)的交流伺服電機(jī)運(yùn)動(dòng)控制系統(tǒng)研究與開(kāi)發(fā)基于單片機(jī)的泵管內(nèi)壁硬度測(cè)試儀的研制基于單片機(jī)的自動(dòng)找平控制系統(tǒng)研究基于C8051F040單片機(jī)的嵌入式系統(tǒng)開(kāi)發(fā)基于單片機(jī)的液壓動(dòng)力系統(tǒng)狀態(tài)監(jiān)測(cè)儀開(kāi)發(fā)模糊Smith智能控制方法的研究及其單片機(jī)實(shí)現(xiàn)一種基于單片機(jī)的軸快流CO〈,2〉激光器的手持控制面板的研制基于雙單片機(jī)沖床數(shù)控系統(tǒng)的研究基于CYGNAL單片機(jī)的在線間歇式濁度儀的研制基于單片機(jī)的噴油泵試驗(yàn)臺(tái)控制器的研制基于單片機(jī)的軟起動(dòng)器的研究和設(shè)計(jì)基于單片機(jī)控制的高速快走絲電火花線切割機(jī)床短循環(huán)走絲方式研究基于單片機(jī)的機(jī)電產(chǎn)品控制系統(tǒng)開(kāi)發(fā)基于PIC單片機(jī)的智能手機(jī)充電器基于單片機(jī)的實(shí)時(shí)內(nèi)核設(shè)計(jì)及其應(yīng)用研究基于單片機(jī)的遠(yuǎn)程抄表系統(tǒng)的設(shè)計(jì)與研究基于單片機(jī)的煙氣二氧化硫濃度檢測(cè)儀的研制基于微型光譜儀的單片機(jī)系統(tǒng)單片機(jī)系統(tǒng)軟件構(gòu)件開(kāi)發(fā)的技術(shù)研究基于單片機(jī)的液體點(diǎn)滴速度自動(dòng)檢測(cè)儀的研制基于單片機(jī)系統(tǒng)的多功能溫度測(cè)量?jī)x的研制基于PIC單片機(jī)的電能采集終端的設(shè)計(jì)和應(yīng)用基于單片機(jī)的光纖光柵解調(diào)儀的研制氣壓式線性摩擦焊機(jī)單片機(jī)控制系統(tǒng)的研制基于單片機(jī)的數(shù)字磁通門(mén)傳感器基于單片機(jī)的旋轉(zhuǎn)變壓器-數(shù)字轉(zhuǎn)換器的研究基于單片機(jī)的光纖Bragg光柵解調(diào)系統(tǒng)的研究單片機(jī)控制的便攜式多功能乳腺治療儀的研制基于C8051F020單片機(jī)的多生理信號(hào)檢測(cè)儀基于單片機(jī)的電機(jī)運(yùn)動(dòng)控制系統(tǒng)設(shè)計(jì)Pico專(zhuān)用單片機(jī)核的可測(cè)性設(shè)計(jì)研究基于MCS-51單片機(jī)的熱量計(jì)基于雙單片機(jī)的智能遙測(cè)微型氣象站MCS-51單片機(jī)構(gòu)建機(jī)器人的實(shí)踐研究基于單片機(jī)的輪軌力檢測(cè)基于單片機(jī)的GPS定位儀的研究與實(shí)現(xiàn)基于單片機(jī)的電液伺服控制系統(tǒng)用于單片機(jī)系統(tǒng)的MMC卡文件系統(tǒng)研制基于單片機(jī)的時(shí)控和計(jì)數(shù)系統(tǒng)性能優(yōu)化的研究基于單片機(jī)和CPLD的粗光柵位移測(cè)量系統(tǒng)研究單片機(jī)控制的后備式方波UPS提升高職學(xué)生單片機(jī)應(yīng)用能力的探究基于單片機(jī)控制的自動(dòng)低頻減載裝置研究基于單片機(jī)控制的水下焊接電源的研究基于單片機(jī)的多通道數(shù)據(jù)采集系統(tǒng)基于uPSD3234單片機(jī)的氚表面污染測(cè)量?jī)x的研制基于單片機(jī)的紅外測(cè)油儀的研究96系列單片機(jī)仿真器研究與設(shè)計(jì)基于單片機(jī)的單晶金剛石刀具刃磨設(shè)備的數(shù)控改造基于單片機(jī)的溫度智能控制系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)基于MSP430單片機(jī)的電梯門(mén)機(jī)控制器的研制基于單片機(jī)的氣體測(cè)漏儀的研究基于三菱M16C/6N系列單片機(jī)的CAN/USB協(xié)議轉(zhuǎn)換器基于單片機(jī)和DSP的變壓器油色譜在線監(jiān)測(cè)技術(shù)研究基于單片機(jī)的膛壁溫度報(bào)警系統(tǒng)設(shè)計(jì)基于AVR單片機(jī)的低壓無(wú)功補(bǔ)償控制器的設(shè)計(jì)基于單片機(jī)船舶電力推進(jìn)電機(jī)監(jiān)測(cè)系統(tǒng)基于單片機(jī)網(wǎng)絡(luò)的振動(dòng)信號(hào)的采集系統(tǒng)基于單片機(jī)的大容量數(shù)據(jù)存儲(chǔ)技術(shù)的應(yīng)用研究基于單片機(jī)的疊圖機(jī)研究與教學(xué)方法實(shí)踐基于單片機(jī)嵌入式Web服務(wù)器技術(shù)的研究及實(shí)現(xiàn)基于AT89S52單片機(jī)的通用數(shù)據(jù)采集系統(tǒng)基于單片機(jī)的多道脈沖幅度分析儀研究機(jī)器人旋轉(zhuǎn)電弧傳感角焊縫跟蹤單片機(jī)控制系統(tǒng)基于單片機(jī)的控制系統(tǒng)在PLC虛擬教學(xué)實(shí)驗(yàn)中的應(yīng)用研究基于單片機(jī)系統(tǒng)的網(wǎng)絡(luò)通信研究與應(yīng)用基于PIC16F877單片機(jī)的莫爾斯碼自動(dòng)譯碼系統(tǒng)設(shè)計(jì)與研究基于單片機(jī)的模糊控制器在工業(yè)電阻爐上的應(yīng)用研究基于雙單片機(jī)沖床數(shù)控系統(tǒng)的研究與開(kāi)發(fā)基于Cygnal單片機(jī)的μC/OS-Ⅱ的研究基于單片機(jī)的一體化智能差示掃描量熱儀系統(tǒng)研究基于TCP/IP協(xié)議的單片機(jī)與Internet互聯(lián)的研究與實(shí)現(xiàn)變頻調(diào)速液壓電梯單片機(jī)控制器的研究基于單片機(jī)γ-免疫計(jì)數(shù)器自動(dòng)換樣功能的研究與實(shí)現(xiàn)基于單片機(jī)的倒立擺控制系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)單片機(jī)嵌入式以太網(wǎng)防盜報(bào)警系統(tǒng)基于51單片機(jī)的嵌入式Internet系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)單片機(jī)監(jiān)測(cè)系統(tǒng)在擠壓機(jī)上的應(yīng)用MSP430單片機(jī)在智能水表系統(tǒng)上的研究與應(yīng)用HYPERLINK"/d