在MATLAB下實(shí)現(xiàn)語音信號(hào)的頻譜分析(論文正文)

1、畢 業(yè) 設(shè) 計(jì)（論 文）基于matlab實(shí)現(xiàn)對(duì)語音信號(hào)的頻譜分析學(xué)生姓名：張彥利學(xué) 號(hào)：所在系部：信息工程學(xué)院專業(yè)班級(jí)：指導(dǎo)教師：日 期：二一貳年五月基于matlab實(shí)現(xiàn)對(duì)語音信號(hào)的頻譜分析摘 要matlab是一個(gè)數(shù)據(jù)分析和處理功能十分強(qiáng)大的工程實(shí)用軟件，運(yùn)用它來進(jìn)行信號(hào)的分析和處理相當(dāng)便捷。本文介紹了在matlab環(huán)境中如何采集語音信號(hào)和語音信號(hào)采集后的頻譜分析處理，并通過實(shí)例分析了語音信號(hào)處理的matlab。并以理想正弦分布的紋理圖像的 matlab仿真和正弦紋理圖像的頻譜特征分析為例，敘述了基于傅立葉能量譜的紋理圖像分析的基本原理和基本過程。編寫了程序，獲得了具有理想正弦分布的空域紋理圖

2、像，并對(duì)其頻譜特征進(jìn)行了分析。并且根據(jù)離散傅里葉變換dft的定義，推導(dǎo)出一種用 dft計(jì)算離散信號(hào)幅值譜的方法，通過 matlab 仿真驗(yàn)證了該方法的有效性。關(guān)鍵詞：matlab；頻譜分析；傅里葉變換；語音信號(hào)；信號(hào)分析analysis of signal spectrum based on matlababstractmatlab is very powerful and practical engineering software which is used in data analysis and processing fraction,and use it to analysis an

3、d process voice signal,which is very convenient.this paper describes the matlab enviironment,how to collect and process speech signal,and through example analysis of the speech signal processing in matlab.then take the matlab simulation of ideal sinusoid texture image and analysis of its spectrum ch

4、aracteristics for example; basic principles and process of analyzing the texture image based on fourier energy spectrum were described. the program was written for obtaining spatial domain texture image with ideal sinusoid distribution. and its spectrum characteristics were analyzed.moreover, accord

5、ing to the definition of dft, a calculation method of amplitude spectrum for periodic signal is derived in this paper. through simulation of amplitude spectrum calculation in matlab, the efficiency of the method is validated to satisfy the need for project. in the end, the problems of spectrum leaka

6、ge and picket fence effect are explained and the corresponding solution is proposed.keywords: matlab; spectrum analysis;fft; speech signal; signal analysis目 錄1. 緒論11.1 課題背景11.2 研究意義11.3 本文研究內(nèi)容12. 頻譜分析及matlab簡介22.1 頻譜分析技術(shù)22.1.1 時(shí)域抽樣定理722.1.2 離散傅立葉變換（dft）822.1.3 快速傅立葉變換（fft）932.1.4 頻譜分析原理1032.2 matlab

7、簡介42.2.1 matlab軟件的發(fā)展42.2.2 matlab組成52.3 本章小結(jié)53. 程序與算例63.1 聲音信號(hào)頻譜分析63.2 圖像信號(hào)頻譜分析73.3 離散信號(hào)/序列93.4 本章小結(jié)11結(jié)論12致謝13參考文獻(xiàn)141. 緒論1.1 課題背景隨著信息時(shí)代和數(shù)字世界的到來，數(shù)字信號(hào)處理己成為當(dāng)今一門極其重要的學(xué)科和技術(shù)領(lǐng)域，數(shù)字信號(hào)處理在通信、語音、圖像、自動(dòng)控制、醫(yī)療和家用電器等眾多領(lǐng)域得到了廣泛的應(yīng)用。任意一個(gè)信號(hào)都具有時(shí)域與頻域特性，信號(hào)的頻譜完全代表了信號(hào)，因而研究信號(hào)的頻譜就等于研究信號(hào)本身。通常從頻域角度對(duì)信號(hào)進(jìn)行分析與處理，容易對(duì)信號(hào)的特性獲得深入的了解。因此，信號(hào)

