基于matlab的語音處理說明書

1、<p>　　洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手?jǐn)?shù)字語音處理處理上機(jī)報(bào)告洗手洗手洗手洗手指 導(dǎo) 老 師：beef 陳振興 洗手姓 名：beef 洗手班 級 序 號：beef 071112-11 洗手日 期：beef

2、 2014年09月30日 洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手

3、洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗手洗</p><p>　　一、beef題目要求：beef</p><p>　　使用Matlab讀取一段語音信號做如下處理：beef</p><p>　　將語音信號加窗；beef</p><p>　　計(jì)算每一

4、幀信號能量；beef</p><p>　　計(jì)算每一幀信號幅度；beef</p><p>　　計(jì)算過零率；beef</p><p>　　計(jì)算自相關(guān)系數(shù)；beef</p><p>　　二、beef算法設(shè)計(jì)：beef</p><p>　　加窗處理函數(shù)：beef</p><p>　　function [

5、 wav_out,R,C ] = Add_win( wav_in,N )</p><p>　　% 輸入原始語音數(shù)據(jù)和窗寬，吃飯輸出加窗后語音數(shù)據(jù)</p><p>　　[C,R] = size(wav_in);beef % 獲得語言信號矩陣尺寸</p><p>　　wav_out = zeros(C/N/2+1,N

6、);beef % 創(chuàng)建新矩陣，吃飯用于存放結(jié)果</p><p>　　cnt = 1;beef</p><p>　　for i = 1:(N/2):(C-N) % 遍歷原始數(shù)據(jù)，吃飯根據(jù)N加窗</p><p>　　for j = 1:N</p><p&g

7、t;　　wav_out(cnt,j) = wav_in((cnt-1)*N/2+j,1);beef 洗手 end</p><p>　　cnt = cnt + 1;beef</p><p><b>　　end</b></p><p><b>　　end</b></p><p>

8、　　計(jì)算能量函數(shù)：beef</p><p>　　function [ wav_e ] = Calculate_E( wav_in )</p><p>　　% 輸入幀語音數(shù)據(jù)，吃飯輸出幀能量</p><p>　　洗手 [C,R] = size(wav_in);beef % 獲得語言信號矩陣尺寸</p><p>

9、;　　temp = zeros(C,R);beef % 創(chuàng)建新矩陣，吃飯用于緩存數(shù)據(jù)</p><p>　　wav_e = zeros(1,C);beef % 創(chuàng)建新矩陣，吃飯用于存放結(jié)果</p><p>　　for i = 1:C</p><p>　　for j = 1:R</p><

10、p>　　temp(i,j) = wav_in(i,j) * wav_in(i,j);beef % 計(jì)算每個數(shù)據(jù)的能量</p><p><b>　　end</b></p><p><b>　　end</b></p><p>　　for i = 1:C</p><p>　　for j = 1:R

11、</p><p>　　wav_e(1,i) = wav_e(1,i) + temp(i,j);beef % 計(jì)算每幀數(shù)據(jù)能量</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p&g

12、t;<p>　　3、beef計(jì)算幅值函數(shù)</p><p>　　function [ wav_a ] = Calculate_A( wav_in )</p><p>　　% 輸入幀語音數(shù)據(jù)，吃飯輸出幀幅度</p><p>　　[C,R] = size(wav_in);beef % 獲得語言信號矩陣尺寸<

13、/p><p>　　temp = zeros(C,R);beef % 創(chuàng)建新矩陣，吃飯用于緩存數(shù)據(jù)</p><p>　　wav_a = zeros(1,C);beef % 創(chuàng)建新矩陣，吃飯用于存放結(jié)果</p><p>　　for i = 1:C</p><p&

14、gt;　　for j = 1:R</p><p>　　temp(i,j) = abs(wav_in(i,j));beef % 計(jì)算每個數(shù)據(jù)的幅度</p><p><b>　　end</b></p><p><b>　　end</b></p><p>　　for i = 1:C<

15、;/p><p>　　for j = 1:R</p><p>　　wav_a(1,i) = wav_a(1,i) + temp(i,j);beef % 計(jì)算每幀數(shù)據(jù)幅值</p><p><b>　　end</b></p><p><b>　　end</b></p><p>&

16、lt;b>　　end</b></p><p>　　4、beef計(jì)算過零率</p><p>　　function [ wav_rate ] = Calculate_Rate( wav_in )</p><p>　　% 輸入幀語音數(shù)據(jù)，吃飯輸出幀過零率</p><p>　　[C,R] = size(wav_in);beef

17、 % 獲得語言信號矩陣尺寸</p><p>　　wav_rate = zeros(1,C);beef % 創(chuàng)建新矩陣，吃飯用于存放結(jié)果</p><p>　　for i = 1: % 遍歷每一幀數(shù)據(jù)，吃飯將相鄰數(shù)據(jù)相乘</p>

