 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4734 | 回復(fù): 7 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[求助]
matlab數(shù)值求解邊界條件微分方程組
|
||
期望高手指點一下我的疑問，指出一個大方向也好，我現(xià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ù)此樓» 猜你喜歡
0854電子信息求調(diào)劑 324
已經(jīng)有3人回復(fù)
-
303求調(diào)劑
已經(jīng)有3人回復(fù)
-
336求調(diào)劑
已經(jīng)有3人回復(fù)
-
306求調(diào)劑
已經(jīng)有9人回復(fù)
-
工科0856求調(diào)劑
已經(jīng)有4人回復(fù)
-
工科材料085601 279求調(diào)劑
已經(jīng)有10人回復(fù)
-
081700 調(diào)劑 267分
已經(jīng)有4人回復(fù)
-
276求調(diào)劑。有半年電池和半年高分子實習(xí)經(jīng)歷
已經(jīng)有9人回復(fù)
-
263求調(diào)劑
已經(jīng)有9人回復(fù)
-
求老師收我
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- matlab如何求解一個非線性微分方程組
已經(jīng)有8人回復(fù)
- 求助，用matlab解微分方程組，希望高手能給予指點
已經(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ù)
3樓2011-12-20 19:08:24
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
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
5樓2011-12-20 21:20:14
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0854電子信息求調(diào)劑 324 +3 | Promise-jyl 2026-03-23 | 3/150 |
|
|---|---|---|---|---|
|
[考研]
 材料與化工考研調(diào)劑
+3
|
孅華 2026-03-22 | 3/150 |
|
|
|
[考研] 070300化學(xué)求調(diào)劑 +5 | 苑豆豆 2026-03-20 | 5/250 |
|
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +5 | vv迷 2026-03-22 | 5/250 |
|
|
|
[考研] 317求調(diào)劑 +12 | 申子申申 2026-03-19 | 18/900 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 310求調(diào)劑 +4 | baibai1314 2026-03-16 | 4/200 |
|
|
|
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑 +4 | 幸運(yùn)的醬醬 2026-03-22 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|