云变换算法
浏览数：26 /
时间：2015年06月08日
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
郑重声明：本站内容如果来自互联网及其他传播媒体，其版权均属原媒体及文章作者所有。转载目的在于传递更多信息及用于网络分享，并不代表本站赞同其观点和对其真实性负责，也不构成任何其他建议。
