qq_45934720的博客

原创 LMBP算法

function[x]=Dool(A,b) n=length(A); A(2:n,1)=A(2:n,1)/A(1,1); for t=2:n-1 A(t,t:n)=A(t,t:n)-A(t,1:t-1)*A(1:t-1,t:n); A(t+1:n,t)=(A(t+1:n,t)-A(t+1:n,1:t-1)*A(1:t-1,t))/A(t,t); end A(n,n)=A(n,n)-A(n,1:n-1)*A(1:n-1,n); L=tril(A,-1)+eye(n); U=triu(A); f...