close?all;?
clear?all;?
[xfsbit]=wavread（‘2.wav‘）;%讀語音數據數據放入了Xfs代表采樣率，bit代表采樣位數?
framelength=240;%設置幀長，%30ms?under?8khz?
framenumber=fix（length（x）/（framelength））;%總的數據幀數?
totaltime=length（x）/fs;%該段語音總的時間長度?
sp1=x（fix（1.1*fs）:（fix（1.1*fs）+framelength+framelength-1））;%1.1*fs=8.8秒處為濁音始點，取兩幀?
sp2=x（0.8*fs:（0.8*fs+framelength+framelength-1））;%0.8*fs=6.4秒處為清音始點，取兩幀?
d=0;%初始化，濁音波峰判斷?
R1=0;%初始化，濁音一個基音周期內幀數?
%計算每幀的能量?
for?i=1:framenumber;?
???E（i）=0;%短時能量初始化?
???Z（i）=0;%短時過零率初始化?
???M（i）=0;%短時平均幅度初始化?
???K（i）=0;%短時零能比初始化?
???for?j=（framelength*（i-1）+1）:framelength*i;?
???????E（i）=E（i）+x（j）*x（j）;??%第i幀短時能量?
???end?
???for?j=（framelength*（i-1）+2）:framelength*i;?
???????Z（i）=Z（i）+abs（sign（x（j）-0.2）-sign（x（j-1）））;%第i幀短時過零率?
???end?
???for?j=（framelength*（i-1）+1）:framelength*i;?
???????M（i）=M（i）+abs（x（j））;?
???end?
?M（i）=M（i）/framelength;%第i幀短時平均幅度?
????

end?
???for?i=1:framelength;?
???????R（i）=0;?
???????for?j=1:framelength;?
		??R（i）=R（i）+sp1（j）*sp1（i+j）;%濁音的短時自相關函數?
??????????end?
??????if?（R（i）>=110.11）?
??????????R（i）=110.11;%修正最高波峰最大值為110.11?
??????????d=d+1;%當波峰為110.11時，自加一?
??????end?
??if?（d==1）%以第一個波峰110.11為起始點，下一個波峰110.11為終點