三阶矩阵的lu分解详细步骤_矩阵的LU分解

矩阵的

LU

分解

一、题目

求一个

4

阶矩阵的

LU

分解。

(



























112

92

74

48

78

66

54

36

46

40

34

24

18

16

14

12

A

)

二、方法

Doolittle

(杜里特尔)分解法

三、程序

jiangLU.M

的程序如下：

function[L,U,flag]=jiangLU(A)

[n,n]=size(A);

L=eye(n);U=zeros(n);flag='OK';

for k=1:n

for j=k:n

z=0;

for q=1:k-1

z=z+L(k,q)*U(q,j);

end

U(k,j)=A(k,j)-z;

end

if abs(U(k,k))

flag='failure';return;

end

for i=k+1:n

z=0;

for q=1:k-1

z=z+L(i,q)*U(q,k);

end

L(i,k)=(A(i,k)-z)/U(k,k);

end

end

四、结果

>> format compact

>> A=[12 14 16 18;24 34 40 46;36 54 66 78;48 74 92 112]

A =

12

14

16

18

24

34

40

46

36

54

66

78