【正文】 語音信號進行調制再發(fā)送出去 . 下面是對一段語音信號的調制仿真圖，載波頻率為 fc=20xx0。 ④ 設計數字濾波器和畫出頻率響應 根據語音信號的特點給出有關濾波器的性能指標： 低通濾波器性能指標， fp=1000Hz， fc=1200 Hz， As=100dB， Ap=1dB； 濾波器的頻 率響應： 圖 314濾波器的頻率響應 ⑤ 進行低通濾波， 比較濾波前后語音信號的波形及頻譜。最后加以仿真。例如，清濁音的幅度差別、語音信號的端點、語音信號在頻域中的共振峰頻率等。 一個簡單的例子，用 MATLAB 對語音信號最基本的分析就是通過編程對其進行語譜圖分析，從語譜圖中可以很清楚的了解語音信號的特點，語音信號的能量，周期等，這些相當于是將我們平時的不可見的聽覺信號轉換成可見的信號圖 。同時，它與信息科學中最活躍的前沿學科保持密切聯系，并且一起發(fā)展。西安工程大學畢業(yè)設計（論文） 30 參考文獻 [1] 王永琦 . MATLAB 與音視頻技術 . 清華大學 出版社 ,20xx. [2] 胡航 . 現代語音信號處理 .電子工業(yè) 出版社 ,20xx. [3]姚天任 . 數字語音處理 . 華中科技大學 出版社， 1999. [4] 胡航．語音信號處理【 M】．哈爾濱：哈爾濱工業(yè)大學出版社， 20xx． [5] 丁玉美．高西全．數字信號處理【 M】．西安電子科技大學出版社， 20xx． [6] 樊昌信．通信原理【 M】北京：國防工業(yè)出版社， 20xx. [7] 張威． MATLAB 基礎與編程入門【 M】．西安電子科技大學出版社， 20xx． [8] 林福宗．多媒體技術基礎【 M】．北京：清華大學出版社， 20xx． [9] 周淵，王炳和，基于 MATLAB 的噪聲信號采集和分析系統(tǒng)的設計【 J】．噪聲控制． 20xx(7) 5254． [10] 何強，何英． MATLAB 擴展編程【 M】．北京：清華大學出版社， 20xx：293296． [11] ]朱學芳，徐建平．計算機語音信號處理與語音識別系統(tǒng) [J]．南京郵電學院學報， 1998， 18(5— 6)： 113— 11 [12] ]楊崇林，李雪耀，孫羽．強噪聲背景下漢語語音端點檢測和音節(jié)分割 [J]．哈爾濱工程大學學報， 1997． 18(5)： 28— 32． [13] 張雄偉，陳量，楊吉斌．現代語音處理技術及應用 [M]．北京：機械工業(yè)出版社． 20xx [14] 韋曉東，等．應用倒譜特征的帶噪語音端點檢測方法 [J]．上海交通大學學報， 20xx， 34(2)： 185— 188． [15] 胡光銳，韋曉東．基于倒譜特征的帶噪語音端點檢測 [J]．電子學報。 感謝他們對我的培養(yǎng)。D:\39。 subplot(2,1,1) plot(signal) title(39。) 矩形窗 clear all x=linspace(0,100,1000)。矩形窗時域波形 39。幅度 39。 W3=20*log10(abs(W2))。) xlabel(39。) 漢明窗 clear all h=hamming(60)。樣點數 39。 H1=H/H(1)。漢明窗幅度特性 39。幅度 /dB39。 framelength=50。 E(i)=0。 xlabel(39。) legend(39。 for i=1:framenumber。 E(i)=E(i)+framesingal(j)^2。) ylabel(39。) clear E framelength=200。 E(i)=0。 xlabel(39。) legend(39。 for i=1:framenumber。 E(i)=E(i)+framesingal(j)^2。) ylabel(39。) 短時平均過零率 clear all [signal,fs,bit]=wavread(39。 title(39。 西安工程大學本科畢業(yè)設 計（論文） 38 for i=1:framenumber。 Z(i)=Z(i)+abs(sign(framesignal(j))sign(framesignal(j1)))。) 短時自相關函數 clear all [signal,fs,bit]=wavread(39。 unvoiceframe=signal(5001:5000+framelength)。) for i=1:framelength。 end end subplot(3,1,2) plot(RectR)。 unvoiceframe=unvoiceframe*h。 HamR(i)=HamR(i)+unvoiceframe(j)*unvoiceframe(i+j)。) 短時平均幅度差函數 clear all [signal,fs,bit]=wavread(39。 voiceframesignal=signal(1001:1000+framelength*2)。) unvoiceframe=signal(5001:5000+framelength*2)。) for i=1:framelength。 UnvoiceR(i)=UnvoiceR(i)+abs(unvoiceframe(j)unvoiceframe(i+j))。濁音短時幅度差函數 39。) 語音信號的綜合仿真分析程序： [y,fs,bits]=wavread(39。 Y=fft(y,40000)。title(39。plot(abs(Y))。 fc=20xx0。%對原語音信號調制 Y1=fft(y1,40000)。title(39。plot(abs(Y1))。% sound(y1)。 subplot(211)。)。調制信號 +原始信號頻譜 39。 fp=1000。fs1=22050。 [b,a]=ellip(n,Ap,As,wn)。 figure(5)。低通濾波后信號波形 39。title(39。 用雙線性變換法設計的低通濾波器的程序如下： fp=1000。fs=22050。 [b,a]=ellip(n,Ap,As,wn)。 學生簽名： 指導教師簽名： 日 期： 。 西安工程大學本科畢業(yè)設 計（論文） 42 西安工程大學 本科畢業(yè)設計（論文） 誠信聲明 稟承學校優(yōu)良傳統(tǒng)學風，保持我校學生一貫誠信風尚，本人鄭重聲明：所呈交畢業(yè)設計（論文）是在指導老師的指導下獨立完成的，無抄襲和剽竊現象。wp=2*fp/fs。As=100。)。 subplot(212)。plot(x)。 x=filter(b,a,y2)。wp=2*fp/fs1。As=100。 % sound(y2)。plot(abs(Y2))。title(39。 Y2=fft(y2,40000)。調制后信號頻譜 39。)。 subplot(211)。fm39。原始信號頻譜 39。)。 subplot(211)。)。 title(39。 end 西安工程大學本科畢業(yè)設 計（論文） 40 end subplot(2,2,3) plot(VoiceR)。 VoiceR(i)=0。 title(39。 title(39。)。 title(39。 西安工程大學本科畢業(yè)設 計（論文） 39 HamR(i)=0。加矩形窗的清音自相關函數 39。 for j=1:framelengthi。 title(39。)。 title(39。 Z(i)=0。) framelength=150。)。) legend(39。 xlabel(39。 E(i)=0。) clear E framelength=400。) ylabel(39。 E(i)=E(i)+framesingal(j)^2。 for i=1:framenumber。) legend(39。 xlabel(39。 E(i)=0。) clear E framelength=100。) ylabel(39。 E(i)=E(i)+framesingal(j)^2。 for i=1:framenumber。D:\39。歸一化頻率 39。 H3=2*[0:511]/1024。幅度 39。漢明窗時域波形 39。) ylabel(39。 subplot(2,2,2) plot(W4,W3(1:512)) 西安工程大學本科畢業(yè)設 計（論文） 35 title(39。 W=fft(w,1024)。樣點數 39。 h(1001:8000)=1。) subplot(2,1,2) plot(x) title(39。 x=double(signal)。 所援引的著作 論文作者在此一并致以感謝。28(10)： 95— 97． [16] 方佩敏 . 音頻功率放大器 [J]. 電子世界， 20xx 年 08 期 . [13] . Small, ClosedBox Loudspeaker Systems, Part I:Analysis,“ J. Audio ., vol. 20, pp. 798 – 808 (1972 Dec.) [14] Philsp Newell and Keith Holland , Loudspeaker For Music Recording 西安工程大學畢業(yè)設計（論文） 31 and Reproduction, Focal Press is an imprint of Elsevier, First edition 20xx. 西安工程大學本科畢業(yè)設 計（論文） 32 西安工程大學本科畢業(yè)設 計（論文） 33 致 謝 在本次畢業(yè)設計中，我要真誠的感謝崔琳老師給予我的指導，在設計過程老師嚴格督促我的畢業(yè)設計速度，及時的幫我解決我在設計中遇到的各種問題和困難，在設計中一直對我嚴格要求，在整個畢業(yè)設計過程中都給我耐心的指導和講解，并主動為我提供各種相關技術資料，在老師的幫助和指導下我順利完成了本次畢業(yè)設計，也讓我在這次畢業(yè)設計中受益匪淺， 在此謹向老師致以深切的謝意！ 感謝 電信學院 的各位老師在我大學四年學習期間所給予的辛勤教育和幫助。西安工程大學畢業(yè)設計（論文） 29 本章小結： 本章介紹了語音信號時域分析的作用及未來的發(fā)展。語音信號處理的研究目的和處理方法多種多樣，一直是數字信號處理技術的重要推動力量，而數字信號處理很大部分內容也 涉及語音信號處理；數字信號處理學科與技術發(fā)展的一部分來源于數字語音處理的研究。例如，根據語音信號的的短時參數，以及不同語音信號的短時參數的性質對一段給定的信號進行有無聲和清濁音的判斷、對語音信號的基因周期進行估值等。詳細的介紹了語音信號時域分析的全過程。 語音信號的時域參數有短時能量、短時過零率、短時自相關函數和短時平均幅度差函數等。 ③ 用經過調制的信號與原信號疊加。把原始的待傳信號托附到高頻振蕩的過程稱為調制 。每個豎直條相當于一個基音，條紋的起點相當于聲門脈沖的起點，條紋之間的距離表示基音周期 。橫杠是與時間軸平行的幾條深黑色帶紋，它們相應于短時譜中的幾個凸出點，也就是共振峰。雖然對采用語譜圖的講活人識別技術的可靠性過存在相當人的懷疑，但目前這一技術已在司法法庭小得到某些認可及采用 [11]。語語圖的實際應用是用于確定出講話人的本性。帶通濾波器有兩種帶寬可供選擇：窄帶為 45Hz，寬帶為 300Hz。語諾圖是一種三維頻譜，它是表示語音頻譜隨時間變化的圖形，其縱軸為頻率，橫軸為時間，任一給定頻率成分在給定時刻的強弱用相應點的灰度或色調的濃談來表示。但是由于語音信號隨時間變化是很緩慢的，因而在一段短時間內 (如 10～ 30ms 之間，即所謂的一幀之內 )可以認為其頻譜是固定不變的，這種頻譜又稱為短時譜。 下圖為一段語音信