代码之家 › 专栏 › 技术社区 › jett chen

如何在sympy中生成矩阵函数?

sympy python

jett chen · 技术社区 · 6 月前

我定义 f(At) = e^(At) 哪里 A = [[0, 1], [-1, 0]]

from sympy import *
from sympy.abc import x,y

t = symbols("t")
A = Matrix([
    [0, 1],
    [-1, 0]
])*t

A1 = A.exp()

A1 输出 Matrix([[exp(I*t)/2 + exp(-I*t)/2, -I*exp(I*t)/2 + I*exp(-I*t)/2], [I*exp(I*t)/2 - I*exp(-I*t)/2, exp(I*t)/2 + exp(-I*t)/2]]) 当我希望输出为 [ [cos(t), sin(t)], [-sin(t), cos(t)] ] .我怎样才能得到那个结果?

1 回复 | 直到 6 月前

Tim Roberts 6 月前

你需要逐一简化结果:

A1 = A.exp()
A1 = A1.applyfunc(simplify)
print(A1)

输出:

Matrix([[cos(t), sin(t)], [-sin(t), cos(t)]])

推荐文章

Google User · Django管理员在`list_display中未显示`creation_date`字段`

5 月前

user29747013 · 如何创建一个新的数据框架,其中包含原始数据框架中列的聚合列?

5 月前

ÎÎÎ½Î· ÎÎ®Î¹Î½Î¿Ï · Python lxml.html语法错误:使用lxml find时XPATH的谓词无效

5 月前

user29715306 · from_users=和chats=电视节目中的差异

5 月前

Redshoe · 当执行numpy.genfromtxt()时,python是否会读取文件的所有行?

5 月前

RASEL MAHMUD · 为什么以及如何在is_even()函数内的IF条件中递归X变量在满足0后递增?[副本]

5 月前

prayner · 更新嵌套字典包含列表中的项

5 月前

Bringo Jr · 我可以在O(n)中解决这个问题吗?

5 月前

Dave · 如何在for循环中修改列表值

5 月前

Shukurullox Komiljonov · 从记录中获得相互和解。使用SQL

5 月前