Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
ilo kog funkcijskog programa <strong>za</strong> pronalaženje broja argumenata koji se koriste u sesiji<br />
koja poziva funkciju.<br />
<strong>MATLAB</strong> (7)<br />
Primer rešavanja jednačine klatna<br />
Razmotrimo klatno koje ima polugu dužine L i kuglu mase m.<br />
Pretpostavimo da je poluga kruta i da ima <strong>za</strong>nemarljivu mase u poređenju sa kuglom.<br />
Takođe, pretpostavimo da nema trenja ili otpora vazduha. Koristeći II Njutnov <strong>za</strong>kon<br />
dobijamo sledeću diferencijalnu jednačinu <strong>za</strong> ugaoni pomeraj, θ(t):<br />
θ ''(t)+(g/L)sin( θ (t))=0<br />
Pretpostavljajući da je kugla u startu pomerena <strong>za</strong> ugao θ 0 i da je puštena da se njiše,<br />
dobijamo sledeće početne uslove:<br />
θ (0)= θ 0 , θ '(0)=0.<br />
Diferencijalnu jednačinu sa ovim početnim uslovima rešavamo koristeći Runge Kutta<br />
Fehlberg metod 4-og reda, tj. funkciju ODE45.<br />
Ukoliko u komandnoj liniji <strong>MATLAB</strong>-a navedete<br />
help ode45<br />
dobićete iscrpno objašnjenje.<br />
Za nas je u ovom slučaju od interesa da objasnimo sledeći oblik funkcije ODE45:<br />
[T, Y] = ODE45 (‘F’, TSPAN, Y0),