24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wangxingye

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1510.5
散金: 6
红花: 1
帖子: 97
在线: 49.7小时
虫号: 3908308
注册: 2015-06-04
性别: GG
专业: 机械工程

[求助] 计算expm(At)的定积分，A是4阶方阵，t是变量已有1人参与

想用MATLAB计算expm(At)的定积分，A是4阶方阵，t是变量。
随便取个简单的例子，比如A=[1 1 2 1; 0 1 2 1; 1 0 1 2; 0 1 1 1]；
t的变化范围是[0, 0.05]，如何用MATLAB计算expm(A*t)的积分值呢？
（Ps，我试着用int('expm(A*t)', 't', 0, 0.05)计算，输出的结果是int(expm(A*s), s == 0..1/20)，
不是一个数值解，不知道是怎么回事，希望高手指点一下！）

回复此楼

» 猜你喜欢

面上没中，邀请各位路过的虫友分析一下分数已经有5人回复
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生已经有3人回复
找导师已经有8人回复
国社科又开始会评了，不知道这次命运如何已经有14人回复
麻烦专家们看看评委们的意见（F口面上）已经有13人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有14人回复
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够已经有17人回复
学科评审组评审是指会评吗？已经有5人回复
小白求助投论文要求的highlights应该如何写已经有4人回复
面上意见出来了已经有23人回复

1楼 2017-06-24 16:58:01

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wangxingye

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1510.5
散金: 6
红花: 1
帖子: 97
在线: 49.7小时
虫号: 3908308
注册: 2015-06-04
性别: GG
专业: 机械工程

引用回帖:

6楼: Originally posted by 赵_山河 at 2017-06-25 18:04:18
=eig(A*t)
耗时发生在求inv(V)上，看expm的表达式：
= eig(X)
expm(X) = V*diag(exp(diag(D)))/V
/V=*inv(V)，而V过于复杂。
要是简单求e对每个元素的指数，用exp(X)才对。...

我按照你的方法试了一下，inv(V)耗时稍长，求exp(diag(D))耗时更长一些，但也还能算出来，
而求最耗时的是发生在积分int上，运行了半个多小时，一直busy，得不出结果，
而我的目的也是求expm(At)的积分，不是简单求e对每个元素的积分。
其实之前参考书上的解法也试了一下，和你的方法有些类似：
[V,D]=eig(A)
B=V*exp(A*t)*inv(V)
B=vga(B,6)
然后，手算了B中每个元素对t的定积分，但是得到的结果和正确答案相差很大，也不知道是哪里出了问题？