数值分析矩阵的三角分解课件.pptVIP

数值分析06矩阵的三角分解;数值分析;Gauss变换阵的性质：;数值分析;数值分析;;Gauss变换阵的作用：;数值分析;;数值分析;;数值分析;数值分析;数值分析; ;数值分析;数值分析;进行到第k步消元时;数值分析;数值分析;其中;;;function A=lud(A) %功能：对方阵A作三角分解A=LU,其中， % L为单位下三角阵，U为上三角阵， %输入：方阵A。 %输出：紧凑存储A=[L\U]. %注意：当A的主元=0时退出Matlab.;;数值分析;数值分析;;;数值分析;;;;3. Cholesky分解(平方根分解);(k);再用L的第k行乘 LT的第j列( j=k+1 ,k+2,……,n )有;Cholesky分解公式;数值分析; ;;数值分析;;;数值分析;;MATLAB实现 [l,u,p]=lu(A);Lupd.m %功能：对方阵A作列主元三角分解PA=LU,其中， % L为单位下三角阵，U为上三角阵，排列阵P % 用向量p表示。 %输入：方阵A。 %输出：紧凑存储LU=[L\U]，以及p。 %注意：当A奇异时退出Matlab. function [LU,p]=lupd(A) %初始化 n=length(A); p=1:n; LU=A;;%判断矩阵的奇异性 if s==0 quit; end %行交换 if ik~=k m=p(k); p(k)=p(ik); p(ik)=m; lk=LU(k,:); LU(k,:)=LU(ik,:); LU(ik,:)=lk; end;;Lupqd.m %功能：对方阵A作全主元三角分解PAQT=LU,其中， % L为单位下三角阵，U为上三角阵，排列阵P % 和Q分别用向量p,q表示。 %输入：方阵A。 %输出：紧凑存储LU=[L\U]，以及p和q。 %注意：当A奇异时退出Matlab. function [LU,p,q]=lupqd(A) %初始化 n=length(A); p=1:n; q=p; LU=A;;%分解过程 for k=1:n %搜索全主元(ik,jk) [xk,I]=max(abs(LU(k:n,k:n))); %列最大值及所在行 [s,j]=max(xk);ik=I(j)+k-1;jk=j+k-1; %判断矩阵的奇异性 if s==0 quit; end %行交换和列交换 if ik~=k m=p(k); p(k)=p(ik); p(ik)=m; lk=LU(k,:); LU(k,:)=LU(ik,:); LU(ik,:)=lk; end;if jk~=k m=q(k); q(k)=q(jk); q(jk)=m; ck=LU(:,k); LU(:,k)=LU(:,jk); LU(:,jk)=ck; end %用消元法计算LU=[L\U] if k==n break; end LU(k+1:n,k)=LU(k+1:n,k)/LU(k,k); LU(k+1:n,k+1:n)=LU(k+1:n,k+1:n)-LU(k+1:n,k)*… LU(k,k+1:n); end ;;二版习题