 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 4716 | 回復(fù): 7 | |||
[求助]
matlab數(shù)值求解邊界條件微分方程組
|
期望高手指點(diǎn)一下我的疑問(wèn)，指出一個(gè)大方向也好，我現(xiàn)在就是處于無(wú)從下手的階段，求點(diǎn)撥 解一個(gè)方程組，文獻(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 請(qǐng)問(wèn)這個(gè)問(wèn)題怎么解決呢？還有邊界條件中u(0)=v(p)是否必須需要等于一個(gè)常量，方程組才能給出數(shù)值解。手工處理du/dx=-exp(-x)+n*m/(m+n)方程，分離變量也沒(méi)弄明白。。。曾經(jīng)學(xué)的高數(shù)都還給老師了。。。 |
回復(fù)此樓» 猜你喜歡
268求調(diào)劑
已經(jīng)有8人回復(fù)
-
26申博
已經(jīng)有4人回復(fù)
-
211本，11408一志愿中科院277分，曾在中科院自動(dòng)化所實(shí)習(xí)
已經(jīng)有7人回復(fù)
-
求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以
已經(jīng)有3人回復(fù)
-
材料專碩306英一數(shù)二
已經(jīng)有7人回復(fù)
-
材料工程專碩274一志愿211求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑
已經(jīng)有4人回復(fù)
-
332求調(diào)劑
已經(jīng)有5人回復(fù)
-
271求調(diào)劑
已經(jīng)有14人回復(fù)
-
275求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- matlab如何求解一個(gè)非線性微分方程組
已經(jīng)有8人回復(fù)
- 求助，用matlab解微分方程組，希望高手能給予指點(diǎn)
已經(jīng)有17人回復(fù)
- 如何用Matlab求時(shí)變微分方程的解？
已經(jīng)有5人回復(fù)
- matlab求解高階方程組
已經(jīng)有3人回復(fù)
- MATLAB求解非線性方程組
已經(jīng)有5人回復(fù)
- 關(guān)于matlab求解非線性指數(shù)方程組出現(xiàn)問(wèn)題
已經(jīng)有3人回復(fù)
- 有償求助：matlab怎么解聯(lián)立方程組？
已經(jīng)有11人回復(fù)
- 求高手解答Matlab求解微分方程的初始值和邊界值問(wèn)題
已經(jīng)有5人回復(fù)
- 幫忙用Mathematica求解一個(gè)微分方程組
已經(jīng)有6人回復(fù)
- MATLAB求解代數(shù)方程組問(wèn)題，請(qǐng)幫忙
已經(jīng)有9人回復(fù)
- Matlab求解偏微分方程組
已經(jīng)有22人回復(fù)
- 請(qǐng)問(wèn)如何求解二元一階微分方程組
已經(jīng)有24人回復(fù)
- 高金求助matlab解微分方程組
已經(jīng)有12人回復(fù)
- 求matlab求解一個(gè)方程組的問(wèn)題
已經(jīng)有6人回復(fù)
- 一個(gè)偏微分方程組（包括邊界條件）如何求解得到自振頻率和振型？
已經(jīng)有4人回復(fù)
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問(wèn)題
已經(jīng)有21人回復(fù)
- 【求助】請(qǐng)教非齊次常微分方程組的解析解法
已經(jīng)有4人回復(fù)
- 【求助】積分微分方程matlab求解
已經(jīng)有6人回復(fù)
- 【求助】matlab怎么求解偏微分方程組啊，先謝謝了
已經(jīng)有13人回復(fù)
- 【求助】求教這種微分方程的詳細(xì)求解過(guò)程
已經(jīng)有14人回復(fù)
- 【分享】KS方程求解器-matlab版
已經(jīng)有10人回復(fù)
wangyuewu08
木蟲 (著名寫手)
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 2794
- 紅花: 4
- 帖子: 1141
- 在線: 279.1小時(shí)
- 蟲號(hào): 1290073
- 注冊(cè): 2011-05-09
- 性別: GG
- 專業(yè): 制造系統(tǒng)與自動(dò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
|
我對(duì)微分方程組不是很了解，把樓主的代碼運(yùn)行了之后，系統(tǒng)出現(xiàn)了下面的錯(cuò)誤： 我試著把邊界條件中q的值改為1，運(yùn)行之后沒(méi)有報(bào)錯(cuò)，得到了兩個(gè)結(jié)構(gòu)體變量。 所以啊，問(wèn)題應(yīng)該出在q這個(gè)變量上，既然是數(shù)值解法，那么邊界條件中不應(yīng)該有變量。 |
2樓2011-12-20 18:38:48
3樓2011-12-20 19:08:24
wangyuewu08
木蟲 (著名寫手)
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 2794
- 紅花: 4
- 帖子: 1141
- 在線: 279.1小時(shí)
- 蟲號(hào): 1290073
- 注冊(cè): 2011-05-09
- 性別: GG
- 專業(yè): 制造系統(tǒng)與自動(dòng)化
4樓2011-12-20 20:10:30
5樓2011-12-20 21:20:14
dbb627
榮譽(yù)版主 (著名寫手)
-
專家經(jīng)驗(yàn): +4 - 仿真EPI: 6
- 應(yīng)助: 289 (大學(xué)生)
- 貴賓: 0.589
- 金幣: 24640.4
- 散金: 551
- 紅花: 61
- 沙發(fā): 1
- 帖子: 1246
- 在線: 1794.8小時(shí)
- 蟲號(hào): 149791
- 注冊(cè): 2005-12-29
- 性別: GG
- 專業(yè): 污染控制化學(xué)
- 管轄: 計(jì)算模擬
6樓2011-12-21 11:17:32
7樓2011-12-21 13:05:27
8樓2013-12-06 20:44:37
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 acnes 的主題更新
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 271求調(diào)劑 +12 | 生如夏花… 2026-03-11 | 14/700 |
|
|---|---|---|---|---|
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[文學(xué)芳草園] 伙伴們，祝我生日快樂(lè)吧 +17 | myrtle 2026-03-10 | 26/1300 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +3 | heming3743 2026-03-16 | 3/150 |
|
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
|
[考研] 招收0805（材料）調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
|
[考研] 求材料調(diào)劑 +5 | 隔壁陳先生 2026-03-12 | 5/250 |
|
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
|
[考研] 301求調(diào)劑 +6 | Liyouyumairs 2026-03-11 | 6/300 |
|
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|