24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lixingxun

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 27.5
帖子: 3
在线: 2.3小时
虫号: 803914
注册: 2009-07-05
专业: 能源化工

[求助] 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];

回复此楼

» 猜你喜欢

请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

高手进来用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人回复

1楼2012-08-29 08:17:09

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lixingxun

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 27.5
帖子: 3
在线: 2.3小时
虫号: 803914
注册: 2009-07-05
专业: 能源化工

引用回帖:

2楼: Originally posted by dingd at 2012-08-29 09:22:17
如可能，请把初始、边界条件、公式等用文本格式描述重新清楚。

谢谢，您的回复！我用文本把我的方程具体信息又描述了一下，希望能得到您的指点

%微分方程：
%y''=(4*r*cos(A)/d-32*(u1*L-(u2-u1)*y)*y')/d^2-(p2-p1)*y'^2)/(p1*L+(p2-p1)*y)
%边界值：y(0)=0； y(180)=0.006709
%初值：y(0)=0; y'(0)=0