【文件属性】:
文件名称: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