Defining derivatives in ode45

8 vues (au cours des 30 derniers jours)
sina
sina le 5 Mai 2019
Commenté : sina le 6 Mai 2019
Hi,
I want to use ode45 to solve a system of ODE as shown below (see the imgae). I am not sure about coding these two equations as I derived them using finite difference method according to the direction. I coded them as:
dydt(1,1) = f(y(1)) + f(y(2));
dydt(2,1) = - (f(y(1)) + f(y(2))); % with -1
The second one was multiplied by -1 to show a countercurrent flow. I am wondeing if this approach is correct. Or I should ignore the negative sign in coding these equations.
dydt(1,1) = f(y(1)) + f(y(2));
dydt(2,1) = (f(y(1)) + f(y(2))); % without -1
Thank you.
q1.png
  2 commentaires
John D'Errico
John D'Errico le 5 Mai 2019
Why would you ignore the sign?
sina
sina le 5 Mai 2019
This is exactly my question! Do I add -1 to show the countercurrent mode? As there is no -1 in original equations.

Connectez-vous pour commenter.

Réponse acceptée

David Goodmanson
David Goodmanson le 6 Mai 2019
Hello sina,
If you turn the diagram around by 180 degrees and then let Ca <--> Cb, you get exactly the same diagram that you started with. That means the equations have to have the same symmetry under Ca <--> Cb. If both equations have a positive sign you do get the same equations. But if one equation has a negative sign, that negative sign switches to the other equation. So positive works, negative doesn't.
Another way of saying the same thing is, starting out with that symmetric diagram, how would you decide which of dCa/dt or dCb/dt gets the negative sign?
  1 commentaire
sina
sina le 6 Mai 2019
Hi David,
Thank you for the comment.

Connectez-vous pour commenter.

Plus de réponses (0)

Catégories

En savoir plus sur Programming dans Help Center et File Exchange

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by