求解线性方程组LU分解法的Matlab程序.docx

《求解线性方程组LU分解法的Matlab程序.docx》由会员分享，可在线阅读，更多相关《求解线性方程组LU分解法的Matlab程序.docx（2页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

精品文档，仅供学习与交流，如有侵权请联系网站删除function L,U,flag=LU_decom(A)n,m=size(A);if n=m error(The rows and columns of matrix A must be equal!); return;endL=eye(n);U=zeros(n);flag=OK;for k=1:n for j=k:n z=0; for q=1:k-1 z=z+L(k,q)*U(q,j); end U(k,j)=A(k,j)-z; end if abs(U(k,k)eps flag=failure; return; end for i=k+1:n z=0; for q=1:k-1 z=z+L(i,q)*U(q,k); end L(i,k)=(A(i,k)-z)/U(k,k); endend【精品文档】第 2 页