投入产出平衡

1、一、国民经济投入产出综合平衡一、国民经济投入产出综合平衡 设有n个经济部门，xi为部门i的总产出，cij为部门j单位产品对部门i产品的消耗，di为外部对部门i的需求，fj为部门j新创造的价值。那么各经济部门总产出应满足下列关系式：消耗平衡方程组xxcfjjijinj1j=1,2,n令 c =（cij），x = (x1, , xn) ,d = (d1, , dn)，f= (f1, , fn)则 x=cx+d令 a = ec，e为单位矩阵，则 ax = dc称为直接消耗矩阵，a称为列昂杰夫（leontief）矩阵。分配平衡方程组xc xdiijjnji1i =1,2,ny = 1,1,1 by表示

2、各部门的总投入，称为投入向量。新创造价值向量 f=x y b=cb表示各部门间的投入产出关系，称为投入产出矩阵。xxxn12二、数学理论复习：线性代数数学理论复习：线性代数1、线性方程组a xa xa xba xa xa xba xaxaxbnnnnmmmnnm11 11221121 1222221 122记为记为 a x = b 其中a =(aij)mn x = (x1, ,xn), b = (b1, , bm)若秩(a) 秩(a,b)，则无解；若秩(a) = 秩(a,b) = n, 存在唯一解；若秩(a) = 秩(a,b) n, 存在无穷多解； 通解是齐次线性方程组 ax=0 的基础解系与

3、 ax=b 的一个特解之和。对于线性方程组 ax = b：ax = 0 称为齐次的线性方程组2、逆矩阵方阵a称为可逆的，如果存在方阵b，使a b = b a = e，记 b = a-1方阵a可逆的充分必要条件:a0求逆矩阵方法： a-1 =a*/|a| 这里a*为a的伴随矩阵 （a e） 行变换（e a-1）3、特征值与特征向量对于方阵a，若存在数和非零向量x 使 a x = x，则称为a的一个特征值，x 为a 的一个对应于特征值的特征向量。特征值计算归结为:特征多项式|a - e|=0的求根。对应于特征值的特征向量是齐次线性方程组 (a - e) x = 0的所有非零解三、使用三、使用mat

4、lab det 方阵的行列式 diag 对角阵inv 方阵的逆 cond 方阵的条件数trace 方阵的迹 orth 正交规范化rank 矩阵的秩 null 求基础解系rref 矩阵的行最简形eig 特征值与特征向量jordan 约当标准形分解norm 矩阵或向量范数1、特殊矩阵生成zeros(m,n) 生成m行n列的零矩阵;ones(m,n) 生成m行n列的元素全为1的阵;eye(n) 生成n阶单位矩阵;当a是矩阵,diag(a)返回a的对角线元素构成的向量;当x是向量,diag(x)返回由x的元素构成的对角矩阵；rand(m,n) 生成m行n列0,1上均匀分布随机数矩阵；linspace(

5、x1,x2,n) 生成x1与x2间的n维等距行向量,即将x1,x2 n-1等分。2、行列式和逆矩阵det(a) 返回方阵a的行列式;inv(a) 返回a的逆矩阵。3、矩阵除法左除法 ab 求解矩阵方程ax=b右除法 b/a 求解矩阵方程xa=b(1) 当a为方阵，ab与inv(a)*b基本一致： (2) 当a不是方阵，除法将自动检测。 若方程组无解,除法给出最小二乘意义上的近似解，即使向量axb的长度达到最小； 若方程组有无穷多解，除法将给出一个具有最多零元素的特解； 若为唯一解，除法将给出解。例1 解下列方程组42312yxyx42312zyxzyxa=1 2;3 -2; b=1;4;x=a

6、b 求得唯一解a=1 2 1;3 -2 1; b=1;4;x=ab 求得一特解242312yxyxyx24212yxyx a=1 2;3 -2:1 -1; b=1;4;2;x=a b 求得一最小二乘近似解a=1 2;-2 -4; b=1;-2;x=ab 不能直接求解a=1 2;-2 -4;0 0; b=1;-2;0;x=ab仍可求一近似特解增加方程 0 x+0y=0例2 线性方程组的通解12211432143214321xxxxxxxxxxxx解 在无穷多解情况下可用三种方法求通解， 用rref化为行最简形以后求解；用除法求出一个特解，再用null求得一个齐次组的基础解系；用符号工具箱中的so

7、lve求解。a=1 -1 1 -1;-1 1 1 -1;2 -2 -1 1; b=1;1;-1; rank(a),rank(a,b) rref(a,b)x0=ab,x=null(a)4、特征值和特征向量d=eig(a) 返回方阵a的特征值构成的列向量；v,d=eig(a) 返回方阵a的特征值构成的对角阵d和每个特征值对应的特征向量按列构成的矩阵v。其中每个特征向量都是模等于1的向量，并且属于同一特征值的线性无关特征向量已正交化。四、实验例题 例3 某地有三个产业，一个煤矿，一个发电厂和一条铁路，开采一元钱的煤，煤矿要支付0.25元的电费及0.25元的运输费; 生产一元钱的电力，发电厂要支付0.

8、65元的煤费，0.05元的电费及0.05元的运输费; 创收一元钱的运输费,铁路要支付0.55元的煤费和0.10元的电费，在某一周内煤矿接到外地金额50000元定货，发电厂接到外地金额25000元定货，外界对地方铁路没有需求。解：这是一个投入产出分析问题。设x1为本周内煤矿总产值，x2为电厂总产值, x3为铁路总产值, 则xxxxxxxxxxxx11232123312300 65055500000 250 05010250000 250 0500(.)( .)( .) 问三个企业间一周内总产值多少才能满足自身及外界需求？三个企业间相互支付多少金额？三个企业各创造多少新价值?直接消耗矩阵c= 02

9、500050000外界需求向量 d =产出向量x = xxx12300650550250050100250050.则原方程为则原方程为 (e-c)x=d 投入产出矩阵为 b=c*diag(x)总投入向量 y= ones(1,3)*b 新创造价值向量 f=x-y例4 （隐性病遗传）染色体遗传中，后代是从父母体的基因对中各继承一个基因，形成自己的基因型。如果所考虑的遗传特征是由两个基因a和a控制，那么就有三种基因型，父 母概率aa-aa aa-aaaa-aaaa-aaaa-aaaa-aaaa11/201/400aa01/211/21/20后代aa0001/41/21上表给出父母基因型的所有可能组合使其后代形成每种基因对的概率。设金鱼某种遗传病染色体的正常基因为a，不正常基因为a, 那么aa，aa，aa分别表示正常金鱼，隐性患者，显性患者。设初始分布为90%正常金鱼，10%的隐性患者，无显性患者。考虑下列两种配种方案对后代该遗传病基因型分布的影响方案一：同类基因结合，均可繁殖；方案二：显性患者不允许繁殖，