《MATLAB課程設(shè)計(jì)》報(bào)告DOC

《MATLAB課程設(shè)計(jì)》報(bào)告DOC《MATLAB課程設(shè)計(jì)》報(bào)告DOC/《MATLAB課程設(shè)計(jì)》報(bào)告DOC《MATLAB課程設(shè)計(jì)》報(bào)告設(shè)計(jì)題目：鑒于MATLAB的語(yǔ)音信號(hào)收集與辦理學(xué)生姓名：學(xué)生學(xué)號(hào)：專業(yè)班級(jí)：光信息科學(xué)與技術(shù)辯論時(shí)間：2013年12月指導(dǎo)教師：馮明庫(kù)廣東技術(shù)師范學(xué)院電子與信息學(xué)院一、設(shè)計(jì)目的及意義1.MATLAB

軟件功能簡(jiǎn)介MATLAB的名稱源自

MatrixLaboratory,1984

年由美國(guó)

Mathworks

企業(yè)推向市場(chǎng)。它是一種科學(xué)計(jì)算軟件，特意以矩陣的形式辦理數(shù)據(jù)。MATLAB將高性能的數(shù)值計(jì)算和可視化集成在一同，并供給了大批的內(nèi)置函數(shù)，進(jìn)而被寬泛的應(yīng)用于科學(xué)計(jì)算、控制系統(tǒng)和信息辦理等領(lǐng)域的剖析、仿真和設(shè)計(jì)工作。MATLAB軟件包含五大通用功能，數(shù)值計(jì)算功能（Nemeric）、符號(hào)運(yùn)算功能（Symbolic）、數(shù)據(jù)可視化功能（Graphic）、數(shù)字圖形文字一致辦理功能（Notebook）和建模擬真可視化功能（Simulink）。此中，符號(hào)運(yùn)算功能的實(shí)現(xiàn)是經(jīng)過(guò)懇求MAPLE內(nèi)核計(jì)算并將結(jié)果返回到MATLAB命令窗口。該軟件有三大特色，一是功能強(qiáng)盛；二是界面友好、語(yǔ)言自然；三是開(kāi)放性強(qiáng)。當(dāng)前，Mathworks企業(yè)已推出30多個(gè)應(yīng)用工具箱。MATLAB在線性代數(shù)、矩陣剖析、數(shù)值及優(yōu)化、數(shù)值統(tǒng)計(jì)和隨機(jī)信號(hào)剖析、電路與系統(tǒng)、系統(tǒng)動(dòng)力學(xué)、次那好和圖像辦理、控制理論剖析和系統(tǒng)設(shè)計(jì)、過(guò)程控制、建模和仿真、通訊系統(tǒng)以及財(cái)政金融等眾多領(lǐng)域的理論研究和工程設(shè)計(jì)中獲得了寬泛應(yīng)用。MATLAB在信號(hào)與系統(tǒng)中的應(yīng)用主要包含符號(hào)運(yùn)算和數(shù)值計(jì)算仿真剖析。因?yàn)樾盘?hào)與系統(tǒng)課程的很多內(nèi)容都是鑒于公式演算，而MATLAB借助符號(hào)數(shù)學(xué)工具箱供給的符號(hào)運(yùn)算功能，能基本知足信號(hào)與系統(tǒng)課程的需求。比如解微分方程、傅里葉正反變換、拉普拉斯正反變換和z正反變換等。MATLAB在信號(hào)與系統(tǒng)中的另一主要應(yīng)用是數(shù)值計(jì)算與仿真剖析，主要包含函數(shù)波形繪制、函數(shù)運(yùn)算、沖擊響應(yīng)與階躍響應(yīng)仿真剖析、信號(hào)的時(shí)域剖析、信號(hào)的頻譜剖析、系統(tǒng)的S域剖析和零極點(diǎn)圖繪制等內(nèi)容。數(shù)值計(jì)算仿真剖析能夠幫助學(xué)生更深入地理解理論知識(shí)，并為未來(lái)使用MATLAB進(jìn)行信號(hào)辦理領(lǐng)域的各樣剖析和實(shí)質(zhì)應(yīng)用打下基礎(chǔ)。此題目的意義本次課程設(shè)計(jì)的課題為《鑒于MATLAB的語(yǔ)音信號(hào)收集與辦理》，學(xué)會(huì)運(yùn)用MATLAB的信號(hào)辦理功能，收集語(yǔ)音信號(hào)，并對(duì)語(yǔ)音信號(hào)進(jìn)行濾涉及變換辦理，察看其時(shí)域和頻域特征，加深對(duì)信號(hào)辦理理論的理解，并為此后嫻熟使用MATLAB進(jìn)行系統(tǒng)的剖析仿真和設(shè)計(jì)確立基礎(chǔ)。此次實(shí)習(xí)課程主假如為了進(jìn)一步熟習(xí)對(duì)matlab軟件的使用，以及學(xué)會(huì)利用matlab對(duì)聲音信號(hào)這類實(shí)質(zhì)問(wèn)題進(jìn)行辦理，將理論應(yīng)用于實(shí)質(zhì)，加深對(duì)它的理解。二、設(shè)計(jì)任務(wù)及指標(biāo)利用MATLAB對(duì)語(yǔ)音信號(hào)進(jìn)行剖析和辦理，收集語(yǔ)音信號(hào)后，利用MATLAB軟件平臺(tái)進(jìn)行頻譜剖析；并對(duì)所收集的語(yǔ)音信號(hào)加入擾亂噪聲，對(duì)加入噪聲的信號(hào)進(jìn)行頻譜剖析，設(shè)計(jì)適合的濾波器濾除噪聲，恢還原信號(hào)。語(yǔ)音信號(hào)的“短時(shí)譜”關(guān)于非安穩(wěn)信號(hào),它是非周期的,頻譜隨時(shí)間連續(xù)變化,所以由傅里葉變換獲得的頻譜沒(méi)法獲知其在各個(gè)時(shí)刻的頻譜特征。假如利用加窗的方法從語(yǔ)音流中拿出此中一個(gè)短斷,再進(jìn)行傅里葉變換,就能夠獲得該語(yǔ)音的短時(shí)譜。。三、設(shè)計(jì)過(guò)程(1).語(yǔ)音信號(hào)的收集及剖析鑒于聲卡進(jìn)行數(shù)字信號(hào)的收集。將話筒插入計(jì)算機(jī)的語(yǔ)音輸入插口上,啟動(dòng)錄音機(jī)。按下錄音按鈕，對(duì)話筒說(shuō)話,說(shuō)完后停止錄音。要保存文件時(shí),利用了計(jì)算機(jī)上的A/D變換器,把模擬的聲音信號(hào)變?yōu)榱耸⒌牧炕说臄?shù)字信號(hào),放音時(shí),它又經(jīng)過(guò)D/A變換器,把保存的數(shù)字?jǐn)?shù)據(jù)恢復(fù)為本來(lái)的模擬的聲音信號(hào)。在Matlab軟件平臺(tái)下能夠利用函數(shù)wavread對(duì)語(yǔ)音信號(hào)進(jìn)行采樣,獲得了聲音數(shù)據(jù)變量x1,同時(shí)把x1的采樣頻次fs=22050Hz和數(shù)據(jù)位Nbits=16Bit放進(jìn)了MATALB的工作空間。圖figure1為原始語(yǔ)音信號(hào)的時(shí)域圖形。從圖中能夠看出在時(shí)域環(huán)境下,信號(hào)體現(xiàn)出6不規(guī)則的信號(hào)峰值。經(jīng)過(guò)freqz函數(shù)繪制原始語(yǔ)音信號(hào)的頻次響應(yīng)圖figure2。而后對(duì)語(yǔ)音信號(hào)進(jìn)行頻譜剖析,在Matlab中能夠利用函數(shù)fft對(duì)信號(hào)行迅速傅里葉變換,獲得信號(hào)的頻譜圖figure3，從圖中能夠看出對(duì)各個(gè)頻點(diǎn)上的隨機(jī)信號(hào)在頻域進(jìn)行抽樣,抽樣頻次為22050Hz。(2).給原始信號(hào)加上一個(gè)高頻噪聲在Matlab中人為設(shè)計(jì)一個(gè)固定頻次5500Hz的噪聲擾亂信號(hào)。噪聲信號(hào)往常為隨機(jī)序列,在本設(shè)計(jì)頂用正弦序列取代,擾亂信號(hào)建立命令函數(shù)為d=[Au*sin(2*pi*5500*t)]',給出的擾亂信號(hào)為一個(gè)正弦信號(hào)，針對(duì)上邊的語(yǔ)音信號(hào),收集了此中一段。再對(duì)噪音信號(hào)進(jìn)行頻譜變換獲得其頻譜圖,從圖中能夠看出擾亂信號(hào)，在4000Hz和6000Hz頻點(diǎn)處有一頂峰,此中5500Hz正是本設(shè)計(jì)所要利用的。(3).設(shè)計(jì)一個(gè)濾波器，濾除高頻噪聲由模擬濾波器變換為數(shù)字濾波器時(shí),采納的是雙線性變換法,它保存的是從模擬到數(shù)字域的系統(tǒng)函數(shù)表示。在濾波器的選用時(shí),因?yàn)樵O(shè)計(jì)方法的重視點(diǎn)不一樣,作出比較是困難的。假如FIR濾波器狀況下,最優(yōu)的設(shè)計(jì)將是橢圓濾波器。用雙線性變換法設(shè)計(jì)低通濾波器。的本設(shè)計(jì)是用雙線性變換法設(shè)計(jì)BW帶阻型濾波器。在MATLAB中,能夠利用函數(shù)butterworth,設(shè)計(jì)FIR濾波器,利用MATLAB中的函數(shù)freqz畫出各濾波器的頻次響應(yīng)。用設(shè)計(jì)好的帶阻濾波器對(duì)含噪語(yǔ)音信號(hào)進(jìn)行濾波,在Matlab中FIR濾波器利用函數(shù)fftfilt對(duì)信號(hào)進(jìn)行濾波,FIR濾波器利用函數(shù)filter對(duì)信號(hào)進(jìn)行濾波。在一個(gè)窗口同時(shí)畫出濾波前后的波形及頻譜。從圖中能夠看出,5500Hz看到的頂峰消逝了,語(yǔ)音信號(hào)與開(kāi)始的同樣,濾波器成功的濾除了擾亂信號(hào)。利用MATLAB對(duì)語(yǔ)音信號(hào)進(jìn)行剖析和辦理，收集語(yǔ)音信號(hào)后，利用MATLAB軟件平臺(tái)進(jìn)行頻譜剖析；并對(duì)所收集的語(yǔ)音信號(hào)加入擾亂噪聲，對(duì)加入噪聲的信號(hào)進(jìn)行頻譜剖析，設(shè)計(jì)適合的濾波器濾除噪聲，恢還原信號(hào)。四、結(jié)論及剖析．用MATLAB對(duì)原始語(yǔ)音信號(hào)進(jìn)行剖析，畫出它的時(shí)域波形和頻譜程序：fs=22050;

