如题,当已知t的值的时候可以求,使用scipy中的expm方法,当t是带数符号怎么办。A为已知的六阶矩阵。
使用是用于符号数学的Python库SymPy来操作,示例如下,你可将其中的A替换成六阶阵,代入不同的t值测试计算结果:
from sympy import symbols,Matrix,zeros,eye,simplify,exp
t=symbols('t')
A = Matrix([[16, 1], [0, 16]])
y=exp(A*t)
print(y.evalf(subs={t:3}))