应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 计算expm(At)的定积分，A是4阶方阵，t是变量
71/1返回列表
查看: 1178  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wangxingye

新虫 (小有名气)

[求助] 计算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)，
不是一个数值解，不知道是怎么回事，  希望高手指点一下！）
回复此楼

» 猜你喜欢
1楼 2017-06-24 16:58:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赵_山河

新虫 (初入文坛)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
wangxingye: 金币+2, ★有帮助 2017-06-25 17:15:12
A=[1 2;3 4];
syms t
eval(int(expm(A*t),t,0,0.05))
赞一下 回复此楼
洪兴掌管一带
2楼2017-06-24 18:02:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangxingye

新虫 (小有名气)
引用回帖:
2楼: Originally posted by 赵_山河 at 2017-06-24 18:02:41
A=;
syms t
eval(int(expm(A*t),t,0,0.05))

感谢您的帮助，我也试过这种方法，2*2的矩阵很快就能运行出结果，但是A如果是4*4的矩阵的话，MATLAB运行很长时间也出不来答案，一直busy，这个怎么解决？
赞一下 回复此楼
3楼2017-06-25 17:18:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赵_山河

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by wangxingye at 2017-06-25 17:18:07
感谢您的帮助，我也试过这种方法，2*2的矩阵很快就能运行出结果，但是A如果是4*4的矩阵的话，MATLAB运行很长时间也出不来答案，一直busy，这个怎么解决？...

换Mathematica。要不你问问我老大陈浩南？
赞一下 回复此楼
洪兴掌管一带
4楼2017-06-25 17:29:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangxingye

新虫 (小有名气)
引用回帖:
4楼: Originally posted by 赵_山河 at 2017-06-25 17:29:43
换Mathematica。要不你问问我老大陈浩南？...

惭愧啊，没用过Mathematica，不知道MATLAB能不能解决这个问题，如果你老大陈浩南能帮忙的话，那就更好了
赞一下 回复此楼
5楼2017-06-25 17:37:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赵_山河

新虫 (初入文坛)

【答案】应助回帖

★ ★
wangxingye: 金币+2, ★有帮助 2017-06-25 19:03:27
引用回帖:
5楼: Originally posted by wangxingye at 2017-06-25 17:37:48
惭愧啊，没用过Mathematica，不知道MATLAB能不能解决这个问题，如果你老大陈浩南能帮忙的话，那就更好了...

[V,D]=eig(A*t)
耗时发生在求inv(V)上，看expm的表达式：
[V,D] = eig(X)
expm(X) = V*diag(exp(diag(D)))/V
/V=*inv(V)，而V过于复杂。
要是简单求e对每个元素的指数，用exp(X)才对。
赞一下 回复此楼
洪兴掌管一带
6楼2017-06-25 18:04:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangxingye

新虫 (小有名气)
引用回帖:
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的定积分，但是得到的结果和正确答案相差很大，也不知道是哪里出了问题？
赞一下 回复此楼
7楼2017-06-25 19:44:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangxingye 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 情人节自我反思：在爱情中有过遗憾吗？ +5 瞬息宇宙 2026-02-15 6/300 2026-02-18 12:51 by 月下雪林
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:54 by lqtl9djx19
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:38 by lqtl9djx19
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 07:55 by lotyj5cz79
[基金申请] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:08 by lotyj5cz79
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 3/150 2026-02-18 06:53 by lotyj5cz79
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-17 23:40 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错，点赞！ +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题，请指教！ 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
关闭
请填处理意见
关闭