f.asv

本算法为自适应遗传算法

function  x=f(v)
global G ymax1 tmax1 ymin1 tmin1 ymax2 tmax2
G_Nu=tf([68.81 0],[82*12 94 1]);
G_Pu=tf([-23.57 -2.19],[124*3 127 1]);
G_NB=tf([1],[83*83 2*83 1]);
G_PB=tf([2.194],[80*80 160 1]);

%传递函数
G_1=[G_Nu G_NB;G_Pu G_PB];
G_2=[v(1) -1*v(3);v(2) v(4)];
G_3=series(G_2,G_1);
G_4=[1 0;0 1];
G=feedback(G_3,G_4);
[y t]=step(G(1,1));
[ymax1 tmax1]=max(y);
[ymin1 tmin1]=min(y(tmax1+1:end));
if tmax1+tmin1==size(t,1)
    tmax2=inf;
else
    [ymax2 tmax2]=max(y(tmax1+tmin1+1:end));
end
[yp tp]=step(G(2,1));
[ypmax1 tpmax1]=max(yp);
[ypmin1 tpmin1]=min(yp(tpmax1+1:end));
if tpmax1+tpmin1==size(tp,1)
    tpmax2=inf;
else
[ypmax2 tpmax2]=max(yp(tpmax1+tpmin1+1:end));
end
if tmax1+tmin1==size(t,1)
x(1)=(tmax2-tmin1)^2;
x(2)=t(tmax1+tmax2+tmin1)-t(tmax1);
x(3)=tp(tpmax1+tpmax2+tpmin1)-tp(tpmax1);