| 查看: 2012 | 回复: 5 | ||
[求助]
Matlab: 用bvp解二阶微分方程（边界值）遇到困难，请高手指点
|
|
遇到的error是??? Error using ==> bvp4c Unable to solve the collocation equations -- a singular Jacobian encountered Error in ==> xBVP1 at 18 sol = bvp4c(@ODEfun,@BCfun,solinit) 下面是我的code function xBVP1 clear all clc global r A d u1 L u2 p1 p2; r=30E-3;%water/oil surfac tension(N/m) A=pi/6; %water/oil pore contac angle d=155E-6;%size (m) u1=6.14E-03; %oil viscosity(Pa.s) L=2.2E-2; %total length (m) u2=8.94E-04;%water viscosity(Pa.s) p1=738.5; % oil density p2=998.2;% water density % ------------------------------------------------------------------ a=0; b=180; solinit = bvpinit(linspace(a,b,10),[0 0]); sol = bvp4c(@ODEfun,@BCfun,solinit) format long y=deval(sol,sol.x); plot(sol.x,sol.y);legend('y1','y2') % ------------------------------------------------------------------ function dydx=ODEfun(x,y) global r A d u1 L u2 p1 p2; dydx=[y(2);(4*r*cos(A)/d-32*(u1*L-(u2-u1)*y(1))*y(2)/d^2-(p2-p1)*y(2)^2)/(p1*L+(p2-p1)*y(1))]; % ------------------------------------------------------------------ function bc=BCfun(ya,yb) bc=[ya(1); yb(1)-0.006709]; |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- 如何求解描述振动的二阶微分方程
已经有21人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- matlab解微分方程
已经有10人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 高金求助matlab解微分方程组
已经有12人回复
- matlab解微分方程组
已经有15人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】一个二阶常微分方程的解法（重奖金币！）
已经有16人回复
- 【求助】如何用Runge-Kutta迭代求解二阶常微分方程组【已解决】
已经有9人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2012-08-29 09:22:17
3楼2012-08-29 18:51:56
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2012-08-30 09:48:16
5楼2012-08-30 21:25:47
Figo1937
银虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 441.8
- 帖子: 7
- 在线: 30.7小时
- 虫号: 3847232
- 注册: 2015-05-04
- 性别: GG
- 专业: 聚合物共混与复合材料
6楼2015-05-25 10:05:19
