应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » Matlab: 用bvp解二阶微分方程（边界值）遇到困难，请高手指点
51/1返回列表
查看: 2017  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lixingxun

铁虫 (初入文坛)

[求助] 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];
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-08-29 08:17:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lixingxun

铁虫 (初入文坛)
引用回帖:
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
赞一下 回复此楼
3楼2012-08-29 18:51:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
如可能，请把初始、边界条件、公式等用文本格式描述重新清楚。
赞一下 回复此楼
2楼2012-08-29 09:22:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

xzhdty: 谢谢参与 2012-08-31 22:50:58
如果有：初值：y(0)=0; y'(0)=0，那就成为初值问题了；
如果：y(0)=0； y(180)=0.006709，就是两点边值问题。

你想得是哪种啊。
赞一下 回复此楼
4楼2012-08-30 09:48:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lixingxun

铁虫 (初入文坛)
引用回帖:
4楼: Originally posted by dingd at 2012-08-30 09:48:16
如果有：初值：y(0)=0; y'(0)=0，那就成为初值问题了；
如果：y(0)=0； y(180)=0.006709，就是两点边值问题。

你想得是哪种啊。

谢谢您的回复，我要用边界值。我指的初值是初始猜想值
赞一下 回复此楼
5楼2012-08-30 21:25:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自科面上基金字体 +6 iwuli 2026-03-12 7/350 2026-03-16 21:18 by sculhf
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 9/450 2026-03-16 17:35 by 沐霖12138
[考研] 化学调剂0703 +8 啊我我的 2026-03-11 8/400 2026-03-16 17:23 by 我的船我的海
[考研] 326求调剂 +4 诺贝尔化学奖觊� 2026-03-15 7/350 2026-03-16 17:11 by 诺贝尔化学奖觊�
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +3 大火山小火山 2026-03-16 5/250 2026-03-16 16:54 by barlinike
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 材料与化工 323 英一+数二+物化，一志愿：哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 321求调剂 +3 CUcat 2026-03-10 3/150 2026-03-14 00:25 by JourneyLucky
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 工科，求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 321求调剂（食品/专硕） +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
信息提示
关闭
请填处理意见
关闭