applying conditions in odesolver
1 visualizzazione (ultimi 30 giorni)
Mostra commenti meno recenti
equation of motion of projectile with drag is dvx/dt = -k*v*vx ; dvy/dt = -k*v*vy - g where k = drag coefficient and g = gravitational acceleration i solved this equation with ode solver (ode45) by giving a constant value of k and g. but in real the value of k and g vary with altitude i have the matrix for 'k' for different height interval .and the formula for variation of g with altitude. how to apply this conditional input to odesolver(ode45). so that odesolver(ode45) take different value of k and g at respective height range during steps .
please reply thanks pawan
0 Commenti
Risposta accettata
Jan
il 12 Set 2011
Simply insert the formula for g in the ODE equation. For k I'd suggest an interpolation, if the values are read from a table: ODE45 expectes a smooth ODE function, but a naive lookup-table leads to discontinuities.
0 Commenti
Più risposte (0)
Vedere anche
Categorie
Scopri di più su Ordinary Differential Equations in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!