 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 2042 | 回復: 5 | |||
[求助]
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]; |
回復此樓» 猜你喜歡
328求調(diào)劑
已經(jīng)有7人回復
-
340求調(diào)劑
已經(jīng)有5人回復
-
一志愿南昌大學324求調(diào)劑
已經(jīng)有4人回復
-
求調(diào)劑推薦 材料 304
已經(jīng)有15人回復
-
291求調(diào)劑
已經(jīng)有19人回復
-
331環(huán)境科學與工程求調(diào)劑
已經(jīng)有3人回復
-
289求調(diào)劑
已經(jīng)有8人回復
-
315分求調(diào)劑
已經(jīng)有7人回復
-
295求調(diào)劑
已經(jīng)有6人回復
-
291求調(diào)劑
已經(jīng)有7人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 高手進來用matlab求解常微分方程，十分感謝！
已經(jīng)有17人回復
- 如何求解描述振動的二階微分方程
已經(jīng)有21人回復
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(shù)值比較
已經(jīng)有4人回復
- matlab解微分方程
已經(jīng)有10人回復
- Matlab求解偏微分方程組
已經(jīng)有22人回復
- matlab數(shù)值求解邊界條件微分方程組
已經(jīng)有7人回復
- 高金求助matlab解微分方程組
已經(jīng)有12人回復
- matlab解微分方程組
已經(jīng)有15人回復
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題
已經(jīng)有21人回復
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已經(jīng)有22人回復
- 【求助】積分微分方程matlab求解
已經(jīng)有6人回復
- 【求助】一個二階常微分方程的解法（重獎金幣�。�
已經(jīng)有16人回復
- 【求助】如何用Runge-Kutta迭代求解二階常微分方程組【已解決】
已經(jīng)有9人回復
dingd
鐵桿木蟲 (職業(yè)作家)
- 應助: 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è)作家)
- 應助: 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
Figo1937
銀蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 441.8
- 帖子: 7
- 在線: 30.7小時
- 蟲號: 3847232
- 注冊: 2015-05-04
- 性別: GG
- 專業(yè): 聚合物共混與復合材料
6樓2015-05-25 10:05:19
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 lixingxun 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿哈爾濱工業(yè)大學材料與化工方向336分 +6 | 辰沐5211314 2026-03-26 | 6/300 |
|
|
|
[考研] 083000學碩274求調(diào)劑 +5 | Li李魚 2026-03-26 | 5/250 |
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
|
[考研] 307求調(diào)劑 +8 | 超級伊昂大王 2026-03-24 | 9/450 |
|
|
|
[考研] 考研化學308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
|
[考研] 22408 359分調(diào)劑 +3 | Qshers 2026-03-27 | 3/150 |
|
|
|
[考研] 334求調(diào)劑 +3 | 雨清天晴 2026-03-21 | 3/150 |
|
|
|
[考研]
 材料與化工考研調(diào)劑
+16
|
孅華 2026-03-22 | 16/800 |
|
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
|
[考研] 286求調(diào)劑 +13 | Faune 2026-03-21 | 13/650 |
|
|
|
[考研] 總分322求生物學/生化與分子/生物信息學相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
|
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
|
[考研] 各位老師您好：本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
|
[考研] 318求調(diào)劑 +5 | plum李子 2026-03-21 | 8/400 |
|
|
|
[考研] 一志愿吉大化學322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|