手把手教你用matlab进行矩阵LDL

1、开始第一步我们打开在电脑桌面找到matlab小程序，然后鼠标右击打开桌面上matlab程序，运行起来。由于不同人电脑的配置不一样，软件打开的速度也有所不同，一般固态硬盘比机械硬盘运行的要快好多。大家稍微等待一下。

2、我们为了便于保存数据，我们点击matlab左上角新建脚本命令，创建新的脚本，创建M文件，也便于程序的保存，我们可以将其保存在电脑的其他盘，以便于我们的寻找和使用，这也是比较常见的方式。

3、在新建脚本命令输入窗口输入程序

>>clear all;

注意书写的方式哦

关闭所有正在运行的matlab程序，以便于我们的程序的运行和使用

也可以让我们得到的结果更加的准确

4、输入程序：

A=[1 2 3 4;4 6 10 2;1 1 0 1;0 0 2 3];

这是一个4*4的矩阵

注意矩阵表达方式

5、输入程序：

[L,D]=ldl(A)

若A为一对称矩阵且其任意一k阶主子阵均不为零，则A有如下惟一的分解形式：

A=LDL^T

求解矩阵的分解函数ldl

调用格式

Y=ldl（X）

Y和X都是矩阵

6、其中L的值为：

L =

    0.6667    1.0000         0         0

    1.0000         0         0         0

    0.1667   -0.2000    0.6667    1.0000

         0         0    1.0000         0

7、其中D值为

D =

    6.0000         0         0         0

         0   -1.6667         0         0

         0         0    3.0000         0

         0         0         0   -1.4333

8、我们验证下是不是分解正确

输入程序

>> L*D*L'

9、运行结果：

>> L*D*L'

ans =

     1     4     1     0

     4     6     1     0

     1     1     0     2

     0     0     2     3

10、和初始值一样，计算正确