→ jatj:你寫錯了吧 應該是dy/dt=f(t)=一個矩陣之類的吧 10/30 10:58
→ leven13:不是 是f本身就是一個矩陣 矩陣的變數 10/30 22:26
推 profyang:所以你左邊的f根右邊的f是同一個未知矩陣?那不就A*e^t... 10/30 22:35
→ profyang:A是矩陣 10/30 22:35
恩 我的敘述方式不對 更改一下
f和g都是3*3的矩陣 這樣的情況下有辦法解嗎
※ 編輯: leven13 來自: 128.211.247.69 (10/30 23:03)
→ jatj:解答就是f=exp(g(t)) 收工 3F解答並不完整 10/31 02:44
推 profyang:我那個解答是他之前df/dt=f的解答 另外他現在df/dt=g的解 10/31 12:19
→ profyang:應該也是f=∫g(t)dt+C 不過他應該是希望matlab ode功能f 10/31 12:21
→ profyang:g都是矩陣他也可以一起算而不是像自己寫for loop那樣一個 10/31 12:21
→ profyang:個元素去積 10/31 12:21
→ jatj:矩陣的ode並不是一個個元素去積 每個元素是相關的... 10/31 13:42
→ jatj:所以我才說不完整 10/31 13:42
→ jatj:類似但不完全是原PO要的 10/31 13:44
推 profyang:我知道阿 可是他原文PO的問題沒有你貼的wiki那個A矩陣 10/31 15:57
→ profyang:他就只是df/dt=g f和g都是矩陣 除非他又PO錯 不然這個的 10/31 15:57
→ profyang:f的解就是把g每個對應元素去對t積分阿 不過這樣應該沒有 10/31 15:58
→ profyang:用matlab ode的功能去解的必要 10/31 15:59
→ leven13:感謝各位解答 我自己也想了一下 現在大概是用P大的解法 10/31 23:11
→ leven13:不過未來會碰到類似像j大那樣的情況 在矩陣相乘的情況下 10/31 23:12
→ leven13:解ode 到時候應該會用到wiki裡面的概念 10/31 23:13