线性代数的约束条件绝对值如何处理？

来源：蜘蛛抓取(WebSpider) 时间：2022-09-21 05:11 标签： 线性代数极大线性无关组怎么求

《数值分析：第六章 线性代数方程组的直接法》由会员分享，可在线阅读，更多相关《数值分析：第六章 线性代数方程组的直接法（54页珍藏版）》请在人人文库网上搜索。

1、第六章,线性方程组 的直接解法,问题驱动：投入产出分析,投入产出分析是20世纪30年代由美国经济学家首先提出的， 它是研究整个经济系统各部门之间“投入”与“产出”关系的线性 模型，一般称为投入产出模型。国民经济各个部门之间存在着 相互依存的关系，每个部门在运转中将其它部门的成品或半成 品经过加工（称为投入）变为自己的产品（称为产出），如何 根据各部门之间的投入-产出关系，确定各部门的产出水平，以 满足社会的需求，是投入产出综合平衡模型研究的问题，试讨论 如下简化问题。,设国民经济仅由农业、制造业和服务业三个部门构成，已知某年它们之间的投入和产出关系、外部需求、初始投 入等如表6.1.1所示（数

2、字表示产值，单位为亿元）。,表6.1.1 国民经济各个部门间的关系,表中第一行数字表示农业总产出为100亿元，其中15亿元农产品用于农业生产本身，20亿元用于制造业，30亿元用于服务业，剩下的35亿元农产品用于满足外部需求。类似地可以解释第二、三行数字。第一列数字中，15亿元如前所述，30亿元是制造业对农业的投入，20亿元是服务业对农业的投入，35亿元的初始投入包括工资、税收、进口等，总投入100亿元和总产出相等。假定每个部门的产出和投入是成正比的，由表6.1.1能够确定这三个部门的投入产出表，如表6.1.2所示。,表6.1.2 投入产出表,表中的第一行，第二列的数字表示生产1个单位产值的制造

3、业产品需要投入0.10个单位的产值的农产品，同样第三行、第一列的数字表示，生产1个单位产值的农产品需要0.20个单位的服务业产值。表6.1.2的数字称为投入系数和消耗系数，如果技术水平没有变化，可以假设投入系数是常数。已知投入系数如表2.1.2所示，若今年对农业、制造业和服务业的外部需求分别为50、150、100亿元，试计算三个部门的总产出分别为多少？,若共有n个部门，记一定时期内第i个部门的总产出为 xi, 其中对第 j个部门的投入为xij ，满足的外部需求为 di ，则,(6.1.1),记第j个部门的单位产出需要第i个部门的投入为aij ，在每个部门的产出与投入成正比的假定下，有,(6.1

4、.2),投入系数即为aij，将(6.1.2)式代入(6.1.1)）式得方程组,用矩阵表示为,因此投入产出模型最终可归结为求解线性方程组的问题，下面 介绍求解线性方程组数值方法。,AX = b,（3.1）,线性方程组数值解法的分类,直接法（适用于中等规模的n阶线性方程组） Gauss消去法及其变形 矩阵的三角分解法,迭代法（适用于高阶线性方程组） Jacobi迭代法 Gauss-Seidel迭代法 逐次超松弛法 共轭斜量法,1 高斯消去法,1三角形方程组的解法-回代法,（3.2）,（3.3）,2顺序高斯消去法,基本思想：通过消元将上述方程组化为上三角形 方程组，再利用回代法进行求解。,记,消元公

5、式,回代公式,对k=1,2,n-1, 若,则可得,顺序Gauss消去法可执行的前提,定理 1 给定线性方程组 ，如果n阶方阵 的所有顺序主子式都不为零，即 则按顺序Gauss消去法所形成的各主元素 均不为零，从而顺序Gauss 消去法可顺利执行。,注：当线性方程组的系数矩阵为对称正定或严格对角占优阵时，按顺序Gauss消去法计算是稳定的。,3、列主元Gauss消去法计算步骤：,1、输入矩阵阶数n，增广矩阵 A(n,n+1);,2、对于,(1) 按列选主元：选取 l 使,(2) 如果 ，交换 A(n,n+1) 的第k行与第l 行元素,(3) 消元计算 ：,3、回代计算,4无回代过程的主元消去法(

6、Gauss-Jordan消去法),第一步：选主元，在第一列中选绝对值最大的元素，设第k行为主元行，将主元行换至第一行，将第一个方程中x1的系数变为1，并从其余n 1个方程中消去x1。,第二步：在第二列后n 1个元素中选主元，将第二个方程中 x2的系数变为1，并从其它n 1个方程中消去x2。,第k步：在第k列后n k个元素中选主元，换行，将第k个方程 xk的系数变为1，从其它n - 1个方程中消去变量xk,消元公式为：,对k = 1, 2, , 按上述步骤进行到第n步后，方程组变为：,即为所求的解.,注：无回代的Gauss消元法实际上就是将 方程组的系数矩阵化为行最简形矩阵。,5无回代消去法的应

7、用,（1）解线性方程组系,设要解的线性方程组系为：,AX = b1, AX = b2, AX = bm,上述方程组系可以写为,AX = B = (b1, , bm),因此X = A-1B 即为线性方程组系的解。,在计算机上只需要增加几组右端常数项的存贮单元， 其结构和解一个方程组时一样。,行,系数,右端,（2）求逆矩阵,设A = (aij)nn是非奇矩阵，A 0，且令,由于 AA-1 = AX = I,因此，求A-1的问题相当于解下列线性方程组,相当于（1）中m = n, B = I 的情形。,（3）求行列式的值,用高斯消去法将 A化成上三角形,解：把系数矩阵、单位矩阵和右端项组成增广矩阵，对

8、增广矩 阵实行Gauss-Jordan消元过程,。,2 解三对角方程组的追赶法,3 矩阵的三角分解法, 高斯消元法的矩阵形式,每一步消去过程相当于左乘初等下三角矩阵Lk,记,其中，,记,其中，,一般的，对,A 的 LU 分解 ( LU factorization ),记,则,其中，,于是, 矩阵三角分解的定义,定义1,设A为n n实矩阵，如果存在下三角矩阵L与 上三角矩阵U，使得 A=LU，则称之为矩阵 A的三角分解；若存在单位下三角矩阵L，对角 矩阵D及单位上三角矩阵R，使得A=LDR,则称之 为矩阵A的LDR分解。,注： （1） L 为单位下三角阵而 U 为一般上三角阵的分解称为Dooli

我要回帖

更多关于 线性代数极大线性无关组怎么求 的文章

 

随机推荐