Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Rešavanje diferencijalnih jednaĉina pomoću ode23 ili ode45<br />
• može se rešiti obiĉna diferencijalna jednaĉina višeg reda<br />
• treba napraviti sistem diferencijalnih jednaĉina 1. reda i zapisati ga u<br />
zasebnu M-funkciju (sistem jednaĉina može biti nelinearan)<br />
zaglavlje funkcije mora biti oblika: function xprim = ime( t, x )<br />
gde je t vreme, a x vektor promenljivih stanja<br />
postoji problem sa prenošenjem parametara u ovakvu M-funkciju<br />
• poziv metode: [ t, x ] = ode23( @ime, [tp, tk], x0 )<br />
ili [ t, x ] = ode45( @ime, [tp, tk], x0 )<br />
Ulazni parametri: ime - ime M-funkcije koja opisuje sistem jednaĉina,<br />
tp - poĉetni trenutak integracije, tk - krajnji trenutak, x0 - vektor<br />
poĉetnih vrednosti promenljivih stanja.<br />
Izlazni parametri: t - vektor vremenskih trenutaka u kojima su<br />
izraĉunata rešenja, a x - matrica kretanja promenljivih stanja poreĊanih<br />
po kolonama. Prva vrsta odgovara x0 i poĉetnom trenutku tp, dok je u<br />
poslednjoj vrsti krajnja vrednost prom. stanja (u tk).<br />
• ode metode imaju promenljivi (adaptivni) korak integracije, tako da<br />
dobijena rešenja nisu ekvidistantna