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.
Самое читаемое:
Конструкторско-технологическое проектирование печатной платы
печатная плата
Проектирование печатных плат (ПП) представляет трудоемкий, но очень
важный процесс. Для того, чтобы обеспечить функционирование электронной
аппаратуры (ЭА) необходимы не только схемотехнические решения, функциональная
точность, надежность, но и учет влияния внешней среды, конструктивных,
эксплуатационных требований, пр ...