克劳特（Crout）（LU）分解法求解线性方程组的matlab实现(共3页).doc

精选优质文档-倾情为你奉上1、克劳特（Crout）（LU）分解法求解线性方程组function x,L,U=Crout(A,b)%Crout分解法求解线性方程组%系数矩阵：AN=size(A);n=N(1);L=zeros(n,n); %下三角矩阵U=eye(n,n); %上三角矩阵L(1:n,1)=A(1:n,1); %L的第一列U(1,1:n)=A(1,1:n)/L(1,1); %U的第一行?for k=2:n for i=k:n L(i,k)=A(i,k)-L(i,1:(k-1)*U(1:(k-1),k); %L的第k列 end for j=(k+1):n U(k,j)=(A(k,j)-L(k,1:(k-1)*U(1:(k-1),j)/(L(k,k); %U的第k行 endend%y=inv(L)*b;%x=inv(U)*y;y=Sol