How to Solve an ODE Problem when one of parameters is dependent to derivative?

Click For Summary

Discussion Overview

The discussion revolves around solving a system of non-linear ordinary differential equations (ODEs) in MATLAB, specifically when one of the parameters is dependent on the derivative itself. Participants explore the implications of this dependency and how to approach solving the equations using MATLAB's ODE solvers.

Discussion Character

  • Technical explanation
  • Debate/contested
  • Mathematical reasoning

Main Points Raised

  • One participant presents a general form of the ODE, indicating that it involves a parameter dependent on the derivative.
  • Another participant questions whether the function f is known, to which the original poster confirms it is known but not reversible.
  • A participant argues that the equation does not contain y or t, suggesting that dy/dt must be constant and can be solved numerically rather than using an ODE solver.
  • The original poster clarifies that the parameters A, B, C, and D are not constant and emphasizes the need for a solution method that accommodates the dependency on dy/dt.
  • A later reply proposes a recursive algorithm to compute the values of A, B, C, and D, and suggests using a numerical equation solver to find dy/dt, followed by incrementing y and t.

Areas of Agreement / Disagreement

Participants express differing views on the nature of the ODE and the appropriate methods for solving it. There is no consensus on the best approach, and the discussion remains unresolved regarding the specific implementation in MATLAB.

Contextual Notes

Participants highlight the importance of the dependency of parameters on the derivative and the implications for the form of the ODE. The discussion includes assumptions about the nature of the functions involved and the potential need for numerical methods rather than standard ODE solvers.

mahdi_zabchek
Messages
4
Reaction score
0
Hello Guys!
I have an ODE problem that I'm solving it by MATLAB ODE solvers!
in fact I have a system of non-linear differential equations in one of these equations I have a parameter that it's value is dependent to derivative! the general form of equation is like this (big letter parameters are known!):

dy/dt = A + B + f(C,D,dy/dt)

how can I solve this problem by ode45 or other MATLAB ODE solvers?
 
Physics news on Phys.org
Is the function f known?
 
hunt_mat said:
Is the function f known?

yes! it is.
but it's not reversible
 
mahdi_zabchek said:
Hello Guys!
I have an ODE problem that I'm solving it by MATLAB ODE solvers!
in fact I have a system of non-linear differential equations in one of these equations I have a parameter that it's value is dependent to derivative! the general form of equation is like this (big letter parameters are known!):

dy/dt = A + B + f(C,D,dy/dt)

how can I solve this problem by ode45 or other MATLAB ODE solvers?

The ODE : dy/dt = A + B + f(C,D,dy/dt) contains no y and no t. As a consequence dy/dt = constant.
Let X= dy/dt . X is solution of the equation X = A + B + f(C, D, X) which is not an ODE.
It doesn't matter if the function is not revertsible. We don't need to know the analytical expression of the solution(s) X. We know that dy/dt = constant (or = several different constants if there are several solutions). Each one can be numerically computed, not using an ODE solver, but using an usual numerical equation solver.
The solution(s) is (are) : y(t) = X*t +c
c is a constant to be determined by the boundary condition.
 
JJacquelin said:
The ODE : dy/dt = A + B + f(C,D,dy/dt) contains no y and no t. As a consequence dy/dt = constant.
Let X= dy/dt . X is solution of the equation X = A + B + f(C, D, X) which is not an ODE.
It doesn't matter if the function is not revertsible. We don't need to know the analytical expression of the solution(s) X. We know that dy/dt = constant (or = several different constants if there are several solutions). Each one can be numerically computed, not using an ODE solver, but using an usual numerical equation solver.
The solution(s) is (are) : y(t) = X*t +c
c is a constant to be determined by the boundary condition.

No! No! it has y and t!
A and B and C and D are NOT constant parameters!
I did't write them because they were not necessary!
in fact You don't need to know what's the equation exactly to answer my question!

My question is simple:

MATLAB ODE solvers solve equations in form of dy/dt=f(t,y) but I want to solve an equation in form of dy/dt=f(t,y,dy/dt) ... How I can do that by MATLAB?
 
mahdi_zabchek said:
No! No! it has y and t!
A and B and C and D are NOT constant parameters!
I did't write them because they were not necessary!
in fact You don't need to know what's the equation exactly to answer my question!

My question is simple:

MATLAB ODE solvers solve equations in form of dy/dt=f(t,y) but I want to solve an equation in form of dy/dt=f(t,y,dy/dt) ... How I can do that by MATLAB?

OK. Sorry for the missunderstanding.
May be, you could use an algorithm of this kind:
Start with given initial values y and t.
Recursive process :
Compute A(y,t), B(y,t), C(y,t) and D(y,t)
Solve X=A+B+f(C,D,X) with a numerical equation solver, introduced as sub-program.
With the computed value X=dy/dt the incrementation of y is done, as well as the incrementation of t.
Then continue the recursive process.
 
JJacquelin said:
OK. Sorry for the missunderstanding.
May be, you could use an algorithm of this kind:
Start with given initial values y and t.
Recursive process :
Compute A(y,t), B(y,t), C(y,t) and D(y,t)
Solve X=A+B+f(C,D,X) with a numerical equation solver, introduced as sub-program.
With the computed value X=dy/dt the incrementation of y is done, as well as the incrementation of t.
Then continue the recursive process.

I'll try it ... thank you so much
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 5 ·
Replies
5
Views
3K
Replies
3
Views
3K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 11 ·
Replies
11
Views
3K