8、的頻譜分析是數(shù)字信號(hào)處理技術(shù)中的一種較為重要的工具1-3。在工程領(lǐng)域中，matlab是一種倍受程序開發(fā)人員青睞的語言，對(duì)于一些需要做大量數(shù)據(jù)運(yùn)算處理的復(fù)雜應(yīng)用以及某些復(fù)雜的頻譜分析算法matlab顯得游刃有余4。1.2 研究意義信號(hào)處理幾乎涉及到所有的工程技術(shù)領(lǐng)域，而頻譜分析又是信號(hào)處理中一個(gè)非常重要的分析手段。一般的頻譜分析都依靠傳統(tǒng)頻譜分析儀來完成，價(jià)格昂貴，體積龐大，不便于工程技術(shù)人員的攜帶。虛擬頻譜分析儀改變了原有頻譜分析儀的整體設(shè)計(jì)思路，用軟件代替了硬件，使工程技術(shù)人員可以用一部筆記本電腦到現(xiàn)場就可輕松完成信號(hào)的采集、處理及頻譜分析5-6。 1.3 本文研究內(nèi)容信號(hào)的頻譜分析就是利用

9、傅里葉分析的方法，求出與時(shí)域描述相對(duì)應(yīng)的頻域描述,從中找出信號(hào)頻譜的變化規(guī)律，以達(dá)到特征提取的目的。不同信號(hào)的傅里葉分析理論與方法，在有關(guān)專業(yè)書中都有介紹。但實(shí)際的待分析信號(hào)一般沒有解析式，直接利用公式進(jìn)行傅里葉分析非常困難. dft是一種時(shí)域和頻域均離散化的傅里葉變換，適合數(shù)值計(jì)算且有快速算法，是分析信號(hào)的有力工具。dft及fft是數(shù)字信號(hào)處理的重要內(nèi)容。dft是fft的基礎(chǔ)，fft是dft的快速算法，在 mat lab中可以利用函數(shù) fft來計(jì)算序列的離散傅里葉變換dft?；诖耸紫冉榻B了matlab軟件；其次給出了基于matlab軟件的dft和 fft頻譜分析的方法，利用 matlab軟

10、件方法，使得設(shè)計(jì)方便、快捷，大大減輕了工作量；最后結(jié)合實(shí)例給出了仿真結(jié)果。本文將重點(diǎn)介紹基于matlab的頻譜分析設(shè)計(jì)，包括：(1) 音頻信號(hào)頻譜分析；(2) 圖像信號(hào)頻譜分析；(3) 離散信號(hào)/序列頻譜分析。2. 頻譜分析及matlab簡介2.1 頻譜分析技術(shù)2.1.1 時(shí)域抽樣定理7時(shí)域抽樣定理給出了連續(xù)信號(hào)抽樣過程中信號(hào)不失真的約束條件：對(duì)于基帶信號(hào)，信號(hào)抽樣頻率大于等于2倍的信號(hào)最高頻率，即。時(shí)域抽樣是把連續(xù)信號(hào)變成適于數(shù)字系統(tǒng)處理的離散信號(hào)。對(duì)連續(xù)信號(hào)以間隔t抽樣，則可得到的離散序列為。圖2-1 連續(xù)信號(hào)抽樣的離散序列若，則信號(hào)與的頻譜之間存在：其中，的頻譜為，的頻譜為?？梢?，信號(hào)時(shí)

