摘要： 求解Ax=d的解x1. LU分解% file: myLU.mfunction [L,U]=myLU(A)%实现对矩阵A的LU分解，L为下三角矩阵A[n,n]=size(A);L=zeros(n,n);U=zeros(n,n);for i=1:n L(i,i)=1;endfor k=1:n for j=k:n U(k,j)=A(k,j)-sum(L(k,1:k-1).*U(1:k-1,j)'); end for i=k+1:n L(i,k)=(A(i,k)-sum(L(i,1:k-1).*U(1:k-1,k)'))/U(k,k); end... 阅读全文