 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4735 | 回復(fù): 7 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
matlab數(shù)值求解邊界條件微分方程組
|
||
期望高手指點(diǎn)一下我的疑問，指出一個大方向也好，我現(xiàn)在就是處于無從下手的階段，求點(diǎn)撥 解一個方程組，文獻(xiàn)上給出了如下的形式： du/dx+dv/dx=0 u=dn/dx v=-dm/dx du/dx=-exp(-x)+n*m/(m+n) 邊界條件和初值條件 u(p)=0,u(0)=v(p),v(0)=0 n(0)=q 其中函數(shù)自變量為x，邊界條件和初值條件中p為常量，q為變量 按照matlab求解邊界條件微分方程bvp4c方法，我進(jìn)行了如下處理： y(1)=n,y(2)=m,y(3)=u,y(4)=v y'(1)=y(3),y'(2)=-y(4),y'(3)=-exp(-x)+y(1)*y(2)/(y(1)+y(2)),y'(4)=exp(-x)-y(1)*y(2)/(y(1)+y(2)) 方程M文件 function dydx=ivode(x,y) dydx=[y(3);-y(4);-exp(-x)+2*y(1)*y(2)/(2*y(2)+y(1));exp(-x)-2*y(1)*y(2)/(2*y(2)+y(1))]; 邊界條件M文件 function res=ivbc(ya,yb) res=[ya(1)-q;ya(4);yb(3);ya(3)-yb(4)]; command: solinit=bvpinit(linspace(0,1,10),[1 1 0 0]); sol=bvp4c(@ivode,@ivbc,solinit); ??? Error using ==> bvp4c Unable to solve the collocation equations -- a singular Jacobian encountered 請問這個問題怎么解決呢？還有邊界條件中u(0)=v(p)是否必須需要等于一個常量，方程組才能給出數(shù)值解。手工處理du/dx=-exp(-x)+n*m/(m+n)方程，分離變量也沒弄明白。。。曾經(jīng)學(xué)的高數(shù)都還給老師了。。。 |
回復(fù)此樓» 猜你喜歡
316求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料與化工考研調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿211 初試270分 求調(diào)劑
已經(jīng)有4人回復(fù)
-
08工學(xué)調(diào)劑
已經(jīng)有11人回復(fù)
-
350求調(diào)劑
已經(jīng)有6人回復(fù)
-
求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601）
已經(jīng)有4人回復(fù)
-
接收2026碩士調(diào)劑(學(xué)碩+專碩)
已經(jīng)有6人回復(fù)
-
081700 調(diào)劑 267分
已經(jīng)有5人回復(fù)
-
328求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- matlab如何求解一個非線性微分方程組
已經(jīng)有8人回復(fù)
- 求助，用matlab解微分方程組，希望高手能給予指點(diǎn)
已經(jīng)有17人回復(fù)
- 如何用Matlab求時變微分方程的解？
已經(jīng)有5人回復(fù)
- matlab求解高階方程組
已經(jīng)有3人回復(fù)
- MATLAB求解非線性方程組
已經(jīng)有5人回復(fù)
- 關(guān)于matlab求解非線性指數(shù)方程組出現(xiàn)問題
已經(jīng)有3人回復(fù)
- 有償求助：matlab怎么解聯(lián)立方程組？
已經(jīng)有11人回復(fù)
- 求高手解答Matlab求解微分方程的初始值和邊界值問題
已經(jīng)有5人回復(fù)
- 幫忙用Mathematica求解一個微分方程組
已經(jīng)有6人回復(fù)
- MATLAB求解代數(shù)方程組問題，請幫忙
已經(jīng)有9人回復(fù)
- Matlab求解偏微分方程組
已經(jīng)有22人回復(fù)
- 請問如何求解二元一階微分方程組
已經(jīng)有24人回復(fù)
- 高金求助matlab解微分方程組
已經(jīng)有12人回復(fù)
- 求matlab求解一個方程組的問題
已經(jīng)有6人回復(fù)
- 一個偏微分方程組（包括邊界條件）如何求解得到自振頻率和振型？
已經(jīng)有4人回復(fù)
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題
已經(jīng)有21人回復(fù)
- 【求助】請教非齊次常微分方程組的解析解法
已經(jīng)有4人回復(fù)
- 【求助】積分微分方程matlab求解
已經(jīng)有6人回復(fù)
- 【求助】matlab怎么求解偏微分方程組啊，先謝謝了
已經(jīng)有13人回復(fù)
- 【求助】求教這種微分方程的詳細(xì)求解過程
已經(jīng)有14人回復(fù)
- 【分享】KS方程求解器-matlab版
已經(jīng)有10人回復(fù)
wangyuewu08
木蟲 (著名寫手)
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 2794
- 紅花: 4
- 帖子: 1141
- 在線: 279.1小時
- 蟲號: 1290073
- 注冊: 2011-05-09
- 性別: GG
- 專業(yè): 制造系統(tǒng)與自動化
4樓2011-12-20 20:10:30
wangyuewu08
木蟲 (著名寫手)
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 2794
- 紅花: 4
- 帖子: 1141
- 在線: 279.1小時
- 蟲號: 1290073
- 注冊: 2011-05-09
- 性別: GG
- 專業(yè): 制造系統(tǒng)與自動化
【答案】應(yīng)助回帖
★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai(金幣+1): 謝謝應(yīng)助！ 2011-12-21 22:25:39
acnes(金幣+20): ★★★很有幫助 嗯，非常感謝積極的回復(fù)！ 2011-12-22 15:56:40
感謝參與，應(yīng)助指數(shù) +1
xiegangmai(金幣+1): 謝謝應(yīng)助！ 2011-12-21 22:25:39
acnes(金幣+20): ★★★很有幫助 嗯，非常感謝積極的回復(fù)！ 2011-12-22 15:56:40
|
我對微分方程組不是很了解，把樓主的代碼運(yùn)行了之后，系統(tǒng)出現(xiàn)了下面的錯誤： 我試著把邊界條件中q的值改為1，運(yùn)行之后沒有報錯，得到了兩個結(jié)構(gòu)體變量。 所以啊，問題應(yīng)該出在q這個變量上，既然是數(shù)值解法，那么邊界條件中不應(yīng)該有變量。 |
2樓2011-12-20 18:38:48
3樓2011-12-20 19:08:24
5樓2011-12-20 21:20:14
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 一志愿070300浙大化學(xué)358分，求調(diào)劑！
+4
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 317求調(diào)劑 +12 | 申子申申 2026-03-19 | 18/900 |
|
|
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑 +4 | 幸運(yùn)的醬醬 2026-03-22 | 5/250 |
|
|
|
[考研] 306求調(diào)劑 +5 | 來好運(yùn)來來來 2026-03-22 | 5/250 |
|
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
|
[考研] 一志愿中國石油大學(xué)（華東） 本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-17 | 3/150 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
|
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|