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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » Matlab: 用bvp解二阶微分方程（边界值）遇到困难，请高手指点
51/1返回列表
查看: 2018  |  回复: 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

铁虫 (初入文坛)
引用回帖:
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 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
如可能，请把初始、边界条件、公式等用文本格式描述重新清楚。
赞一下 回复此楼
2楼2012-08-29 09:22:17
已阅   回复此楼   关注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的回帖

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的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 328求调剂，英语六级551，有科研经历 +3 生物工程调剂 2026-03-16 4/200 2026-03-16 20:13 by Wangjingyue
[考研] 化学调剂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 诺贝尔化学奖觊�
[基金申请] 今年的国基金是打分制吗？ 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 一志愿华中师范071000，325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[硕博家园] 深圳大学硕士招生（2026秋，传感器方向，仅录取第一志愿） +4 xujiaoszu 2026-03-11 8/400 2026-03-16 09:45 by xujiaoszu
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[基金申请] 面上和青基一样限30页不合理 +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 一志愿华中农业大学071010，总分三百二,求调剂 +3 困困困困坤坤 2026-03-10 3/150 2026-03-14 00:35 by JourneyLucky
[考研] 321求调剂 +3 CUcat 2026-03-10 3/150 2026-03-14 00:25 by JourneyLucky
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 285化工学硕求调剂（081700） +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 08食品或轻工求调剂，本科发表3篇sci一区top论文，一志愿南师大食品科学与工程 +3 我是一个兵， 2026-03-10 3/150 2026-03-13 10:21 by Yuyi.
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考研] 一志愿：武汉理工，材料工程，英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
关闭
请填处理意见
关闭