| 查看: 908 | 回复: 2 | ||
[求助]
非线性微分方程组的边值问题 已有1人参与
|
|
我现在是一个六阶的非线性微分方程，但是只有三个初始值，其他都是边界条件。我现在的难处就是想把边界条件转换成初始条件。 clear clc close all global T; for T = 1:20000 solinit = bvpinit(linspace(0.00000001,20,200),[10^(-5) 0 0 0.2 0 0.9]);% 解的初始估计 bcfun = @(ya,yb)[ya(1)-10^(-5); yb(1)-2.25;ya(2);yb(2)-10;ya(3);yb(3)-18]; sol = bvp4c(@ringbal1,bcfun,solinit); s=linspace(0.00000001,20,2000); y=deval(sol,s); end X=10.*y(:,1).*cos(y(:,2)); Z=-10.*y(:,3); plot3(T,X,Z); 程序报错为： ??? Error using ==> bvp4c at 252 Unable to solve the collocation equations -- a singular Jacobian encountered Error in ==> dabafa_ at 8 sol = bvp4c(@ringbal1,bcfun,solinit); 求大神指点，急急急 |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有6人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有6人回复
-
现代”学阀”该如何界定
已经有14人回复
-
求合成方法
已经有8人回复
-
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？
已经有4人回复
-
有机合成以后会不会被AI改变？做科研的虫友怎么看
已经有8人回复
-
我的奶奶
已经有3人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2019-01-11 17:31:39
liuichini
铁杆木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 7743.3
- 红花: 13
- 帖子: 492
- 在线: 104.6小时
- 虫号: 3322728
- 注册: 2014-07-14
- 专业: 固体力学
3楼2019-01-31 13:28:04