 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1481 | 回復: 7 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
priest木蟲 (正式寫手)
|
[求助]
幫忙：求解一組微分代數(shù)方程，獎勵100金幣 已有1人參與
|
||
| 請各位幫忙求解一組微分代數(shù)方程，14個方程，其中7個微分方程，7個代數(shù)方程。我的思路是這樣的：把代數(shù)方程對時間求導，轉(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,但是我仔細檢查了代數(shù)方程的微分形式，例如，dx1-2*dx3 = x1+4*x2,就是所有代數(shù)方程的變量對時間的一階導數(shù)＝f(x, dx/dt), 全部都是一階導數(shù)。所以，我認為這組微分代數(shù)方程的index為1，不應(yīng)該有這個錯誤。（2）如果x1為狀態(tài)變量， M(1,1) = c1/x1, matlab 代碼里應(yīng)該如何表示，x1作為數(shù)組還是標量，是要用c1./x1 還是c1/x1？或者，我可以在方程兩邊同時乘以x1?（3） 如果x1=c1*sinh(x2-x3), c1=exp(-c2*t), c1只是作為輸入函數(shù)，這種情況下，x1對時間進行微分，c1也要對對時間求導么，如果這樣，c1也變成了一個狀態(tài)變量。還是把c1當作常熟處理，dx1/dt=c1*(sinh(x2-x3))’。多謝大家的幫助！ |
回復此樓» 猜你喜歡
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復
-
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有228人回復
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復
-
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復
-
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！
已經(jīng)有1人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復
-
計算與實驗科學工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學換能器分會S8
已經(jīng)有0人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 哪位大神幫忙看一下MATLAB求解微分方程的問題？
已經(jīng)有3人回復
- 求高手幫我解一階常微分方程組，非常感謝！！
已經(jīng)有5人回復
- 求解一個二階偏微分方程~~有結(jié)果，但是不知道中間過程（完美解決者答謝100￥）~
已經(jīng)有4人回復
- 微分方程組求參數(shù)問題，求高人指教，感謝
已經(jīng)有11人回復
- 求組！求解一個微分方程組的符號解，并畫出圖像�。�！
已經(jīng)有7人回復
- 幫忙求一個微分方程的數(shù)值解或者 真實解
已經(jīng)有6人回復
- 一階常微分方程組，matlab進行數(shù)值求解
已經(jīng)有3人回復
- MATLAB求解代數(shù)方程組問題，請幫忙
已經(jīng)有9人回復
- 一階微分方程求解，并優(yōu)化參數(shù)值
已經(jīng)有14人回復
- 一個微分方程組，求解
已經(jīng)有18人回復
- LABVIEW四階龍格庫塔法求解一階微分方程組
已經(jīng)有12人回復
- 求解一個含字母系數(shù)的微分方程
已經(jīng)有5人回復
- 誰那里可以運行maple，幫我解一個微分方程組
已經(jīng)有3人回復
- 【求助】20金幣求一個微積分方程
已經(jīng)有15人回復
- 【求助】求解以二階微分方程，贈金幣��！
已經(jīng)有13人回復
feixiaolin
榮譽版主 (文壇精英)
-
專家經(jīng)驗: +518 - 應(yīng)助: 942 (博后)
- 貴賓: 1.275
- 金幣: 3880
- 散金: 58785
- 紅花: 532
- 沙發(fā): 11
- 帖子: 24215
- 在線: 2601.8小時
- 蟲號: 2139575
- 注冊: 2012-11-21
- 專業(yè): 光學信息獲取與處理
- 管轄: 數(shù)學
3樓2014-09-27 08:16:11
priest
木蟲 (正式寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 1743.8
- 散金: 20
- 紅花: 2
- 帖子: 658
- 在線: 463.7小時
- 蟲號: 442901
- 注冊: 2007-10-27
- 性別: GG
- 專業(yè): 催化化學
4樓2014-09-27 17:49:03
feixiaolin
榮譽版主 (文壇精英)
-
專家經(jīng)驗: +518 - 應(yīng)助: 942 (博后)
- 貴賓: 1.275
- 金幣: 3880
- 散金: 58785
- 紅花: 532
- 沙發(fā): 11
- 帖子: 24215
- 在線: 2601.8小時
- 蟲號: 2139575
- 注冊: 2012-11-21
- 專業(yè): 光學信息獲取與處理
- 管轄: 數(shù)學
5樓2014-09-27 17:56:06
priest
木蟲 (正式寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 1743.8
- 散金: 20
- 紅花: 2
- 帖子: 658
- 在線: 463.7小時
- 蟲號: 442901
- 注冊: 2007-10-27
- 性別: GG
- 專業(yè): 催化化學
|
(∂(Fx(14)x(1)))/∂t=(-3x(10)-x(8))*ρ_c ∂(Fx(14)x(2))/∂t=(2x(10)-x(9))*ρ_c x(9)=k_K x 〖(12)sinh〗(((x(5)-x(7)))⁄〖2v〗_T ) x(8)=k_Na x(13)x(11)sinh〖(((x(5)-x(6)))⁄〖2v〗_T ) 〗 x(10)=k_NaK tanh(((x(5)+2x(7)-3x(6)-v_ATP ))⁄(2v_T )) (dx(12))/dt=1/τ_K cosh(2 ((x(5)-v_x ))⁄v_T ) {0.5*[1+tanh(2 ((x(5)-v_x ))⁄v_T ) ]-x(12)} (dx(13))/dt=1/τ_Na cosh(2 ((x(5)-h))⁄v_T ) {0.5*[1+tanh(2 ((x(5)-v_h ))⁄v_T ) ]-x(13)} x(11)=0.5*(1+tanh(2(x(5)-v_m )/v_T ) ) x(5)=v_T ln〖(x(3))⁄(x(1))〗 x(6)=v_T ln〖(x(4))⁄(x(2))〗 x(5)=RT/F ln((P_K x(4)+P_Na x(2)+P_Cl 〖Cl〗_i)/(P_K x(3)+P_Na x(1)+P_Cl 〖Cl〗_e )) (∂((1-x(14))x(4)))/∂t=∇D_2 ∇x(4)-∂(x(14)x(2))/∂t (∂((1-x(14))x(3)))/∂t=∇D_1 ∇x(3)-∂(x(14)x(1))/∂t (∂x(14))/∂t=L_p [(x(1)+x(2))-(x(3)+x(4)+Ci)] 全部14個方程如上，x是狀態(tài)變量，其他都是參數(shù)（常數(shù)）。盡管x（3）和x（4）隨空間變換，但可以通過有限差分進行離散。多謝啦！ |
6樓2014-09-28 06:43:48
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856材料化工調(diào)劑 總分330 +14 | zhubinhao 2026-03-27 | 14/700 |
|
|---|---|---|---|---|
|
|
[考研] 085701求調(diào)劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
|
[考研] 275求調(diào)劑 +10 | Micky11223 2026-03-25 | 14/700 |
|
|
|
[考研] 材料與化工272求調(diào)劑 +9 | 阿斯蒂芬2004 2026-03-28 | 9/450 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
|
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
|
[考研] 359求調(diào)劑 +4 | 王了個楠 2026-03-25 | 4/200 |
|
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
|
[考研] 【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】 +4 | Higraduate 2026-03-22 | 8/400 |
|
|
|
[考研] 340求調(diào)劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
|
[考研] 化學調(diào)劑一志愿上海交通大學336分-本科上海211 +4 | 小魚愛有機 2026-03-25 | 4/200 |
|
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
|
[考研] 機械學碩總分317求調(diào)劑�。。�！ +4 | Acaciad 2026-03-25 | 4/200 |
|
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +6 | 研研，接電話 2026-03-24 | 7/350 |
|
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|