1 / 11
文档名称:

数学实验矩阵运算.docx

格式:docx   大小:112KB   页数:11页
下载后只包含 1 个 DOCX 格式的文档,没有任何的图纸或源代码,查看文件列表

如果您已付费下载过本站文档,您可以点这里二次下载

分享

预览

数学实验矩阵运算.docx

上传人:AIOPIO 2021/7/29 文件大小:112 KB

下载得到文件列表

数学实验矩阵运算.docx

文档介绍

文档介绍:数学实验报告
学  院:            ﻩ
班 级:  ﻩ
学  号:   ﻩ
姓  名:       ﻩ
完成日期:  
ﻬ实验四  矩阵的运算
(一)投入产出分析

1。理解投入产出分析中的基本概念和模型;
,理解矩阵乘法、逆矩阵等的含义。
二.问题描述
设国民经济由农业、制造业和服务业三个部门构成,已知某年它们之间的投入产出关系、部需求、初始投入等如表1-1所示
表1—1国民经济三产部门之间的投入产出表
产出
农业
制造业
服务业
外部需求
总产出


农业
制造业
3
服务业
20
60
0
70
150
初始投入
35
110
75
总投入
100
200
150
根据表回答下列问题:
(1)如果农业、制造业、服务业外部需求为50,150,100,问三个部门总产出分别为多少?
(2)如果三个部门的外部需求分别增加一个单位,问他们的总产出分别为多少?
三。实验过程
(1)的求解
(1)求直接消耗矩阵A
根据直接消耗的计算公式
aij=xij/xj
和各部门中间需求;
xnan
运行如下代码可得直接消耗系数表。
X=[15 20 30;30 10 45;20 60 0];
X_colsum=[100 200 150];
X_rep=repmat(X_colsum,3,1)
A=X。/ X_rep
运行结果为:
A =
0。1500    0.1000   0。2000
         0。3000
0。2000       0
(2)求解
根据公式
X=(I—A)-1y
在运行如下代码
y=[50;150;100];
n=size(y,1);
W=eye(n)—A;
X=W\y
运行结果为
X =
  

208。1377
即三个部门的总产出分别为139.2801,, 208.1377亿元。
2。问题2求解
设外部需求由y增加至y+Δy,则产出x的增量为
    Δx=(I—A)-1(y+Δy)— (I—A)—1y=(I—A)-1Δy
利用问题(1)求得的I-A矩阵,再运行如下的MATLAB代码可得问题的结果:
dx=inv(W)
运行结果:
dx =
    0。3443
 0.5634   1.2676   0。4930
  0.4382    0。4304  1.2167
根据上述结果可知,当农业的外部需求增加1个单位时,农业、制造业、,0。5634,0。4382个单位;当制造业的外部需求增加
1个单位时,农业、制造业、服务业的总产出分别增加0.2504,1。2676,0。4304个单位;当服务业的外部需求增加1个单位时,农业、制造业、,0。4930,1.2167个单位。

投入产出的理论依据,主要是矩阵运算和逆矩阵,投入和产出分析在编制和修订宏观计划、开展紧急预测和项目效果预测、研究价格水平及其变动影响,研究产业结构及其关联程度,投资对消费的带动分析等方面对会有很重要的应用。
(二)Hill密码的加密、解密与破译
一、实验目的
1. 复****线性代数, 矩阵, 线性空间与线性变换等概念和运算.  
2。 熟悉Hill密码体制的加密, 解密和破译过程。 
3。 Hill密码体制的加密, 解密和破译过程的MATLAB编程实践。
二、问题描述
对明文为“Mr Hill made this code”,利运用不同的密匙矩阵加密矩阵,实现Hill加密与解密过程。
三。实验过程
(1)模型建立于问题分析
简单起见,为了使明文参与矩阵运算,建立如下表格2—1:

B

Y
Z
空格  
句号
1

25
26
0
27
表2-1 26个英文字母、空格、句号与数字之间的对应关系
设明文为M=(m1, m2,…,mi)T,密匙矩阵为可逆的l x l的方阵,则经过线性变换得到密文C=EK(M)=(c1c2…cI)T,其中
C1=(k11m1+ k12m2+…+ k1lml)mod28,
C2=(k21m1+ k22m2+…+ k2lml)mod28,
……