EEMD程序求助

黎明的夕阳

哪位大神用过论坛中发的EEMD程序吗？我在用的过程中出现各种错误，请大神指点一下，不甚感激
EEMD分析全球温度异常的源程序.zip (4.93 KB, 下载次数: 30)

2013-7-7 09:32 上传

点击文件名下载附件
function allmode=eemd(Y,Nstd,NE)
xsize=length(Y);
dd=1:1:xsize;
Ystd=std(Y);
Y=Y/Ystd;

TNM=fix(log2(xsize))-1;
TNM2=TNM+2;
for kk=1:1:TNM2,
    for ii=1:1:xsize,
        allmode(ii,kk)=0.0;
    end
end

for iii=1:1:NE,
    for i=1:xsize,
        temp=randn(1,1)*Nstd;
        X1(i)=Y(i)+temp;
    end

    for jj=1:1:xsize,
        mode(jj,1) = Y(jj);
    end

    xorigin = X1;
    xend = xorigin;

    nmode = 1;
    while nmode <= TNM,
        xstart = xend;
        iter = 1;

        while iter<=10,
            [spmax, spmin, flag]=extrema(xstart);
            upper= spline(spmax(:,1),spmax(:,2),dd);
            lower= spline(spmin(:,1),spmin(:,2),dd);
            mean_ul = (upper + lower)/2;
            xstart = xstart - mean_ul;
            iter = iter +1;
        end
        xend = xend - xstart;

               nmode=nmode+1;

        for jj=1:1:xsize,
            mode(jj,nmode) = xstart(jj);
        end
    end

    for jj=1:1:xsize,
        mode(jj,nmode+1)=xend(jj);
    end

    allmode=allmode+mode;

end

allmode=allmode/NE;
allmode=allmode*Ystd;

会出现function  未定义啊，然后NE或者Y没定义，不知道怎么弄，如能解惑，奉上金钱和贡献，谢谢各位了

geonuist

可以问群里的hillside，他是这方面的专家，为人十分热情！

zihan6688

function是程序定义的，你可以在这儿下载
http://rcada.ncu.edu.tw/research1_clip_program.htm

wode_q_x

zihan6688 发表于 2013-7-7 11:43
function是程序定义的，你可以在这儿下载
http://rcada.ncu.edu.tw/research1_clip_program.htm

感谢!!!!!今后遇到问题还要多多请教!

xfx

金钱不够              

xfx

金钱不够                                   

zoiujing456

又有好东西