這個是三角分解法的一個標準程式function[ L,U] =lu_decompose(A,b)A=[-4 1 1 1;1 -4 1 1;1 1 -4 1;1 1 1 -4];b=[1;1;1;1];n=length(b);L=eye(n);U=zeros(n,n);%對A進行LU分解for i=1:n U(1,i)=A(1,i); if i==1 L(i,1)=1; else L(i,1)=A(i,1)/U(1,1); endendfor i=2:n for j=i:n sum=0; for k=1:i-1 sum=sum+L(i,k)*U(k,j); end U(i,j)=A(i,j)-sum; if j~=n sum=0; for k=1:i-1 sum=sum+L(j+1,k)*U(k,i); end L(j+1,i)=(A(j+1,i)-sum)/U(i,i); end endenddisp(U)
這個是三角分解法的一個標準程式function[ L,U] =lu_decompose(A,b)A=[-4 1 1 1;1 -4 1 1;1 1 -4 1;1 1 1 -4];b=[1;1;1;1];n=length(b);L=eye(n);U=zeros(n,n);%對A進行LU分解for i=1:n U(1,i)=A(1,i); if i==1 L(i,1)=1; else L(i,1)=A(i,1)/U(1,1); endendfor i=2:n for j=i:n sum=0; for k=1:i-1 sum=sum+L(i,k)*U(k,j); end U(i,j)=A(i,j)-sum; if j~=n sum=0; for k=1:i-1 sum=sum+L(j+1,k)*U(k,i); end L(j+1,i)=(A(j+1,i)-sum)/U(i,i); end endenddisp(U)