syms a b c d
syms A S D

A = [a, b; c, d]
A =
[ a, b]
[ c, d]

[V,D] = eig(A)
V =
[ -(-1/2*a+1/2*d-1/2*(a^2-2*a*d+d^2+4*b*c)^(1/2))/c, -(- 1/2*a+1/2*d+1/2*(a^2-2*a*d+d^2+4*b*c)^(1/2))/c]
[ 1, 1]
D =
[ 1/2*a+1/2*d+1/2*(a^2-2*a*d+d^2+4*b*c)^(1/2), 0]
[ 0, 1/2*a+1/2*d-1/2*(a^2-2*a*d+d^2+4*b*c)^(1/2)]

Snum = [-1,1;1,0]'
Snum =
-1 1
1 0

Dnum = [2, 0; 0, 3]
Dnum =
2 0
0 3

Anum = Snum * Dnum * inv(Snum)
Anum =
3 1
0 2

[Vnum,Dnum] = eig(Anum)
Vnum =
1 -0.70711
0 0.70711
Dnum =
3 0
0 2

4楼: Originally posted by 我怀念的TQ at 2016-09-12 15:52:58
将一个矩阵变成对角阵？？？是P‘AP=B，B是对角阵，这个意思么？这个矩阵论或线性代数都有啊，你可以查一下