Solving Boundry value problem with fsolve

11 次查看(过去 30 天)
  2 个评论
Alex Sha
Alex Sha 2021-4-9
Hi, Justin, try to solve your BVP problems by using 1stOpt:
Code (Case-1):
ODEStep = 0.05;
Constant B=-0.056, Peh=2, Pem=2, r=17.6;
ConstStr Rct=3.36*c*exp(r-r/T);
SubjectTo c'[0]=(c[0]-1)*Pem, T'[0]=(T[0]-1)*Peh;
Variable x,c',T';
ODEFunction c''=(c'+Rct)*Pem;
T''=(T'+B*Rct)*Peh;
Data;
0 NAN NAN
1 0 0
Result (Case-1):
Root of Mean Square Error (RMSE): 3.98105846789366E-10
Sum of Squared Residual: 6.33953060991512E-19
Parameter Best Estimate
-------------------- -------------
c' Initial Value -1.15593943548473
t' Initial Value 0.0647326083792722
c Initial Value 0.422030282257633
t Initial Value 1.03236630418964
SubjectTo:
t'[initialvalue]-((t[initialvalue]-1)*2): -2.77555756156289E-16
c'[initialvalue]-((c[initialvalue]-1)*2): 3.88578058618805E-15
chart of x-c
x-t
chart of x-t
x-c
chart of x-dc/dx
chart of x-dt/dx

请先登录,再进行评论。

回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Linear Programming and Mixed-Integer Linear Programming 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by