线性代数应用案例1.doc

行列式的应用
案例1 大学生在饮食方面存在很多问题，多数大学生不重视吃早餐，日常饮食也没有规律，为了身体的健康就需要注意日常饮食中的营养。大学生每天的配餐中需要摄入一定的蛋白质、脂肪和碳水化合物，下表给出了这三种食物提供的营养以及大学生的正常所需营养（它们的质量以适当的单位计量）。
营养
单位食物所含的营养
所需营养
食物1
食物2
食物3
蛋白质
36
51
13
33
脂肪
0
7
1.1
3
碳水化合物
52
34
74
45
试根据这个问题建立一个线性方程组，并通过求解方程组来确定每天需要摄入的上述三种食物的量。

解：设分别为三种食物的摄入量，则由表中的数据可以列出下列方程组

利用matlab可以求得
x =
0.27722318361443
0.39192086163701
0.23323088049177
案例2 一个土建师、一个电气师、一个机械师组成一个技术服务社。假设在一段时间内，每个人收入1元人民币需要支付给其他两人的服务费用以及每个人的实际收入如下表所示，问这段时间内，每人的总收入是多少？（总收入=实际收入+支付服务费）
服务者
被服务者
实际收入
土建师
电气师
机械师
土建师
0
0.2
0.3
500
电气师
0.1
0
0.4
700
机械师
0.3
0.4
0
600

解：设土建师、电气师、机械师的总收入分别是元，根据题意，建立方程组

利用matlab可以求得
x =
1.0e+003 *
1.25648414985591
1.44812680115274
1.55619596541787

案例3 医院营养师为病人配制的一份菜肴由蔬菜、鱼和肉松组成，这份菜肴需含1200cal热量，30g蛋白质和300mg维生素c，已知三种食物每100g中的有关营养的含量如下表，试求所配菜肴中每种食物的数量。

蔬菜
鱼
肉松
热量/cal
60
300
600
蛋白质/g
3
9
6
维生素c/mg
90
60
30

解：设所配菜肴中蔬菜、鱼和肉松的数量分别为百克，根据题意，建立方程组

利用matlab可以求得
x =
1.52173913043478
2.39130434782609
0.65217391304348

矩阵的应用
案例1 矩阵概念的引入
（1）线性方程组

的系数按原来的位置构成一数表

该数表决定着上述方程组是否有解，以及如果有解，解是什么等问题，因而研究这个数表就很重要。
A
B
C
D
（2）某航空公司在A,B,C,D四城市之间开辟了若干航线，下图所示表述了四城市间的航班图，若从A到B有航班，则用带箭头的线连接A和B。

为了便于研究，表中√为1，空白为0，得到下列数表：

列表表示到站

A
B
C
D
行标表示发站
A

√
√

0
1
1
0
B
√

√
√
1
0
1
1
C
√
√

√
1
1
0
1
D

√

0
1
0
0

(3)某中学学生身高体重的测量，得到如下一份统计如下表

体重（kg）
人数
身高（m）

40
50
60
70
1.5
60
80
70
20
1.6
30
120
150
90
1.7
10
15
80
150
1.8
0
2
5
10

此表反映身高与体重这种关系时也可将上面表格写成一个简化的4行4列的矩形数表，

40（kg）
50（kg）
60（kg）
70（kg）
1.5
60
80
70
20
1.6
30
120
150
90
1.7
10
15
80
150
1.8
0
2
5
10
如果只反映1.5米与体重的关系，则可以用（60 80 70 20）；如果只反映60kg与身高的关系，则可以用。
案例5 矩阵概念的应用——逻辑判断问题
甲、乙、丙、丁四人各从图书馆借来一本小说，他们约定读完后互相交换，这四本书的厚度以及他们四人的阅读速度差不多，因此，四人总是同时交换书，经三次交换后，他们四人读完了这四本书，现已知：
（1）乙读的最后一本书是甲读的第二本书；
（2）丙读的第一本书是丁读的最后一本书。
问四人的阅读顺序是怎样的？
1
2
3
4
甲 乙 丙 丁
解：设甲、乙、丙、丁最后读的书的代号依次为A,B,C,D,则根据题设条件可以列出初始矩阵

下面我们来分析矩阵中各位置的书名代号。已知每个人都读完了所有的书，所以并第二次读的书不可能是C,D。又甲第二次读的书是B，所以丙第二次读的书也不可能是B，从而丙第二次读的书是A，同理可依次推出丙第三次读的书是B，丁第二次读的书是C，丁第三次读的书是A，丁第一次读的书是B，乙第二次读的书是D，甲第一次读的书是C，乙第一次读的书是A，乙第三次读的书是C，甲第三次读的书是D。故四人阅读的顺序可用矩阵表示如下：

