資源簡介
在網絡上找了半天找不到協方差法的程序,一怒之下自己按照那又臭又長的公式編寫了這段函數代碼。因為這個方法較為冷門,代碼中注釋并不多,有問題可以私信我。如果喜歡的話,別忘了回頭評個分哦~
代碼片段和文件信息
function?peakLoc?=?FDPK(dataFWHMMLMRTHR)
%?協方差法尋峰,輸入范圍、半高寬和閾值,輸出區(qū)域內峰
%?半高寬FWHM,下界ML,上界MR,閾值THR(2~5)
????k=floor(FWHM);%常數K
????
????c=@(j)exp(-2.773*(j/FWHM^2));%峰形函數
????g=@(mj)exp(-2*(j/FWHM)^4)/data(m+j);%g函數
????SP=max(k+1ML);%起始尋峰點
????EP=min(max(size(data))-kMR);%終止尋峰點
????rm=[];
????
????for?m=SP:EP
????????G=0;
????????for?j=-k:k
????????????G=G+g(mj);?
????????end
????????GCY=0;
????????for?j=-k:k
????????????GCY=GCY+g
- 上一篇:倒頻譜分析
- 下一篇:地震波S變換matlab程序
評論
共有 條評論