program Krit;crt;n=5; Tem=0.131; Tya=0.05; Typr=0.076; Tkz=0.148;massiv= array [1 n,1 n] of real;a,b:massiv; Kv,diag:real; i,j: integer;Vivod;textcolor(9); writeln(' Коэффициенты определителя Гурвица :'); textcolor(7); writeln(' a1=',1); writeln(' a2=',Tya+Typr+Tkz:9:8); writeln(' a3=',Tem*Tya+Typr*Tya+Tkz*Tya+Tkz*Typr:9:8); writeln(' a4=',Typr*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.
Самое читаемое:
Испытания схем увязки САУТ-ЦМ и МПЦ в лаборатории
Главной задачей решаемой на железнодорожном транспорте всеми его
службами, хозяйствами, техническими и организационными мероприятиями является
обеспечение непрерывного, нормального (штатного) протекания основного
технологического процесса (ТП) движения поездов. Движение поездов - это
комплексный технологический процесс, состоящий из ...