案例6 矩阵乘法的应用
某企业某年出口到三个国家的两种货物的数量及两种货物的单位价格、重量、体积如下表所示：
国
家
数
量
货

美国
德国
日本

3000
1500
2000

1400
1300
800

单位价格（万元）
单位重量（吨）
单位体积（）

0.5
0.04
0.2

0.4
0.06
0.4

利用矩阵乘法计算该企业出口到三个国家的货物总价值、总重量、总体积各为多少？

解：设矩阵

则矩阵
总价值 总重量 总体积

美国
德国
日本
=
案例7 逆矩阵的应用
一个城市有三个重要的企业：一个煤矿，一个发电厂和一条地方铁路。开采一块钱的煤，煤矿必须支付0.25元的运输费。而生产一块钱的电力，发电厂需支付煤矿0.65元的燃料费，自己亦需支付0.05元的电费来驱动辅助设备及支付0.05元的运输费。而提供一块钱的运输费铁路需支付煤矿0.55元的燃料费，0.10元的电费驱动它的辅助设备。某个星期内，煤矿从外面接到50000元煤的订货，发电厂从外面接到25000元电力的订货，外界对地方铁路没有要求。问这三个企业在那一个星期的生产总值各为多少时才能精确地满足它们本身的要求和外界的要求？

解：各企业产出一元钱的产品所需费用为
企
业
产
品
费
用

煤矿
发电厂
铁路
燃料费（元）
0
0.65
0.55
电力费（元）
0
0.05
0.10
运输费（元）
0.25
0.05
0

对于一个星期的周期，设表示煤矿的总产值，表示电厂的总产值，表示铁路的总产值。
煤矿的总消耗为
电厂的总消耗为
铁路的总消耗为
则

联立三个方程并整理得方程组

上述方程组可化为 ，其中
，

利用matlab求解，可知 ，所以方程组有唯一解，其解为

所以煤矿总产值为80423元，发电厂总产值为28583元，铁路总产值为21535元。

案例8 求解线性方程组
（1）假设你是一个建筑师，某小区要建设一栋公寓，现在有一个模块构造计划方案需要你来设计，根据基本建筑面积每个楼层可以有三种设置户型的方案，如下表所示。如果要设计出含有136套一居室，74套两居室，66套三居室，是否可行？设计方案是否唯一？
方案
一居室（套）
两居室（套）
三居室（套）
A
8
7
3
B
8
4
4
C
9
3
5

解：设公寓的每层采用同一种方案，有层采用方案A，有层采用方案B，有层采用方案C，根据题意，可得

利用matlab计算方程组的系数矩阵A、增广矩阵的秩：
，
所以方程组有无穷多个解。
利用matlab将增广矩阵化为行简化阶梯型矩阵：

矩阵对应的方程组为，
取，则方程组的全部解为

又由题意可知，都为正整数，则方程组有唯一解。
所以设计方案可行且唯一，设计方案为：6层采用方案A，2层采用方案B，8层采用方案C。

（2） 在一个原始部落中，农田耕作记为F，农具及工具的制作记为M，织物的编织记为C。人们之间的贸易是实物交易系统（见下图）。由图中可以看出，农夫将每年的收获留下一半，分别拿出四分之一给工匠和织布者；工匠平均分配他们制作的用具给每个组。织布者则留下四份之一的衣物为自己，四分之一给工匠，二分之一给农夫。
随着社会的发展，实物交易形式需要改为货币交易。假设没有资本和负债，那么如何对每类产品定价才能公正地体现原有的实物交易系统？
F
M
C

也可以用下表表示：
组名
F
M
C
F

M

C

解：令为农作物的价值，为工具的价值，为织物价值。那么从上表第一列，农夫生产的价值应该等于他们交换到的产品的价值，即

同理可以得到工匠和纺织者产品价值的方程
，
从而得到下列方程组：

利用matlab将系数矩阵化为行简化阶梯型矩阵，为
A=
令，写成方程组，为

写成向量形式为

所以当农作物价值、工具价值与织物价值的定价之比为时，才能公正地体现原有的实物交易系统。

（3）某药厂生产3种中成药，每件中成药的生产要经过3个车间加工。3个车间每周的工时、每件中成药在各车间需要的工时数如下表所示，问3中中成药每周的产量各是多少？

