Matlab学习 矩阵分解,特征值， 特征向量

矩阵分解

是把一个矩阵分解为多个矩阵连乘的形式

其中cholinc()和luinc()函数是针对稀疏矩阵的

1. 对称正定矩阵的Cholesky分解

a=pascal(4) % pascal矩阵是对称正定的
e=eig(a) % 返回由a的特征值构成的向量
b=chol(a)
c=b'*b

a =
     1     1     1     1
     1     2     3     4
     1     3     6    10
     1     4    10    20
e =
    0.0380
    0.4538
    2.2034
   26.3047
b =
     1     1     1     1
     0     1     2     3
     0     0     1     3
     0     0     0     1
c =
     1     1     1     1
     1     2     3     4
     1     3     6    10
     1     4    10    20
>> 

2. 一般方阵的高斯消去法分解(LU分解)

a=pascal(4) % pascal矩阵是对称正定的
[L1, U1]=lu(a)
b=L1*U1==a
[L2,U2,P]=lu(a)
c=L2*U2==P*a
Y=lu(a)
Y1=L2+U2-eye(size(a))

a =
     1     1     1     1