线性代的MATLAB求解及线性议程组的应用.pptVIP

第三章 线性代数问题的MATLAB求解 踊脏且溃由托戚吭喷曰楔民潮扳狐之继畏篇韭索街窃熔找率服察惜乔腻枉线性代数的MATLAB求解及线性议程组的应用线性代数的MATLAB求解及线性议程组的应用 主要内容 一. 矩阵及其运算 二. 矩阵的初等变换与线性方程组 三. 矩阵的对角化 觅帆盈兹勒囚网河破顽难志窥料涅饮度除鞠屉泞壁坎豺膝勋座傈婶往滔仗线性代数的MATLAB求解及线性议程组的应用线性代数的MATLAB求解及线性议程组的应用 一. 矩阵及其运算 3.1.1矩阵的算术运算 MATLAB中矩阵的基本运算有+,-,*,/（右除）,\（左除）,^（乘方） 1、矩阵的线性运算 两个矩阵A,B,A+B和A-B的运算规则是：若两个矩阵的维数相同，则可以执行加减运算，相应的元素相加减，如果尾数不相同则出错。 2、矩阵的乘法 设两个矩阵A,B分别为m*n和n*p的矩阵，则C=A*B为m*p的矩阵 3、矩阵的除法（ mldivide \和mrdivide / ） 如果两个矩阵A,B是非奇异的，则A\B和A/B运算都可以实现。 A\B等价于inv(A)*B,A/B等价于A*inv(B) 对含有标量的运算，两种除法运算的结果相同3/4和4/3都是0.75 4、矩阵的乘方^ 誊叹鸭阀他赶铝猿迂烫虫装矮教侮偿匙澡络靶议愧霄媒忠铂咸军捎诲皱货线性代数的MATLAB求解及线性议程组的应用线性代数的MATLAB求解及线性议程组的应用 5、点运算 两个矩阵进行点运算是指他们的对应元素进行相关运算，要求两矩阵的维参数相同 6、矩阵的转置’ 3.1.2方阵的行列式 方阵A的行列式对应的求值函数为det(A) 3.1.3矩阵的逆矩阵 对于一个可逆的方阵A，求解函数为inv(A) 3.1.4应用实例-投入产出模型 问题描述：国民经济各部门之间存在相互依存的关系，每个部门在运转中将其他部门的半成品（称为投入）经过加工变为自己的产品（产出），如果根据各部门间的投入产出关系，确定各部门的产出水平以满足社会需求。 影吏怜彤拼染见开趣姆铁很给带良努靡帝馁翼削谱靴舍冻傣甫辐躁诵膘撕线性代数的MATLAB求解及线性议程组的应用线性代数的MATLAB求解及线性议程组的应用 设国民经济由农业、制造业、服务业三个部门构成，关系如下表（3-1）：单位(亿元) 投入/产出 农业 制造业 服务业 初始投入 总投入 农业 15 30 20 35 100 制造业 20 10 60 110 200 服务业 30 45 0 75 150 外部需求 35 115 7 总产出 100 200 150 三个部门的投入产出表如下（表3-2）：单位(亿元) 投入/产出 农业 制造业 服务业 农业 0.15 0.30 0.20 制造业 0.10 0.05 0.30 服务业 0.20 0.30 0 如第一行第二列的数字0.10表示生产1个单位产值的制造业产品需要投入0.10个单位产值的农产品，由表3-1中20亿元农产品投入制造业可以产出200亿元制造业总产值，20/200=0.1。该表中的数字称为投入系数或消耗系数。 饰央谚恰丰领苦习簧胃盾匀酱体月晓辜章格移全陀让哦紫纲笔常侨论彬愁线性代数的MATLAB求解及线性议程组的应用线性代数的MATLAB求解及线性议程组的应用 需要解决的问题： 1、如果今年对农业制造业和服务业的外部需求分别为50亿元、150亿元、100亿元，问这个三个部门的总产出分别是多少？ 2、如果三个部门的外部需求分别增加1个单位，问他们的总产出应分别增加多少？ 3、投入产出分析称为可行的，如果对于任意给定的、非负的外部需求，都能得到非负的总产品。为了可行，投入系数应该满足什么条件？ 模型建立：设有n个部门，记一定时期内第i个部门的总产出为xi,其中对第j个部门的投入是xij,外部需求为di,则： （1） 表3-1的每一行都满足(1)式，投入系数记为aij，即为第j个部门的单位产出所需要的第i个部门的投入： （2） 运仟蜕稚悄失用悬饰颓当蛛孵侮糠藉宾莎鲜弥脐边悉档任罩掂隋秒责时崩线性代数的MATLAB求解及线性议程组的应用线性代数的MATLAB求解及线性议程组的应用 由(1)和(2)式得到 （3） 记投入系数矩阵 产出向量 需求向量 (3)式可写为： （4） 或 （5） 当投入系数A和外部需求d给定后，求线性方程组(5)即可得到各部门的总产出x。 问题(1)的解答： 编写MATLAB程序：a=[0.15 0.1 0.2;0.3 0.05 0.3;0.2 0.3 0]; d=[50 150 100]’; b=eye(3)-a;x=b\d 摈饼瓦范鞍崖都隘卧饿魂橇个漱能耘外沈篆鸣间泣牙谴翔倡艾狈结晤安堂线性代数的MATLAB求解及线性议程组的应用线性代数的MATLAB求解及线性议程组的应用 问题(2)