【问题标题】:how to find the integral of a matrix exponential in python如何在python中找到矩阵指数的积分
【发布时间】:2013-01-31 16:33:37
【问题描述】:

我有一个矩阵形式,比如e^(Ax),其中A 是一个方阵。如何将它从给定值 a 集成到另一个值 b 以便输出是相应的数组?

【问题讨论】:

标签: python matrix numpy integration exponential


【解决方案1】:

如果 A 具有正确的属性,您可以通过计算其特征向量和特征值将其转换为对角线形式 A0。在对角线形式中,解是sol = [exp(A0*b) - exp(A0*a)] * inv(A0),其中A0 是具有特征值的对角矩阵,inv(A0) 仅包含其对角线中特征值的逆矩阵。最后,通过将解与左侧特征值和右侧特征值的转置相乘来转换解:transpose(eigvecs) * sol * eigvecs

【讨论】:

  • 非常感谢 Aradi 提供的解决问题的好方法和 tcaswell 的参考。但是,我怎样才能在 python 中实现这个,例如 scipy.integrate.quad 或类似的东西?我知道如何为单个函数而不是矩阵执行此操作。形成循环是唯一的解决方案吗?再次感谢。哦,a 和 b 是标量,是 x 的积分极限。
猜你喜欢
  • 1970-01-01
  • 1970-01-01
  • 2015-04-25
  • 1970-01-01
  • 1970-01-01
  • 2019-03-05
  • 2012-11-12
  • 2014-09-10
  • 2018-04-19
相关资源
最近更新 更多