11、域抽樣導(dǎo)致信號(hào)頻譜的周期化。(rad/s)為抽樣角頻率，為抽樣頻率。數(shù)字角頻率與模擬角頻率的關(guān)系為：=t。2.1.2 離散傅立葉變換（dft）8有限長序列的離散傅立葉變換（dft）為逆變換為2.1.3 快速傅立葉變換（fft）9在各種信號(hào)序列中，有限長序列占重要地位。對(duì)有限長序列可以利用離散傅立葉變換(dft)進(jìn)行分析。dft不但可以很好的反映序列的頻譜特性，而且易于用快速算法(fft)在計(jì)算機(jī)上進(jìn)行分析。有限長序列的dft是其z變換在單位圓上的等距離采樣，或者說是序列傅立葉的等距離采樣，因此可以用于序列的譜分析。fft是dft的一種快速算法，它是對(duì)變換式進(jìn)行一次次分解，使其成為若干小數(shù)據(jù)點(diǎn)的

12、組合，從而減少運(yùn)算量。 matlab為計(jì)算數(shù)據(jù)的離散快速傅立葉變換，提供了一系列豐富的數(shù)學(xué)函數(shù)，主要有fft、ifft、fft2 、ifft2, fftn、ifftn和fftshift、ifftshift等。當(dāng)所處理的數(shù)據(jù)的長度為2的冪次時(shí)，采用基-2算法進(jìn)行計(jì)算，計(jì)算速度會(huì)顯著增加。所以，要盡可能使所要處理的數(shù)據(jù)長度為2的冪次或者用添零的方式來添補(bǔ)數(shù)據(jù)使之成為2的冪次。fft函數(shù)調(diào)用方式：y=fft(x)；yfft(x,n)；yfft(x,dim)或yfft(x,n,dim)。函數(shù)ifft的參數(shù)應(yīng)用與函數(shù)fft完全相同。2.1.4 頻譜分析原理10時(shí)域分析只能反映信號(hào)的幅值隨時(shí)間的變化情況，

13、除單頻率分量的簡單波形外，很難明確提示信號(hào)的頻率組成和各頻率分量大小，而頻譜分析能很好的解決此問題。由于從頻域能獲得的主要是頻率信息，所以本節(jié)主要介紹頻率(周期)的估計(jì)與頻譜圖的生成。1、頻率、周期的估計(jì)對(duì)于y(kf)，如果當(dāng)kf=f時(shí)，y(kf)取最大值，則f為頻率的估計(jì)值，由于采樣間隔的誤差，f也存在誤差，其誤差最大為f/ 2。周期t=1/f。從原理上可以看出，如果在標(biāo)準(zhǔn)信號(hào)中混有噪聲，用上述方法仍能夠精確地估計(jì)出原標(biāo)準(zhǔn)信號(hào)的頻率和周期，這個(gè)將在下一章做出驗(yàn)證2、頻譜圖為了直觀地表示信號(hào)的頻率特性，工程上常常將fourier變換的結(jié)果用圖形的方式表示，即頻譜圖。以頻率f為橫坐標(biāo)，|y(f)

14、|為縱坐標(biāo)，可以得到幅值譜；以頻率f為橫坐標(biāo)，argy(f)為縱坐標(biāo)，可以得到相位譜；以頻率f為橫坐標(biāo)，rey(f)為縱坐標(biāo)，可以得到實(shí)頻譜；以頻率f為橫坐標(biāo)，imy(f)為縱坐標(biāo)，可以得到虛頻譜。根據(jù)采樣定理，只有頻率不超過fs/2的信號(hào)才能被正確采集，即fourier變換的結(jié)果中頻率大于fs/2的部分是不正確的部分，故不在頻譜圖中顯示。即橫坐標(biāo)f0,fs/22.2 matlab簡介2.2.1 matlab軟件的發(fā)展matlab 軟件11是由美國 mathworks 公司推出的用于數(shù)值計(jì)算和圖形處理的科學(xué)計(jì)算系統(tǒng)環(huán)境。matlab 是英文 matrix laboratory(矩陣實(shí)驗(yàn)室)的縮