中成药1
中成药2
中成药3
车间工时（时/周）
车间1
1
1
2
40
车间2
3
2
3
75
车间3
1
1
1
28

解：设3种中成药每周的产量分别为，则由题意得

利用matlab将方程组的增广矩阵化为行简化阶梯型矩阵，得

由此可以得出
所以三种中成药每周的产量分别为7件，9件，12件。

案例9解线性方程组应用—人口迁移模型
在生态学、经济学和工程学等许多领域中经常需要对随时间变化的动态系统进行数学建模，此类系统中的某些量常按离散时间间隔来测量，这样就产生了与时间间隔相应的向量序列其中表示第n次测量时系统状态的有关信息，而常被称为初始向量。
如果存在矩阵A，并给定初始向量，使得即

则上述方程为一个线性差分方程或者递归方程。

（1）已知某城市2009年的城市人口为5000000人，农村人口为7800000人。假设每年大约有5%的城市人口迁移到农村（95%仍然留在城市），12%的农村人口迁移到城市（88%仍然留在农村），如下图所示，忽略其他因素对人口规模的影响。计算2011年的人口分布。

城市

农村
0.95
0.12
0.05
0.88

解：由题意可得迁移矩阵为

设2009年的初始人口为，2010年和2011年的人口分别为，则

即2011年的人口分布情况是：城市人口为6255380，农村人口为6544620.

(2)在某个地区，每年约有4%的城市人口移居到周围的农村，大约5%的农村人口移居到城市中。在2009年，城市中有400000居民，农村有600000居民。建立一个差分方程来描述这种情况，用表示2009年的初始人口，然后估计两年之后，即2011年城市和农村的人口数量（忽略其他因素对人口规模的影响）

（3）某公司有一个车队，大约有450辆车，分布在三个地点。一个地点租出去的车可以归还到三个地点中的任意一个，但租出的车不许当日归还。下面的矩阵给出了汽车归还到每个地点的不同比率。假设星期一在机场有304辆车，东部办公区有48辆车，西部办公区有98辆车，那么在星期三时，车辆的大致分布式怎么样？

解：设星期一机场、东部和系部的车辆为，星期二和星期三三个地方的车辆分别为，由题意可得，迁移矩阵为

则

所以，星期三时，机场有310辆车，东部办公区有48辆车，系部办公区有92辆车。

案例10 解线性方程组应用—网络流模型
网络流模型广泛应用于交通、运输、通信、电力分配、城市规划、任务分派以及计算机辅助设计等众多领域。当科学家、工程师和经济学家研究某种网络中的流量问题是，线性方程组就自然而然地产生了，例如：城市规划设计人员和交通工程师监控城市道路网络内的交通流量，电气工程师计算电路中流经的电流，经济学家分析产品通过批发商和零售商网络从生产者到消费者的分配等。大多数网络流模型中的方程组都包含了数百甚至上千个未知量和线性方程。
一个网络由一个点集以及连接部分或全部点的直线或弧线构成。网络中的点称作联结点（或节点），网络中的连接线称作分支。每一分支中的流量方向已经指定，并且流量（或流速）已知或者已标为变量。
网络流的基本假设是网络流中流入与流出的总量相等，并且每个联接点流入和流出的总量也相等。例如：下图所示分别说明了流量从一个或两个分支流入联结点，和分别表示从其他分支流出的流量，表示从其他分支流入的流量。因为流量在每个联结点守恒，所以有。
网络分析要解决的问题就是：在部分信息（如网络的输入量）已知的情况下，确定每一分支中的流量。
60
60
60

80

下图的网络给出了在下午两点钟，某市区部分单行道的交通流量（以每刻钟通过的汽车数量来度量）。试确定网络的流量模式。
40
50
30
30
20
10

B
A
C
D

解：根据网络流模型的基本假设，在节点（交叉口）A,B,C,D处，我们可以得到下列方程：

此外，该网络的总流入等于网络的总流出，即

联立以上方程的方程组：

取，则网络的流量模式表示为

线性规划问题
案例1、.生产计划问题
（1）假设某厂计划生产甲、乙两种产品，现库存主要原料有A类3600kg，B类2000kg,C类3000kg.每件甲产品需用材料A类9kg,B类4kg,C类3kg。每件乙产品需用材料A类4kg,B类5kg,C类10kg。甲单位产品的利润70元，乙单位产品的利润120元。问如何安排生产，才能使该厂所获的利润最大。
解：建立模型：
设生产甲、乙产品的件数分别为，为该厂所获总利润，则

