第6章矩阵与线性方程.docVIP

第6章 矩阵与线性方程 本章主要涉及线性代数中向量和矩阵的运算、线性方程组的解法等知识，介绍线性代数再经济方面的应用，以及MATLAB有关线性代数运算的命令。 6.1 引例： 某城镇有三个重要产业，一个煤矿，一个发电厂和一条地方铁路。开采一元钱的煤，煤矿要支付0.25元的电费和0.35元的运输费；生产一元钱的电，发电厂要支付0.40元的煤费，0.05元的电费及0.10元的运输费；提供一元钱的运输服务，铁路要支付0.45元的煤费的0.10元的电费和0.10元的运输费，在某一周内煤矿接到外地金额50000元定货，发电厂接到外地金额25000元定货，外界对地方铁路需求为30000元。问： （1）三个企业间一周内总产值多少才能满足自身及外界需求？ （2）三个企业间相互支付多少金额？三个企业各创造多少新价值？ （3）如果煤矿需要增加总产值10000元，它对各个企业的产品或服务的完全需求分别将是多少？ （4）假定三企业的外部需求仍是用于城镇的各种消费和积累，其中用于消费的产品价值分别为35000元、18000元和20000元，而假定三个企业的新创造价值又包括支付劳动报酬（工资等）和纯收入，其中支付劳动报酬分别为25488元、10146元和14258元，试分析各企业产品使用情况的比例关系；以及该星期系统的经济效益； （5）若在以后的三周内，企业外部需求的增长速度分别是15%、3%和12%；那么各企业的总产值将增长多少？ 投入产出综合平衡分析： 在一个国家或区域的经济系统中，各部门（或企业）既有投入又有产出。生产的产品满足系统内部各部门和系统外的需求，同时也消耗系统内各部门内的产品。应如何组织生产呢？ 俄裔美国经济学家W.Leontief于20世纪30年代首先提出并成功地建立了研究国民经济投入产出地数学模型，他数次主持制定了美国的国民经济投入产出列表，且由此对国民经济各部门的结构和各种比例关系进行了定量分析。这一方法即投入产出法，投入产出法以其重要的应用价值迅速为世界各国经济学界和决策部门所采纳。W.Leontief因此于1973年获得了Nobel经济学奖。 数学模型：假设有n个经济部门，xi为部门i的总产出；aij为部门j单位产品对部门i产品的消耗； yi为外部对部门i的需求；zj为部门j 新创造的价值。则 （6.1）称为分配平衡方程组，每一个等式以价值形式说明了对每一部门： 中间产品（作为系统内个部门的消耗）+最终产品（外部需求）=总产品 （6.2）称为消耗平衡方程组，说明对每一部门： 对系统内各企业产品的消耗+新创价值=总价值 令 （6.1）化为 X=AX+Y (6.3) 令C=E-A, E—单位矩阵，(7.3)化为 CX=Z (6.4) 经济学上称：A—直接消耗矩阵，C—Leontief矩阵，令 B为投入产出矩阵，D为总投入向量，F为新创价值向量。 6.2 线性代数基本知识复习 线性方程组就是n元一次方程。实际的问题，如大型的土建结构、机械结构、大型的输电网络、管道网络、投入产出分析都可以归结为线性方程组，有些数值计算的方法也会导致线性方程组的求解，如数据拟合、非线性方程组和偏微分方程组数值解问题等等。n往往达到几百、几千甚至上万。 6.2.1 线性方程组 n个未知量m个方程的线性方程组的一般形式为 令 则得矩阵形式 (6.7) 若b=0，即 Ax=0 (6.8) 则称方程组为齐次的。 方程组(6.7)可能有唯一解、有无穷多解，也可能无解，这取决于系数矩阵A及增广矩阵(A,b)的秩。若秩(A)=秩(A,b)=n，存在唯一解。其理论上可用Cramer法则求出，但计算量太大。若秩(A)=秩(A,b)n，存在无穷多解，其通解可表示为对应齐次方程组(6.8)的一个基础解系与(6.7)式的一个特解的叠加；若秩(A)≠秩(A,b)，则无解，这时一般寻求最小二乘近似解，即求x使向量Ax-b长度最小。 6.2.2 逆矩阵 方阵A称为可逆的，如果存在方阵B，使 AB=BA=E 这里E表示单位阵。B称为A的逆矩阵，记。A可逆的充分必要条件是。且 (6.9) 这里A*为A的伴随矩阵。(6.7)可表示为 由于这个公式涉及大量的行列式计算，数值计算不采用。 6.2.3 特征值与特征向量 对于方阵A，若存在数(和非零向量x使 称(为A的一个特征值，x为A的一个对应于特征值(的特征向量。特征值的计算归结为特征多项式的求根。对于n阶实数方阵，特征多项式在复数范围内总有n个根。对应