线性代数方程组的直接法 (2)幻灯片.ppt
《线性代数方程组的直接法 (2)幻灯片.ppt》由会员分享,可在线阅读,更多相关《线性代数方程组的直接法 (2)幻灯片.ppt(54页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、线性代数方程组的直接法第1页,共54页,编辑于2022年,星期一问题驱动:投入产出分析问题驱动:投入产出分析 投入产出分析是投入产出分析是2020世纪世纪3030年代由美国经济学家首先提出的,年代由美国经济学家首先提出的,它是研究整个经济系统各部门之间它是研究整个经济系统各部门之间“投入投入”与与“产出产出”关系的线性关系的线性模型,一般称为投入产出模型。国民经济各个部门之间存在着模型,一般称为投入产出模型。国民经济各个部门之间存在着相互依存的关系,每个部门在运转中将其它部门的成品或半成相互依存的关系,每个部门在运转中将其它部门的成品或半成品经过加工(称为投入)变为自己的产品(称为产出),如何
2、品经过加工(称为投入)变为自己的产品(称为产出),如何根据各部门之间的投入根据各部门之间的投入-产出关系,确定各部门的产出水平,以产出关系,确定各部门的产出水平,以满足社会的需求,是投入产出综合平衡模型研究的问题,试讨论满足社会的需求,是投入产出综合平衡模型研究的问题,试讨论如下简化问题。如下简化问题。设国民经济仅由农业、制造业和服务业三个部门构成,已知设国民经济仅由农业、制造业和服务业三个部门构成,已知某年它们之间的投入和产出关系、外部需求、初始投某年它们之间的投入和产出关系、外部需求、初始投入等如表入等如表6.1.16.1.1所示(数字表示产值,单位为亿元)所示(数字表示产值,单位为亿元)
3、。第2页,共54页,编辑于2022年,星期一表表6.1.1 6.1.1 国民经济各个部门间的关系国民经济各个部门间的关系表中第一行数字表示农业总产出为表中第一行数字表示农业总产出为100100亿元,其中亿元,其中1515亿元农产品用亿元农产品用于农业生产本身,于农业生产本身,2020亿元用于制造业,亿元用于制造业,3030亿元用于服务业,剩下的亿元用于服务业,剩下的3535亿元农产品用于满足外部需求。类似地可以解释第二、三行数字。亿元农产品用于满足外部需求。类似地可以解释第二、三行数字。第一列数字中,第一列数字中,1515亿元如前所述,亿元如前所述,3030亿元是制造业对农业的投入,亿元是制造
4、业对农业的投入,2020亿元是服务业对农业的投入,亿元是服务业对农业的投入,3535亿元的初始投入包括工资、税收、亿元的初始投入包括工资、税收、进口等,总投入进口等,总投入100100亿元和总产出相等。假定每个部门的产出和投亿元和总产出相等。假定每个部门的产出和投入是成正比的,由表入是成正比的,由表6.1.16.1.1能够确定这三个部门的投入产出表,如能够确定这三个部门的投入产出表,如表表6.1.26.1.2所示。所示。第3页,共54页,编辑于2022年,星期一表表6.1.2 6.1.2 投入产出表投入产出表表中的第一行,第二列的数字表示生产表中的第一行,第二列的数字表示生产1 1个单位产值的
5、制造业个单位产值的制造业产品需要投入产品需要投入0.100.10个单位的产值的农产品,同样第三行、第个单位的产值的农产品,同样第三行、第一列的数字表示,生产一列的数字表示,生产1 1个单位产值的农产品需要个单位产值的农产品需要0.200.20个单位个单位的服务业产值。表的服务业产值。表6.1.26.1.2的数字称为投入系数和消耗系数,如的数字称为投入系数和消耗系数,如果技术水平没有变化,可以假设投入系数是常数。已知投入果技术水平没有变化,可以假设投入系数是常数。已知投入系数如表系数如表2.1.22.1.2所示,若今年对农业、制造业和服务业的外部所示,若今年对农业、制造业和服务业的外部需求分别为
6、需求分别为5050、150150、100100亿元,试计算三个部门的总产出分亿元,试计算三个部门的总产出分别为多少?别为多少?第4页,共54页,编辑于2022年,星期一若共有若共有n个部门,记一定时期内第个部门,记一定时期内第i个部门的总产出为个部门的总产出为 xi,其中对第其中对第 j个部门的投入为个部门的投入为xij,满足的外部需求为,满足的外部需求为 di,则,则 (6.1.1)记第记第j个部门的单位产出需要第个部门的单位产出需要第i个部门的投入为个部门的投入为aij,在每个部门的产出,在每个部门的产出与投入成正比的假定下,有与投入成正比的假定下,有(6.1.2)投入系数即为投入系数即为
7、aij,将,将(6.1.2)式代入式代入(6.1.1))式得方程组)式得方程组 用矩阵表示为用矩阵表示为或(6.1.3)因此投入产出模型最终可归结为求解线性方程组的问题,下面因此投入产出模型最终可归结为求解线性方程组的问题,下面介绍求解线性方程组数值方法。介绍求解线性方程组数值方法。第5页,共54页,编辑于2022年,星期一AXAX=b b(3.1)第6页,共54页,编辑于2022年,星期一线性方程组数值解法的分类线性方程组数值解法的分类 直接法直接法(适用于中等规模的(适用于中等规模的n n阶线性方程组)阶线性方程组)GaussGauss消去法及其变形消去法及其变形 矩阵的三角分解法矩阵的三
8、角分解法 迭代法迭代法(适用于高阶线性方程组)(适用于高阶线性方程组)JacobiJacobi迭代法迭代法 Gauss-SeidelGauss-Seidel迭代法迭代法 逐次超松弛法逐次超松弛法 共轭斜量法共轭斜量法第7页,共54页,编辑于2022年,星期一1 1 高斯消去法高斯消去法1 1三角形方程组的解法三角形方程组的解法-回代法回代法(3.2)(3.3)第8页,共54页,编辑于2022年,星期一2 2顺序高斯消去法顺序高斯消去法 基本思想:基本思想:通过消元将上述方程组化为上三角形通过消元将上述方程组化为上三角形 方程组,再利用回代法进行求解。方程组,再利用回代法进行求解。第9页,共54
9、页,编辑于2022年,星期一记记第10页,共54页,编辑于2022年,星期一消元公式消元公式消元公式消元公式回代公式回代公式回代公式回代公式对对k=1,2,k=1,2,n-1,n-1,若若则可得则可得第11页,共54页,编辑于2022年,星期一顺序顺序GaussGauss消去法可执行的前提消去法可执行的前提定理定理 1 1 给定线性方程组给定线性方程组 ,如果,如果n n阶方阵阶方阵 的所有顺序主子式都不为零,即的所有顺序主子式都不为零,即 则按顺序则按顺序GaussGauss消去法所形成的各主元素消去法所形成的各主元素 均不为零,从而顺序均不为零,从而顺序Gauss 消去法可顺利执行。消去法
10、可顺利执行。注:注:当线性方程组的系数矩阵为当线性方程组的系数矩阵为对称正定对称正定或或严格对角占优阵严格对角占优阵时,时,按顺序按顺序GaussGauss消去法计算是稳定的。消去法计算是稳定的。第12页,共54页,编辑于2022年,星期一3、列主元、列主元Gauss消去法计算步骤:消去法计算步骤:1、输入矩阵阶数输入矩阵阶数n,增广矩阵增广矩阵 A(n,n+1);2、对于对于(1)按列选主元:选取按列选主元:选取 l 使使(2)如果如果 ,交换,交换 A(n,n+1)的第的第k行与第行与第l 行元素行元素(3)消元计算消元计算 :3、回代计算回代计算第13页,共54页,编辑于2022年,星期
11、一4 4无回代过程的主元消去法无回代过程的主元消去法(Gauss-Jordan(Gauss-Jordan消去法消去法)第一步:选主元,在第一列中选绝对值最大的元素,设第第一步:选主元,在第一列中选绝对值最大的元素,设第k行为行为主元行,将主元行换至第一行,将第一个方程中主元行,将主元行换至第一行,将第一个方程中x1的系数变为的系数变为1,并从其余并从其余n 1n 1个方程中消去个方程中消去x1。第二步:在第二列后第二步:在第二列后n 1个元素中选主元,将第二个方程中个元素中选主元,将第二个方程中 x2的系数变为的系数变为1,并从其它,并从其它n n 1 1个方程中消去个方程中消去x2。第第k步
12、:在第步:在第k列后列后n k个元素中选主元,换行,将第个元素中选主元,换行,将第k个方程个方程 xk的系数变为的系数变为1,从其它,从其它n-1n-1个方程中消去变量个方程中消去变量xk第14页,共54页,编辑于2022年,星期一消元公式为:消元公式为:对对k=1,2,按上述步骤进行到第按上述步骤进行到第n步后,方程组变为:步后,方程组变为:即为所求的解即为所求的解.第15页,共54页,编辑于2022年,星期一注:注:无回代的无回代的GaussGauss消元法实际上就是将消元法实际上就是将 方方程组的系数矩阵化为行最简形矩阵。程组的系数矩阵化为行最简形矩阵。第16页,共54页,编辑于2022
13、年,星期一5无回代消去法的应用无回代消去法的应用(1)解线性方程组系解线性方程组系设要解的线性方程组系为:设要解的线性方程组系为:AX=b1,AX=b2,AX=bm上述方程组系可以写为上述方程组系可以写为AX=B=(b1,bm)第17页,共54页,编辑于2022年,星期一因此因此X=A-1B 即为线性方程组系的解。即为线性方程组系的解。在计算机上只需要增加几组右端常数项的存贮单元,在计算机上只需要增加几组右端常数项的存贮单元,其结构和解一个方程组时一样。其结构和解一个方程组时一样。行行系数系数右端右端第18页,共54页,编辑于2022年,星期一(2)求逆矩阵求逆矩阵设设A=(aij)n n是非
14、奇矩阵是非奇矩阵,A 0,且令且令由于由于 AA-1=AX=I因此,求因此,求A-1的问题相当于解下列线性方程组的问题相当于解下列线性方程组相当于相当于(1)中中m=n,B=I 的情形。的情形。第19页,共54页,编辑于2022年,星期一(3)求行列式的值求行列式的值用高斯消去法将用高斯消去法将 A 化成上三角形化成上三角形第20页,共54页,编辑于2022年,星期一例例 用用Gauss-Jordan消去法解方程组消去法解方程组,并求出并求出 其中其中 解:把系数矩阵、单位矩阵和右端项组成增广矩阵,对增广矩解:把系数矩阵、单位矩阵和右端项组成增广矩阵,对增广矩 阵实行阵实行Gauss-Jord
15、an消元过程消元过程。第21页,共54页,编辑于2022年,星期一故故,系数矩阵的逆为系数矩阵的逆为 第22页,共54页,编辑于2022年,星期一2 2 解三对角方程组的追赶法解三对角方程组的追赶法第23页,共54页,编辑于2022年,星期一第24页,共54页,编辑于2022年,星期一33 矩阵的矩阵的三角分解法三角分解法 高斯消元法的矩阵形式高斯消元法的矩阵形式 每一步消去过程相当于左乘初等下三角矩阵每一步消去过程相当于左乘初等下三角矩阵L Lk k记记其中,其中,第25页,共54页,编辑于2022年,星期一记记其中,其中,一般的,对一般的,对第26页,共54页,编辑于2022年,星期一A
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 线性代数方程组的直接法 2幻灯片 线性代数 方程组 直接 幻灯片
限制150内