15、寫。它的第一版（dos 版本 1.0）發(fā)行于1984年，經(jīng)過10余年的不斷改進(jìn)，現(xiàn)今已推出它的 windows 98/nt 版本（6.1 版）。新的版本集中了日常數(shù)字處理中的各種功能，包括高效的數(shù)值計(jì)算、矩陣運(yùn)算、信號(hào)處理和圖形生成等功能。在 matlab 環(huán)境下，用戶可以集成地進(jìn)行程序設(shè)計(jì)、數(shù)值計(jì)算、圖形繪制、輸入輸出、文件管理等各項(xiàng)操作。matlab 提供了一個(gè)人機(jī)交互的數(shù)學(xué)系統(tǒng)環(huán)境，該系統(tǒng)的基本數(shù)據(jù)結(jié)構(gòu)是矩陣，在生成矩陣對(duì)象時(shí)，不要求作明確的維數(shù)說明。與利用 c 語言或 fortran 語言作數(shù)值計(jì)算的程序設(shè)計(jì)相比，利用 matlab 可以節(jié)省大量的編程時(shí)間。在美國的一些大學(xué)里，matl

16、ab 正在成為對(duì)數(shù)值計(jì)算、算法預(yù)設(shè)計(jì)與驗(yàn)證，以及一些特殊的短陣計(jì)算應(yīng)用，如自動(dòng)控制理論、統(tǒng)計(jì)、數(shù)字信號(hào)處理（時(shí)間序列分拆）等。matlab 系統(tǒng)最初是由 cievemoler用fortran語言設(shè)計(jì)的，有關(guān)矩陣的算法來自linpack和eispack課題的研究成果；現(xiàn)在的 matlab 程序是 mathworks 公司用 c 語言開發(fā)的，第一版由 steve bangert 主持開發(fā)編譯解釋程序，steve kleiman完成圖形功能的設(shè)計(jì)，john little 和 cleve moler 主持開發(fā)各類數(shù)學(xué)分析的子模塊，撰寫用戶指南和大部分 m 文件。自從第 1 版發(fā)行以來，已有眾多的科技工

17、作者加入到 atlab 的開發(fā)隊(duì)伍中，并為形成今天的matlab 系統(tǒng)做出了巨大的貢獻(xiàn)，matlab 以商品形式出現(xiàn)后，僅短短幾年，就以其良好的開放性和運(yùn)行可靠性，使原先控制領(lǐng)域里的封閉式軟件包（如英國的 umist，瑞典的 lund，德國的 keddc）紛紛淘汰，而改以matlab為平臺(tái)加以重建。在時(shí)間進(jìn)入20世紀(jì)九十年代的時(shí)候，matlab已經(jīng)成為國際控制界公認(rèn)的標(biāo)準(zhǔn)計(jì)算軟件。到九十年代初期，在國際上30 幾個(gè)數(shù)學(xué)類科技應(yīng)用軟件中，matlab 在數(shù)值計(jì)算方面獨(dú)占鰲頭。2.2.2 matlab組成matlab 系統(tǒng)由五個(gè)主要部分組成，下面分別加以介紹：1. matlab 語言體系。matl

18、ab 是高層次的矩陣/數(shù)組語言，具有條件控制、函數(shù)調(diào)用、數(shù)據(jù)結(jié)構(gòu)、輸入輸出、面向?qū)ο蟮瘸绦蛘Z言特性。利用它既可以進(jìn)行小規(guī)模編程，完成算法設(shè)計(jì)和算法實(shí)驗(yàn)的基本任務(wù)，也可以進(jìn)行大規(guī)模編程，開發(fā)復(fù)雜的應(yīng)用程序。2. matlab 工作環(huán)境。這是對(duì) matlab 提高給用戶使用的管理功能的總稱。包括管理工作空間中的變量輸入輸出的方式和方法，以及開發(fā)、調(diào)試、管理 m 文件的各種工具。3. 圖形句相系統(tǒng)。 這是 matlab 圖形系統(tǒng)的基礎(chǔ)，包括完成 2d和 3d 數(shù)據(jù)圖示、圖象處理、動(dòng)畫生成、圖形顯示等功能的高層 matlab命令，也包括用戶對(duì)圖形圖象等對(duì)象進(jìn)行野性控制的低層 matlab 命令，以及開

