Streeter Pheps Equations
显示 更早的评论
I need to do a program in Matlab with two first-order differential equations with the Streeter Phelps model. Please help me.
BOD: dL/dt=-(k1+ks)L+B
DO: d(D)/dt=k1L-k2D-P
回答（2 个）
Matt Kindig
2012-5-22
What have you done so far? To get you started, look at the documentation for ode45 (and the other ODE solvers). The documentation gives really clear examples of how to set up your equations into the form required by the solvers. For first-order systems like this, it is almost trivial.
doc ode45
Walter Roberson
2012-5-22
I have renamed your "D" to "Df" in order to prevent confusion with differentiation.
The Maple syntax would be
dsolve( [diff(L(t), t) = -(k1+ks)*L(t)+B, diff(Df(t), t) = k1*L(t)-k2*Df(t)-P ])
with result
Df(t) = (-k1 * k2 * C2 * (k1+ks) * exp(-t * (-k2+k1+ks)) + ((-P * ks + (B-P) * k1) * exp(k2*t) + C1 * k2 * (k1+ks)) * (-k2+k1+ks)) * exp(-k2*t) / ((k1+ks) * k2 * (-k2+k1+ks))
L(t) = B / (k1+ks) + exp(-(k1+ks) * t) * C2
In the above, C1 and C2 are constants of integration
2 个评论
Carlos Gomes
2012-5-22
Walter Roberson
2012-5-22
http://www.mathworks.com/matlabcentral/answers/8026-best-way-s-to-master-matlab
类别
在 帮助中心 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息
2012-5-22
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
