数字语音处理b第二次实验.pdfVIP

数字语音处理 B 2016 –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; end frameCount=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 (nframeLength1) frameSet=frameSet.*w(:,ones(frameCount,1)); end timeAxis=(((1:frameCount)-1)*inc+frameLength/2)/fs; % 求出每帧对应的时间 主函数： clc clear all close 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 ． 语音单帧频谱图 clc clear all close all % load ah ; fs=10000; x2=enframe(data, fs, 20, 0.5,