| 查看: 870 | 回复: 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); 求大神指点，急急急 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有123人回复
-
今天基金会出结果吗？20260819
已经有13人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
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
