资源描述：
《数字语音处理实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、实验一、用MATLAB实现语音信号的时域分析学院：信息与通信工程学院专业：通信工程班级：通信144学号：2014136410姓名：刘新雨指导教师：崔艳秋1.实验目的观察并验证语音信号的时域特性，理解并掌握典型的语音信号时域分析方法和时域特征，为深入学习语音信号处理的各种应用奠定基础。2.MATLAB程序代码（1）由麦克风采集语音数据，将采集的数据存成WAV文件（采样率为8000Hz），存在本人的文件夹中。所用程序代码为：clear;closeall;Fs=11025;y=wavrecord(5*Fs,Fs,'double');w
2、avwrite(y,'D:\lxy');soundview(y,Fs);（2）读取WAV文件，显示语音信号的波形。所用程序代码为：clear;closeall;x=wavread('D:\lxy.wav');figure;plot(x);axis([0,size(x,1),-0.35,0.3]);title('语音信号时域波形');xlabel('样点数');ylabel('幅度');（3）读取WAV文件，计算并显示语音的短时能量（窗函数为矩形窗，帧长为50,200,400,600）。所用程序代码为：clear;closea
3、ll;x=wavread('d:\k.wav');x=double(x);f1=enframe(x,50,50);energy1=sum(abs(f1),2);subplot(2,2,1);plot(energy1);title('语音信号的短时能量');legend('帧长LEN=50');f2=enframe(x,200,200);energy2=sum(abs(f2),2);subplot(2,2,2);plot(energy2);title('语音信号的短时能量');legend('帧长LEN=200');f3=enf
4、rame(x,400,400);energy3=sum(abs(f3),2);subplot(2,2,3);plot(energy3);title('语音信号的短时能量');legend('帧长LEN=400');9f4=enframe(x,600,600);energy4=sum(abs(f4),2);subplot(2,2,4);plot(energy4);title('语音信号的短时能量');legend('帧长LEN=600');（4）读取WAV文件，计算并显示语音的短时过零率（窗函数为矩形窗，帧长为256，帧移为128
5、），所用程序代码为：clear;closeall;x=wavread('D:\lxy.wav')；x=double(x);LEN=300;INC=150;f=enframe(x,LEN,INC);%计算短时过零率z=zeros(size(f,1),1);difs=0.01;fori=1:size(f,1)s=f(i,:);forj=1:(length(s)-1)ifs(j)*s(j+1)<0&abs(s(j)-s(j+1))>difs;z(i)=z(i)+1;endendendsubplot(2,1,1);plot(x);ax
6、is([0,size(x,1),-0.35,0.3]);title('语音信号时域波形');xlabel('样点数');ylabel('幅度');subplot(2,1,2);plot(z);title('语音信号的短时过零率');xlabel('帧数');ylabel('短时过零率');（5）读取WAV文件，计算浊音与清音的短时自相关函数所用的程序代码为：clear;closeall;x=wavread('D:\lxy.wav');x=double(x);LEN=256;INC=128;f=enframe(x,LEN,INC
7、);ff=f(15,:);N=LEN;R=zeros(1,N);fork=1:NR(k)=sum(ff(k:N).*ff(1:N-k+1));9endfork=1:NR1(k)=R(k)/R(1);endsubplot(2,1,1);plot(ff);axis([0,N,-0.5,0.75])title('加矩形窗的语音帧');xlabel('样点数');ylabel('幅度');subplot(2,1,2);plot(R1);axis([0,N,-1,1])title('加矩形窗的短时自相关函数');xlabel('k');y
8、label('R(k)');clear;closeall;x=wavread('d:\k.wav');x=double(x);LEN=256;INC=128;f=enframe(x,LEN,INC);ff=f(124,:);ff=ff'.*hamming(le
