文件名称:LU分解法解线性方程组
文件大小:96KB
文件格式:DOC
更新时间:2013-07-19 09:40:44
LU
MTLAB用LU分解法解线性方程function [L,U]=do(A) [n,n]=size(A); U=zeros(n,n); L=eye(n); U(1,:)=A(1,:); L(2:n,1)=A(2:n,1)/U(1,1); for i=2:n for j=i:n U(i,j)=A(i,j)-L(i,1:i-1)*U(1:i-1,j); if j~=i L(j,i)=A(j,i)-L(j,1:i-1)*U(1:i-1,i)/U(i,i); end end end