| 查看: 729 | 回复: 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); 求大神指点，急急急 |
回复此楼» 猜你喜欢
硝基苯如何除去
已经有3人回复
-
职称评审没过，求安慰
已经有37人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
回收溶剂求助
已经有6人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
申请26博士
已经有5人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
求助文献
已经有3人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 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