| 查看: 2226 | 回复: 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]; |
回复此楼» 猜你喜欢
可靠消息，8月21号(明天)放榜
已经有10人回复
-
只有每年这种时候来逛逛小木虫
已经有20人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有125人回复
-
今天基金会出结果吗？20260819
已经有15人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高手进来用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
