Курсовая работа: Численные методы интегрирования и оптимизации сложных систем

J2=int2(Kp2, Kd2, Ki1);

if J2>J1

Kd2=Kd1-h;

J2=int2(Kp2,Kd2,Ki1);

if J2>J1

Kd2=Kd1;

end

end

Ki2=Ki1+h;

J2=int2(Kp2, Kd2, Ki2,h);

if J2>J1

Ki2=Ki1-h;

J2=int2(Kp2,Kd2,Ki2,h);

if J2>J1

Ki2=Ki1;

end

end

h=fibon(a,b,h);

while J2<J1

Kp=Kp1+2*(Kp2-Kp1); Kd=Kd1+2*(Kd2-Kd1); Ki=Ki1+2*(Ki2-Ki1);

J1=J2;

J2=int2(Kp,Kd,Ki,h);

Kp1=Kp2;Kp2=Kp; Kd1=Kd2;Kd2=Kd; Ki1=Ki2;Ki2=Ki;

end

end

disp(Kp)

disp(Kd)

disp(Ki)

int2.m