24小时热门版块排行榜    

查看: 1552  |  回复: 11

wiu2

新虫 (初入文坛)

[求助] 急求e^(At)的定积分matlab程序

e^(At)的定积分matlab程序,其中t[0,2],A为矩阵,怎么求出它的数值解?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

信彼南山

木虫 (著名写手)

引用回帖:
10楼: Originally posted by wiu2 at 2013-08-18 18:43:08
所以我才问如何离散化,得到积分的数值解啊!!!...

我发现我的表达能力似乎有问题啊。
这个不能把A与t当矩阵来相乘
你要是非按矩阵算的话,那应该

X(t)是n×1×m的
t是1×m的
A是n×n×1的,A*t是n×n×m的
X(0)是n×1×m的

这块我自己已经糊涂了

11楼2013-08-18 19:39:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

信彼南山

木虫 (著名写手)

做展开呗
2楼2013-08-17 12:38:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wiu2

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 信彼南山 at 2013-08-17 12:38:45
做展开呗

展开之后积分里面它还是一个矩阵啊,求数值解,将时间t离散化后,相乘后还是会出现维数不匹配的问题啊~怎么办??
3楼2013-08-17 14:43:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)

★ ★
xiegangmai: 金币+2, 鼓励讨论交流! 2013-08-24 21:51:58
你的公式应该是
CODE:
X(t)=e^(A*t)*X(0)+int

的形式吧

或者考虑用拉普拉斯变换呗。
4楼2013-08-17 18:28:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wiu2

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by 信彼南山 at 2013-08-17 18:28:36
你的公式应该是

X(t)=e^(A*t)*X(0)+int

的形式吧

或者考虑用拉普拉斯变换呗。

额。。。可它还是解析形式,,我想要能不能直接算出数值解,给定A后,就接得到一个数值矩阵~
5楼2013-08-17 20:46:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)

X(t)是n×1的
A是n×n的,e^(A*t)是n×n的
X(0)是n×1的

这有什么问题么
6楼2013-08-17 22:10:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)

你的意思不是说A跟t相乘的时候维数不匹配吧????
7楼2013-08-17 22:12:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wiu2

新虫 (初入文坛)

引用回帖:
7楼: Originally posted by 信彼南山 at 2013-08-17 22:12:59
你的意思不是说A跟t相乘的时候维数不匹配吧????

t离散化之后并不是nx1的啊。。比如说A为3x3,t不可能也只化为3步吧,取步长为0.1的话,t离散为[0:0.1:2]...维数是不匹配的啊!!
8楼2013-08-18 08:36:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)

引用回帖:
8楼: Originally posted by wiu2 at 2013-08-18 08:36:19
t离散化之后并不是nx1的啊。。比如说A为3x3,t不可能也只化为3步吧,取步长为0.1的话,t离散为...维数是不匹配的啊!!...

果然是这样。
你这数学跟谁学的啊。
A跟t是不能这么相乘的
9楼2013-08-18 11:19:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wiu2

新虫 (初入文坛)

引用回帖:
9楼: Originally posted by 信彼南山 at 2013-08-18 11:19:18
果然是这样。
你这数学跟谁学的啊。
A跟t是不能这么相乘的...

所以我才问如何离散化,得到积分的数值解啊!!!
10楼2013-08-18 18:43:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wiu2 的主题更新
信息提示
请填处理意见