24.02.2013 Views

A1 Matlab

A1 Matlab

A1 Matlab

SHOW MORE
SHOW LESS

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

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

Saved successfully!

Ooh no, something went wrong!