Shooting method with a constraint

6 次查看（过去 30 天）
sina
sina 2018-9-21
评论： Torsten 2018-9-24
Hi,
I want to solve a system of ODE. It's a BVP. I changed it to IVP and used Shooting method to solve this problem. For the method I defined
X_i_new = X_i - fX_i.*((X_i - X_i_1)/(fX_i - fX_i_1));
X_i and X_i_1 are matrixes like [0.018 0.0258 0.9562] and [0.0150 0.0259 0.9591]. The sum of X_i elements is unity, and also for X_i_1.
The question is X_i_new calculates the new X for the 3rd loop. But the sum of X_i_new does not equal unity. I don't know how to put this criterion in the shooting method.
I am not sure but I changed the shooting method like:
X_i_new = [X_i_new(1) X_i_new(2) 1-(X_i_new(1)+X_i_new(2))];
was that the correct change?
Thank you, Regards

请先登录，再回答此问题。

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by