19、發(fā) gui 應(yīng)用程序的各種工具。4. matlab數(shù)學(xué)函數(shù)庫。這是對(duì) matlab 使用的各種數(shù)學(xué)算法的總稱。包括各種初等函數(shù)的算法，也包括矩陣運(yùn)算、矩陣分析等高層次數(shù)學(xué)算法。5. matlab應(yīng)用程序接口（api）。這是matlab為用戶提供的一個(gè)函數(shù)庫，使得用戶能夠在matlab環(huán)境中使用c程序或 fortran 程序，包括從matlab中調(diào)用程序（動(dòng)態(tài)連接），讀寫 mat 文件的功能?？梢钥闯?matlab 是一個(gè)功能十分強(qiáng)大的系統(tǒng)，是集數(shù)值計(jì)算、圖形管理、程序開發(fā)為一體的環(huán)境。除此之外，matlab 還具有很強(qiáng)的功能擴(kuò)展能力，與它的系統(tǒng)一起，可以配備各種各樣的工具箱。2.3 本章小結(jié)本

20、章詳細(xì)介紹了頻譜分析的基本原理，并對(duì)matlab軟件的發(fā)展和組成展開介紹。3. 程序與算例3.1 聲音信號(hào)頻譜分析matlab語言是一種數(shù)據(jù)分析和處理功能十分強(qiáng)大的計(jì)算機(jī)應(yīng)用軟件， 它可以將聲音文件變換為離散的數(shù)據(jù)文件， 然后利用其強(qiáng)大的矩陣運(yùn)算能力處理數(shù)據(jù)， 如數(shù)字濾波、傅里葉變換、時(shí)域和頻域分析、聲音回放以及各種分析圖的呈現(xiàn)等等。 下面以語音信號(hào)的波形圖、頻譜圖和倒譜圖分析為例來說明matlab在語音信號(hào)處理中的具體實(shí)現(xiàn)方法。程序代碼：s=wavread(1.wav)subplot(3,1,1);plot(s);title(波形圖)n=size(s);d=100;l=512;k=l-d;t

21、=fix(n/k);frame_length=512;r=(rectwin(frame_length);for a=1:t n1=(l-d)*(a-1)+1; n2=(l-d)*(a-1)+frame_length; y=s(n1:n2); sf=fft(y.*r,512); x1(n1:n2)=sf(1:frame_length);endfftsize=8000;y=zeros(fftsize,1);y=20*log10(abs(x1);subplot(3,1,2);plot(y(1:4000),title(頻譜圖);t1=length(x1);t1=fix(t1/l);for i=1:t1

22、z=ifft(x1(i-1)*l+1:i*l),512);z=real(z./r);y(i-1)*k+1:i*k)=z(1:k);end subplot(3,1,3);plot(y),title(倒譜圖)上程序?qū)⒄Z音文件1.wav進(jìn)行頻譜分析，分析結(jié)果如圖3-1所示。圖3-1 聲音信號(hào)的頻譜分析結(jié)果3.2 圖像信號(hào)頻譜分析紋理圖像的頻譜可以通過離散傅里葉變換（dft）得到。用表示一幅空域紋理圖像，用表示該圖像的頻譜，圖像的大小為mn，則和質(zhì)檢可以通過dft計(jì)算，計(jì)算公式如下：其中能量譜可采用公式：基于傅立葉能量譜的紋理圖像分析的前提是假設(shè)紋理有不同的正弦波組成。理想正弦分布的紋理圖像，是最為典

