云变换算法

function [cent1,band1,he1,peak1,i]=ybh(err)

 yzi=xlsread('sheet1.xls','B2:B261'); 

zd=max(yzi);

zx=max(yzi);

erro=err;

i=1;

n=0;

table=tabulate(yzi);

for m=1:size(table,1)

    if table(m,2)~=0

        n=n+1;

        x(n)=table(m,1);%第一列所有的值

        %y(n)=table(m,3)/100;%每一行数出现的频率数

        y(n)=table(m,3);%每一行数出现的频率数

        y1(n)=table(m,2);%每一行数出现的频率数

    end

end

 plot(x,y1,'r');

   hold on;

varz=var(yzi);%求方差

%数据规整

while max(y)>=erro

    [my,I]=max(y);

    core=x(I);

    delta=0.1;

    %初始的delta可以设置

    z=y;

    while min(z)>=(erro+0.005)

        delta=delta+0.01;

        %步长可以设置

        z=y-my*exp(-((x-core).^2)/(2*delta^2));

    end

    delta=delta-0.01;

    y=y-my*exp(-((x-core).^2/(2*delta^2)));

    peak1(i)=my;

    cent1(i)=core;

    band1(i)=delta;

    i=i+1;

end

for k=1:(i-1)

    he1(k)=sqrt(varz-band1(k)^2);

end

k0=i-1;

for gg=1:i-1

    for  jsx=cent1(gg)-4*band1(gg):0.001:cent1(gg)+4*band1(gg)

         y=peak1(gg)*exp(-((jsx-cent1(gg)).^2/(2*band1(gg)^2)));

        plot(jsx,y);

        hold on;

    end

end