| 查看: 1121 | 回复: 0 | |||
[交流]
Matlab求解微分方程bvp4c
|
|
solinit = bvpinit(linspace(0,14,10),[1;0]); sol=bvp4c(@ODEfun,@BCfun,solinit); x = linspace(0,14,127); y = deval(sol,x); plot(x,y(2,) function dydx=ODEfun(x,y) a=pi*14.701^4/64; b=pi*(24^4-14.701^4)/64; c=3.54365802196686; d=5.94058826800217; e=2/14.701/pi; f=14.701/2; dydx(1,1)=y(2); dydx(2,1)=(1/a+1/b)/(c+d)*f^2/e*(cos(x/(f*e)))^2*y(1)-(tan(x/(f*e)))/f/e*y(2); function res=BCfun(ya,yb) res=[ya(1); yb(1)-1]; 哪位大神给看看我这个微分方程的数值解是多少。用bvp4c解出来的总感觉不对 |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
200