矩陣LU分解程序實現（Matlab）

n=4;%確定需要LU分解的矩陣維數
%A=zeros(n,n);
L=eye(n,n);P=eye(n,n);U=zeros(n,n);%初始化矩陣
tempU=zeros(1,n);tempP=zeros(1,n);%初始化中間變量矩陣
A=[1 2 -3 4;4 8 12 -8;2 3 2 1;-3 -1 1 -4];%需要LU分解矩陣賦值
for p=1:n                %將A矩陣賦值給U
    for q=1:n
        U(p,q)=A(p,q);
    end
end
jt=1;kt=0;
for i=1:n-1           
    jt=jt+1;
    kt=kt+1;
    ii=U(i,i);    
    if ii==0                    %主元為零，進行行變換
        for m=i:n
            if U(m,i)~=0
                tempU=U(i,:);
                U(i,:)=U(m,:);
                U(m,:)=tempU;
                ii=U(i,i);
                %%
                tempP=P(i,:);  %行變換結果存儲在P中
                P(i,:)=P(m,:);
                P(m,:)=tempP;
                break;
            end
        end
        %disp(ii);
    end
    disp(ii);
    for j=jt:n                    %%兩重循環，完成高斯消元
          perj=U(j,i)/ii;
          L(j,i)=perj;
          for k=kt:n
                U(j,k)=U(j,k)-perj*U(i,k);
          end
     end
 end

savefile='LUdapart';
save(savefile)