求矩阵的特征值和特征向量
代码如下
#
import sympy
A=Matrix([[0,2,1],[2,6,1],[1,1,4]])
#特征值
lamda=symbols('lamda')
lamdaA=lamda*eye(3)-A
det_lamdaA=lamdaA.det()
A24_eig=list(solveset(Eq(det_lamdaA,0),lamda))
A24_eig
为什么在输出的时候不是数学表达式
输出如下:
[-2*sqrt(46)*cos(atan(9*sqrt(1167)/53)/3)/3 + 10/3, -sqrt(138)*sin(atan(9*sqrt(1167)/53)/3)/6 - 46*re(1/((-1/2 - sqrt(3)*I/2)*(53 + 9*sqrt(1167)*I)**(1/3)))/3 + sqrt(46)*cos(atan(9*sqrt(1167)/53)/3)/6 + 10/3 + I*(-46*im(1/((-1/2 - sqrt(3)*I/2)*(53 + 9*sqrt(1167)*I)**(1/3)))/3 + sqrt(46)*sin(atan(9*sqrt(1167)/53)/3)/6 + sqrt(138)*cos(atan(9*sqrt(1167)/53)/3)/6), sqrt(138)*sin(atan(9*sqrt(1167)/53)/3)/6 + sqrt(46)*cos(atan(9*sqrt(1167)/53)/3)/6 - 46*re(1/((-1/2 + sqrt(3)*I/2)*(53 + 9*sqrt(1167)*I)**(1/3)))/3 + 10/3 + I*(-sqrt(138)*cos(atan(9*sqrt(1167)/53)/3)/6 + sqrt(46)*sin(atan(9*sqrt(1167)/53)/3)/6 - 46*im(1/((-1/2 + sqrt(3)*I/2)*(53 + 9*sqrt(1167)*I)**(1/3)))/3)]