| 查看: 768 | 回复: 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); 求大神指点，急急急 |
回复此楼» 猜你喜欢
278求调剂
已经有9人回复
-
297求调剂
已经有12人回复
-
0703化学297求调剂
已经有3人回复
-
材料学硕301分求调剂
已经有6人回复
-
280求调剂
已经有11人回复
-
311求调剂
已经有3人回复
-
求调剂
已经有3人回复
-
313求调剂
已经有4人回复
-
317求调剂
已经有15人回复
-
307求调剂
已经有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
