 |
 |
 |
| 24小時熱門版塊排行榜 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2033 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[求助]
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]; |
回復(fù)此樓» 猜你喜歡
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有8人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有5人回復(fù)
-
一志愿華中科技大學(xué)，080502，354分求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有6人回復(fù)
-
085410人工智能專碩317求調(diào)劑（0854都可以）
已經(jīng)有3人回復(fù)
-
330求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 高手進(jìn)來用matlab求解常微分方程，十分感謝！
已經(jīng)有17人回復(fù)
- 如何求解描述振動的二階微分方程
已經(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求解微分方程組遇到的一個問題
已經(jīng)有21人回復(fù)
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已經(jīng)有22人回復(fù)
- 【求助】積分微分方程matlab求解
已經(jīng)有6人回復(fù)
- 【求助】一個二階常微分方程的解法（重獎金幣！）
已經(jīng)有16人回復(fù)
- 【求助】如何用Runge-Kutta迭代求解二階常微分方程組【已解決】
已經(jīng)有9人回復(fù)
3樓2012-08-29 18:51:56
dingd
鐵桿木蟲 (職業(yè)作家)
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
2樓2012-08-29 09:22:17
dingd
鐵桿木蟲 (職業(yè)作家)
- 應(yīng)助: 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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +8 | llllkkkhh 2026-03-18 | 8/400 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +7 | 手機(jī)用戶 2026-03-17 | 7/350 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
|
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|
|
|
[考研] 321求調(diào)劑（食品/專碩） +3 | xc321 2026-03-12 | 6/300 |
|