C = [b0 b1 0 0 0]
% Перевірка
syms s
W_s = collect (simplify (C * (s. * eye (5)-A) ^ (-1) * B), s)
pretty (W_s)
В
% 2. Паралельна декомпозиція
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%
b1 = b1/a5;
b0 = b0/a5;
В В
s1 = 0;
s2 = -6615/487;
s3 = -1022/1747 + 4016/1451 * i;
s4 = -1022/1747 - 4016/1451 * i;
s5 = -415/606;
В
alfa = real (s3);
beta = imag (s3);
В
syms s A B C D E
В
W_s_etal = collect (((b1 * s + b0)/((s-s1) * (s-s2) * ((s + alfa) ^ 2 + beta ^ 2) * (s-s5))), s)
% pretty (W_s_etal)
В
Slag_1 = simplify (collect (A * (s-s2) * ((s + alfa) ^ 2 + beta ^ 2) * (s-s5), s));
Slag_2 = simplify (collect (B * (s-s1) * ((s + alfa) ^ 2 + beta ^ 2) * (s-s5), s));
Slag_3 = simplify (collect (C * (s-s1) * ((s + alfa) ^ 2 + beta ^ 2) * (s-s2), s));
Slag_4 = simplify (collect ((D * s + E) * (s-s1) * (s-s2) * (s-s5), s));
В
Chislit_W_s = collect (Slag_1 + Slag_2 + Slag_3 + Slag_4, s);
В
% Рішення системи лінійних рівнянь
В
MS = Double ([1 1 1 1 0;
6753029497/515578134 -513659/1058682 10560977/850789 4210795/295122 1;
+77456808434995506239663107/+126764366837761533378822144 +1874500571398143988939141/+260296441145300889894912 -3300780600401725219142291/418364246989311991349248 915075/98374 4210795/295122;
+26189071674868424275768861465/253528733675523066757644288 2853037197681682345182805/+520592882290601779789824 +45476725452203201718998205/+418364246989311991349248 0 915075/98374;
+6290947020888109571128085025/+84509577891841022252548096 0 0 0 0])
В
PCH = [0, 0, 0; b1; b0];
В
Koeff = MS ^ (-1) * PCH
В
% Перевірка
MS * [Koeff (1); Koeff (2); Koeff (3); Koeff (4); Koeff (5)];
В
Slag_1 = simplify (collect (Koeff (1) * (s-s2) * ((s + alfa) ^ 2 + beta ^ 2) * (s-s5), s));
Slag_2 = simplify (collect (Koeff (2) * (s-s1) * ((s + alfa) ^ 2 + beta ^ 2) * (s-s5), s));
Slag_3 = simplify (collect (Koeff (3) * (s-s1) * ((s + alfa) ^ 2 + beta ^ 2) * (s-s2), s));
Slag_4 = simplify (collect ((Koeff (4) * s + Koeff (5)) * (s-s1) * (s-s2) * (s-s5), s));
В
Chislit_W_s = collect ((Slag_1 + Slag_2 + Slag_3 + Slag_4), s);
Znamena_W_s = collect ((s-s1) * (s-s2) * ((s + alfa) ^ 2 + beta ^ 2) * (s-s5), s);
В
W_s =
pretty (W_s)
% Розрахунок матриць стану
A = [s1 0 0 0 0;
0 s2 0 0 0;
0 0 0 1 0;
0 0 - (alfa ^ 2 + beta ^ 2) -2 * alfa 0;
0 0 0 0 s5]
В
B = [Koeff (1); Koeff (2); 0; 1; Koeff (3)]
В
C = [1 січня Koeff (5) Koeff (4) 1]
В
% Перевірка
W_s = collect (simplify (C * (s. * eye (5)-A) ^ (-1) * B), s)
pretty (W_s)
В
% ВСЕ Підрахували ВІРНО! br/>
SimplexMetod2.m
function Simple...