%

語(yǔ)音信號(hào)采樣頻次為

22050x1=wavread('C:\Documents

and

桌面

聲音文件.wav');%讀取語(yǔ)音信號(hào)的數(shù)據(jù)，賦給變量x1，sound(x1,22050);%播放語(yǔ)音信號(hào)y1=fft(x1,1024);%對(duì)信號(hào)做1024點(diǎn)FFT變換f=fs*(0:511)/1024;%將0到511，步長(zhǎng)為1的序列的值與fs相乘并除以1024的值，賦值給ffigure(1)%創(chuàng)立圖形窗1plot(x1)%做原始語(yǔ)音信號(hào)的時(shí)域圖形title('原始語(yǔ)音信號(hào)');xlabel('timen');

%x

軸的名字是“

timen

”ylabel('fuzhin');

%y

軸的名字是“

fuzhin

”figure(2)freqz(x1)

%

繪制原始語(yǔ)音信號(hào)的頻次響應(yīng)圖title('頻次響應(yīng)圖')figure(3)subplot(2,1,1);plot(abs(y1(1:512)))

%

%

創(chuàng)立兩行一列畫圖區(qū)間的第1個(gè)畫圖區(qū)間做原始語(yǔ)音信號(hào)的FFT頻譜圖title('原始語(yǔ)音信號(hào)FFT頻譜')subplot(2,1,2);plot(f,abs(y1(1:512)));%abs是絕對(duì)值，plot是直角坐標(biāo)下線性刻度曲線title('原始語(yǔ)音信號(hào)頻譜')xlabel('Hz');ylabel('幅值');時(shí)域波形和頻譜：圖1原始語(yǔ)音信號(hào)圖2語(yǔ)音信號(hào)頻次響應(yīng)圖圖3原始語(yǔ)音信號(hào)FFT與信號(hào)頻譜給原始的語(yǔ)音信號(hào)加上一個(gè)高頻余弦噪聲，頻次為5500hz。對(duì)加噪后的語(yǔ)音進(jìn)行剖析，并畫出其信號(hào)時(shí)域和頻譜圖。程序：fs=22050;x1=wavread('C:\Documentsand桌面聲音文件.wav');%讀取語(yǔ)音信號(hào)的數(shù)據(jù)，賦給變量x1f=fs*(0:511)/1024;%將0到511，步長(zhǎng)為1的序列的值與fs相乘并除以1024的值，賦值給ft=0:1/fs:(length(x1)-1)/fs;%將0到x1的長(zhǎng)度減1后的值除以fs的值，且步長(zhǎng)為1/fs的值，的序列的值，給予tAu=;%噪聲幅值d=[Au*sin(2*pi*5500*t)]';%所加的噪聲是正弦信號(hào)x2=x1+d;%將正弦信號(hào)噪聲加在語(yǔ)音信號(hào)上sound(x2,22050);%播放語(yǔ)音信號(hào)y1=fft(x1,1024);%對(duì)信號(hào)y1做1024點(diǎn)FFT變換y2=fft(x2,1024);%對(duì)信號(hào)y2做1024點(diǎn)FFT變換figure(1);%創(chuàng)立圖形窗1plot(t,x2);%做加噪后的信號(hào)時(shí)域圖形（繪出t對(duì)x2的線性圖）title('加噪后的信號(hào)');xlabel('timen');%x軸的名字是“timen”ylabel('fuzhin');%y軸的名字是“fuzhin”figure(2)%創(chuàng)立圖形窗2subplot(2,1,1);%創(chuàng)立兩行一列畫圖區(qū)間的第1個(gè)畫圖區(qū)間plot(f,abs(y1(1:512)));%做原始語(yǔ)音信號(hào)的頻譜圖title('原始語(yǔ)音信號(hào)頻譜');xlabel('Hz');ylabel('fuzhi');subplot(2,1,2);plot(f,abs(y2(1:512)));

%x%y%%

軸的名字是“timen”軸的名字是“fuzhi”創(chuàng)立兩行一列畫圖區(qū)間的第2個(gè)畫圖區(qū)間做加噪后的語(yǔ)音信號(hào)的頻譜圖（abs是絕對(duì)值，plot是直角坐標(biāo)下線性刻度曲線）title('加噪后的信號(hào)頻譜');xlabel('Hz');ylabel('fuzhi');

%x%y

軸的名字是“軸的名字是“

timenfuzhi

””時(shí)域波形和頻譜：圖4加噪后的信號(hào)時(shí)域圖圖5原始信號(hào)與加噪后信號(hào)頻譜對(duì)照與原始信號(hào)對(duì)照，差別：先原始信號(hào)沒(méi)加噪音從前0到2000有幅值，在4000到6000之間沒(méi)有幅值，可是在加了噪音以后4000到6000之間出現(xiàn)最大幅值12，高出正常值。如圖5原始信號(hào)與加噪后信號(hào)頻譜對(duì)照?qǐng)D所示。設(shè)計(jì)適合的濾波器，濾除高頻噪聲，繪出濾波后的信號(hào)頻域和時(shí)域波形：程序：%用雙線性變換法設(shè)計(jì)巴特沃思低通濾波器wp=*pi;%通帶截止頻次ws=*pi;%阻帶截止頻次Rp=1;%通帶最大衰減（db）Rs=15;%阻帶最大衰減（db）Fs=22050;Ts=1/Fs;wp1=2/Ts*tan(wp/2);%將模擬指標(biāo)變換成數(shù)字指標(biāo)ws1=2/Ts*tan(ws/2);%將模擬指標(biāo)變換成數(shù)字指標(biāo)[N,Wn]=buttord(wp1,ws1,Rp,Rs,'s');%選擇濾波器的最小階數(shù)（估量獲得Butterworth低通濾波器的最小階數(shù)N和3dB截止頻次Wc）[Z,P,K]=buttap(N);%創(chuàng)立Butterworth低通濾波器原型[Bap,Aap]=zp2tf(Z,P,K);%將零極點(diǎn)增益變換為廣泛分子，分母[b,a]=lp2lp(Bap,Aap,Wn);%將廣泛的分子和分母變換為以Wn為截止頻次[bz,az]=bilinear(b,a,Fs);%用雙線性變換法實(shí)現(xiàn)模擬濾波器到數(shù)字濾波器的變換（模擬變換為數(shù)字）[H,W]=freqz(bz,az);%求頻次響應(yīng)figure(9)plot(W*Fs/(2*pi),abs(H))%繪制Butterworth低通濾波器頻次響應(yīng)曲線grid%增添圖格xlabel('頻次／Hz')ylabel('頻次響應(yīng)幅度')title('Butterworth')f1=filter(bz,az,z);%利用上邊已做好的濾波器（已bz和az為特色的濾波器）對(duì)x2信號(hào)進(jìn)行濾波，賦值給f1figure(10)subplot(2,1,1)%創(chuàng)立兩行一列畫圖區(qū)間的第1畫圖區(qū)間plot(t,z)%做加噪后的信號(hào)時(shí)域圖形title('濾波前的時(shí)域波形');subplot(2,1,2)%創(chuàng)立兩行一列畫圖區(qū)間的第2畫圖區(qū)間plot(t,f1);title('濾波后的時(shí)域波形');sound(f1,22050);%播放語(yǔ)音信號(hào)p=length(f1);F0=fft(f1,p);f=0:fs/p:fs*(s-1)/p;figure(11)y2=fft(z,1024);%對(duì)信號(hào)y2做1024點(diǎn)FFT變換subplot(2,1,1);%創(chuàng)立兩行一列畫圖區(qū)間的第1畫圖區(qū)間plot(f,abs(y4));title('濾波前的頻譜')xlabel('Hz');%x軸的名字是“Hz”ylabel('fuzhi');%y軸的名字是“fuzhi”subplot(2,1,2)F1=plot(f,abs(F0));title('濾波后的頻譜')xlabel('Hz');ylabel('fuzhi');時(shí)域波形和頻譜：圖6低通濾波器頻次響應(yīng)曲線圖7濾波前后時(shí)域波形圖圖8濾波前后頻譜圖與原噪音信號(hào)對(duì)照，差別：經(jīng)過(guò)對(duì)照剖析可知，濾波后的輸出波形和原始語(yǔ)音加噪聲信號(hào)的圖形發(fā)生了一些變化，在3000-6000Hz之間的信號(hào)消逝，出現(xiàn)0-1000和7000-8000之間的信號(hào)。濾波后的輸出波形顯然在對(duì)應(yīng)時(shí)間幅度比原語(yǔ)音加噪聲信號(hào)的要小，并且濾波的成效也與濾波器的選擇相關(guān)，能夠看出濾波器的性能差別以及參數(shù)的調(diào)理會(huì)對(duì)濾波器產(chǎn)生必定的影響，并且經(jīng)過(guò)回放能夠發(fā)現(xiàn)濾波前后的聲音有變化.低通濾波后,已很靠近本來(lái)的聲音，人耳幾乎鑒別不出。從頻譜圖中我們還能夠看出聲音的能量信號(hào)主要集中在低頻部分，說(shuō)明高頻語(yǔ)音信號(hào)被濾出，濾波器達(dá)到既定要求。五、設(shè)計(jì)領(lǐng)會(huì)正所謂“紙上得來(lái)終覺(jué)淺，覺(jué)知此事要躬行?！睂W(xué)習(xí)任何知識(shí)，僅從理論上去求知，而不去實(shí)踐、探究是不夠的。所以在學(xué)期末到臨之際，我們迎來(lái)了MATLAB課程設(shè)計(jì)。經(jīng)過(guò)為期一周半的MATLAB課程設(shè)計(jì)，我對(duì)MATLAB這個(gè)仿真軟件有了更進(jìn)一步的認(rèn)識(shí)和認(rèn)識(shí)。在這一周半時(shí)間里，我經(jīng)過(guò)自己探究，查閱資料，并且在老師的指導(dǎo)下達(dá)成了：語(yǔ)音信號(hào)的收集及剖析；給原始信號(hào)加上一個(gè)高頻噪聲；設(shè)計(jì)一個(gè)濾波器，濾除高頻噪聲；并最后將課程設(shè)計(jì)報(bào)告總結(jié)完成。在整個(gè)設(shè)計(jì)過(guò)程中我懂得了很多東西，也培育了獨(dú)立思慮和設(shè)計(jì)的能力，建立了對(duì)知識(shí)應(yīng)用的信心，相信會(huì)對(duì)此后的學(xué)習(xí)工作和生活有特別大的幫助，并且提升了自己的著手實(shí)踐操作能力，使自己充分領(lǐng)會(huì)到了在設(shè)計(jì)過(guò)程中的成功愉悅。固然這個(gè)設(shè)計(jì)做的不怎么好，可是在設(shè)計(jì)過(guò)程中所學(xué)到的東西是此次課程設(shè)計(jì)的最大收獲和財(cái)產(chǎn)，使我終生得益。在沒(méi)有做課程設(shè)計(jì)從前，感覺(jué)課程設(shè)計(jì)不過(guò)對(duì)知識(shí)的純真總結(jié)，可是經(jīng)過(guò)此次課程設(shè)計(jì)發(fā)現(xiàn)自己的見(jiàn)解有點(diǎn)太片面，課程設(shè)計(jì)不單是對(duì)前面所學(xué)知識(shí)的一種查驗(yàn)，也是對(duì)自己能力的一種提升，經(jīng)過(guò)此次課程設(shè)計(jì)使自己理解了本來(lái)的那點(diǎn)知識(shí)是特別短缺的，要學(xué)習(xí)的東西還好多，經(jīng)過(guò)此次課程設(shè)計(jì)，理解學(xué)習(xí)是一個(gè)長(zhǎng)久累積的過(guò)程，在此后的工作和生活中都應(yīng)當(dāng)不停的學(xué)習(xí)，努力提升自己的知識(shí)和綜合素質(zhì)。希望此后像這樣的課程設(shè)計(jì)在多一點(diǎn)。參照文件：姚東.MATLAB及在電子信息課程中的應(yīng)用[M].第二版.北京：北京電子工業(yè)第一版社,2000張智星.MATLAB程序設(shè)計(jì)與應(yīng)用.北京：清華大學(xué)第一版社，2002陳懷琛.MATLAB及其在理工課程中的應(yīng)用指南.西安：西安電子科技大學(xué)第一版社，肖偉、劉忠.MATLAB程序設(shè)計(jì)與應(yīng)用[M].北京：清華大學(xué)第一版社2005程序清單：fs=22050;%語(yǔ)音信號(hào)采樣頻次為22050x1=wavread('C:\Documentsand桌面聲音文件.wav');%讀取語(yǔ)音信號(hào)的數(shù)據(jù)，賦給變量x1，sound(x1,22050);%播放語(yǔ)音信號(hào)y1=fft(x1,1024);%對(duì)信號(hào)做1024點(diǎn)FFT變換f=fs*(0:511)/1024;%將0到511，步長(zhǎng)為1的序列的值與fs相乘并除以1024的值，賦值給ffigure(1)%創(chuàng)立圖形窗1plot(x1)%做原始語(yǔ)音信號(hào)的時(shí)域圖形title('原始語(yǔ)音信號(hào)');xlabel('timen');

