协方差法寻峰函数

在网络上找了半天找不到协方差法的程序，一怒之下自己按照那又臭又长的公式编写了这段函数代码。因为这个方法较为冷门，代码中注释并不多，有问题可以私信我。如果喜欢的话，别忘了回头评个分哦~

function peakLoc = FDPK（dataFWHMMLMRTHR）
% 协方差法寻峰，输入范围、半高宽和阈值，输出区域内峰
% 半高宽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