findselfitm1.m

ｍａｔｌａｂ程序

function [w,f,x]=findselfitm1(t);
v=100;cn=0;i=1;
while cn<10
    cn=cn+1;
  for k=0.001:0.001:0.3
      [c,f]=selfitm2(t/100,1000,k); 
     if f<v,
        v=f;
        w=c;
     end; 
  end;
  t(10,:)=t(9,:).*w(1,:)+w(2,:).*t(8,:); 
  t(10,:)=t(10,:)*100/sum(t(10,:));
  t(1,:)=[];
end;
x=t;
f=v;