How to solve the differential equation numerically
56 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
My equation is this:
(-Vs/500)-((6e-10)*diff(Vs,t))+(7/250)+(7.226*Vs*((Vs/1000)+((3e-10)*diff(Vs,t))-7/500)^2)-(3/25000)*(Vs^2)-((4.336e-9)*diff(Vs,t))-(0.014*Vs)-((1.08e5)*((Vs/1000)+((3e-10)*diff(Vs,t))-7/500))+0.2023=0
Vs(0)=7.9
I tried using dsolve but there is no explicit solution. I need solved vector with numerical values for Vs with constant rise in 't'.
Can anyone please help me?
5 commentaires
J. Alex Lee
le 9 Juin 2020
By "no explicit solution" I think you mean "no analytical solution".
If diff(Vs(t),t) means and that is the only differential term, your equation still looks "explicit in ", which means you can put your equation in the form
rather than how you have it. I'm not going to go through your equation term by term to verify, you should do so yourself.
Then, you should be able to follow Asvin's advice to look at examples of how to solve that class of problems using ode45.
If your problem is truly "implicit" in , then look at https://www.mathworks.com/help/matlab/math/choose-an-ode-solver.html under "Types of ODEs" near the top
Réponse acceptée
Asvin Kumar
le 9 Juin 2020
As others have mentioned in the comments, ODE seems best suited for this problem. Here’s an example from which you can adapt a solution. Solve a Nonstiff Equation using ode45
2 commentaires
Asvin Kumar
le 9 Juin 2020
I see what you're talking about. ode15i is the recommended way to go about these problems. Have a look at this example on solving Weissinger Implicit ODE to get started. This should give you an idea of how you can use the function.
Plus de réponses (1)
Voir également
Catégories
En savoir plus sur Ordinary Differential Equations dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!