矩阵的分解精选PPT.ppt

《矩阵的分解精选PPT.ppt》由会员分享，可在线阅读，更多相关《矩阵的分解精选PPT.ppt（43页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、关于矩阵的分解第1页，讲稿共43张，创作于星期二定理3.12 L是单位下三角矩阵U一个上三角矩阵 Gauss消元法的消元过程实际上是对线性代数方程组进行一消元法的消元过程实际上是对线性代数方程组进行一系列系列初等行变换初等行变换的过程。由线性代数知识知，线性代数方程组的过程。由线性代数知识知，线性代数方程组的初等变换相当于对其增广矩阵实行的初等变换相当于对其增广矩阵实行初等行变换初等行变换，也即相当，也即相当于增广矩阵于增广矩阵左边乘以一个初等矩阵左边乘以一个初等矩阵。第2页，讲稿共43张，创作于星期二第3页，讲稿共43张，创作于星期二也可以直接用比较法导出矩阵A的LU分解的计算公式。上式可记

2、为比较第1行比较第r行第4页，讲稿共43张，创作于星期二同样，由比较第r列第5页，讲稿共43张，创作于星期二综合以上分析，有因此可以推导出U的第一行L的第一列-(1)-(2)第6页，讲稿共43张，创作于星期二思考U的第r行L的第r列-(3)-(4)称上述(1)(4)式所表示的分解过程为矩阵A的Doolittle分解第7页，讲稿共43张，创作于星期二function l,u=lu_Doolittle1(A)%求可逆矩阵的LU分解%A为可逆矩阵，l为单位下三角矩阵，u为上三角矩阵n=length(A);u=zeros(n);l=eye(n);u(1,:)=A(1,:);l(2:n,1)=A(2:n

3、,1)/u(1,1);for k=2:n for j=k:n u(k,j)=A(k,j)-l(k,1:k-1)*u(1:k-1,j);end u(k,k:n)=A(k,k:n)-l(k,1:k-1)*u(1:k-1,k:n);for i=k+1:n l(i,k)=(A(i,k)-l(i,1:k-1)*u(1:k-1,k)/u(k,k);end l(k+1:n,k)=(A(k+1:n,k)-l(k+1:n,1:k-1)*u(1:k-1,k)/u(k,k);end第8页，讲稿共43张，创作于星期二对于线性方程组系数矩阵非奇异，经过Doolittle分解后线性方程组可化为下面两个三角形方程组第9页，

4、讲稿共43张，创作于星期二第10页，讲稿共43张，创作于星期二上述解线性方程组的方法称为直接三角分解法的 Doolittle分解例例3.2.1 用Doolittle分解求解方程组解解下面再用Doolittle分解方法求解第11页，讲稿共43张，创作于星期二第12页，讲稿共43张，创作于星期二Doolittle分解在计算机上实现是比较容易的但如果按上述流程运算仍需要较大的存储空间:第13页，讲稿共43张，创作于星期二因此可按下列方法存储数据:第14页，讲稿共43张，创作于星期二直接三角分解的Doolittle分解可以用以下过程表示:存储单元(位置)第15页，讲稿共43张，创作于星期二Doolit

5、tle分解的紧凑格式第16页，讲稿共43张，创作于星期二Doolittle分解的结果与Gauss消元法所得结果完全一样，但却避免了中间过程。第17页，讲稿共43张，创作于星期二定理3.2.3 设矩阵A非奇异，当且仅当矩阵A的所有顺序主子式全非零时，其Doolittle分解式存在，且分解是惟一的。下面给出Doolittle分解存在惟一的一个充要条件第18页，讲稿共43张，创作于星期二例例3.2.2 用紧凑格式的Doolittle分解求解方程组解解第19页，讲稿共43张，创作于星期二第20页，讲稿共43张，创作于星期二所以第21页，讲稿共43张，创作于星期二例例3.2.3 用Doolittle分解

6、求解方程组解解直接利用Doolittle分解的紧凑格式算得第22页，讲稿共43张，创作于星期二第23页，讲稿共43张，创作于星期二 列选主元列选主元Doolittle分解分解在Doolittle分解(包括紧凑格式)中，会反复用到公式仍有可能是小主元做除数为此，也要考虑在算法中加入选取列主元第24页，讲稿共43张，创作于星期二 Crout 分解分解L为下三角矩阵，U为单位上三角矩阵第25页，讲稿共43张，创作于星期二三、三、Cholesky分解与平方根法分解与平方根法 对称正定矩阵的三角分解（对称正定矩阵的三角分解（Cholesky分解）分解）第26页，讲稿共43张，创作于星期二第27页，讲稿共

7、43张，创作于星期二因此可以证明这种分解是唯一的可以证明这种分解是唯一的设存在另外的一个分解设存在另外的一个分解则则单位下三角单位下三角上三角上三角所以：所以：第28页，讲稿共43张，创作于星期二又因为：又因为：即即所以：所以：即即则：则：令：令：第29页，讲稿共43张，创作于星期二综合以上分析,则有为了方便我们记：为了方便我们记：定理定理3.2.3 (Cholesky分解)且该分解式唯一这种关于对称正定矩阵的分解称为Cholesky分解第30页，讲稿共43张，创作于星期二第31页，讲稿共43张，创作于星期二-(6)-(7)-(8)第32页，讲稿共43张，创作于星期二第33页，讲稿共43张，创

8、作于星期二 对称正定线性方程组的解法对称正定线性方程组的解法线性方程组-(10)-(11)因而线性方程组(10)可化为两个三角形方程组-(12)-(13)第34页，讲稿共43张，创作于星期二例例3.2.7用平方根法解对称正定方程组解解第35页，讲稿共43张，创作于星期二第36页，讲稿共43张，创作于星期二即所以原方程组的解为第37页，讲稿共43张，创作于星期二 平方根法的数值稳定性平方根法的数值稳定性用平方根法求解对称正定方程组时不需选取主元由可知因此平方根法是数值稳定的第38页，讲稿共43张，创作于星期二3.2.6 解三对角方程组的追赶法解三对角方程组的追赶法第39页，讲稿共43张，创作于星期二如果矩阵A满足第40页，讲稿共43张，创作于星期二以下先以gauss消元法导出三对角方程组的解法第41页，讲稿共43张，创作于星期二设经过经过n-1消元以后消元以后第42页，讲稿共43张，创作于星期二感感谢谢大大家家观观看看第43页，讲稿共43张，创作于星期二