* h4_Tt + ks5 * h5_Tt) ^ 2, t, 0, T)), 50)
% d_v_2 = double (d_v_2)
gamma_v_2 = 1/d_v_2
% gamma_v_2 = double (gamma_v_2)
u = vpa (expand (simplify (gamma_v_2 * (ks1 * h1_Tt + ks2 * h2_Tt + ks3 * h3_Tt + ... p> ks4 * h4_Tt + ks5 * h5_Tt))), 50)
% u = vpa (u, 6)
u_0 = subs (u, t, 0)
u_T = subs (u, t, T)
ezplot (u, [0 T], 1)
hl = legend ('u (t)');
set (hl, 'FontName', 'Courier');
title ('u (t)');
xlabel ('t')
grid on
% -------------------------------------------------- ----------------------%
% ---------------------------------------------- --------------------------%
% Знаходження X
% Обчислення матричної експоненти
MatrEx = simplify (vpa (ilaplace (inv (s * eye (5) - A)), 50));
В
syms t tay
X_svob = MatrEx * X_0;
X_vinyg = int ((subs (MatrEx, t, t - tay)) * B * (subs (u, t, tay)), tay, 0, t);
X_real = X_svob + X_vinyg;
В
save Sostoyaniya X_real u
В
X_real = vpa (expand (simplify (X_real)), 50)
X_real_0 = double (subs (X_real, t, 0))
X_real_T = double (subs (X_real, t, T))
% Похибка X
delta_X_T = double (vpa (X_T - X_real_T, 50))
delta_X_0 = double (vpa (X_0 - X_real_0, 50))
В
% Знаходження Y
for i = 1: poryadok - 1
Y_real (i) = B_ (i, :) * X_real;
end
Y_real = vpa (expand (simplify (Y_real ')), 50)
Y_real_0 = double (subs (Y_real, t, 0))
Y_real_T = double (subs (Y_real, t, T))
% Похибка Y
delta_Y_T = double (vpa (Y_T - Y_real_T, 50))
delta_Y_0 = double (vpa (Y_0 - Y_real_0, 50))
% -------------------------------------------------- ----------------------%
% -------------------------------------------------- ----------------------%
% Обчислення max значень для задачі АКОР
h = 0.01;
tic
tt = 0: h: T;
for i = 1: poryadok
X_max (i) = max (abs (subs (X_real (i), t, tt)));
end
U_max = max (abs (subs (u, t, tt)));
toc
save Sostoyaniya X_max U_max
% -------------------------------------------------- ----------------------%
% -------------------------------------------------- ----------------------%
% Побудова результатів X (t)
ezplot (X_real (1), [0 T], 2)
title ('x_1 (t)');
grid on
В
ezplot (X_real (2), [0 T], 3)
title ('x_2 (t)');
grid on
В
ezplot (X_real (3), [0 T], 4)
title ('x_3 (t)');
grid on
В
ezplot (X_real (4), [0 T], 5)
title ('x_4 (t)');
grid on
В
ezplot (X_real (5), [0 T], 6)
title ('x_5 (t)');
grid on
В
% Побудова результатів Y (t)
ezplot (Y_re...
