findselfit.m

ｍａｔｌａｂ程序

function [w,f,x]=findselfit(t);
v=10;cn=0;i=1;
while cn<10
    cn=cn+1;
  for k=0.001:0.001:0.5
      [c,e,f]=selfit1(t(:,i)'/100,5000,k); 
    if f<v,
        v(i)=f;
        w(i,:)=c;
    end; 
  end;
  t(11,:)=w(i,1)*t(10,i)+w(i,2)*t(9,i); 
  t(1)=[];
end;
x=t;
f=v;