= b_0/a_5;
b1 = b_1/a_5;
В
a5 = a_5/a_5;
a4 = a_4/a_5;
a3 = a_3/a_5;
a2 = a_2/a_5;
a1 = a_1/a_5;
a0 = a_0/a_5;
% -------------------------------------------------- ----------------------%
% Порядок системи
poryadok = 5;
% Початкові і кінцеві умови щодо вектора Y
Y_0 = [3 2 1 5] ';
Y_T = [0 -1 0 3] ';
% Кінцевий час переходу
T = 3;
% Матриця переходу від Н.У. Y до Н.У. X
B_ = [B0 b1 0 0 0;
0 b0 b1 0 0;
0 0 b0 b1 0;
0 0 0 b0 b1];
% -------------------------------------------------- ----------------------%
% -------------------------------------------------- ----------------------%
% Початкові умови для впорядкованої системи
X_0 = B_ '* inv (B_ * B_') * Y_0
X_T = B_ '* inv (B_ * B_') * Y_T
% -------------------------------------------------- ----------------------%
% -------------------------------------------------- ----------------------%
% Представлення системи в просторі станів
A = [0 1 0 0 0;
0 0 1 0 0;
0 0 0 1 0
0 0 0 0 1;
-a0 -A1-a2-a3-a4]
B = [0, 0, 0, 0, 1]
C = [B0 b1 0 0 0]
% -------------------------------------------------- ----------------------%
% -------------------------------------------------- ----------------------%
% Обчислення матричної експоненти
syms s t
MatrEx = simplify (vpa (ilaplace (inv (s * eye (5) - A)), 50))
% -------------------------------------------------- ----------------------%
В
RETURN = 1;
В
while RETURN == 1
disp ('L - Проблема моментів у просторі вхід-вихід: 1 ')
disp ('L - Проблема моментів в просторі станів: 2 ')
reply = Input ('Виберете метод рішення [1 або 2]:', 's');
В
switch reply
case '1 'p> disp ('L - Проблема моментів у просторі вхід-вихід ')
% ------------------------ L - проблема моментів ---------------------- -----%
% ---------------------- В просторі вхід-вихід ----------------------- -%
% ---------------------------------------------- --------------------------%
% Передавальна функція
W_obj_s = 1/(a5 * s ^ 5 + a4 * s ^ 4 + a3 * s ^ 3 + a2 * s ^ 2 + a1 * s + a0);
% Полюса передавальної функції
polyusa_TF = roots ([a5 a4 a3 a2 a1 a0]);
% ИПФ
K_t = simplify (vpa (ilaplace (1/(a5 * s ^ 5 + a4 * s ^ 4 + a3 * s ^ 3 + ...
a2 * s ^ 2 + a1 * s + a0)), 50))
% K_t = vpa (K_t, 6)
% -------------------------------------------------- ----------------------%
% Складання матриці Вронського
for i = 1: poryadok
Matrix_Vron (i, 1) = diff (exp (polyusa_TF (1) * t), t, i - 1);
Matrix_...
