數字語音處理B第二次實驗

1、精選優(yōu)質文檔-傾情為你奉上數字語音處理B2016 2017學年第2學期 電信1405 龍為花第二次實驗一、實驗目的1. 掌握短時傅里葉變換的定義，并理解其兩種解釋2. 理解窗函數對短時傅里葉譜的影響，并能根據需要選擇適當的窗函數3. 理解同態(tài)信號處理的原理4. 掌握倒譜的定義及計算方法5. 對以上各項能用Matlab代碼實現, 并能畫出相應譜圖 二、實驗內容1 畫出ah第20幀單幀語音時域圖，窗長20ms, 幀移10ms, 同時加矩形窗和hamming窗分幀函數：function frameSet, timeAxis,frameCount = enframe(data, fs, timeper

2、Frame, overlapRate, winF) winfunc = str2func(winF);ndata=length(data(:);frameLength=fs*timeperFrame*0.001;nframeLength=length(winfunc(frameLength);%取窗長overlap=frameLength*overlapRate;inc=frameLength-overlap;if(nframeLength=1) len=frameLength;else len=nframeLength;endframeCount=fix(ndata-overlap)/inc

3、); %計算幀數frameSet=zeros(len,frameCount);startIndexf =(0:(frameCount-1)*inc;%每幀在數據data中開始的位置startIndexs=(1:len);startIndexsT=startIndexs(:);frameSet(:)=data(startIndexsT(:,ones(frameCount,1)+startIndexf(ones(1,len),:);w=winfunc(frameLength);if (nframeLength1) frameSet=frameSet.*w(:,ones(frameCount,1);

4、endtimeAxis=(1:frameCount)-1)*inc+frameLength/2)/fs;% 求出每幀對應的時間主函數：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(矩形窗信號,漢明窗信號);xlabel(時間（采樣點數）); ylabel(幅度);suptitle(

5、電信 1405 龍為花 29)圖1：語音ah 第20幀時域波形圖與加hammming窗 10k采樣2 語音單幀頻譜圖clcclear allclose all%-load ah;fs=10000;x2=enframe(data, fs, 20, 0.5, boxcar);x1=enframe(data, fs, 20, 0.5, hamming);u1=x1(:,20); % 取得一幀數據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

6、(1);xlabel(Hz); ylabel(db);plot(f,y1(1:512),r,LineWidth,1);%做原始語音信號的FFT頻譜legend(漢明窗的頻譜);figure(2);plot(f,y2(1:512),b,LineWidth,1); xlabel(Hz); ylabel(db);legend(矩形窗的頻譜);suptitle(電信 1405 龍為花 29) 圖2：ah第20幀單幀語音頻譜圖，窗長20ms, 幀移10ms，矩形窗和hamming窗，10k采樣3 語音頻域功率譜圖clcclear allclose all%-load ah;fs=10000;x2=enf

7、rame(data, fs, 20, 0.5, boxcar);x1=enframe(data, fs, 20, 0.5, hamming);u1=x1(:,20); % 取得一幀數據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);%做原始語音信號的頻譜hold on;plot(f,y2(1:512),r,LineWidth,1.5);%做原始語音信號的功

8、率譜 xlabel(Hz); ylabel(db);suptitle(電信 1405 龍為花 29) legend(漢明窗的頻譜,漢明窗功率譜);z1=20*log10(fft(u2,N);z2=20*log10(fft(u2,N).2);figure(2);plot(f,z1(1:512),b,LineWidth,1);%做原始語音信號的頻譜hold on;plot(f,z2(1:512),r,LineWidth,1);%做原始語音信號的功率譜 xlabel(Hz); ylabel(db);suptitle(電信 1405 龍為花 29) legend(矩形窗的頻譜,矩形窗功率譜);圖3：A

9、h 第20幀時域波形圖，10k采樣, 幀長20ms, hamming窗圖4：Ah 第20幀時域波形圖，10k采樣, 幀長20ms, 矩形窗4 語譜圖（1）clc clear close all x,fs=audioread(E:1.wav); m1=42230; m2=; 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)/inc

10、len); 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(時間 /s);ylabel(頻率 /Hz); colormap(jet); suptitle(電信 1405 龍為花 29)圖5：語音WHUT語譜圖，10k采樣（2）clc clear close all x,fs=audioread(E:1.wav); m1=422

11、30; m2=; x=x(m1:m2); win=0.04;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; imagesc(frameTime,freq,2

12、0*log10(abs(Yn)+eps); axis xy; xlabel(時間 /s);ylabel(頻率 /Hz); colormap(jet); suptitle(電信 1405 龍為花 29)圖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(data)-winle

13、n+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(時間 /s);ylabel(頻率 /Hz); colormap(jet); suptitle(電信 1405 龍為花 29)圖7（2）clc clear close all load(should.mat); fs=10000; win=0.

14、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; frameTime=(1:n)*inclen; imagesc(frameTime,freq,20*log10(abs(Yn)+eps);

15、axis xy; xlabel(時間 /s);ylabel(頻率 /Hz); colormap(jet); suptitle(電信 1405 龍為花 29)圖85 語音倒譜圖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); % 取得一幀數據u2=x2(:,15);y1=rceps(u1);y2=rceps(u2);figure(1);xlabel(Hz); ylabel(幅度);plot(y1,r);%做原始語音信號