[問題]如何用ode解矩陣

看板MATLAB作者 (楊麗花)時間12年前 (2013/10/30 08:59), 編輯推噓3(3推 0噓 17→)
留言20則, 3人參與, 最新討論串1/1
如果今天一個ode方程式 df/dt=g(t) 而方程式f和g是個3*3的矩陣 這樣matlab有辦法解嗎 目前只用過ode解過純量解 如果是3*3的矩陣該怎麼做 把單獨每一行(或每一列)單獨抓出來跑ode嗎 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 128.211.247.69
→
10/30 10:58, , 1F
你寫錯了吧 應該是dy/dt=f(t)=一個矩陣之類的吧
10/30 10:58, 1F
→
10/30 22:26, , 2F
不是 是f本身就是一個矩陣 矩陣的變數
10/30 22:26, 2F
推
10/30 22:35, , 3F
所以你左邊的f根右邊的f是同一個未知矩陣?那不就A*e^t...
10/30 22:35, 3F
→
10/30 22:35, , 4F
A是矩陣
10/30 22:35, 4F
恩 我的敘述方式不對 更改一下 f和g都是3*3的矩陣 這樣的情況下有辦法解嗎 ※ 編輯: leven13 來自: 128.211.247.69 (10/30 23:03)
→
10/31 02:44, , 5F
解答就是f=exp(g(t)) 收工 3F解答並不完整
10/31 02:44, 5F
推
10/31 12:19, , 6F
我那個解答是他之前df/dt=f的解答 另外他現在df/dt=g的解
10/31 12:19, 6F
→
10/31 12:21, , 7F
應該也是f=∫g(t)dt+C 不過他應該是希望matlab ode功能f
10/31 12:21, 7F
→
10/31 12:21, , 8F
g都是矩陣他也可以一起算而不是像自己寫for loop那樣一個
10/31 12:21, 8F
→
10/31 12:21, , 9F
個元素去積
10/31 12:21, 9F
→
10/31 13:42, , 10F
矩陣的ode並不是一個個元素去積 每個元素是相關的...
10/31 13:42, 10F
→
10/31 13:42, , 11F
所以我才說不完整
10/31 13:42, 11F
→
10/31 13:44, , 13F
類似但不完全是原PO要的
10/31 13:44, 13F
推
10/31 15:57, , 14F
我知道阿 可是他原文PO的問題沒有你貼的wiki那個A矩陣
10/31 15:57, 14F
→
10/31 15:57, , 15F
他就只是df/dt=g f和g都是矩陣 除非他又PO錯 不然這個的
10/31 15:57, 15F
→
10/31 15:58, , 16F
f的解就是把g每個對應元素去對t積分阿 不過這樣應該沒有
10/31 15:58, 16F
→
10/31 15:59, , 17F
用matlab ode的功能去解的必要
10/31 15:59, 17F
→
10/31 23:11, , 18F
感謝各位解答 我自己也想了一下 現在大概是用P大的解法
10/31 23:11, 18F
→
10/31 23:12, , 19F
不過未來會碰到類似像j大那樣的情況 在矩陣相乘的情況下
10/31 23:12, 19F
→
10/31 23:13, , 20F
解ode 到時候應該會用到wiki裡面的概念
10/31 23:13, 20F
‣ 更多 leven13 的文章
文章代碼(AID): #1IS5dh3b (MATLAB)