From 24b4bddcfec835760ed9496604a1ee677b080d1f Mon Sep 17 00:00:00 2001
From: pu3mybsfx <3034189236@qq.com>
Date: Wed, 24 May 2023 19:15:18 +0800
Subject: [PATCH] =?UTF-8?q?=E7=AC=AC=E5=9B=9B=E9=A2=98LU=E5=87=BD=E6=95=B0?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 lu_decompose.m | 32 ++++++++++++++++++++++++++++++++
 1 file changed, 32 insertions(+)
 create mode 100644 lu_decompose.m

diff --git a/lu_decompose.m b/lu_decompose.m
new file mode 100644
index 0000000..8123d22
--- /dev/null
+++ b/lu_decompose.m
@@ -0,0 +1,32 @@
+function [L,U] = lu_decompose(A)
+%   lu decompose
+%   L:下三角矩阵
+%	U:上三角矩阵
+%	A:输入矩阵
+[m,n]=size(A);
+
+L=eye(m);
+L(:,1)=A(:,1)/A(1,1);%L第一列赋值
+
+U=zeros(m,n);
+U(1,:)=A(1,:);%U第一行赋值
+
+for i=2:m
+    for j=2:n
+        if i<=j
+            U(i,j)=A(i,j)-sum(L(i,1:i-1).*U(1:i-1,j)');%递推表达式(1-6)
+        else
+            if U(j,j)==0
+                L(i,j)=0;
+            else
+                L(i,j)=(A(i,j)-sum(L(i,1:j-1).*U(1:j-1,j)'))/U(j,j);%递推表达式(1-7)
+            end
+        end
+    end
+end
+
+
+
+
+
+