【正文】 。while(pos+tt3)=length(x)y=x(pos:pos+tt31)。plot(t,x)。 %每次移動半個窗 end for w=1:j Q(w)=Zm(230*(w1)+1)/(2*N)。頻率39。采樣點39。短時平均過零率39。短時能量39。F:\xinhaochuli\39。if(i==2) legend(39。其中著重對Matlab軟件進(jìn)行深入學(xué)習(xí)。同時利用sound函數(shù)對語音信號進(jìn)行回放在聽到“信號處理—濾波器”這句話時，伴隨著著明顯的噪聲。被成為可視語言。因此，我們用短時自相關(guān)函數(shù)來測定語音的相似特性。2）可用于區(qū)分濁音變?yōu)榍逡艋蚯逡糇優(yōu)闈嵋舻臅r間（根據(jù)值的變化趨勢）。任何語音信號的分析和處理必須建立在“ 短時” 的基礎(chǔ)上, 即進(jìn)行“ 短時分析”。對語音信號進(jìn)行時頻域分析后，進(jìn)行加白噪聲處理并進(jìn)行了相關(guān)分析，設(shè)計濾波器并運用所設(shè)計的濾波器對加噪信號進(jìn)行濾波, 繪制濾波后信號的時域波形和頻譜。 定義短時平均過零率： （）其中為符號函數(shù)。反之，清音接近于隨機(jī)噪聲，清音的短時自相關(guān)函數(shù)不具有周期性，也沒有明顯突起的峰值，且隨著延時k的增大迅速減小。語譜圖上因其不同的灰度，形成不同的紋路，稱之為“聲紋”。當(dāng)用sound函數(shù)回放時，噪聲消失達(dá)到濾波效果，但是聲音與原始信號比較相對沉悶。希望今后能有跟多的實踐機(jī)會，從中不斷提升自己。)。xlabel(39。 else b(i+1)= 1。)。n=length(a)。 else b(i) = 1。)。波形圖39。%找出 2ms(=500Hz)和 20m(=50Hz)之間的最大值fx=[fx fs/(tt+fxval1)]。F:\xinhaochuli\39。x=filter(b,a,y)。 )。 )。 )。plot(w/ pi,X(1:L))。plot (w/ pi,Y(1:L))。plot (w/ pi,Y(1:L))。y(:,1)=y(:,1) mean(y(:,1))。倒譜圖39。fx=[]。subplot(2,1,1)。 k=k+N/2。xlabel(39。xlabel(39。 %短時平均過零率 end subplot(3,1,3),plot(Q),ylabel(39。),ylabel(39。endend程序二： ；；；；；短時平均過零率；；；；；；；；；；；a=wavread(39。 % 求短時能量函數(shù)En subplot(6,1,1),plot(En)。第一天把查閱的相關(guān)資料進(jìn)行了系統(tǒng)的分析和理解。：，表示白噪聲基本覆蓋了原始信號的頻譜。與離散傅立葉變換逼近傅立葉變換一樣，如令ω=2πk/N，則得離散的短時傅立葉吧如下： () 語譜圖語譜圖反映了語音信號的動態(tài)頻率特性，在語音分析中具有重要的實用價值。濁音的時間波形呈現(xiàn)出一定的周期性，波形之間相似性較好；清音的時間波形呈現(xiàn)出隨機(jī)噪聲的特性，樣點間的相似性較差。值大對應(yīng)于濁音段，值小對應(yīng)于清音段。語音信號可以認(rèn)為在短時間內(nèi)(一般認(rèn)為在 10~30ms 的短時間內(nèi))近似不變， 因而可以將其看作是一個準(zhǔn)穩(wěn)態(tài)過程, 即語音信號具有短時平穩(wěn)性。整體設(shè)計框圖如下圖所示：第一章 語音信號時域分析 語音信號的時域分析可直接對語音信號進(jìn)行時域波形分析，在此只只針對語音信號的短時能量、短時平均過零率、短時自相關(guān)函數(shù)進(jìn)行討論。在矩形窗條件下，可以簡化為： （） 矩形窗（N=320）條件下的短時平均過零率。1）基音頻率的估計a) 可利用時域分析判定某一語音有效的清音和濁音段。聲紋因人而異，因此可以在司法、安全等場合得到應(yīng)用。若要達(dá)到較高的濾波要求可對濾波器的參數(shù)進(jìn)行高精度調(diào)試。再次感謝老師的悉心教導(dǎo)！參考文獻(xiàn)[1] 程佩青. 數(shù)字信號處理教程[M]. 北京: 清華大學(xué)出版社,2002.[2] 劉敏,魏玲. Matlab 通信仿真與應(yīng)用[M]. 北京: 國防工業(yè)出版社, 2001.[3] 樓順天, 劉小東, 的系統(tǒng)分析與設(shè)計[J].: 西安電子科技大學(xué) , 2005.[4] 謝平，王娜，：機(jī)械工業(yè)出版社，2008.[5] . 北京：機(jī)械工業(yè)出版社，2003.附 錄 第 頁 共18頁程序一：；；；；；；短時能量函數(shù)；；；；；；；；；；；；[x,fs,bits]=wavread(39。elseif(i==4