首页 > 代码库 > 云变换算法
云变换算法
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
云变换算法
声明：以上内容来自用户投稿及互联网公开渠道收集整理发布，本网站不拥有所有权，未作人工编辑处理，也不承担相关法律责任，若内容有误或涉及侵权可进行投诉： 投诉/举报 工作人员会在5个工作日内联系你，一经查实，本站将立刻删除涉嫌侵权内容。