[线性代数实践教师班第三讲.pptVIP

线性代数实践（第三讲） 第7章 矩阵运算法解方程 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 7.1 矩阵运算的规则 在MATLAB入门中已讲过的，不再重复。 由于其乘法不符合交换律，有些公式不能乱用； 单列向量与单行向量的左右两种乘法要加区别，而且往往有特别的用途。 例如向量长度（范数）的计算； 例如二维坐标网格的生成； X=ones(21,1)*[-10:10],Y= [-10:10]’*ones(1,21) 矩阵的乘幂An, eA和(I-A) -1的级数展开，都要求A是方阵。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 矩阵乘法不满足交换律 有许多我们习惯的公式，其中隐含地包含了交换律，这些公式在矩阵运算中也不能直接使用。比如： 正确的做法是展开时不交换次序 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 平面上网格坐标系的产生 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 用列矩阵乘行矩阵生成网格坐标 这两个矩阵都是21行21列的，都有441个元素，如何快捷地输入呢？这时可以用到列乘行的乘法运算。可用下面的语句： h10:10; lhlength(h) % 输入均分行向量 % 用全么列乘均分行生成X Xones(lh,1)*h % 用均分列乘全么行生成Y Yh‘*ones(1,lh) Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 7.2 初等变换乘子矩阵的生成 行交换E1gen(n,i,j)：使n行矩阵中的第i,j两行交换 function E=E1gen(n,i,j) n=size(A); E=eye(n); E(i,i)=0; E(j,j)=0; E(i,j)=1; E(j,i)=1; 乘子矩阵E2gen (n,i,k),使n行矩阵中的第i行乘以k function E=E2gen(n,i,k) n=size(A); E=eye(n); E(i,i)=k; E3gen(n,i,j,c)使n行矩阵中的第i行乘以k加到第j行上 function E=E3gen(n,i,j,k) n=size(A); E=eye(n); E(j,i)=k; Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 初等变换乘子矩阵示例 E=E1gen(8,4,6) E2=E2gen(8,4,6) E3=E3gen(8,4,6,5) 例如E3=E3gen(3,1,3,4) Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 例7.2.4 求消元所需的乘子矩阵 要消去下列矩阵的A(2,1),求乘子矩阵E3 在第二行加以第一行乘A(2,1)/A(1,1)3， 故令B E3gen（A,1,2,3） Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 行阶梯生成等价于矩阵左乘 因此，整个行阶梯形式U的生成过程，可以看作把原矩阵左乘以一系列的初等变换矩阵E1和E3。把这些初等矩阵的连乘积写成Ex，设其逆为L: 从而有 L*UA