| 查看: 1086 | 回复: 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解出来的总感觉不对 |
回复此楼» 猜你喜欢
一个有机合成实验室都需要哪些设备？
已经有7人回复
-
2026年国自然面上资助率
已经有20人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有22人回复
-
微信指数没变化，科研之友没阅读
已经有18人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有8人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有15人回复
-
你们的时间戳变了吗
已经有4人回复
-
产物和副产物价值比较
已经有5人回复
-
时间戳他又来了
已经有17人回复