《投入产出平衡》PPT课件.ppt

实验一,投入产出平衡,矩阵和线性方程组,一、国民经济投入产出综合平衡,设有n个经济部门，xi为部门i的总产出，cij为部门j单位产品对部门i产品的消耗，di为外部对部门i的需求，fj为部门j新创造的价值。那么各经济部门总产出应满足下列关系式：,消耗平衡方程组,j=1,2,n,令C=（cij），X=(x1,xn),D=(d1,dn)，F=(f1,fn)则X=CX+D令A=EC，E为单位矩阵，则AX=D,C称为直接消耗矩阵，A称为列昂杰夫（Leontief）矩阵。,分配平衡方程组,i=1,2,n,Y=1,1,1B,Y表示各部门的总投入，称为投入向量。,新创造价值向量F=XY,B=C,B表示各部门间的投入产出关系，称为投入产出矩阵。,二、数学理论复习：线性代数,1、线性方程组,记为Ax=b其中A=(aij)mnx=(x1,xn),b=(b1,bm),若秩(A)秩(A,b)，则无解；若秩(A)=秩(A,b)=n,存在唯一解；若秩(A)=秩(A,b)n,存在无穷多解；通解是齐次线性方程组Ax=0的基础解系与Ax=b的一个特解之和。,对于线性方程组Ax=b：,Ax=0称为齐次的线性方程组,2、逆矩阵,方阵A称为可逆的，如果存在方阵B，使AB=BA=E，记B=A-1方阵A可逆的充分必要条件:A0,求逆矩阵方法：,A-1=A*/|A|这里A*为A的伴随矩阵,（AE）行变换,（EA-1）,3、特征值与特征向量,对于方阵A，若存在数和非零向量x使Ax=x，则称为A的一个特征值，x为A的一个对应于特征值的特征向量。,特征值计算归结为:特征多项式|A-E|=0的求根。对应于特征值的特征向量是齐次线性方程组(A-E)x=0的所有非零解,三、使用MATLAB,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(x1,x2,n)生成x1与x2间的n维等距行向量,即将x1,x2n-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解下列方程组,A=12;3-2;B=1;4;x=AB求得唯一解,A=121;3-21;B=1;4;x=AB求得一特解,A=12;3-2:1-1;B=1;4;2;x=AB求得一最小二乘近似解,A=12;-2-4;B=1;-2;x=AB不能直接求解,A=12;-2-4;00;B=1;-2;0;x=AB仍可求一近似特解,增加方程0 x+0y=0,例2线性方程组的通解,解在无穷多解情况下可用三种方法求通解，用rref化为行最简形以后求解；用除法求出一个特解，再用null求得一个齐次组的基础解系；用符号工具箱中的solve求解。,a=1-11-1;-111-1;2-2-11;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.65元的煤费，0.05元的电费及0.05元的运输费;创收一元钱的运输费,铁路要支付0.55元的煤费和0.10元的电费，在某一周内煤矿接到外地金额50000元定货，发电厂接到外地金额25000元定货，外界对地方铁路没有需求。,解：这是一个投入产出分析问题。设x1为本周内煤矿总产值，x2为电厂总产值,x3为铁路总产值,则,问三个企业间一周内总产值多少才能满足自身及外界需求？三个企业间相互支付多少金额？三个企业各创造多少新价值?,直接消耗矩阵C=,外界需求向量D=,产出向量X=,则原方程为(E-C)X=D,投入产出矩阵为B=C*diag(X)总投入向量Y=ones(1,3)*B新创造价值向量F=X-Y,例4（隐性病遗传）染色体遗传中，后代是从父母体的基因对中各继承一个基因，形成自己的基因型。如果所考虑的遗传特征是由两个基因A和a控制，那么就有三种基因型，,上表给出父母基因型的所有可能组合使其后代形成每种基因对的概率。,设金鱼某种遗传病染色体的正常基因为A，不正常基因为a,那么AA，Aa，aa分别表示正常金鱼，隐性患者，显性患者。设初始分布为90%正常金鱼，10%的隐性患者，无显性患者。考虑下列两种配种方案对后代该遗传病基因型分布的影响,方案一：同类基因