摘要:
%直接三角分解法function my_LU(a, b) n = length(a); l = zeros(n, n);%初始化 u = zeros(n, n); for i=1:n l(i,i) = 1; end u(1,1:n) = a(1,1:n); l(2:n, 1) = a(2:n, 1) ./ u(1,1); for r=2:n for i=r:n u(r, i) = a(r, i) - sum(l(r,1:r-1) .* (u(1:r-1,i))'); end for i=r+1:n if... 阅读全文