%x

軸的名字是“

timen

”ylabel('fuzhin');

%y

軸的名字是“

fuzhin

”figure(2)freqz(x1)%繪制原始語(yǔ)音信號(hào)的頻次響應(yīng)圖title('頻次響應(yīng)圖')figure(3)subplot(2,1,1);

%

創(chuàng)立兩行一列畫圖區(qū)間的第

1個(gè)畫圖區(qū)間plot(abs(y1(1:512)))

%

做原始語(yǔ)音信號(hào)的

FFT頻譜圖title('原始語(yǔ)音信號(hào)FFT頻譜')subplot(2,1,2);plot(f,abs(y1(1:512)));%abs是絕對(duì)值，plot是直角坐標(biāo)下線性刻度曲線title('原始語(yǔ)音信號(hào)頻譜')xlabel('Hz');ylabel('幅值');fs=22050;x1=wavread('C:\Documentsand桌面聲音文件.wav');%讀取語(yǔ)音信號(hào)的數(shù)據(jù)，賦給變量x1f=fs*(0:511)/1024;%將0到511，步長(zhǎng)為1的序列的值與fs相乘并除以1024的值，賦值給ft=0:1/fs:(length(x1)-1)/fs;%將0到x1的長(zhǎng)度減1后的值除以fs的值，且步長(zhǎng)為1/fs的值，的序列的值，給予tAu=;%噪聲幅值d=[Au*sin(2*pi*5500*t)]';%所加的噪聲是正弦信號(hào)x2=x1+d;%將正弦信號(hào)噪聲加在語(yǔ)音信號(hào)上sound(x2,22050);%播放語(yǔ)音信號(hào)y1=fft(x1,1024);%對(duì)信號(hào)y1做1024點(diǎn)FFT變換y2=fft(x2,1024);%對(duì)信號(hào)y2做1024點(diǎn)FFT變換figure(1);%創(chuàng)立圖形窗1plot(t,x2);%做加噪后的信號(hào)時(shí)域圖形（繪出t對(duì)x2的線性圖）title('加噪后的信號(hào)');xlabel('timen');%x軸的名字是“timen”ylabel('fuzhin');%y軸的名字是“fuzhin”figure(2)%創(chuàng)立圖形窗2subplot(2,1,1);%創(chuàng)立兩行一列畫圖區(qū)間的第1個(gè)畫圖區(qū)間plot(f,abs(y1(1:512)));%做原始語(yǔ)音信號(hào)的頻譜圖title('原始語(yǔ)音信號(hào)頻譜');xlabel('Hz');ylabel('fuzhi');subplot(2,1,2);plot(f,abs(y2(1:512)));

%x%y%%

軸的名字是“timen”軸的名字是“fuzhi”創(chuàng)立兩行一列畫圖區(qū)間的第2個(gè)畫圖區(qū)間做加噪后的語(yǔ)音信號(hào)的頻譜圖（abs是絕對(duì)值，plot是直角坐標(biāo)下線性刻度曲線）title('加噪后的信號(hào)頻譜');xlabel('Hz');ylabel('fuzhi');

%x%y

軸的名字是“軸的名字是“

timenfuzhi

””%用雙線性變換法設(shè)計(jì)巴特沃思低通濾波器wp=*pi;%通帶截止頻次ws=*pi;%阻帶截止頻次Rp=1;%通帶最大衰減（db）Rs=15;%阻帶最大衰減（db）Fs=22050;Ts=1/Fs;wp1=2/Ts*tan(wp/2);%將模擬指標(biāo)變換成數(shù)字指標(biāo)ws1=2/Ts*tan(ws/2);%將模擬指標(biāo)變換成數(shù)字指標(biāo)[N,Wn]=buttord(wp1,ws1,Rp,Rs,'s');%選擇濾波