倒谱分析2.docVIP

实验三、用MATLAB实现语音信号的倒谱分析 学院：信息与通信工程学院 专业：通信工程 班级：通信112 学号：******* 姓名: *** 指导教师：*** 实验目的 观察语音信号的倒谱系数，理解并掌握语音信号倒谱分析的原理和语音信号倒谱系数的特点，为深入学习语音信号处理的各种应用奠定基础。 实验内容 （1）由麦克风采集语音数据，将采集的数据存成WAV文件（要求采样率为8000Hz），存在本人的文件夹中。 （2）读取WAV文件，分别计算并显示一帧浊音和一帧清音的原始语音信号、加窗信号、倒谱（要求窗函数为汉宁窗、帧长为256，帧移为128），观察并分析浊音和清音倒谱的差异，并利用浊音的倒谱近似的估计出浊音的基音周期。 （3）读取WAV文件，选取一帧浊音并显示其波形，计算并显示不同窗函数情况下这帧浊音的加窗信号、倒谱（要求窗函数分别为矩形窗和汉宁窗、帧长为256，帧移为128），观察并分析不同的窗函数对倒谱分析的影响。 3. （1）由麦克风采集语音数据，将采集的数据存成WAV文件（要求采样率为8000Hz），存在本人的文件夹中。 供参考的程序代码： clear; close all; Fs=8000; y=wavrecord(5*Fs,Fs,double); wavwrite(y,e:\\my.wav); soundview(y,Fs); （2读取WAV文件，分别计算并显示一帧浊音和一帧清音的原始语音信号、加窗信号、倒谱（要求窗函数为汉宁窗、帧长为256，帧移为128），观察并分析浊音和清音倒谱的差异，并利用浊音的倒谱近似的估计出浊音的基音周期。 供参考的程序代码 clear; close all; x= wavread(e:\\my.wav); % 从文件my.wav读取语音 x = double(x); LEN = 256; INC= 128; f = enframe(x, LEN, INC); % 分帧 ff=f(36,:); % 选取一帧浊音信号 ff1=ff.*hamming(length(ff)); % 加汉宁窗 % 计算倒谱 r=fft(ff1); LogMag=log(abs(r)); Ceps=ifft(LogMag); LEN1=LEN/2; Ceps1=zeros(LEN,1); Ceps1(1: LEN1-1)= flipud(Ceps(2:LEN1)); Ceps1(LEN1:LEN)= Ceps(1:LEN1+1); subplot(3,1,1); plot(ff); % 绘制截取的一帧语音信号 axis([0,LEN,-0.5,0.5]); title(截取的一帧语音信号); xlabel(样点数); ylabel(幅度); subplot(3,2,1); plot(ff1); % 绘制一帧加窗后的语音信号 axis([0,LEN,-0.5,0.5]); title(加汉宁窗的语音帧); xlabel(样点数); ylabel(幅度); subplot(3,2,5); Q=[-LEN1+1:1:LEN1]; plot(Q,Ceps1); % 绘制倒谱 axis([-LEN1+1,LEN1,-2,0.7]); xlabel(倒频/s ); ylabel(倒频谱幅度); title (加汉宁窗的倒谱); x= wavread(e:\\my.wav); % 从文件my.wav读取语音 x = double(x); LEN = 256; INC= 128; f = enframe(x, LEN, INC); % 分帧 ff=f(36,:); % 选取一帧浊音信号 ff1=ff.*hamming(length(ff)); % 加汉宁窗 % 计算倒谱 r=fft(ff1); LogMag=log(abs(r)); Ceps=ifft(LogMag); LEN1=LEN/2; Ceps1=zeros(LEN,1); Ceps1(1: LEN1-1)= flipud(Ceps(2:LEN1)); Ceps1(LEN1:LEN)= Ceps(1:LEN1+1); subplot(3,2,2); plot(ff);