首页 > 代码库 > 云变换算法

云变换算法

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

云变换算法