 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4737 | 回復(fù): 7 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[求助]
matlab數(shù)值求解邊界條件微分方程組
|
||
期望高手指點一下我的疑問，指出一個大方向也好，我現(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方法，我進行了如下處理： 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ù)此樓» 猜你喜歡
生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338
已經(jīng)有6人回復(fù)
-
321求調(diào)劑
已經(jīng)有7人回復(fù)
-
343求調(diào)劑
已經(jīng)有4人回復(fù)
-
305求調(diào)劑
已經(jīng)有3人回復(fù)
-
本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑
已經(jīng)有6人回復(fù)
-
332求調(diào)劑
已經(jīng)有3人回復(fù)
-
學(xué)碩274求調(diào)劑
已經(jīng)有8人回復(fù)
-
壓國家一區(qū)線，求導(dǎo)師收留，有恩必謝！
已經(jīng)有7人回復(fù)
-
081200-11408-276學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
22408 359分調(diào)劑
已經(jīng)有4人回復(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ù)
- 【求助】求教這種微分方程的詳細求解過程
已經(jīng)有14人回復(fù)
- 【分享】KS方程求解器-matlab版
已經(jīng)有10人回復(fù)
7樓2011-12-21 13:05:27
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
|
我對微分方程組不是很了解，把樓主的代碼運行了之后，系統(tǒng)出現(xiàn)了下面的錯誤： 我試著把邊界條件中q的值改為1，運行之后沒有報錯，得到了兩個結(jié)構(gòu)體變量。 所以啊，問題應(yīng)該出在q這個變量上，既然是數(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小時
- 蟲號: 1290073
- 注冊: 2011-05-09
- 性別: GG
- 專業(yè): 制造系統(tǒng)與自動化
4樓2011-12-20 20:10:30
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 343求調(diào)劑 +4 | 愛羈絆 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 312，生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
|
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +5 | 芝士雪baoo 2026-03-28 | 6/300 |
|
|
|
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
|
[考研] 085600，材料與化工321分求調(diào)劑 +9 | 大饞小子 2026-03-28 | 9/450 |
|
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
|
[考研] 277跪求調(diào)劑 +5 | 1915668 2026-03-27 | 9/450 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
|
[考研] 085600，材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
|
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
|
[考研] 303求調(diào)劑 +6 | 藍山月 2026-03-25 | 6/300 |
|
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達問津 2026-03-23 | 5/250 |
|
|
|
[考研]
 一志愿070300浙大化學(xué)358分，求調(diào)劑！
+4
|
酥酥魚.. 2026-03-21 | 4/200 |
|