Solving Boundry value problem with fsolve
11 次查看(过去 30 天)
显示 更早的评论

2 个评论
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

chart of x-t

chart of x-dc/dx

chart of x-dt/dx

回答(0 个)
另请参阅
类别
在 Help Center 和 File 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!