Question

如何使用matlab命令进行LU分解

Answer 1

1、双击matlab软件图标，打开matlab软件，可以看到matlab软件的界面。

2、使用函数pascal(5)创建5x5的矩阵A。

3、使用函数chol(A)对矩阵A进行Cholesky分解。

4、使用函数表达式[B,C]=lu(A)对矩阵进行LU分解，也成为高斯消去法。其中B是下三角矩阵，C是上三角矩阵。

5、使用函数magic(4)创建一个4x4的矩阵A。使用函数表达式[Q,R]=qr(A)，对矩阵A进行QR分解，其中Q是正交矩阵。

Answer 2

matlab有多种LU分解程序下面算一种：

function [L,U]=myLU(A)

%实现对矩阵A的LU分解，L为下三角矩阵

A

[n,n]=size(A);

L=zeros(n,n);

U=zeros(n,n);

for i=1:n

L(i,i)=1;

end

for k=1:n

for j=k:n

U(k,j)=A(k,j)-sum(L(k,1:k-1).*U(1:k-1,j)');

end

for i=k+1:n

L(i,k)=(A(i,k)-sum(L(i,1:k-1).*U(1:k-1,k)'))/U(k,k);

end

end

用法，在控制台输入

A=[1 2 3 -4;-3 -4 -12 13;2 10 0 -3;4 14 9 -13];

Answer 3

matlab自带函数lu就可以了~
用法:
[L,U] = lu(A);
得到的结果就满足 A=L*U.