f.m

来自「本算法为自适应遗传算法」· M 代码 · 共 30 行

function  x=f1(v)
global G
global tn1 tn3 tn5 tp1 tp3 tp5
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) 0;0 v(2)];
G_3=series(G_2,G_1);
G_4=[1 0;0 1];
G=feedback(G_3,G_4);
dc=dcgain(G);
dc1=dc(1,1);
dc2=dc(2,1);
t=0:1:4000;
x=[];
if stable_check(G(1,1))==1
    yn=step(G(1,1),t);
    yp=step(G(2,1),t);
    
    x(1)=(tn5-2*tn3+tn1)^2+(tp5-2*tp3+tp1)^2;
    x(2)=(t(tn5)-t(tn1))/2;
    x(3)=(t(tp5)-t(tp1))/2;
else
    x=[inf inf inf];
end