findselfitm.m

来自「ｍａｔｌａｂ程序」· M 代码 · 共 17 行

function [w,f,x]=findselfitm(t);
v=100;cn=0;i=1;%w=[0.5 0.5 0.5 0.5 0.5 0.5;0.5 0.5 0.5 0.5 0.5 0.5];
while cn<10
    cn=cn+1;
  for k=0.01:0.01:0.3
      [c,f]=selfit1(t/100,5000,k); 
     if f<v,
        v=f;
        w=c;
     end; 
  end;
  t(11,:)=t(10,:).*w(1,:)+w(2,:).*t(9,:); 
  %t(11,:)=t(11,:)*100/sum(t(11,:));
  t(1,:)=[];
end;
x=t;
f=v;