作者leven13 (杨丽花)
看板MATLAB
标题[问题]如何用ode解矩阵
时间Wed Oct 30 08:59:21 2013
如果今天一个ode方程式
df/dt=g(t)
而方程式f和g是个3*3的矩阵
这样matlab有办法解吗
目前只用过ode解过纯量解
如果是3*3的矩阵该怎麽做
把单独每一行(或每一列)单独抓出来跑ode吗
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 128.211.247.69
1F：→ jatj:你写错了吧 应该是dy/dt=f(t)=一个矩阵之类的吧 10/30 10:58
2F：→ leven13:不是 是f本身就是一个矩阵 矩阵的变数 10/30 22:26
3F：推 profyang:所以你左边的f根右边的f是同一个未知矩阵?那不就A*e^t... 10/30 22:35
4F：→ profyang:A是矩阵 10/30 22:35
恩 我的叙述方式不对 更改一下
f和g都是3*3的矩阵 这样的情况下有办法解吗
※ 编辑: leven13 来自: 128.211.247.69 (10/30 23:03)
5F：→ jatj:解答就是f=exp(g(t)) 收工 3F解答并不完整 10/31 02:44
6F：推 profyang:我那个解答是他之前df/dt=f的解答 另外他现在df/dt=g的解 10/31 12:19
7F：→ profyang:应该也是f=∫g(t)dt+C 不过他应该是希望matlab ode功能f 10/31 12:21
8F：→ profyang:g都是矩阵他也可以一起算而不是像自己写for loop那样一个 10/31 12:21
9F：→ profyang:个元素去积 10/31 12:21
10F：→ jatj:矩阵的ode并不是一个个元素去积 每个元素是相关的... 10/31 13:42
11F：→ jatj:所以我才说不完整 10/31 13:42
13F：→ jatj:类似但不完全是原PO要的 10/31 13:44
14F：推 profyang:我知道阿 可是他原文PO的问题没有你贴的wiki那个A矩阵 10/31 15:57
15F：→ profyang:他就只是df/dt=g f和g都是矩阵 除非他又PO错 不然这个的 10/31 15:57
16F：→ profyang:f的解就是把g每个对应元素去对t积分阿 不过这样应该没有 10/31 15:58
17F：→ profyang:用matlab ode的功能去解的必要 10/31 15:59
18F：→ leven13:感谢各位解答 我自己也想了一下 现在大概是用P大的解法 10/31 23:11
19F：→ leven13:不过未来会碰到类似像j大那样的情况 在矩阵相乘的情况下 10/31 23:12
20F：→ leven13:解ode 到时候应该会用到wiki里面的概念 10/31 23:13