matlab自定义矩阵分解,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

1 2 3 4

1 3 6 10

1 4 10 20

L1 =

1.0000 0 0 0

1.0000 0.3333 1.0000 0

1.0000 0.6667 1.0000 1.0000

1.0000 1.0000 0 0

U1 =

1.0000 1.0000 1.0000 1.0000

0 3.0000 9.0000 19.0000

0 0 -1.0000 -3.3333

0 0 0 -0.3333

b =

4×4 logical 数组

1 1 1 1

1 1 1 1

1 1 1 1

1 1 1 1

L2 =

1.0000 0 0 0

1.0000 1.0000 0 0

1.0000 0.3333 1.0000 0

1.0000 0.6667 1.0000 1.0000

U2 =

1.0000 1.0000 1.0000 1.0000

0 3.0000 9.0000 19.0000

0 0 -1.0000 -3.3333

0 0 0 -0.3333

P =

1 0 0 0

0 0 0 1

0 1 0 0

0 0 1 0

c =

4×4 logical 数组

1 1 1 1

1 1 1 1

1 1 1 1

1 1 1 1

Y =

1.0000 1.0000 1.0000 1.0000

1.0000 3.0000 9.0000 19.0000

1.0000 0.3333 -1.0000 -3.3333

1.0000 0.6667 1.0000 -0.3333

Y1 =

1.0000 1.0000 1.0000 1.0000

1.0000 3.0000 9.0000 19.0000

1.0000 0.3333 -1.0000 -3.3333

1.0000 0.6667 1.0000 -0.3333

>>

3. 舒尔分解

4. 矩形矩阵的正交分解

标签：10,9.0000,1.0000,矩阵,分解,0.3333,Matlab

来源： https://blog.csdn.net/qq_36607894/article/details/99889144