23、型的紋理圖像之一，下面討論理想正弦分布的紋理圖像的仿真及其頻譜特征分析。編寫下面的程序獲得具有理想正弦分布的空域紋理圖像，其中a為正弦紋理的幅值，uf0、vf0分別為x軸（垂直方向）、y軸（水平方向）方向的模擬頻率，m、n分別為x軸、y軸的采樣的點(diǎn)數(shù)，tsu、tsv分別為x、y軸的采樣間隔，為了便于觀察和處理，取tsu=1/m、tsv= 1/n，即x軸、y軸的采樣頻率分別為m和n，這樣在空域中得到了01范圍的紋理圖像（如考慮成時(shí)域抽樣信號(hào)的話，相當(dāng)于在01s間的抽樣）。 a=1;uf0=0;vf0=25;m=200;n=200;tsu = 1/m; tsv = 1/n; r = 0:m-1;

24、c = 0:n-1; c,r = meshgrid(c,r); g = a * sin(2* pi* uf0 * r * tsu + 2 * pi * vf0 * c * tsv); f = mat2gray(g); figure imshow(f)程序中輸出變量g返回的是實(shí)際的理想正弦函數(shù)的取值，f返回的是g平移后的結(jié)果（取值限定在01范圍內(nèi)）。 在matlab軟件中仿真得到的一組具有理想正弦分布的紋理圖像：圖3-2中各紋理的方向或周期是不同的。其中(a)和(b)中紋理呈豎直方向分布（即水平方向變化），(a) 的y軸方向的頻率為25hz，(b)的y軸頻率為5hz，即(b)圖是周期為(a)圖的

