Linearna i nelinearna reaktivna TR - LEDA
Linearna i nelinearna reaktivna TR - LEDA
Linearna i nelinearna reaktivna TR - LEDA
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Analiza kola<br />
Analiza kola<br />
Analiza greške diskretizacije<br />
Intuitivno je jasno (a znanja iz numeričke matematike<br />
to potvrđuje) da diskretizacija unesi određenu grešku,<br />
i da može da se očekuje da greška bude manja ako je<br />
korak diskretizacije manji i ako je promena sporija.<br />
Želimo da utvrdimo<br />
-koliko iznosi greška i<br />
-od čega zavisi.<br />
Analiza greške diskretizacije<br />
Neka je x(t n+1 ) tačna vrednost<br />
a x n+1 izračunata vrednost pomenljive x.<br />
Tada je lokalna greška zaokruživanja<br />
(Local trncation Error, LTE)<br />
ε Tx = x(t n+1 ) - x n+1<br />
29<br />
30<br />
Analiza kola<br />
Analiza greške diskretizacije<br />
Razvojem funkcije x(t) u Tajlorov red u okolini<br />
tačke t=t n+1 dobija se<br />
x(t) = x(t<br />
za t = t<br />
x(t ) = x(t<br />
h = t<br />
x(t ) = x(t<br />
x(t<br />
n<br />
n<br />
n+<br />
1<br />
n<br />
n+<br />
1<br />
n+<br />
1<br />
n<br />
) + (t − t<br />
n+<br />
1<br />
− t ,<br />
n+<br />
1<br />
) + (t<br />
) − hx&<br />
) = x(t ) + hx&<br />
n<br />
n<br />
n+<br />
1<br />
− t<br />
t=<br />
t<br />
t=<br />
t<br />
)x&<br />
n+<br />
1<br />
n+<br />
1<br />
n+<br />
1<br />
+<br />
t=<br />
t<br />
)x&<br />
−<br />
n+<br />
1<br />
t=<br />
t<br />
1<br />
h<br />
2<br />
1<br />
h<br />
2<br />
n+<br />
1<br />
2<br />
2<br />
+<br />
&& x<br />
&& x<br />
1<br />
(t<br />
2<br />
+<br />
t=<br />
t<br />
t=<br />
t<br />
− t<br />
1<br />
(t<br />
2<br />
n+<br />
1<br />
n+<br />
1<br />
n+<br />
1<br />
n<br />
− t<br />
+ ...<br />
− ...<br />
2<br />
) && x<br />
n+<br />
1<br />
t=<br />
t<br />
n+<br />
1<br />
2<br />
) && x<br />
+ ...<br />
t=<br />
t<br />
n+<br />
1<br />
+ ...<br />
31<br />
x(t &<br />
Na osnovu<br />
n+<br />
1<br />
x(t<br />
) =<br />
t<br />
x<br />
Analiza kola<br />
Analiza greške diskretizacije<br />
n+<br />
1<br />
n+<br />
1<br />
) − x(t<br />
− t<br />
n<br />
n<br />
sledi da je približna vrednost promenljive x u trenutku t=t n+1<br />
n+<br />
1 n<br />
= x + hx&<br />
ε<br />
Tx<br />
) x(t<br />
=<br />
t=<br />
t n + 1<br />
n+<br />
1<br />
n+<br />
1<br />
n+ 1)<br />
− x = −<br />
) − x(t<br />
h<br />
= x(t<br />
h 2<br />
& x<br />
t=<br />
n<br />
t n + 1<br />
) x<br />
=<br />
n+<br />
1<br />
− x<br />
h<br />
Ako se pretpostavi da je u t=t n , poznato tačno rešenje i da je<br />
x(t n )=x n , tada je<br />
n<br />
32