Solving Coupled Differential Equations (in Matlab)

Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
1 reply · 18K views
ksmanis
Messages
1
Reaction score
0
Hi,

I have the following problem to solve. My system can be defined as a set of coupled differential equations as described below:

dy/dt = -a*y(t) + b*x(t) + (m-y(t))*r;
dx/dt = a*y(t) - b*x(t) + (n-x(t))*r;

where m,n and r are constants. I have the values of y(t) and x(t) for different values of t. The above dy/dt and dx/dt are based on theoretical interpretation (and I do not have their values). For the given set of x(t), y(t), I would like to find the coefficients 'a' and 'b' fitting the above system.

Appreciate your help in solving this problem in any package (preferably matlab)

Kind Regards
Subbu
 
Physics news on Phys.org
The ODE is linear with constant coefficients, so can be solved analytically to obtain [itex](x,y)(t)[/itex] expressly. One could then attempt to minimize the square error [itex]\sum_i \|(x,y)(t_i) - (x_i,y_i)\|^2[/itex] in the usual way.