weiwang_math的博客

今天数值分析上机课， LU 分解 Doolittle 分解，写个小程序， 由于MATLAB对向量和矩阵乘法的优势，可简化其中一些步骤 % LU 分解 Doolittle 分解 A = [1,2,1,-2;2,5,3,-2;-2,-2,3,5;1,3,2,3]; b = [2;8;4;9];[m,n]=size(A); u = zeros(n); l=eye(n); % LU分解 u(

% Jacobi方法和Gauss-Seidel方法 clear all A = [8,-1,1;2,10,-1; 1 1 -5]; b = [1;4;3]; n = length(b); itermax = 100; iter = 0; tol = 1e-4 x = zeros(n,1); while iter