Předložený studijní materiál slouží jako základní studijní materiál distanční formy
studia předmětu Elektrotechnika 2, který navazuje na předmět Elektrotechnika 1 a spolu s ním
vytváří nezbytně nutné teoretické základy společné pro všechny elektrotechnické obory, které
jsou potřebné pro studium předmětů specializací v dalších ročnících studia.
Autor: Doc. Ing. Jiří Sedláček, CSc. Prof. Ing. Juraj Valsa, CSc.
Strana 183 z 186
Vámi hledaný text obsahuje tato stránku dokumentu který není autorem určen k veřejnému šíření.
7071067; cb=0.928*(x[1]+x[2]);
dx[2]:=-59523.
program runge;
const
tmin pocatek intervalu reseni *)
tmax 500e-6; konec intervalu reseni *)
nt 101; celkovy pocet bodu reseni *)
nn pocet rovnic *)
type
vek array[1.Elektrotechnika 183
Diferenciální rovnice mají tvar
( 2
221
21
1
11
2
21
11
1 1
,
1
CC
C
CC
C
u
CRR
RR
u
CRdt
du
uu
CRdt
du +
−−=+−=
Parametry obvodu jsou 1680 0,235µF, kΩ, nF.nn] real;
const
x vek (1,-3); pocatecni podminky *)
var
h,t real;
i,j,k integer;
procedure RKG(nn:integer;h:real;var t:real;var x:vek);
(* algoritmus Runge-Kutta-Gill radu *)
const
ca=1.1*x[2];
end;
begin
(* vlastni algoritmus *)
difr;
for i:=1 do
begin q1:=0.5*h;
difr;
for i:=1 do
begin x[i]:=x[i]+h*dx[i]/6-q[i]/3; end;
.8*x[1]-67857.5*h;
difr;
for i:=1 do
begin q1:=h*dx[i]; q2:=q1-q[i]; x[i]:=x[i]+cb*q2;
q[i]:=q[i]+3*cb*q2-cb*q1; end;
difr;
for i:=1 do
begin q1:=h*dx[i]; q2:=q1-q[i]; x[i]:=x[i]+ca*q2;
q[i]:=q[i]+3*ca*q2-ca*q1; end;
t:=t+0..5*h*dx[i]; x[i]:=x[i]+q1; q[i]:=2*q1; end;
t:=t+0. rovnic normalnim tvaru *)
begin
dx[1]:=-2532.2928933;
var
dx,q vek;
q1,q2 real;
i integer;
procedure difr;
(* vypocet levych stran dif