数字语音处理B第二次试验

数字语音处理 B2016– 2017 学年第 2 学期电信 1405 龙为花第二次实验一、实验目的1. 掌握短时傅里叶变换的定义，并理解其两种解释2. 理解窗函数对短时傅里叶谱的影响，并能根据需要选择适当的窗函数3. 理解同态信号处理的原理4. 掌握倒谱的定义及计算方法5. 对以上各项能用 Matlab 代码实现, 并能画出相应谱图二、实验内容1． 画出 ah 第 20 帧单帧语音时域图，窗长 20ms, 帧移 10ms, 同时加矩形窗和 hamming 窗分帧函数：function [frameSet, timeAxis,frameCount] = enframe(data, fs,timeperFrame, 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);%计算帧数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 (nframeLength>1) frameSet=frameSet.*w(:,ones(frameCount,1));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('电信 1405 龙为花 0121409341029')图 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(1);xlabel('Hz');ylabel('db');plot...