matlab实现矩阵LU分解

最新推荐文章于 2021-03-16 17:23:33 发布

damant

最新推荐文章于 2021-03-16 17:23:33 发布

阅读量3.5k

点赞数 2

分类专栏： matlab 矩阵

本文链接：https://blog.csdn.net/damant/article/details/53621294

版权

matlab 同时被 2 个专栏收录

7 篇文章 1 订阅

订阅专栏

矩阵

3 篇文章 0 订阅

订阅专栏

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%matrix LU Factorize
%authour: j.cai
%mail: jcai@mail.com
%time: 2016-09-29
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [L,U,P]= LU(A)
    disp("LU factorization");
    L=zeros(size(A));
    U=zeros(size(A));
    P=zeros(size(A));
    if size(A,1)~=size(A,2)
        disp("不是方阵");
        return 
    end
    n=size(A,1);
    seqence=linspace(1,n,n);
    A=[A,seqence']; %加上行号 
    for c=1:n %部分主元法 找主元
        [~,pivotnum]=max(abs(A(c:end,c)));
        pivotnum=pivotnum+c-1;
        if pivotnum ~= c %交换行 
            row_temp=A(c,:);
            A(c,:)=A(pivotnum,:);
            A(pivotnum,:)=row_temp;
        end
        for r=(c+1):n %消元 
            if abs(A(c,c))