find_remodify_pts.m

免疫算法实现多峰、多极值函数平面曲面拟合

function [idx,pts_x,pts_y]=find_remodify_pts(result2,t_x,t_y,threshold)
% find out the samples that need to be remodified according to the
% threshold set by clients
% -------------------------------------
count=find(result2(:,1));
nodenum=length(count);

% calculate the output of the network
centers=result2(:,2:size(result2,2)-1);
radius=result2(:,size(result2,2));
cou1=find(radius<0.01);
radius(cou1)=0.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);
result1=weights*p1;
error=abs(t_y-result1);

% judge the errors between the network output and the standard output
idx=find(error>threshold);
pts_x=t_x(:,idx);
pts_y=t_y(:,idx);