资源描述：
《数值分析课件06》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、第五节高斯变换阵与矩阵的三角分解一、Gauss变换阵定义Gauss变换阵为数值分析数值分析数值分析数值分析Gauss变换阵的性质：数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析Gauss变换阵的作用：数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析顺序高斯消元的基本思想：将矩阵A的下三角部分消为零，即二、矩阵的三角分解1.顺序高斯消元与LU分解的等价性数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析进行到第k步消元时数值分析数值分析数值分析数值分析数值分析数值分析其中
2、数值分析数值分析数值分析数值分析消元过程全部完成后，原来的二维数组中存放的元素实际上是一个新的矩阵，记为数值分析数值分析functionA=lud(A)%功能：对方阵A作三角分解A=LU,其中，%L为单位下三角阵，U为上三角阵，%输入：方阵A。%输出：紧凑存储A=[LU].%注意：当A的主元=0时退出Matlab.LU分解的MATLAB程序[n,n]=size(A);%确定A的维数fork=1:n-1fori=k+1:nifA(k,k)==0quit;endA(i,k)=A(i,k)/A(k,k);A(i,k+1:n)=A(i,k+1:n)-A(i,k)*A(k
3、,k+1:n);endend数值分析数值分析2.矩阵三角分解的基本定理数值分析数值分析数值分析数值分析数值分析数值分析2.矩阵三角分解的基本定理数值分析数值分析主要结论数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析3.Cholesky分解(平方根分解)(1)（用直接三角分解法）数值分析数值分析(k)若已求出了L的前k-1列元素，则为求第k列元素，先用L的第k行乘LT的第k列数值分析数值分析再用L的第k行乘LT的第j列(j=k+1,k+2,……,n)有数值分析数值分析Cholesky分解公式数值分析数值分析数值分析数值分析4.列主元
4、LU分解数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析数值分析MATLAB实现[l,u,p]=lu(A)数值分析数值分析Lupd.m%功能：对方阵A作列主元三角分解PA=LU,其中，%L为单位下三角阵，U为上三角阵，排列阵P%用向量p表示。%输入：方阵A。%输出：紧凑存储LU=[LU]，以及p。%注意：当A奇异时退出Matlab.function[LU,p]=lupd(A)%初始化n=length(A);p=1:n;LU=A;%分解过程fork=1:n%搜索列主元ik[s,i]=max(abs(LU(k:
5、n,k)));ik=i+k-1;数值分析数值分析%判断矩阵的奇异性ifs==0quit;end%行交换ifik~=km=p(k);p(k)=p(ik);p(ik)=m;lk=LU(k,:);LU(k,:)=LU(ik,:);LU(ik,:)=lk;end%用消元法计算LU=[LU]ifk==nbreak;endLU(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数值分析数值分析5.全主元LU分解数值分析数值分析Lupqd.m%功能：
6、对方阵A作全主元三角分解PAQT=LU,其中，%L为单位下三角阵，U为上三角阵，排列阵P%和Q分别用向量p,q表示。%输入：方阵A。%输出：紧凑存储LU=[LU]，以及p和q。%注意：当A奇异时退出Matlab.function[LU,p,q]=lupqd(A)%初始化n=length(A);p=1:n;q=p;LU=A;数值分析数值分析%分解过程fork=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;%判断矩阵的奇异性ifs==0
7、quit;end%行交换和列交换ifik~=km=p(k);p(k)=p(ik);p(ik)=m;lk=LU(k,:);LU(k,:)=LU(ik,:);LU(ik,:)=lk;end数值分析数值分析ifjk~=km=q(k);q(k)=q(jk);q(jk)=m;ck=LU(:,k);LU(:,k)=LU(:,jk);LU(:,jk)=ck;end%用消元法计算LU=[LU]ifk==nbreak;endLU(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
8、(k,k+
