my_lu.m

来自「包括使用修正Gram-Schmit算法实现QR分解」· M 代码 · 共 27 行

function [L,U]=my_lu(A)
% 自编函数, 用来计算矩阵A的LU分解
%   [L,U]=my_lu(A);
%   A为待分解的矩阵,A=LU,
%   其中L为单位下三角矩阵,U为阶梯型矩阵;

%   版权所有: gogomx, 2005.5.
%   Copyright (c) 2005 by gogomx.


[m,n]=size(A);
A_tmp=A;

L=eye(m);U=[];
for j=1:n
    if j==1
        v(j:n)=A_tmp(j:n,j);
    else 
        z=inv(L(1:j-1,1:j-1))*A_tmp(1:j-1,j);
        U(1:j-1,j)=z;
        v(j:n)=A_tmp(j:n,j)-L(j:n,1:j-1)*z;
    end
    if j<n
        L(j+1:n,j)=(v(j+1:n)/v(j))';
    end
    U(j,j)=v(j);
end