function?y1=cut（s_address）

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?function:?cut（）?把靜音段裁剪掉
%?input?:???音頻文件地址
%?output:???裁剪之后的音頻
%?author:???yuhansgg（Shi?Gaige）
%?time:?????2017.4.9
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
y=audioread（s_address）;
h=hamming（320）;

%?計(jì)算短時(shí)平均能量SAE（short?average?energe）
%信號(hào)的平方在與窗函數(shù)相卷
%?E（n）＝[x（m）]^2*h（n-m）m從負(fù)無窮到正無窮求和，h（n-m）為漢明窗
e=conv（y.*yh）;????%?y.*2對(duì)y中各元素平方；conv（uv）?求u與v的卷積

%?對(duì)語音信號(hào)進(jìn)行切割，當(dāng)SAE小于能量大值的1/100時(shí)，認(rèn)為是起點(diǎn)或終點(diǎn)

mx=max（e）;
n=length（e）;
y（n）=0;?%?將原始語音信號(hào)矩陣擴(kuò)充至n維
for?i=1:n
????if?e（i）????????e（i）=0;
????else?e（i）=1;????%?e中非0的數(shù)用1來代替
????end
end
y1=y.*e;
y1（find（y1==0））=[];?%?把0元素剔除
fs=16000;
audiowrite（s_addressy1fs）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件?????????845??2017-04-09?23:25??cut.m
?????文件????????2804??2017-04-09?23:25??PitchDetect.m
?????文件????????1068??2018-01-03?20:07??pitchUI.m
?????文件?????????245??2018-01-03?20:09??ReadMe.txt
?????目錄???????????0??2018-01-03?20:07??test_record\
?????文件?????????698??2017-04-09?23:21??test_record.m
?????文件???????50750??2018-01-03?20:07??test_record\1.wav