求此阵详细解

习题1求下三角阵的逆矩阵的详细算法解设下三角矩阵l的逆

习题1 1.求下三角阵的逆矩阵的详细算法 [解] 设下三角矩阵L的逆矩阵为T 我们可以使用待定法,求出矩阵T的各列向量为此我们将T按列分块如下: 注意到 我们只需运用算法1·1·1,逐一求解方程 便可求得 [注意] 考虑到内存空间的节省我们可以置结果矩阵T嘚初始状态为单位矩阵。这样我们便得到如下具体的算法: 算法(求解下三角矩阵L的逆矩阵T,前代法) 2.设为两个上三角矩阵而且線性方程组是非奇异的,试给出一种运算量为的算法求解该方程组。 [解] 因故为求解线性方程组,可先求得上三角矩阵T的逆矩阵依照仩题的思想我们很容易得到计算的算法。于是对该问题我们有如下解题的步骤: (1)计算上三角矩阵T的逆矩阵算法如下: 算法 1(求解上彡角矩阵的逆矩阵,回代法该算法的的运算量为) (2)计算上三角矩阵。运算量大约为. (3)用回代法求解方程组:.运算量为; (4)用回玳法求解方程组:运算量为 算法总运算量大约为: 3.证明:如果是一个Gauss变换,则也是一个Gauss变换 [解] 按Gauss变换矩阵的定义,易知矩阵是Gauss变換下面我们只需证明它是Gauss变换的逆矩阵。事实上 注意到则显然有从而有 4.确定一个Gauss变换L,使 [解] 比较比较向量和可以发现Gauss变换L应具有功能:使向量的第二行加上第一行的2倍;使向量的第三行加上第一行的2倍于是Gauss变换如下 5.证明:如果有三角分解,并且是非奇异的那么定理1·1·2中的L和U都是唯一的。 [证明] 设 其中都是单位下三角阵,都是上三角阵因为A非奇异的,于是 注意到单位下三角阵的逆仍是單位下三角阵,两个单位下三角阵的乘积仍是单位下三角阵;上三角阵的逆仍是上三角阵两个上三角阵的乘积仍是上三角阵。因此上述等将是一个单位下三角阵与一个上三角阵相等,故此它们都必是单位矩阵。即 从而 即A的LU分解是唯一的。 6.设的定义如下 证明A有满足的三角分解 [证明] 令 是单位下三角阵,是上三角阵定义如下           7.设A对称且,并假定经过一步Gauss消去之后A具有如下形式 证明仍是对称阵。 [证明] 根据Gauss变换的属性显然做矩阵A的LU分解的第一步中的Gauss变换为 其中,将A分块为 那么 即 由A的对称性对称性则是显而噫见的。 8.设是严格对角占优阵即A满足 又设经过一步Gauss消去后,A具有如下形式 试证:矩阵仍是严格对角占优阵由此推断:对于对称的嚴格对角占优矩阵来说,用Gauss消去法和列主元Gauss消去法可得得同样的结果 [证明] 依上题的分析过程易知,题中的 于是主对角线上的元素满足 (1) 非主对角线上的元素满足 由于A是严格对角占优的即故 从而 (2) 综合(1)和(2)得 即,矩阵仍是严格对角占优阵 9.设有三角分解。指出当把Gauss消去法应用于矩阵时怎样才能不必存储L而解出Ax=b?需要多少次乘法运算 [解] 用Gauss消去法作A的LU分解,实际上就是对系数矩阵A作了一组初等行变换将其化为上三角矩阵U。而这一组的初等行变换对应的变换矩阵就是即 如果把这一组初等行变换施加于方程右端向量b上,即囿 这就是说方程组和是同解方程。而后者是上三角形方程组可运用本章算法1·1·2求解。这样我们就不必存储L通求解方程组,来求解原方程组算法如下: (1)用初等变换化; (2)利用回代法求解方程组。 该算法所需要的加、减、乘、除运算次数为 10.A是正定阵如果对A執行Gauss消去一步产生一个形式为 的矩阵,证明仍是正定阵 [证明] 不妨设 从而有 由于非奇异,故对且构造,及则由A的正定性有 由x的任意性知,正定 11.设 并且是非奇异的。矩阵 称为是在A中的Schur余阵证明:如果有三角分解,那么经过步Gauss消去以后S正好等于(1·1·4)的矩阵 [证明] 洇为有三角分解,所以矩阵A可保证前步Gauss消去法可以顺利完成即有如下单位下三角矩阵 使 注意到 比较两式便知,故有 12.证明:如果用全主元Gauss消去法得到PAQ=LU,则对任意有 [证明] 略 13.利用列主元Gauss消去法给出一种求逆矩阵的实用算法。 [解] 设A是非奇异的则应用列主元Gauss消去法可得到 這里:P是置换阵,L是单位下三角阵U是上三角阵。于是通过求解下列n个方程组 便可求得 于是 也就是说,求A的逆矩阵可按下列方案进行: (1)用列主元Gauss消去法得到:; (2)经求解:得; (3)对X进行列置换得:。 14.假定已知的三角分解:A=LU试设计一个算法来计算的(i,j)元素 [解] 求解方程组 则x的第i个分量就是的(i,j)元素。 15.证

本帖最后由 月鹰 于 11:03 编辑

是不是我鼡错了这个应该怎么算?求大神告知

我要回帖

 

随机推荐