程序如下：
f=[-70 -120]';
A=[9 4;4 5;3 10];
b=[3600 2000 3000]';
[x,maxf]=linprog(f,A,b);
maxf=-maxf;
结果为：x =
200.0000
240.0000
maxf =
4.2800e+004

（2）某工厂生产A,B两种产品，已知生产A产品每公斤需耗煤9吨，耗电400度，用工3个工作日；生产B产品每公斤需耗煤4吨，耗电500度，用工10个劳动日。A产品每公斤利润700元，B产品每公斤利润1200元，因客观条件限制，该厂只能得到煤360吨，电20000度，劳动力300个，问该厂如何安排生产才能使总利润最大？
解：建立模型
设分别表示生产A,B两种产品的数量，则

案例2、投资问题
某公司有一批资金用于4个工程项目的投资，其投资各项目时所得的净收益（投入资金的百分比）如下表所示。由于某种原因，决定用于项目A的投资不大于其他各项投资之和，而用于项目B和C的投资要大于项目D的投资。试确定该公司收益最大的投资分配方案。
工程项目
A
B
C
D
收益（%）
15
10
8
12
解：建立模型：
设分别表示用于项目A,B,C,D的投资百分数，则

程序如下：
f=[-0.15 -0.1 -0.08 -0.12]';
A=[1 -1 -1 -1;0 -1 -1 1];
b=[0 0]';
Aeq=[1 1 1 1];
beq=[1];
ub=[];
lb=zeros(4,1);
[x,maxf]=linprog(f,A,b,Aeq,beq,lb,ub);
maxf=-maxf;
结果如下：
x =
0.5000
0.2500
0.0000
0.2500
maxf =
0.1300

案例3、运输问题
有A,B,C三个食品加工厂，负责供给甲、乙、丙、丁4个市场。3个工厂每天生产食品箱数如表1所示；4个市场每天的需求量如表2所示；从各厂运到各市场的运输费（元/箱）如表3所示。求在基本满足供需平衡的约束条件下使总运输费用最小。
表1
工厂
A
B
C
生产数
60
40
50
表2
市场
甲
乙
丙
丁
需求量
20
35
33
34
表3
收点
发点
市场
甲
乙
丙
丁
工厂
A
2
1
3
2
B
1
3
2
1
C
3
4
1
1
解：建立数学模型：
设为由工厂运到市场的费用，是由工厂运到市场的箱数。是工厂的产量，是市场的需求量，则
，

程序如下：
B=[2 1 3 2;1 3 2 1;3 4 1 1];
f=B(:);%将矩阵B转化为列向量
A=[1 0 0 1 0 0 1 0 0 1 0 0
0 1 0 0 1 0 0 1 0 0 1 0
0 0 1 0 0 1 0 0 1 0 0 1];
b=[60 40 50]';
Aeq=[1 1 1 0 0 0 0 0 0 0 0 0
0 0 0 1 1 1 0 0 0 0 0 0
0 0 0 0 0 0 1 1 1 0 0 0
0 0 0 0 0 0 0 0 0 1 1 1];
beq=[20 35 33 34]';
lb=zeros(12,1);
ub=[];
[x,minf]=linprog(f,A,b,Aeq,beq,lb,ub); 结果如下：
x =
0.0000
20.0000
0.0000
35.0000
0.0000
0.0000
0.0000
0.0000
33.0000
0.0000
18.4682
15.5318
minf =
122.0000
案例4、合理下料问题
某工厂有一批5米的钢管（数量充足），为制造零件的需要，要将它们截成140厘米，95厘米，65厘米的管料，并要求这三种管料按照2：4：1的比例配套，问如何下料才能使残料最少？
解：合理下料问题：
分析：截取的方法很多，但是残料较大的可以排除

1
2
3
4
5
6
7
8
140
3
2
2
1
1
0
0
0
95
0
2
0
3
1
5
3
1
65
1
0
3
1
4
0
3
6
残料
15
30
25
10
5
25
20
15
设表示第种截法截得的钢管数量

案例5、经济配料问题
某饲养场有5中饲料，已知各种饲料的单位价格和每百公斤饲料的蛋白质、矿物质、维生素含量如表，又知该饲养场每日至少需蛋白质70单位，矿物质3单位，维生素10单位，问如何混合调配5种饲料，能使总成本最低？
饲料种类
有关成分
饲料单价
（百公斤）
蛋白质
矿物质
维生素
1
0.30
0.1
0.05
2
2
2.20
0.05
0.10
7
3
1.00
0.02
0.02
4
4
0.60
0.20
0.20
3
5
1.80
0.05
0.08
5
解：设用表示第种饲料的用量，