实验一、用 MATLAB 实现语音信号的时域分析 1 1. 实验目的 观察并验证语音信号的时域特性，理解并掌握典型的语音信号时域分析方法和时域特征，为深入学习语音信号处理的各种应用奠定基础。 2. MATLAB 程序代码 （1）由麦克风采集语音数据，将采集的数据存成WAV 文件（采样率为8000Hz），存在本人的文件夹中。所用程序代码为： clear; close all; Fs=8000; y=wavrecord(5*Fs,Fs,'double'); wavwrite(y,'f:\\a'); soundview(y,Fs,'name') （2）读取 WAV 文件，显示语音信号的波形。所用程序代码为： x = wavread('f:\\a.wav'); figure; plot(x); axis([0,size(x,1),-0.35,0.3]); %坐标轴范围 title(' 语音信号时域波形 ');xlabel('样点数 '); ylabel('幅度');（ 3）读取 WAV 文件，计算并显示语音的短时能量（窗函数为矩形窗，帧长为20）。所用程序代码为： x = wavread('f:\\a.wav'); x = double(x); f = enframe(x,20, 100); energy = sum(abs(f), 2); subplot(2,2,1);plot(x); axis([0,size(x,1),-0.7,0.9]); title('语音信号时域波形');xlabel('样点数'); ylabel('幅度'); subplot(2,2,2);plot(energy); title('语音信号的短时能量'); legend('帧长 LEN = 20'); f = enframe(x,20, 200); energy = sum(abs(f), 2); subplot(2,2,3);plot(energy); title('语音信号的短时能量'); legend('帧长 LEN = 100'); f = enframe(x,20, 400); energy = sum(abs(f), 2); subplot(2,2,4);plot(energy); title('语音信号的短时能量'); legend('帧长 LEN = 200'); 2 （4）读取WAV 文件，计算并显示语音的短时过零率（窗函数为矩形窗，帧长为256，帧移为128），所用程序代码为： x = wavread('f:\\a.wav'); x = double(x); LEN =100; INC=100; f = enframe(x, LEN, INC); %分帧 % 计算短时过零率 z = zeros(size(f,1),1); difs =0.01; for i=1:size(f,1) s=f(i,:); for j=1:(length(s)-1) if s(j)* s(j+1)<0&abs(s(j)-s(j+1))>difs; z(i)= z(i)+1; end end end subplot(2,1,1); plot(x); axis([0,size(x,1),-0.35,0.3]); title('语音信号时域波形'); xlabel('样点数'); ylabel('幅度'); subplot(2,1,2);plot(z); title('语音信号的短时过零率');xlabel('帧数'); ylabel('短时过零率'); （5）读取WAV 文件，当窗函数...
