| 查看: 754 | 回复: 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); 求大神指点，急急急 |
回复此楼» 猜你喜欢
有院领导为了换新车，用横向课题经费买了俩车
已经有9人回复
-
酰胺脱乙酰基
已经有13人回复
-
博士延得我，科研能力直往上蹿
已经有8人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有4人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有10人回复
-
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有5人回复
-
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优
已经有6人回复
-
面上项目申报
已经有3人回复
-
遇见不省心的家人很难过
已经有22人回复
liuichini
铁杆木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 7743.3
- 红花: 13
- 帖子: 492
- 在线: 104.6小时
- 虫号: 3322728
- 注册: 2014-07-14
- 专业: 固体力学
3楼2019-01-31 13:28:04
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2019-01-11 17:31:39