数值分析矩阵分解法.pptxVIP

会计学;;;第3页/共69页; ;第5页/共69页;第6页/共69页;进行到第k步消元时;第8页/共69页;第9页/共69页;其中;;;function A=lud(A) %功能：对方阵A作三角分解A=LU,其中， % L为单位下三角阵，U为上三角阵， %输入：方阵A。 %输出：紧凑存储A=[L\U]. %注意：当A的主元=0时退出Matlab.;function A=lud(A) %功能：对方阵A作三角分解A=LU,其中， % L为单位下三角阵，U为上三角阵， %输入：方阵A。 %输出：紧凑存储A=[L\U].;function A=lud(A) %功能：对方阵A作三角分解A=LU,其中， % L为单位下三角阵，U为上三角阵， %输入：方阵A。 %输出：紧凑存储A=[L\U]. [n,n]=size(A); % 确定A的维数 for k=1:n-1 if A(k,k) ==0 break; end A(k+1:n,k) =A(k+1:n,k)/ A(k,k); A(k+1:n ,k+1:n)= A(k+1:n ,k+1:n)- A(k+1:n ,k) *A(k,k+1:n); e