second order partial differential element in an ODE
Afficher commentaires plus anciens
Hi, Does anyone know who to write a script for the below equation?
z' = 1/L*[Ω Pctr *Vs *cos(Ωt)−(R +2∂L/∂u*v)*z−(1/c+ (∂^2 L)/(∂u^2 )v^2+ ∂L/∂u *v')*y].
The main concern is the second order partial differential element.
Thanks a lot.
4 commentaires
Walter Roberson
le 3 Sep 2011
Additionally, L, u, v, y and z are all functions of t.
Davis sent me a copy of the paper, but I do not know how to deal with the diff(L(t),u(t)) part.
Bjorn Gustavsson
le 3 Sep 2011
There must be something more relating z to some of u, v and y (L?)? Or some of them to make this a second order ODE.
Walter Roberson
le 3 Sep 2011
There is; Davis has shown only one of the equations here.
L(t) has a known exponential formula, but u(t) is defined through another differential equation.
There are two particular difficulties:
1) I do not know how to differentiate a function with respect to another function, diff(L(t),u(t))
2) because u(t) is being worked out numerically by an ode*() call (as u' is one of the equations), the formula for u(t) is not initially known, which makes differentiating with respect to u(t) suspect.
Bjorn Gustavsson
le 3 Sep 2011
Ok. I would (simplistically) start trying to make functions for dL/du and d^2L/du^2 - which should work if the OP just want to solve the ODE for one given L (or if L, dL/du, d^2L/du^2 is called with additional parameters). This might work if L doesn't show up in the LHS of the ODE - which would turn it into something else than an ODE...
Réponses (0)
Catégories
En savoir plus sur Ordinary Differential Equations dans Centre d'aide et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!