%學(xué)生實(shí)驗(yàn)--語音信號(hào)時(shí)域波形\短時(shí)幅值和短時(shí)過零率
clear;
y1_sound=wavread（‘J:\文檔\003\hanshen001.wav‘）;%打開聲音信號(hào)
%sound（y1_sound）;
y_length=length（y1_sound）;???????????????????????????????????%聲音信號(hào)的長(zhǎng)度
Ny=80;????????????????????????????????????????????????????%短時(shí)幀長(zhǎng)，10ms
if?y_length>Ny
frame_all=floor（（y_length-（Ny/2））/（Ny/2））;????%frame_moving?equal?to?a?half?of?frame_length?get?number?of?all?frameto?use?“enframe“?too.
????y2_sound=y1_sound（1:（frame_all+1）*（Ny/2））;
????y3=reshape（y2_sound（Ny/2）[]）;
????y4=[y3（:1:frame_all）;y3（:2:（frame_all+1））];
else
????frame_all=1;?
????y4=y1_sound;
end
%%%%%%%加方窗????
for?i=1:（frame_all）?????????????????????????????????????
????m1_y（i）=sum（abs（y4（:i）））;????????????????????????????%計(jì)算短時(shí)幅值
????m2_y（i）=sum（y4（:i）.*y4（:i））;????????????????????????????????%計(jì)算短時(shí)能量
????
????zero_y（i）=sum（abs（sign（y4（2:80i））-sign（y4（1:79i））））/2;???%計(jì)算短時(shí)過零率
end
%%%%%%%加哈明窗
win_h=hamming（Ny）;??????%win_y=y_sound.*win_h;
for?i=