sysidenobj.m

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

function tobj=sysidenobj(Ti,pnum,popsize,lchrom,ts,rin,rout)
global yout time
zh=round(Ti(pnum));
tobj=0;
if pnum==5
    a=Ti([1 2]);
    b=Ti([3 4]);
    differ='yout(k+zh)=-den(2)*y_1-den(3)*y_2+num(2)*rin(k)+num(3)*r_1;';
else
    a=Ti([1 2 3]);
    b=Ti([4 5 6]);
     differ='yout(k+zh)=-den(2)*y_1-den(3)*y_2-den(4)*y_3+num(2)*rin(k)+num(3)*r_1++num(4)*r_2;';
end
sys=tf(a,[1 b]);
dsys=c2d(sys,ts,'z');
[num,den]=tfdata(dsys,'v');
P=100;
y_1=0;y_2=0;y_3=0;r_1=rin(1);r_2=0;r_3=0;yout(1:zh)=zeros(1,zh);
for k=1:1:P
   time(k)=k;      
   eval(differ);
   error(k)=yout(k)-rout(k);
   
   y_3=y_2;y_2=y_1;y_1=yout(k+zh);r_2=r_1;r_1=rin(k);
end
tobj=sum(abs(error));
tobj=tobj+1.0e-6;
tobj=100/tobj;