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是对角阵，这个意思么？这个矩阵论或线性代数都有啊，你可以查一下

3楼: Originally posted by FMStation at 2016-09-12 10:25:48
http://www.ece.utah.edu/eceCTools/LinearAlgebra/EigenValVec/Diagonalization/LinAlgEigenDiagMatlabEx.pdf


syms a b c d
syms A S D

A =
A =



= eig(A)
V =


D =



Snum = '
Snum ...

2楼: Originally posted by yogort at 2016-09-11 16:20:36
rref或rrefmovie

5楼: Originally posted by guanmm at 2016-09-12 16:13:45
不懂，学习了