nlsysidenobj.m

这是一个遗传算法的优化软件

function tobj=nlsysidenobj(Ti,popsize,lchrom,ts,rin,rout)
global yout time
tobj=0;
P=100;y_1=0;y_2=0;r_1=0;r_2=0;
for k=1:1:P
   time(k)=k;      
   for i=1:1:6
         hu(i)=logsig(Ti(i)*rin(k)+Ti(6+i)*y_1+Ti(12+i)*y_2+Ti(18+i)*r_1+Ti(24+i)*r_2);
   end;
   r_1=rin(k);r_2=r_1;
   for i=1:1:6
         hy(i)=logsig(hu*Ti(30+i:6:60+i)');
   end;
   yout(k)=hy*Ti(73:78)';
   error(k)=yout(k)-rout(k);
   y_2=y_1;
   y_1=yout(k);   
end
tobj=sum(abs(error));
tobj=tobj+1.0e-6;
tobj=10/tobj;