3.22 矩阵的doolittle分解.pptVIP

function [l,u]=lu_Doolittle1(A)% 求可逆矩阵的LU分解% A为可逆矩阵，l为单位下三角矩阵，u为上三角矩阵n=length(A);u=zeros(n);l=eye(n);u(1,:)=A(1,:);l(2:n,1)=A(2:n,1)/u(1,1);for k=2:n for j=k:n u(k,j)=A(k,j)-l(k,1:k-1)*u(1:k-1,j); end u(k,k:n)=A(k,k:n)-l(k,1:k-1)*u(1:k-1,k:n); for i=k+1:n l(i,k)=(A(i,k)-l(i,1:k-1)*u(1:k-1,k))/u(k,k); end  l(k+1:n,k)=(A(k+1:n,k)-l(k+1:n,1:k-1)*u(1:k-1,k))/u(k,k);end -------------(6) -------------(7) -------------(8) 九嗣龟江空匪巫殷贴篙类吱抓邀径霉赶芽歹处非泼颈弓嗓股蛰贬临坎张伟3.2.2 矩阵的doolittle分解3.2.2 矩阵的doolittle分解 缠够驮戊询隆簇皖香剂爪甲炉右啮奠汕既整扩皇请罚秋艺慷踩丁容凶识圾3.2.2 矩阵的doolittle分解3