ypr * Tem * Tya + Tkz * Tem * Tya + Tkz * Typr * Tya: 9:8); writeln ('a5 = ', Tkz * Typr * Tem * Tya: 9:8);
end;
Обнулення елементів стоять під головною діагоналлю
procedure Obnylenie (m: integer; a: massiv; var b: massiv); i, j: integer;; for i: = m +1 to n do for j: = m +1 to n do b [i, j]: = a [i, j]-a ​​[i, m]/a [m, m] * a [m, j];;
Обчислення визначників
procedure Diag_opred (b: massiv; var sum: real); i, j: integer; clrscr; Vivod; writeln; textcolor (9); writeln ('Обчислені визначники:'); sum: = 1; for i: = 1 to n do for j: = 1 to n do if i = j then begin sum: = sum * b [i, i]; textcolor (7); writeln ('', sum: 7: 9); end;;
Формування визначника Гурвіца
procedure Matrix (Kv: real; var a: massiv); clrscr; for i: = 1 to n-3 do begin a [2 * i, i +2]: = 1; a [ 2 * i, i +1]: = Tem * Tya + Typr * Tya + Tkz * Tya + Tkz * Typr; a [2 * i, i]: = Tkz * Typr * Tem * Tya; end; for i: = 1 to n-2 do begin a [2 * i-1, i +2]: = Kv; a [2 * i-1, i +1]: = Tya + Typr + Tkz; a [2 * i-1 , i]: = Typr * Tem * Tya + Tkz * Tem * Tya + Tkz * Typr * Tya; end;;
Основне тіло програми
begin Kv: = 0; repeat begin Kv: = Kv +0.05; Matrix (Kv, a); end; b: = a; for i: = 1 to n-1 do Obnylenie (i , b, b); begin Diag_opred (b, diag); end; until diag <= 0; begin writeln; textcolor (9); writeln ('Критичний коефіцієнт дорівнює:'); textcolor (7); writeln ('Kv = ', Kv: 5:2); readkey; end;
end.
3.1 Опис змінних
У програмі використані такі змінні: a-масив чисел, з елементів якого формуємо визначник Гурвіца; n-розмірність масиву; i-номер рядка; j-номер стовпчика; b-масив чисел. За допомогою його, змінюючи номери рядків і стовпців наводимо матрицю до трикутного вигляду; Tem-електро-механічна постійна часу двигуна; Tya-постійна якірного ланцюга двигуна; Typr-постійна часу обмотки управління ему; Tkz-постійна часу короткозамкненою обмотки ЕМУ; m-деяка мінлива приймаюча значення від 1 до 5; sum-лічильник, який використовується для обчислення діагональних визначників; Kv-дискретність чи крок; diag-один з обчислених діагональних визначників;
3.2 Опис програми
У процедурі Vivod здійснюється розрахунок коефіцієнтів a1, a2, a3, a4, a5 визначника Гурвіца і висновок значень на екран ; У процедурі Obnylenie задаючи деякої змінної m значення від 1 до 5 ми в циклі змінюючи номери рядків і стовпців наводимо матрицю розмірністю 5 '