行列式的应用
案例1 大学生在饮食方面存在很多问题,多数大学生不重视吃早餐,日常饮食也没有规
律,为了身体的健康就需要注意日常饮食中的营养。大学生每天的配餐中需要摄入一定的蛋白质、脂肪和碳水化合物,下表给出了这三种食物提供的营养以及大学生的正常所需营养(它们的质量以适当的单位计量)。
试根据这个问题建立一个线性方程组,并通过求解方程组来确定每天需要摄入的上述三种食物的量。
解:设123,,x x x 分别为三种食物的摄入量,则由表中的数据可以列出下列方程组
123231
2336511333
7 1.1352347445
x x x x x x x x ++=??
+=?
?++=? 利用matlab 可以求得 x =
0.27722318361443 0.39192086163701 0.23323088049177
案例2 一个土建师、一个电气师、一个机械师组成一个技术服务社。假设在一段时间,
每个人收入1元人民币需要支付给其他两人的服务费用以及每个人的实际收入如下表所示,问这段时间,每人的总收入是多少?(总收入=实际收入+支付服务费)
解:设土建师、电气师、机械师的总收入分别是123,,x x x 元,根据题意,建立方程
组
1232133
120.20.35000.10.47000.30.4600
x x x x x x x x x --=??
--=??--=? 利用matlab 可以求得 x =
1.0e+003 *
1.25648414985591 1.44812680115274 1.55619596541787
案例3
医院营养师为病人配制的一份菜肴由蔬菜、鱼和肉松组成,这份菜肴需含
1200cal 热量,30g 蛋白质和300mg 维生素c ,已知三种食物每100g 中的有关营养的含量如下表,试求所配菜肴中每种食物的数量。
解:设所配菜肴中蔬菜、鱼和肉松的数量分别为123,,x x x 百克,根据题意,建立方程组
12312312360300600120039630906030300
x x x x x x x x x ++=??
++=?
?++=?
利用matlab 可以求得 x =
1.52173913043478
2.39130434782609 0.65217391304348
矩阵的应用
案例1 矩阵概念的引入 (1)线性方程组
11112211211222221122n n n n n n nn n n
a x a x a x
b a x a x a x b a x a x a x b ++
+=??+++=????++
+=?
的系数(,1,2,
,),(1,2,
,)i j j a i j n b j n ==按原来的位置构成一数表
11121121222212
n n n n nn
n a a a b a
a a
b a a a b ?????????
???
该数表决定着上述方程组是否有解,以及如果有解,解是什么等问题,因而研究这个数表就很重要。
(2)某航空公司在A,B,C,D 四城市之间开辟了若干航线,下图所示表述了四城市间的航
班图,若从A 到B 有航班,则用带箭头的线连接A 和B
。
为了便于研究,表中√为1,空白为0,得到下列数表:
A
B
C
D
(3)某中学学生身高体重的测量,得到如下一份统计如下表
此表反映身高与体重这种关系时也可将上面表格写成一个简化的4行4列的矩形数表,
如果只反映1.5米与体重的关系,则可以用(60 80 70 20);如果只反映60kg 与身高
的关系,则可以用70150805??
?
? ? ???
。 案例5 矩阵概念的应用——逻辑判断问题
甲、乙、丙、丁四人各从图书馆借来一本小说,他们约定读完后互相交换,这四本书的厚度以及他们四人的阅读速度差不多,因此,四人总是同时交换书,经三次交换后,他们四人读完了这四本书,现已知:
(1)乙读的最后一本书是甲读的第二本书;
(2)丙读的第一本书是丁读的最后一本书。 问四人的阅读顺序是怎样的?
解:设甲、乙、丙、丁最后读的书的代号依次为A,B,C,D,则根据题设条件可以列出初
始矩阵
??
? ? ? ???
D B A B C D 下面我们来分析矩阵中各位置的书名代号。已知每个人都读完了所有的书,所以并第二次读的书不可能是C,D 。又甲第二次读的书是B ,所以丙第二次读的书也不可能是B ,从而丙第二次读的书是A ,同理可依次推出丙第三次读的书是B ,丁第二次读的书是C ,丁第三次读的书是A ,丁第一次读的书是B ,乙第二次读的书是D ,甲第一次读的书是C ,乙第一次读的书是A ,乙第三次读的书是C ,甲第三次读的书是D 。故四人阅读的顺序可用矩阵表示如下:
1
2 3 4
甲 乙 丙 丁
C A B
D A D C B A ??
? ?
?
???
D
B C A B C D
案例6 矩阵乘法的应用
某企业某年出口到三个国家的两种货物的数量及两种货物的单位价格、重量、体积如下表所示:
利用矩阵乘法计算该企业出口到三个国家的货物总价值、总重量、总体积各为多少?
解:设矩阵
300014000.50.040.215001300,0.40.060.42000800A B ??
?? ?
== ? ??? ?
??
则矩阵
C AB ==2060 204 11601270 138 8201320 128 720?? ?
? ???
总价值 总重量 总体积
美国 德国
日本
案例7 逆矩阵的应用
一个城市有三个重要的企业:一个煤矿,一个发电厂和一条地方铁路。开采一块钱的煤,煤矿必须支付0.25元的运输费。而生产一块钱的电力,发电厂需支付煤矿0.65元的燃料费,自己亦需支付0.05元的电费来驱动辅助设备及支付0.05元的运输费。而提供一块钱的运输费铁路需支付煤矿0.55元的燃料费,0.10元的电费驱动它的辅助设备。某个星期,煤矿从外面接到50000元煤的订货,发电厂从外面接到25000元电力的订货,外界对地方铁路没有要求。问这三个企业在那一个星期的生产总值各为多少时才能精确地满足它们本身的要求和外界的要求?
解:各企业产出一元钱的产品所需费用为
对于一个星期的周期,设1x 表示煤矿的总产值,2x 表示电厂的总产值,3x 表示铁路的总产值。
煤矿的总消耗为 12300.650.55x x x ++ 电厂的总消耗为 12300.050.10x x x ++ 铁路的总消耗为 1230.250.050x x x ++ 则
1123(00.650.55)50000x x x x -++=
2123(00.050.10)25000x x x x -++= 3123(0.250.050)0x x x x -++=
联立三个方程并整理得方程组
123231230.650.55500000.950.10250000.250.050
x x x x x x x x --=??
-=?
?--+=?
上述方程组可化为 AX b = ,其中
10.650.550
0.950.100.250.051A --?? ?=- ? ?--??, 12350000,250000x X x b x ???? ? ?
== ? ? ? ?????
利用matlab 求解,可知 det()0.7981250A =≠,所以方程组有唯一解,其解为
1123804232858321535x X x A b x -???? ? ?=== ? ? ? ?????
所以煤矿总产值为80423元,发电厂总产值为28583元,铁路总产值为21535元。
案例8 求解线性方程组
(1)假设你是一个建筑师,某小区要建设一栋公寓,现在有一个模块构造计划方
案需要你来设计,根据基本建筑面积每个楼层可以有三种设置户型的方案,如下表所示。如果要设计出含有136套一居室,74套两居室,66套三居室,是否可行?设计方案是否唯一?
解:设公寓的每层采用同一种方案,有1x 层采用方案A ,有2x 层采用方案B ,有3
x 层采用方案C ,根据题意,可得
1231231
238891367437434566
x x x x x x x x x ++=??
++=??++=? 利用matlab 计算方程组的系数矩阵A 、增广矩阵()A A b =的秩:
()()23r A r A ==<,
所以方程组有无穷多个解。
利用matlab 将增广矩阵化为行简化阶梯型矩阵:
1102213
()011580000A A b ??
- ? ?
?=→ ? ? ? ???
矩阵对应的方程组为1323122
3158x x x x ?
=+????=-??
,
取3(x c c =为正整数),则方程组的全部解为
1
2
3122
13158x c x c x c ?
=+??
?
=-??
=???
又由题意可知,123,,x x x 都为正整数,则方程组有唯一解1236,2,8x x x ===。
所以设计方案可行且唯一,设计方案为:6层采用方案A,2层采用方案B,8层采用方案C。
(2)在一个原始部落中,农田耕作记为F,农具及工具的制作记为M,织物的编织记为C。人们之间的贸易是实物交易系统(见下图)。由图中可以看出,农夫将每年的收获留下一半,分别拿出四分之一给工匠和织布者;工匠平均分配他们制作的用具给每个组。织布者则留下四份之一的衣物为自己,四分之一给工匠,二分之一给农夫。
随着社会的发展,实物交易形式需要改为货币交易。假设没有资本和负债,那么如何对每类产品定价才能公正地体现原有的实物交易系统?
也可以用下表表示:
解:令1x 为农作物的价值,2x 为工具的价值,3x 为织物价值。那么从上表第一列,
农夫生产的价值应该等于他们交换到的产品的价值,即
1123111
232
x x x x =
++ 同理可以得到工匠和纺织者产品价值的方程
2213111344x x x x =++, 3312111443
x x x x =++
从而得到下列方程组:
12312
312311102321
2104
341
130434x x x x x x x x x ?--=??
?-+-=???--+=??
利用matlab 将系数矩阵化为行简化阶梯型矩阵,为
A=5103011000?
?-????-????????
令3x c =,写成方程组,为
1253x c
x c
?
=???=? 写成向量形式为
1235311x x c x ??????
????=????
??????????
所以当农作物价值、工具价值与织物价值的定价之比为1235
:::1:13
x x x =时,才能公正地体现原有的实物交易系统。
(3)某药厂生产3种中成药,每件中成药的生产要经过3个车间加工。3个车间
每周的工时、每件中成药在各车间需要的工时数如下表所示,问3中中成药每周的产量各是多少?
解:设3种中成药每周的产量分别为123,,x x x ,则由题意得
123123123240
3237528
x x x x x x x x x ++=??
++=??++=?
利用matlab 将方程组的增广矩阵化为行简化阶梯型矩阵,得
112401007(,)3237501091112800112A b ????
????=→????
????????
由此可以得出 123
7
912x x x =??
=??=?
所以三种中成药每周的产量分别为7件,9件,12件。
案例9解线性方程组应用—人口迁移模型
在生态学、经济学和工程学等许多领域中经常需要对随时间变化的动态系统进行数学建模,此类系统中的某些量常按离散时间间隔来测量,这样就产生了与时间间隔相应的向量序列012,,,
,x x x 其中n x 表示第n 次测量时系统状态的有关信息,而0x 常被称为初始向量。
如果存在矩阵A ,并给定初始向量0x ,使得1021,,
,x Ax x Ax ==即
1(0,1,2,)n n x Ax n +==
则上述方程为一个线性差分方程或者递归方程。
(1)已知某城市2009年的城市人口为5000000人,农村人口为7800000人。
假设每年大约有5%的城市人口迁移到农村(95%仍然留在城市),12%的农村人口迁移到城市(88%仍然留在农村),如下图所示,忽略其他因素对人口规模的影响。计算2011年的人口分布。
解:由题意可得迁移矩阵为
0.950.120.050.88M ??
=??
??
设2009年的初始人口为0x ,2010年和2011年的人口分别为12,x x ,则
100.950.12500000056860000.050.8878000007114000x Mx ??????
===??????
?????? 210.950.12568600062553800.050.8871140006544620x Mx ??????
===????????????
即2011年的人口分布情况是:城市人口为6255380,农村人口为6544620.
(2)在某个地区,每年约有4%的城市人口移居到周围的农村,大约5%的农村人口移居到
城市中。在2009年,城市中有400000居民,农村有600000居民。建立一个差分方程来描述这种情况,用0x 表示2009年的初始人口,然后估计两年之后,即2011年城市和农村的人口数量(忽略其他因素对人口规模的影响)
(3)某公司有一个车队,大约有450辆车,分布在三个地点。一个地点租出去的车可
以归还到三个地点中的任意一个,但租出的车不许当日归还。下面的矩阵给出了汽车归还到
每个地点的不同比率。假设星期一在机场有304辆车,东部办公区有48辆车,西部办公区有98辆车,那么在星期三时,车辆的大致分布式怎么样?
0.970.050.1
0.9
0.050.030.05 0.85车辆出租地
机场东部西部归还到机场
东部
西部
解:设星期一机场、东部和系部的车辆为0x ,星期二和星期三三个地方的车辆分别为
12,x x ,由题意可得,迁移矩阵为
0.970.050.100.90.050.030.050.85M ??
??=??
????
则
100.970.050.130430700.90.0548480.030.050.859895x Mx ??????
??????==≈??????
????????????
210.970.050.130731000.90.0548480.030.050.859592x Mx ??????
??????==≈??????
????????????
所以,星期三时,机场有310辆车,东部办公区有48辆车,系部办公区有92辆车。
案例10 解线性方程组应用—网络流模型
网络流模型广泛应用于交通、运输、通信、电力分配、城市规划、任务分派以及计算机辅助设计等众多领域。当科学家、工程师和经济学家研究某种网络中的流量问题是,线性方程组就自然而然地产生了,例如:城市规划设计人员和交通工程师监控城市道路网络的交通流量,电气工程师计算电路中流经的电流,经济学家分析产品通过批发商和零售商网络从生产者到消费者的分配等。大多数网络流模型中的方程组都包含了数百甚至上千个未知量和线
性方程。
一个网络由一个点集以及连接部分或全部点的直线或弧线构成。网络中的点称作联结点(或节点),网络中的连接线称作分支。每一分支中的流量方向已经指定,并且流量(或流速)已知或者已标为变量。
网络流的基本假设是网络流中流入与流出的总量相等,并且每个联接点流入和流出的总量也相等。例如:下图所示分别说明了流量从一个或两个分支流入联结点,12,x x 和3x 分别表示从其他分支流出的流量,45,x x 表示从其他分支流入的流量。因为流量在每个联结点守恒,所以有1245360,80x x x x x +=+=+。
网络分析要解决的问题就是:在部分信息(如网络的输入量)已知的情况下,确定每一分支中的流量。
下图的网络给出了在下午两点钟,某市区部分单行道的交通流量(以每刻钟通过的汽车数量来度量)。试确定网络的流量模式。
1
x 4x
3x
5x
80
解:根据网络流模型的基本假设,在节点(交叉口)A,B,C,D 处,我们可以得到下列方程:
12234:2030;:30;A x x B x x x +=++=+ 4551:40;
:5010.C x x D x x =++=+
此外,该网络的总流入等于网络的总流出,即
33203050304010,20x x ++=+++?=
联立以上方程的方程组:
12234451531030404020
x x x x x x x x x x -=?
?--=-??
-=??-=??=?
取5x c =,则网络的流量模式表示为
1234540,30,20,40,x c x c x x c x c =+=+==+=
40
50
30
30
20
10
B
A C
D
2x
3x
4x
5x
1x
线性规划问题
案例1、.生产计划问题
(1)假设某厂计划生产甲、乙两种产品,现库存主要原料有A 类3600kg ,B 类
2000kg,C 类3000kg.每件甲产品需用材料A 类9kg,B 类4kg,C 类3kg 。每件乙产品需用材料A 类4kg,B 类5kg,C 类10kg 。甲单位产品的利润70元,乙单位产品的利润120元。问如何安排生产,才能使该厂所获的利润最大。
解:建立模型:
设生产甲、乙产品的件数分别为12,x x ,f 为该厂所获总利润,则
1212121212max 70120.943600
4520003103000
,0
f x x s t x x x x x x x x =++≤+≤+≤≥ 程序如下: 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 =