25、5倍的豎直紋理，所以在整個(gè)圖像范圍內(nèi)（相當(dāng)于一秒內(nèi)）(a)和(b)分別有25和5個(gè)周期。 (c)圖紋理是x和y軸方向頻率均為25hz。 (a) (b) (c)圖3-2 理想正弦紋理圖像然后利用fft算法對(duì)以上面三個(gè)圖像信號(hào)進(jìn)行頻譜分析。程序代碼為：i=imread(1.tif)i = rgb2gray( i );imshow(i);ffti=fft2(i);sffti=fftshift(ffti);rr=real(sffti);ii=imag(sffti);a=sqrt(rr.2+ii.2);a=(a-min(min(a)/(max(max(a)-min(min(a)*225;figure;i

26、mshow(a);圖3-3為對(duì)應(yīng)的圖中各圖的傅立葉頻譜能量圖（在matlab中執(zhí)行了fft（快速傅立葉變換）后，使用了fftshift函數(shù)調(diào)整，以使頻譜圖像的原點(diǎn)從起始點(diǎn)（0，0），移到圖像的中心點(diǎn)（m/2,n/2）。從圖3-3中可以看出：豎直方向理想單一頻率的正弦分布紋理的頻譜能量集中在水平方向的三個(gè)點(diǎn)上，而水平方向單一頻率的正弦分布紋理的頻譜能量集中在豎直方向的三個(gè)點(diǎn)上；對(duì)于單一頻率的理想正弦分布的紋理圖像，其紋理方向和其頻譜能量圖的頻點(diǎn)連線方向相垂直，如圖 (c)紋理為圖 (b)的紋理順時(shí)針傾斜45度，其頻譜能量圖中頻點(diǎn)也相應(yīng)順時(shí)針傾斜45度；垂直方向紋理和水平方向紋理的疊加圖像的頻譜能

27、量圖表現(xiàn)為水平和豎直方向的五個(gè)點(diǎn)上；紋理周期越小，對(duì)應(yīng)的頻點(diǎn)到頻譜中心（中心點(diǎn)，即頻率為0）的距離越大，如圖3-3 (a)紋理的頻點(diǎn)到中心的距離為圖3-3(b)紋理的五倍。(a) (b) (c)圖 理想正弦紋理圖像的頻譜能量圖3.3 離散信號(hào)/序列以，采樣頻率100hz，采樣100個(gè)點(diǎn)，形成離散信號(hào)，利用matlab進(jìn)行頻譜分析。程序如下：n=100;fs=100;t=(0:(n-1)/fs; xn=1+2*sin(2*pi*10*t)+3*sin(2*pi*15*t); xk=abs(fft(xn,n);f=(0:n/2)*fs/n; xk(1)=xk(1)/n;xk(2:(n/2+1)=x

28、k(2:(n/2+1)*2/n; stem(f,xk(1:(n/2+1);axis(-1 n/2 0 5); grid on;xlabel(f(hz);ylabel(|x( f )|);程序運(yùn)行結(jié)果：圖3-4 n=100的幅值頻譜圖從圖3-4可以看出，dft法分析的結(jié)果和實(shí)際信號(hào)吻合得很好，說明該方法確實(shí)有效。但是，需要注意的是，離散傅里葉變換在頻域是離散的，即限制在基頻整數(shù)倍上，只能得到信號(hào) k 次諧波成分的幅值譜，而對(duì)于非k次諧波成分的頻譜則無法檢測出來，并且由于柵欄效應(yīng)和頻譜泄漏，對(duì)其它k次諧波的幅值也有影響。如將上述程序中的n=100改為n=50，結(jié)果如圖3-5所示。圖3-5 n=50

29、的幅值頻譜圖從圖3-5可以看出，由于基頻頻率為100/50=2hz，x(t)中的15hz分量未檢測出來，其它頻率分量的幅值也出現(xiàn)了誤差，這種柵欄效應(yīng)也是 dft 應(yīng)用中不可避免的問題之一。要減小柵欄效應(yīng)和頻譜泄漏，可以減小采樣頻率，增加采樣點(diǎn)數(shù)，以減小基頻值，使譜線變密，這樣原來漏掉的某些頻譜成分就可以檢測出來，但注意采樣頻率必須滿足采樣定理的要求，增加采樣點(diǎn)數(shù)也會(huì)增加系統(tǒng)計(jì)算 dft 的開銷。3.4 本章小結(jié)本章利用matlab分別對(duì)聲音信號(hào)、圖像信號(hào)和離散信號(hào)/序列開展頻譜分析，得到頻譜分析結(jié)果，并開展結(jié)果分析。結(jié)論本文主要是基于高級(jí)面向?qū)ο箝_發(fā)語言matlab的基本特征，以及matlab

30、強(qiáng)大的工具箱功能，實(shí)現(xiàn)信號(hào)的預(yù)處理和頻譜分析。（1）通過實(shí)例說明，基于matlab可以很好的達(dá)到對(duì)語音信號(hào)的頻譜的分析處理。（2）采用頻譜分析方法可以提取出紋理基本特征，如方向特征和周期特征。本文以理想正弦分布的紋理圖像的matlab仿真和正弦紋理圖像的頻譜特征分析為例，敘述了基于傅立葉能量譜的紋理圖像分析的基本原理和基本過程。（3）通過離散傅里葉變換對(duì)和周期信號(hào)傅立葉級(jí)數(shù)的定義，清楚地推導(dǎo)出 x(k)與信號(hào)幅值的關(guān)系，提出了一種有效的計(jì)算離散信號(hào)/序列幅值譜的方法。致謝本論文從擬定題目到定稿，從對(duì)論文設(shè)計(jì)原理一無所知到模糊了解再到逐漸清晰，此過程歷時(shí)數(shù)月，我也經(jīng)歷了千辛萬苦，但最終有一個(gè)結(jié)果，我感到很欣慰。通過本次設(shè)計(jì)，我感覺自己成長了很多，對(duì)自己所學(xué)的知識(shí)有了更深刻的了