Vron (i, 2) = diff (exp (polyusa_TF (2) * t), t, i - 1);
Matrix_Vron (i, 3) = diff (exp (real (polyusa_TF (3)) * t) * ...
cos (imag (polyusa_TF (3)) * t), t, i - 1);
Matrix_Vron (i, 4) = diff (exp (real (polyusa_TF (4)) * t) * ...
sin (imag (polyusa_TF (4)) * t), t, i - 1);
Matrix_Vron (i, 5) = diff (exp (polyusa_TF (5) * t), t, i - 1);
end
% Матриця Вронського при t = 0;
Matrix_Vron_t_0 = double (subs (Matrix_Vron, t, 0));
% Матриця Вронського при t = T;
T = 3;
Matrix_Vron_t_T = double (subs (Matrix_Vron, t, T));
% vpa (Matrix_Vron_t_0, 6)
% -------------------------------------------------- ----------------------%
% -------------------------------------------------- ----------------------%
% Визначення невідомих коефіцієнтів C
C_ = inv (Matrix_Vron_t_0) * X_0;
% -------------------------------------------------- ----------------------%
% -------------------------------------------------- ----------------------%
% Знаходження моментних функцій
K_Tt_1 = Subs (K_t, t, T - t);
В
K_Tt = diff (K_t);
K_Tt_2 = subs (K_Tt, t, T - t);
В
K_Ttt = diff (K_Tt);
K_Tt_3 = subs (K_Ttt, t, T - t);
В
K_Tttt = diff (K_Ttt);
K_Tt_4 = subs (K_Tttt, t, T - t);
В
K_Ttttt = diff (K_Tttt);
K_Tt_5 = subs (K_Ttttt, t, T - t);
В
h1_Tt = K_Tt_1
h2_Tt = K_Tt_2
h3_Tt = K_Tt_3
h4_Tt = K_Tt_4
h5_Tt = K_Tt_5
% -------------------------------------------------- ----------------------%
% -------------------------------------------------- ----------------------%
% Знаходження моментів
for i = 1: poryadok
Matrix_a (i) = X_T (i) - C_ '* Matrix_Vron_t_T (i, :)';
end
Matrix_a = Matrix_a '
% -------------------------------------------------- ----------------------%
% -------------------------------------------------- ----------------------%
% -------------------------------------------------- ----------------------%
RETURN = 2;
case '2 'p> disp ('L - Проблема моментів у просторі станів ')
% ------------------------ L - проблема моментів ---------------------- -----%
% ---------------------- У просторі станів ------------------------- -%
% ---------------------------------------------- --------------------------%
Matr_Ex_T = subs (MatrEx, t, T);
% -------------------------------------------------- ----------------------%
% -------------------------------------------------- ----------------------%
% Знаходження моментів
for i = 1: poryadok
Matrix_a (i) = X_T (i) - Matr_Ex_T (i, :) * X_0;
end
Matrix_a = Matrix_a '
% -------------------------------------------------...
