數(shù)字語(yǔ)音處理B第二次實(shí)驗(yàn)

1、數(shù)字語(yǔ)音處理B2016 2017學(xué)年第2學(xué)期 電信1405 龍為花第二次實(shí)驗(yàn)一、實(shí)驗(yàn)?zāi)康?. 掌握短時(shí)傅里葉變換的定義，并理解其兩種解釋2. 理解窗函數(shù)對(duì)短時(shí)傅里葉譜的影響，并能根據(jù)需要選擇適當(dāng)?shù)拇昂瘮?shù)3. 理解同態(tài)信號(hào)處理的原理4. 掌握倒譜的定義及計(jì)算方法5. 對(duì)以上各項(xiàng)能用Matlab代碼實(shí)現(xiàn), 并能畫(huà)出相應(yīng)譜圖 二、實(shí)驗(yàn)內(nèi)容1 畫(huà)出ah第20幀單幀語(yǔ)音時(shí)域圖，窗長(zhǎng)20ms, 幀移10ms, 同時(shí)加矩形窗和hamming窗分幀函數(shù)：function frameSet, timeAxis,frameCount = enframe(data, fs, timeperFrame, overla

2、pRate, winF) winfunc = str2func(winF);ndata=length(data(:);frameLength=fs*timeperFrame*0.001;nframeLength=length(winfunc(frameLength);%取窗長(zhǎng)overlap=frameLength*overlapRate;inc=frameLength-overlap;if(nframeLength=1) len=frameLength;else len=nframeLength;endframeCount=fix(ndata-overlap)/inc); %計(jì)算幀數(shù)frame

3、Set=zeros(len,frameCount);startIndexf =(0:(frameCount-1)*inc;%每幀在數(shù)據(jù)data中開(kāi)始的位置startIndexs=(1:len);startIndexsT=startIndexs(:);frameSet(:)=data(startIndexsT(:,ones(frameCount,1)+startIndexf(ones(1,len),:);w=winfunc(frameLength);if (nframeLength>1) frameSet=frameSet.*w(:,ones(frameCount,1);endtimeAx

4、is=(1:frameCount)-1)*inc+frameLength/2)/fs;% 求出每幀對(duì)應(yīng)的時(shí)間主函數(shù)：clcclear allclose all%-load ah; x1=enframe(data, fs, 20, 0.5, 'boxcar');x2=enframe(data, fs, 20, 0.5, 'hamming');plot(x1(:,20),'g','LineWidth',2);hold on;plot(x2(:,20),'b','LineWidth',2);legend(

5、'矩形窗信號(hào)','漢明窗信號(hào)');xlabel('時(shí)間（采樣點(diǎn)數(shù)）'); ylabel('幅度');suptitle('電信 1405 龍為花 0121409341029')圖1：語(yǔ)音ah 第20幀時(shí)域波形圖與加hammming窗 10k采樣2 語(yǔ)音單幀頻譜圖clcclear allclose all%-load ah;fs=10000;x2=enframe(data, fs, 20, 0.5, 'boxcar');x1=enframe(data, fs, 20, 0.5, 'hamming

6、');u1=x1(:,20); % 取得一幀數(shù)據(jù)u2=x2(:,20);N=1024;t=(0:N-1)/fs; y1=20*log10(fft(u1,N);y2=20*log10(fft(u2,N);f = (0: 511)*fs/1024;figure(1);xlabel('Hz'); ylabel('db');plot(f,y1(1:512),'r','LineWidth',1);%做原始語(yǔ)音信號(hào)的FFT頻譜legend('漢明窗的頻譜');figure(2);plot(f,y2(1:512),

7、9;b','LineWidth',1); xlabel('Hz'); ylabel('db');legend('矩形窗的頻譜');suptitle('電信 1405 龍為花 0121409341029') 圖2：ah第20幀單幀語(yǔ)音頻譜圖，窗長(zhǎng)20ms, 幀移10ms，矩形窗和hamming窗，10k采樣3 語(yǔ)音頻域功率譜圖clcclear allclose all%-load ah;fs=10000;x2=enframe(data, fs, 20, 0.5, 'boxcar');x1=en

8、frame(data, fs, 20, 0.5, 'hamming');u1=x1(:,20); % 取得一幀數(shù)據(jù)u2=x2(:,20);N=1024;t=(0:N-1)/fs; y1=20*log10(fft(u1,N);y2=20*log10(fft(u1,N).2);f = (0: 511)*fs/1024;figure(1);plot(f,y1(1:512),'b','LineWidth',1.5);%做原始語(yǔ)音信號(hào)的頻譜hold on;plot(f,y2(1:512),'r','LineWidth',1.

9、5);%做原始語(yǔ)音信號(hào)的功率譜 xlabel('Hz'); ylabel('db');suptitle('電信 1405 龍為花 0121409341029') legend('漢明窗的頻譜','漢明窗功率譜');z1=20*log10(fft(u2,N);z2=20*log10(fft(u2,N).2);figure(2);plot(f,z1(1:512),'b','LineWidth',1);%做原始語(yǔ)音信號(hào)的頻譜hold on;plot(f,z2(1:512),'r&#

10、39;,'LineWidth',1);%做原始語(yǔ)音信號(hào)的功率譜 xlabel('Hz'); ylabel('db');suptitle('電信 1405 龍為花 0121409341029') legend('矩形窗的頻譜','矩形窗功率譜');圖3：Ah 第20幀時(shí)域波形圖，10k采樣, 幀長(zhǎng)20ms, hamming窗圖4：Ah 第20幀時(shí)域波形圖，10k采樣, 幀長(zhǎng)20ms, 矩形窗4 語(yǔ)譜圖（1）clc clear close all x,fs=audioread('E:1.wav&

11、#39;); m1=42230; m2=138000; x=x(m1:m2); win=0.01;inc=win/2; winlen=win.*fs; inclen=inc.*fs; amp1 = enframe(filter(1 -0.9375, 1, x),hamming(winlen), inclen)' n=fix(length(x)-winlen+inclen)/inclen); w=winlen/2+1; n2=1:w; freq=(n2-1)*fs/winlen; Y=fft(amp1); Yn=Y(n2,:);clf; frameTime=(1:n)*inclen; i

12、magesc(frameTime,freq,20*log10(abs(Yn)+eps); axis xy; xlabel('時(shí)間 /s)');ylabel('頻率 /Hz'); colormap(jet); suptitle('電信 1405 龍為花 0121409341029')圖5：語(yǔ)音WHUT語(yǔ)譜圖，10k采樣（2）clc clear close all x,fs=audioread('E:1.wav'); m1=42230; m2=138000; x=x(m1:m2); win=0.04;inc=win/2; winlen

13、=win.*fs; inclen=inc.*fs; amp1 = enframe(filter(1 -0.9375, 1, x),hamming(winlen), inclen)' n=fix(length(x)-winlen+inclen)/inclen); w=winlen/2+1; n2=1:w; freq=(n2-1)*fs/winlen; Y=fft(amp1); Yn=Y(n2,:);clf; frameTime=(1:n)*inclen; imagesc(frameTime,freq,20*log10(abs(Yn)+eps); axis xy; xlabel('

14、時(shí)間 /s');ylabel('頻率 /Hz'); colormap(jet); suptitle('電信 1405 龍為花 0121409341029')圖64.2（1）clc clear close all load('should.mat'); fs=10000; win=0.01;inc=win/2; winlen=win.*fs; inclen=inc.*fs; amp1 = enframe(filter(1 -0.9375, 1, data),hamming(winlen), inclen)' n=fix(length

15、(data)-winlen+inclen)/inclen); w=winlen/2+1; n2=1:w; freq=(n2-1)*fs/winlen; Y=fft(amp1); Yn=Y(n2,:);clf; frameTime=(1:n)*inclen; imagesc(frameTime,freq,20*log10(abs(Yn)+eps); axis xy; xlabel('時(shí)間 /s');ylabel('頻率 /Hz'); colormap(jet); suptitle('電信 1405 龍為花 0121409341029')圖7（2）c

16、lc clear close all load('should.mat'); fs=10000; win=0.04;inc=win/2; winlen=win.*fs; inclen=inc.*fs; amp1 = enframe(filter(1 -0.9375, 1, data),hamming(winlen),inclen)' n=fix(length(data)-winlen+inclen)/inclen); w=winlen/2+1; n2=1:w; freq=(n2-1)*fs/winlen; Y=fft(amp1); Yn=Y(n2,:);clf; fra

17、meTime=(1:n)*inclen; imagesc(frameTime,freq,20*log10(abs(Yn)+eps); axis xy; xlabel('時(shí)間 /s');ylabel('頻率 /Hz'); colormap(jet); suptitle('電信 1405 龍為花 0121409341029')圖85 語(yǔ)音倒譜圖clcclear allclose all%-load ah;fs=10000;x2=enframe(data, fs, 40, 0.5, 'boxcar');x1=enframe(data, fs, 40, 0.5, 'hamming');u1=x1(:,15); % 取得一幀數(shù)據(jù)u2=x2(:,15);y1=rceps(u1);y2=rceps(u2);figure(1);xlabel('Hz'); ylabel('幅度');plot(y1,'r');%做原始語(yǔ)音信號(hào)的FFT頻譜hold on;plot(u1);legend('漢明窗的倒譜','漢明窗的時(shí)域