数学建模基础知识:矩阵和线性方程组_第1页
数学建模基础知识:矩阵和线性方程组_第2页
数学建模基础知识:矩阵和线性方程组_第3页
数学建模基础知识:矩阵和线性方程组_第4页
数学建模基础知识:矩阵和线性方程组_第5页
已阅读5页,还剩17页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

实验一投入产出平衡——矩阵和线性方程组一、国民经济投入产出综合平衡设有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=E-C,E为单位矩阵,则AX=DC称为直接消耗矩阵,A称为列昂杰夫(Leontief)矩阵。分配平衡方程组i=1,2,…,nY=[1,1,…,1]BY表示各部门的总投入,称为投入向量。新创造价值向量F=X–Y'B=CB表示各部门间的投入产出关系,称为投入产出矩阵。二、数学理论复习:线性代数1、线性方程组记为Ax=b其中A=(aij)m×nx=

(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可逆的充分必要条件:

A

0求逆矩阵方法:

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,x2]n-1等分。2、行列式和逆矩阵det(A)返回方阵A的行列式;inv(A)返回A的逆矩阵。3、矩阵除法左除法A\B求解矩阵方程AX=B右除法B/A求解矩阵方程XA=B(1)当A为方阵,A\B与inv(A)*B基本一致:

(2)当A不是方阵,除法将自动检测。若方程组无解,

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论