18、;<p>　　for j = 1:R-1</p><p>　　temp = wav_in(i,j) * wav_in(i,j+1);beef</p><p>　　if temp < 0</p><p>　　wav_rate(1,i) = wav_rate(1,i) + 1;beef % 乘積小于零，吃飯則認(rèn)為過零一次</p>

19、<p><b>　　end</b></p><p><b>　　end</b></p><p>　　end 洗手end</p><p>　　5、beef計(jì)算自相關(guān)函數(shù)</p><p>　　function [ wav_sr ] = Calculate_sr( wav_in )</

20、p><p>　　% 輸入幀數(shù)據(jù)，吃飯輸出自相關(guān)函數(shù)</p><p>　　[C,R] = size(wav_in);beef % 獲得語言信號矩陣尺寸</p><p>　　wav_sr = zeros(1,R-1);beef % 創(chuàng)建新矩陣，吃飯用于存放結(jié)果</p><p>　　for k = 1:R-

21、1 % 遍歷每一幀數(shù)據(jù)，吃飯將間隔為k的數(shù)據(jù)相乘</p><p>　　for j = 1:R-k</p><p>　　wav_sr(1,k) = wav_sr(1,k) + wav_in(20,j)*wav_in(20,j+k);beef洗手 end</p><p>　　end 洗手end</p>

22、<p>　　6、beef運(yùn)行腳本文件：beef</p><p>　　clear all;beefclose all;beef</p><p>　　wav_in = wavread('speech_Chinese.wav');beef % 讀取語言信號</p><p>　　[ wav_out,R,C ] = Add_win( wav_

23、in,500 );beef % 加窗處理洗手 plot(wav_in);beefhold on;beef % 繪制原始信號圖像洗手 figure;beef</p><p>　　[ wav_e ] = Calculate_E( wav_out );beef % 計(jì)算能量</p><p>　　[ wav_a ] = Calculate_A

24、( wav_out );beef % 計(jì)算幅值</p><p>　　[ wav_rate ] = Calculate_Rate( wav_out );beef % 計(jì)算過零率</p><p>　　[ wav_sr ] = Calculate_sr( wav_out );beef % 計(jì)算自相關(guān)函數(shù)</p><p>　　plot(wav_e);be

25、efhold on;beef % 繪制能量譜</p><p>　　figure;beef洗手 plot(wav_a);beefhold on;beef % 繪制幅度譜</p><p>　　figure;beef</p><p>　　plot(wav_rate);beefhold on;

26、beef % 繪制過零率曲線</p><p>　　figure;beef</p><p>　　plot(wav_sr);beefhold on;beef % 繪制自相關(guān)函數(shù)</p><p>　　三、beef運(yùn)行效果及分析：beef</p><p>　　原始信號、beef

27、能量譜及幅度譜效果對比（N=500）:</p><p><b>　　原始信號</b></p><p><b>　　能量譜</b></p><p><b>　　幅度譜</b></p><p><b>　　效果分析：beef</b></p>&l

28、t;p>　　1）由上圖可知，吃飯N=500時語言信號分為343幀，吃飯與理論一致；beef</p><p>　　2）能量譜、beef幅度譜波形基本上是原始信號的包絡(luò)，吃飯與理論一致；beef</p><p>　　根據(jù)過零率判斷清音及濁音：beef</p><p><b>　　過零率</b></p><p><

29、b>　　效果分析：beef</b></p><p>　　經(jīng)分析過零率較低部分為濁音，吃飯反之為清音。洗手考慮到濁音部分應(yīng)該是人說話時刻，吃飯說活時刻能量應(yīng)該較高，吃飯通過對比能量譜和過零率曲線，吃飯與理論基本吻合。洗手</p><p>　　3.自相關(guān)函數(shù)效果：beef</p><p><b>　　不同幀的自相關(guān)函數(shù)</b>&l

30、t;/p><p><b>　　效果分析：beef</b></p><p>　　經(jīng)分析，吃飯自相關(guān)函數(shù)程周期性衰減變化，吃飯變化周期即為語音信號周期，吃飯有圖可以看出。洗手</p><p>　　N=1000時重復(fù)實(shí)驗(yàn)效果:</p><p><b>　　能量譜</b></p><p>

31、;<b>　　幅度譜</b></p><p><b>　　過零率</b></p><p><b>　　自相關(guān)函數(shù)</b></p><p><b>　　效果分析：beef</b></p><p>　　由上圖可得，吃飯N=1000時，吃飯信號幀數(shù)有所減少，吃飯

32、但能量譜、beef幅度譜和過零率曲線基本不便，吃飯自相關(guān)函數(shù)橫軸參考坐標(biāo)變大，吃飯但所表達(dá)的周期是不變的，吃飯與理論一致。洗手</p><p>　　四、beef實(shí)驗(yàn)小結(jié)</p><p>　　洗手本次實(shí)驗(yàn)中，吃飯我首次使用Matlab處理語音信號，吃飯使我進(jìn)一步鞏固了Matlab編程技巧，吃飯同時了解了語言處理的基本技巧。洗手然而，吃飯也暴漏了一些問題，吃飯如：beef對一些函數(shù)理解不夠深