Курсовая работа: Математическое моделирование в задачах расчета и проектирования систем автоматического управления

g54=g45;

g55=(Xr5(K(1),K(2),K(3),K(4),K(5)+2*h,T)-2*Xr5(K(1),K(2),K(3),K(4),K(5)+h,T)+a)/h^2;

G=[g11, g12, g13, g14, g15; g21, g22, g23, g24, g25; g31, g32 ,g33, g34, g35; g41, g42 ,g43, g44, g45; g51, g52 ,g53, g54, g55;];

%G1=G.^-1;

G1=inv(G);

%построение градиента

gr1=(Xr5(K(1)+h,K(2),K(3),K(4),K(5), T)-a)/h;

gr2=(Xr5(K(1),K(2)+h,K(3),K(4),K(5), T)-a)/h;

gr3=(Xr5(K(1),K(2),K(3)+h,K(4),K(5), T)-a)/h;

gr4=(Xr5(K(1),K(2),K(3),K(4)+h,K(5), T)-a)/h;

gr5=(Xr5(K(1),K(2),K(3),K(4),K(5)+h, T)-a)/h;

grad=[gr1 gr2 gr3 gr4 gr5]';

if(Xr5(K(1),K(2),K(3),K(4),K(5),T)>Xr5(K0(1),K0(2),K0(3),K0(4),K0(5),T))

L=lambdamin(K0,G10,grad0);

end

K0=K;

G10=G1;

grad0=grad;

K=K0+L*G1*grad;

J=Xr5(K(1),K(2),K(3),K(4),K(5),T);

end

disp(K0);

disp(J0);

МетодЗолотогоСечения

function L=lambdamin(K,G1,grad)

Xzs=(-1+sqrt(5))/2; %золотое сечение

a=0;

b=1;

while (abs(b-a) >0.01)