 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2032 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
Matlab: 用bvp解二階微分方程（邊界值）遇到困難，請高手指點(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）總分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ù)
-
一志愿中海洋材料工程專碩330分求調(diào)劑
已經(jīng)有5人回復(fù)
-
304求調(diào)劑
已經(jīng)有5人回復(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ù)
5樓2012-08-30 21:25:47
dingd
鐵桿木蟲 (職業(yè)作家)
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
2樓2012-08-29 09:22:17
3樓2012-08-29 18:51:56
dingd
鐵桿木蟲 (職業(yè)作家)
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
4樓2012-08-30 09:48:16
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求材料調(diào)劑 +8 | 隔壁陳先生 2026-03-12 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +4 | Ymlll 2026-03-18 | 6/300 |
|
|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +9 | 秋有木北 2026-03-14 | 9/450 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研]
 機(jī)械專碩調(diào)劑
+3
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 招收0805（材料）調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|
|
|
[考研] 283求調(diào)劑，材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|