| 5 | 1/1 | 返回列表 |
| 查看: 1484 | 回復(fù): 7 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
priest木蟲 (正式寫手)
|
[求助]
幫忙:求解一組微分代數(shù)方程,獎勵100金幣 已有1人參與
|
||
| 請各位幫忙求解一組微分代數(shù)方程,14個方程,其中7個微分方程,7個代數(shù)方程。我的思路是這樣的:把代數(shù)方程對時間求導(dǎo),轉(zhuǎn)化成微分方程,這樣方程組寫成 M(t,x)*dx/dt = F(t,x),可以用matlab ode15s求解,關(guān)鍵是寫出正確的mass matrix M。我遇到的問題就是:(1)matlab會顯示 This DAE appears to be of index greater than 1,但是我仔細(xì)檢查了代數(shù)方程的微分形式,例如,dx1-2*dx3 = x1+4*x2,就是所有代數(shù)方程的變量對時間的一階導(dǎo)數(shù)=f(x, dx/dt), 全部都是一階導(dǎo)數(shù)。所以,我認(rèn)為這組微分代數(shù)方程的index為1,不應(yīng)該有這個錯誤。(2)如果x1為狀態(tài)變量, M(1,1) = c1/x1, matlab 代碼里應(yīng)該如何表示,x1作為數(shù)組還是標(biāo)量,是要用c1./x1 還是c1/x1?或者,我可以在方程兩邊同時乘以x1?(3) 如果x1=c1*sinh(x2-x3), c1=exp(-c2*t), c1只是作為輸入函數(shù),這種情況下,x1對時間進(jìn)行微分,c1也要對對時間求導(dǎo)么,如果這樣,c1也變成了一個狀態(tài)變量。還是把c1當(dāng)作常熟處理,dx1/dt=c1*(sinh(x2-x3))’。多謝大家的幫助! |

榮譽版主 (文壇精英)
![]() |
專家經(jīng)驗: +518 |
木蟲 (正式寫手)
榮譽版主 (文壇精英)
![]() |
專家經(jīng)驗: +518 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 335求調(diào)劑 +3 | 321* 2026-03-31 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿085600中科院寧波所276分求調(diào)劑 +6 | 材料學(xué)257求調(diào)劑 2026-03-28 | 7/350 |
|
|
[考研] 本科211總分289,08工學(xué)真心求調(diào)劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
[考研] 279求調(diào)劑 +12 | j的立方 2026-03-29 | 12/600 |
|
|
[考研] 359求調(diào)劑 +5 | 王了個楠 2026-03-25 | 5/250 |
|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專碩 +10 | Koxui 2026-03-28 | 10/500 |
|
|
[考研] 一志愿鄭州大學(xué),080500學(xué)碩,總分317分求調(diào)劑 +10 | 舉個栗子oi 2026-03-24 | 11/550 |
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +18 | Mr. Z 2026-03-25 | 19/950 |
|
|
[考研] 284求調(diào)劑 +14 | junqihahaha 2026-03-26 | 15/750 |
|
|
[考研] 085701求調(diào)劑初試286分 +5 | secret0328 2026-03-28 | 5/250 |
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|
|
[考研] 348求調(diào)劑 +6 | 小懶蟲不懶了 2026-03-28 | 6/300 |
|
|
[考研]
|
nnnnnnn5 2026-03-25 | 11/550 |
|
|
[考研] 321求調(diào)劑 +7 | 璞玉~~ 2026-03-25 | 8/400 |
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么! 2026-03-26 | 7/350 |
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
[考博] 26申博 +3 | 加油沖。 2026-03-26 | 3/150 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|