注意这里面的WAV文件是双声道的,窗函数选的是汉宁窗
clc;
clear all;
close all;
[filename,pathname]=uigetfile('*.wav','请选择语音文件:');
[y,fs,Nbits]=wavread([pathname,filename]); % y里面放的是采样的数值;fs是采样频率;bits表示采样位数
N=length(y); %信号长度
sound(y,fs,Nbits);
framelength=200; %给出帧长framelength,
frameshift=80; %给出帧移frameshift
hanning_window=hanning(framelength); %给出汉宁窗
Y=enframe(y,framelength,frameshift); %分帧
figure(1);plot(Y),xlabel('帧');ylabel('幅值');title('Hanning-分帧加窗后矩阵');
time=(0:N-1)/fs; %计算出信号的时间刻度,即坐标
figure(2);
plot(time,y);title('语音文件的时间波形');ylabel('幅值');xlabel('时间/s');
n=size(Y);
for i=1:n(1)
frame_energy(i)=Y(i,:)*Y(i,:)';
end
figure(3);plot(frame_energy),xlabel('