%男女聲識別程序1
%方法為通過基音頻率的估計，判斷男女聲
%男聲：50～200Hz??女聲：200～450
%窗口長度為900

function?out2=amdf1（filename）
N=900;
%尋找濁音段在給定的語音信號段內尋找到合適的濁音部分
%短時能量函數
function?out=the（v）
??????n=floor（（length（v）-1）/N）;
????????for?i=0:（n-1）
????????????p=0;
????????????for?j=1:N
??????????????ss=p+v（i*N+j）*v（i*N+j）;
??????????????p=ss;
????????????end
????????????sf（i+1）=p;
????????end
????????out=sf;????????
????end
%;濾波函數
????function?out=voicefilter（v）
????[fFsbit]=wavread（filename）;
????[ba]?=?ellip（40.0580[60?N]*2/Fs）;
????[Hw]?=?freqz（ba512）;
????sf?=?filter（bav）;
????out=sf;
????end
%;中心削波
????function??out=center（v）
????????len=length（v）;
????????k=0;
????????for?i=1:len
?????????????if?abs（v（i））>k??k=abs（v（i））;
?????????????end
????????end
????????p=0.5*k;
????????for?i=1:len
????????????if?v（i）>p??sf（i）=v（i）-p;
????????????else?if?v（i）<-p??sf（i）=v（i）+p;
????????????????else??sf（i）=0;
????????????????end
????????????end
????????end
????????out=sf;
????end
%;三電平削波
????function??out=three（v）
????????for?i=1:length（v）
????????????sf（i）=sign（v（i））;
????????end
????????out=sf;
????end

%?;主程序
????????f=wavread（filename）;
????????len=length（f）;

????????k=0;
????????for?i=1:len
????????????if?abs（f（i））>k??k=abs（f（i））;
????????????end
????????end
????????for?i=1:len
????????????l=f（i）;
????????	f（i）=l/k;
????????end
????????o=voicefilter（f）;
????????len=length（o）;
????????nn=the（o）;
????????energy=0;
????????for?i=1:length（nn）
????????????????re=energy;
????????????????energy=re+nn（i）;
????????end
????????average=energy/length（nn）;
???????woman=0;
???????man=0;
????????for?i=1:length（nn）-1
????????????????if?nn（i）>average
????????????????????x=o（（i-1）*N+1:i*N）;
????????????????????y=center（x）;
????????????????????y1=three（y）;???????
????????????????????p=0;
????????????????????for?k=1:length（y）
????????????????????????for?i=1:length（y）-k
????????????????????????????c=p+y（i）*y1（k+i）;
????????????????????????????p=c;
????????????????????????end
????????????????????????r（k）=p;
????????????????????????p=0;
????????????????????end
????????????????????rr=0;
????%????????????????plot（r）;
????????????????????
????????????????????
????????????????????pl=0;
????????????????????for?j=45:length（r）
????????????????????????if?r（j）>rr?pl=j;rr=r（j）;
????????????????????????end
????????????????????end
????????????????????m=man;
????????????????????w=woman;
????????????????????if?pl<200?man=m+1;?????????%if?pl>110?man=m+1;
????????????????????else?woman=w+1;
????????????????????end
????????????????end
????????end
????????????num（1）=man;
????????????num（2）=woman;
%????????????plot（num）;
????????if?man>woman?out2=1;?
????????else?out2=0;????
????????end

end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1995??2013-01-05?09:59??amdf2.m

?????文件????????120??2007-12-26?20:38??fenpin.m

?????文件???????4912??2013-01-05?13:28??jiemian.fig

?????文件???????9150??2013-01-05?13:48??jiemian.m

?????文件????????459??2007-12-25?13:40??simple.m

?????文件???????2874??2013-01-05?10:00??amdf1.m

-----------?---------??----------?-----??----

????????????????19510????????????????????6