24小时热门版块排行榜

返回列表

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

wangxingye

新虫 (小有名气)

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

[求助] MATLAB计算expm(At)的定积分已有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)，
不是一个数值解，不知道是怎么回事，希望高手指点一下！）

回复此楼

» 猜你喜欢

散金币祈福已经有187人回复
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有158人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
深圳大学应用技术学院招聘凝聚态物理博士后已经有115人回复
深圳大学2027级光电信息工程/物理学研究生推免已经有70人回复
求助中医院临床中医中标的青年标书学习，谢谢已经有4人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有257人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有45人回复

1楼 2017-06-23 22:56:24

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

wangxingye

新虫 (小有名气)

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

引用回帖:

6楼: Originally posted by huab1984666 at 2017-06-27 20:03:02
A=；
= eig(A*t);
expm_At=V*diag(exp(diag(D)))/V;
ZT=V.*diag(exp(diag(D)))./V;
zx=int(ZT,t,0,5)
zx1=double(zx);
____
1.50089258542985 0 0 0
0 1.50089258542985 0 0
0 0 77.7866509151555 0
0 0 ...

不好意思，别的事耽搁了，没能及时回复您！
还有几点不太明白，
一是expm_At=V*diag(exp(diag(D)))/V; 这句执行了一个多小时，MATLAB一直是busy，
得不到结果，这与我的软件是2012有关吗？还是别的其他原因？
二是为什么要按位运算求ZT，对ZT求积分，不应该是对expm_At求积分吗？
三是（2）中inlin=@(t) expm(A.*t);为什么也是用的按位乘法，是对于变量t就要用按位运算，而不能直接A*t吗？
我试了一下A*t运算报错。
还望您能指导一下！