function?[AFPfs]=an_sinsoid（filename）
%正弦模型的分析函數(shù)
%輸入：
%filename——輸入文件名（wav文件）
%輸出：
%A——幅度矩陣
%F——頻率矩陣
%P——相位矩陣
%Version?0.1?，Written?by?徐寧?07/13/2007
%y050922@njupt.edu.cn
%***********************************************
[speechfs]=wavread（filename）;
%以20ms為一幀
framesize=fs*0.02;
%構(gòu)造窗函數(shù)
win=hamming（framesize）;
%歸一化窗函數(shù)
?alpha=sum（win）;
win=win‘/alpha;
????ii=1;
????L=length（speech）;
????frame=floor（L/framesize）;
????A=cell（1frame）;
????F=cell（1frame）;
????P=cell（1frame）;
????while?framesize*ii<=L
????????framespeech=speech（framesize*（ii-1）+1:?framesize*ii）;
????????framespeech=framespeech.*win;
????????%求DFT
????????dftspec1=fft（framespeech1024）;
????????%求對數(shù)幅度譜
????????dftspec2=10*log10（abs（dftspec1（1:512）））;
????????k=findpeaks（dftspec2）;
????????A{ii}=abs（dftspec1（k））;
????????F{ii}=k*fs/1024;
????????P{ii}=angle（dftspec1（k））;
????????ii=ii+1;
????end
end

function?[kv]=findpeaks（xmw）
if?nargin<2
????m=‘?‘;
end
nx=length（x）;
if?any（m==‘v‘）
????x=-x（:）;????????%?invert?x?if?searching?for?valleys
else
????x=x（:）;????????%?force?to?be?a?column?vector
end
dx=x（2:end）-x（1:end-1）;
r=find（dx>0）;
f=find（dx<0）;

if?length（r）>0?&?length（f）>0????%?we?must?have?at?least?one?rise?and?one?fall
????dr=r;
????dr（2:end）=r（2:end）-r（1:end-1）;
????rc=repmat（1nx1）;
????rc（r+1）=1-dr;
????rc（1）=0;
????rs=cumsum（rc）;?%?=?time?since?the?last?rise
????
????df=f;
????df（2:end）=f（2:end）-f（1:end-1）;
????fc=repmat（1nx1）;
????fc（f+1）=1-df;
????fc（1）=0;
????fs=cumsum（fc）;?%?=?time?since?the?last?fall
????
????rp=repmat（-1nx1）;
????rp（[1;?r+1]）=[dr-1;?nx-r（end）-1];
????rq=cumsum（rp）;??%?=?time?to?the?next?rise
????
????fp=repmat（-1nx1）;
????fp（[1;?f+1]）=[df-1;?nx-f（end）-1];
????fq=cumsum（fp）;?%?=?time?to?the?next?fall
????
????k=find（（rs????v=x（k）;
????
????if?any（m==‘q‘）?????????%?do?quadratic?interpolation
????????b=0.5*（x（k+1）-x（k-1））;
????????a=x（k）-b-x（k-1）;
????????j=（a>0）;????????????%?j=0?on?a?plateau
????????v（j）=x（k（j））+0.25*b（j）.^2./a（j）;
????????k（j）=k（j）+0.5*b（j）./a（j）;
????????k（~j）=k（~j）+（fq（k（~j））-rs（k（~j）））/2;????%?add?0.5?to?k?if?plateau?has?an?even?width
????end
????
????%?now?purge?nearby?peaks
????
????if?nargin>2
????????j=find（k（2:end）-k（1:end-1）<=w）;
????????while?any（j）
????????????j=j+（v（j）>=v（j+1））;
????????????k（j）=[];
????????????v（j）=[];
????????????j=find（k（2:end）-k（1:end-1）<=w）;
????????end
????end
else
????k=[];
????v=[];
end
if?any（m==‘v‘）
????v=-v;????%?invert?peaks?if?searching?for?valleys
end
if?~nargout
????if?any（m==‘v‘）
????????x=-x;????%?re-invert?x?if?searching?for?valleys
????????ch=‘v‘;
????else
????????ch=‘^‘;
????end
????plot（1:nxx‘-‘kvch）;
end
end

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????2907??2010-07-27?15:02??an_sinsoid.m

?????文件????????461??2007-07-18?18:10??sy_sinsoid.m

?????文件????????127??2010-07-27?15:00??test.m

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

?????????????????3495????????????????????3