21.01.2015 Views

Linearna i nelinearna reaktivna TR - LEDA

Linearna i nelinearna reaktivna TR - LEDA

Linearna i nelinearna reaktivna TR - LEDA

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!