数值线性代数.docxVIP

1.求下三角阵的逆矩阵的详细算法。 ［解］设下三角矩阵L的逆矩阵为T 我们可以使用待定法，求出矩阵丁的各列向量。为此我们将T按列分块如下: 注意到 我们只需运用算法1?1?1,逐一求解方程 ￡7］二驾， 攵=12???，想 便可求得丁=尸 ［注意］考虑到内存空间的节省，我们可以置结果矩阵T的初始状态为单位矩阵。 这样，我们便得到如下具体的算法： 算法（求解下三角矩阵L的逆矩阵T,前代法） 预置步置T = I for J = 1 ： ? for 止=1:珂 一 1 T也顶=欢虹扑 T优 +1:冬 j） = T优 +1:冬 j） — T{k, + l end =Tg）f end 3 .证明：如果玲=1忒是一个Gauss变换，则牛 E感也是一个Gauss 变换。 ?p ［解］按Gauss变换矩阵的定义，易知矩阵‘ +匕乳是Gauss变换。下面我们只需 证明它是Gauss变换。1一就 的逆矩阵。事实上 （/ t感）5 或）= 注意到A =帼…，口，上旧…志J，则显然有，怎二0?从而有￡? = E感 4 .确定一个Gauss变换L,使 2 2 L 3 — 7 _4_ _8_ ［解］比较比较向量松34）和（2工驴可以发现Gauss变换L应具有功能:使向量仅*）的 第二行加上第一行的2倍;使向量旧34尸的第三行加上第一行的2倍。于是Gauss 变换如下 -1 - Z= 2 1 _2 0 1_ 5 .证明：如果心皿有三角分解，并且是非奇异的，那么定理1?1?2中的L和 U都是唯一的。 ［证明］设妇西=站 ,其中孔扇都是单位下三角阵，叫,％ 都是上三角 阵。因为a』E奇异的，于是￡;￡1 = %叫一1 注意到，单位下三角阵的逆仍是单位下三角阵，两个单位下三角阵的乘积仍是单 位下三角阵；上三角阵的逆仍是上三角阵，两个上三角阵的乘积仍是上三角阵。 因此，上述等将是一个单位下三角阵与一个上三角阵相等，故此，它们都必是单 位矩阵。即宴1 =叩？ = ￡ 从而 l1 = l2ju1 = u2. 即入的LU分解是唯一的。 7 .设A对称且心°，并假定经过一步Gauss消去之后，A具有如下形式 ^ii _o 证明瓦仍是对称阵 ［证明］根据Gauss变换的属性，显然做矩阵A的LU分解的第一步中的Gauss变换为 1 r 1 o- 0… 1 an 其中伐I , 其中伐I ,将A分块为 由入的对称性 鸵对称性则是显而易见的。 8.设网=［。山曾!是严格对角占优阵，即A满足 KI ＞习翎|’ k= 12…m 又设经过一步Gauss消去后，A具有如下形式 白11 af _o A_ 试证：矩阵禹仍是严格对角占优阵。由此推断：对于对称的严格对角占优矩阵 来说，用Gauss消去法和列主元Gauss消去法可得得同样的结果。 ［证明］依上题的分析过程易知，题中的 A =［瓦］其中：瓦=.-组色fJ = 2,…m. 的1 于是瓦主对角线上的元素满足 TOC \o 1-5 \h \z r 22 | 22 ■■ r 22 | I ? 7 ? =11 险 1 (1) 非主对角线上的元素满足 力得二支饱L1KI+档刘知■ L JU2 JU2 队 11 Z PH JU2 知 i JUz .^2 | 稣优=LZ…齐 由于A是严格对角占优的，即 敏 故 22 I知 aH —同1 V JU2 3 M 221知1却1J -屁. JU2 从而 加 PE 、 (2) 综合(1)和(2)得 |祠 立国亦 ，=乙… * 即，矩阵瓦仍是严格对角占优阵。 9.设AeE有三角分解。指出当把Gauss消去法应用于心3+D矩阵［司时， 怎样才能不必存储L而解出Ax=b ?需要多少次乘法运算？ ［解］用Gauss消去法作A的LU分解，实际上就是对系数矩阵入作了一组初等行变 换，将其化 为上三角矩阵U。而这一组的初等行变换对应的变换矩阵就是a：即 * = U 如果把这一组初等行变换施加于方程右端向量b上，即有 这就是说，方程组欧=8和也=￡一与是同解方程。而后者是上三角形方程组， 可运用本章算法1?1?2求解。这样我们就不必存储L，通求解方程组gSb，来 求解原方程组Ax = b。算法如下： (1)用初等变换化皿切知叫厂划, jor 止=1:理一 1 HW+1:珞舫=H（k+l:您先）（妃先） +1:理，止 +1 :也）=A［ +1 :也,止 +1:理）一 A（ +1:理,k）A（k,止 +1 :扪 8（用+1 :打）=8（用+ 1:冷-+ 1:华先）8侬） end （2）利用回代法求解方程组gE。 该算法所需要的加、减、乘、除运算次数为 M-l s , % Q 旗（弘+ 2枝）+炉=_疽+2按—―理 jui 3 3 也11 a 10. A是正定阵，如果对A执行Gauss消去一步产生一个形式为L ° 土 的矩