find_change.m

来自「免疫算法实现多峰、多极值函数平面曲面拟合」· M 代码 · 共 30 行

function index=find_change(result2,t_x,t_y,t_x1,t_y1)
count=find(result2(:,1));
centers=result2(count,2:(size(result2,2)-1));
radius=result2(count,size(result2,2));
cou1=find(radius<.01);
radius(cou1)=.01;

r1=1./(radius*sqrt(2));
r2=r1*ones(1,size(t_x,2));
p1=radbas(dist(centers,t_x).*r2);

weights=t_y*pinv(p1);
result=weights*p1;
error1=t_y-result;

count=find(result2(:,1));
centers=result2(count,2:(size(result2,2)-1));
radius=result2(count,size(result2,2));
cou1=find(radius<.01);
radius(cou1)=.01;

r1=1./(radius*sqrt(2));
r2=r1*ones(1,size(t_x1,2));
p1=radbas(dist(centers,t_x1).*r2);

weights=t_y1*pinv(p1);
result=weights*p1;
error2=t_y1-result;
errorchange=abs(abs(error1)-abs(error2));
index=find(errorchange>0.0025);