版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

lixingxun

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 27.5
帖子: 3
在線: 2.3小時(shí)
蟲號(hào): 803914
注冊(cè): 2009-07-05
專業(yè): 能源化工

[求助] Matlab: 用bvp解二階微分方程（邊界值）遇到困難，請(qǐng)高手指點(diǎn)

遇到的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];

回復(fù)此樓

» 猜你喜歡

085600材料與化工調(diào)劑 324分已經(jīng)有5人回復(fù)
344求調(diào)劑已經(jīng)有6人回復(fù)
266求調(diào)劑已經(jīng)有9人回復(fù)
化學(xué)工程321分求調(diào)劑已經(jīng)有18人回復(fù)
314求調(diào)劑已經(jīng)有8人回復(fù)
08工科 320總分求調(diào)劑已經(jīng)有5人回復(fù)
085601專碩，總分342求調(diào)劑，地區(qū)不限已經(jīng)有5人回復(fù)
311求調(diào)劑已經(jīng)有6人回復(fù)
化工學(xué)碩306求調(diào)劑已經(jīng)有10人回復(fù)
286求調(diào)劑已經(jīng)有8人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

高手進(jìn)來(lái)用matlab求解常微分方程，十分感謝！已經(jīng)有17人回復(fù)
如何求解描述振動(dòng)的二階微分方程已經(jīng)有21人回復(fù)
matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(shù)值比較已經(jīng)有4人回復(fù)
matlab解微分方程已經(jīng)有10人回復(fù)
Matlab求解偏微分方程組已經(jīng)有22人回復(fù)
matlab數(shù)值求解邊界條件微分方程組已經(jīng)有7人回復(fù)
高金求助matlab解微分方程組已經(jīng)有12人回復(fù)
matlab解微分方程組已經(jīng)有15人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問(wèn)題已經(jīng)有21人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)
【求助】積分微分方程matlab求解已經(jīng)有6人回復(fù)
【求助】一個(gè)二階常微分方程的解法（重獎(jiǎng)金幣�。� 已經(jīng)有16人回復(fù)
【求助】如何用Runge-Kutta迭代求解二階常微分方程組【已解決】已經(jīng)有9人回復(fù)

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

如可能，請(qǐng)把初始、邊界條件、公式等用文本格式描述重新清楚。

贊一下

回復(fù)此樓

2樓2012-08-29 09:22:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lixingxun

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 27.5
帖子: 3
在線: 2.3小時(shí)
蟲號(hào): 803914
注冊(cè): 2009-07-05
專業(yè): 能源化工

引用回帖:

2樓: Originally posted by dingd at 2012-08-29 09:22:17
如可能，請(qǐng)把初始、邊界條件、公式等用文本格式描述重新清楚。

謝謝，您的回復(fù)！我用文本把我的方程具體信息又描述了一下，希望能得到您的指點(diǎn)

%微分方程：
%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

贊一下

回復(fù)此樓

3樓2012-08-29 18:51:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

xzhdty: 謝謝參與 2012-08-31 22:50:58

如果有：初值：y(0)=0; y'(0)=0，那就成為初值問(wèn)題了；
如果：y(0)=0； y(180)=0.006709，就是兩點(diǎn)邊值問(wèn)題。

你想得是哪種啊。

贊一下

回復(fù)此樓

4樓2012-08-30 09:48:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lixingxun

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 27.5
帖子: 3
在線: 2.3小時(shí)
蟲號(hào): 803914
注冊(cè): 2009-07-05
專業(yè): 能源化工

引用回帖:

4樓: Originally posted by dingd at 2012-08-30 09:48:16
如果有：初值：y(0)=0; y'(0)=0，那就成為初值問(wèn)題了；
如果：y(0)=0； y(180)=0.006709，就是兩點(diǎn)邊值問(wèn)題。

你想得是哪種啊。

謝謝您的回復(fù)，我要用邊界值。我指的初值是初始猜想值

贊一下

回復(fù)此樓

5樓2012-08-30 21:25:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Figo1937

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 441.8
帖子: 7
在線: 30.7小時(shí)
蟲號(hào): 3847232
注冊(cè): 2015-05-04
性別: GG
專業(yè): 聚合物共混與復(fù)合材料

請(qǐng)問(wèn)樓主用bvp遇到雅克比奇異矩陣的問(wèn)題解決了嗎？我現(xiàn)在也遇到這個(gè)問(wèn)題，不知怎么解決

贊一下

回復(fù)此樓

風(fēng)起于青萍之末

6樓2015-05-25 10:05:19

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 lixingxun 的主題更新

返回列表

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4	脫穎而出 2026-03-16	9/450	2026-03-18 14:13 by 脫穎而出
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 303求調(diào)劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開(kāi)始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 283求調(diào)劑 +3	聽(tīng)風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 265求調(diào)劑 +4	威化餅07 2026-03-12	4/200	2026-03-14 17:23 by userper
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 308 085701 四六級(jí)已過(guò)求調(diào)劑 +7	溫喬喬喬喬 2026-03-12	14/700	2026-03-14 10:49 by JourneyLucky
[基金申請(qǐng)] 有必要更換申報(bào)口嗎 20+3	fannyamoy 2026-03-11	3/150	2026-03-14 00:52 by zhanghaozhu
[考研] 0703化學(xué)求調(diào)劑 +7	綠豆芹菜湯 2026-03-12	7/350	2026-03-13 17:25 by njzyff
[考研] 材料專碩350 求調(diào)劑 +4	王金科 2026-03-12	4/200	2026-03-13 16:02 by ruiyingmiao
[論文投稿] 投稿問(wèn)題 5+4	星光燦爛xt 2026-03-12	6/300	2026-03-13 14:17 by god_tian
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655