数学模型课后详细复习资料

数学模型作业

六道题

作业一

1.P56.8一垂钓俱乐部鼓励垂钓者将钓上的鱼放生，打算按照放生的鱼的重量给予奖励，俱乐部只准备了一把软尺用于测量，请你设计按照测量的长度估计鱼的重量的方法。假定鱼池中只有一种鲈鱼，并且得到8条鱼的如下数据（胸围指鱼身的最大周长）：

解：

要求鱼的体重，我们利用质量计算公式：M=ρV。我们假定鱼池中是同一种鱼，于是可以近似地考虑其密度是相同的。至于鱼的体积问题，由于是同一种类，可以假定这种鱼在体型上是一致的。我们假设鱼的体积和鱼身长的立方成正比。即：V=k

1

L3，因此，模型为：

33

111

M V k l K L

ρρ

===……………………………

模型一

利用Eviews软件，用最小二乘法估计模型中的参数K

1

，如下图1所示：

图1

从图1结果可以得到参数K

1

=0.014591，所以模型为：

3

1

M0.014591 L

=

上述模型存在缺陷，因为它把肥鱼和瘦鱼同等看待。因此，有必要改进模型。如果只假定鱼的横截面是相似的，假设横截面积与鱼身最大周长的平方成正比，

即：V=k

2

d2L，因此，模型为：

身长

/cm

36.8 31.8 43.8 36.8 32.1 45.1 35.9 32.1

质量

/g

765 482 1162 737 482 1389 652 454

胸围

/cm

24.8 21.3 27.9 24.8 21.6 31.8 22.9 21.6

1

2

5

3 4 7

6

22222M V k d K d L L ρρ===………………………………

模型二

利用Eviews 软件，用最小二乘法估计模型中的参数K 2，如下图2所示：

图2

从图2可以得到参数K 2=0. 032248，所以模型为：

22M 0.032248d L

=

将实际数据与模型结果比较如表1所示：

实际数

据M

765 482 1162 737 482 1389 652 454

模型一M 1 727.165 469.214 1226.061 727.165 482.629 1338.502 675.108 482.619 模型二M 2 729.877 465.248 1099.465 729.877 482.960 1470.719 607.106 483.960

2.P131.2 一家出版社准备在某市建立两个销售代理点，向7个区的大学生售书，每个区的大学生数量（单位：千人）已经表示在图上。每个销售代理点只能向本区和一个相邻区的大学生售书，这两个代理点应该建在何处，才能使所能供应的大学生的数量最大？建立该问题的整数线性规划模型并求解。

解：

将大学生数量为34、29、42、21、56、18、71的区分别标号为1、2、3、4、5、6、7区，画出如下区域区之间的相邻关系：

记r 为第i 区的大学生人数，用0-1变量x ij =1表示（i ，j ）区的大学生由

一个代售点供应图书（i

i j ij

i.j

Max r r x s.t.2

1,{0,1}

i j i j ij

j

j

i j x x x

i

x =+≤+≤?∈∑∑∑∑相邻

（）

即：

12132325344546566747121323242534454647566712131223242513233424455646Max 63*x 76*x 71*x 85*x 63*x 77*x 39x *x 74*x 89*x 92*x s.t.x x x x x x x x x x x 2 x x 1

x x x x 1 x x x 1 x x x 1 x =+++++++++++++++++++≤+≤+++≤++≤++≤5667ij ij x x 1 x 0x 1

++≤==或将上述建立的模型输入LINGO ，如下： modle:

max=63*x12+76*x13+71*x23+85*x25+63*x34+77*x45+39x*x46+74*x56+89*x67+92*x47 s.t. x12+x13+x23+x24+x25+x34+x45+x46+x47+x56+x67<=2; x12+x13<=1;

x12+x23+x24+x25<=1; x13+x23+x34<=1;

x24+x45+x56<=1; x46+x56+x67<=1

@gin(x12); @gin(x13); @gin(x23); @gin(x25); @gin(x34); @gin(x45); @gin(x46);@gin(x47); @gin(x67); End 运行，得到的输出如下：

Local optirnal solution found at iteration Objective value: Vauable Value Reduced Cost

x12 0.000000 0000000 x13 0.000000 0000000 x23 0.000000 0000000 x24 0.000000 0000000 x25 1.000000 0000000 x34 0.000000 0000000 x45 0.000000 0000000 x46 0.000000 0000000 x47 1.000000 0000000 x56 0.000000 0.000000 x67 0.000000 0000000

从上述结果可以得到：最优解 2547x x 1==(其他的均为0)，最优值为177人. 即：第2、5区的大学生由一个销售代理点供应图书，代理点在2区或者5区，第4、7区区的大学生由另一个销售代理点供应图书，代理点在4区或者7区。

作业二

3.P181.14 在鱼塘中投放n 0尾鱼苗，随着时间的增长，尾数将减少而每尾的重量将增加。

(1)设尾数n(t) 的(相对)减少率为常数;由于喂养引起的每尾鱼重量的增加率与鱼表面积成正比，由于消耗引起的每尾鱼重量的减少率与重量本身成正比。分别建立尾数和每尾鱼重的微分方程，并求解。

(2)用控制网眼的办法不捕小鱼，到时刻T 才开始捕捞，捕捞能力用尾数的相对减少量|?/n| 表示，记作E ，即单位时间捕获量是En(t)。问如何选择T 和E ，使从T 开始的捕获量最大。 解： （1）

鱼塘的初始鱼苗为n 0尾，且随着时间的增长，尾数将减少。设尾数n(t) 的(相对)减少率为为k 1，因此由题意建立微分方程为：

,(0)(0)dn

kn k dt n n =->= 求解得：

0()kt n t n e -=

在鱼塘里，由于喂养引起的每尾鱼重量的增加率与鱼表面积成正比，即：

S αI(t)=

在鱼塘里，由于消耗引起的每尾鱼重量的减少率与重量本身成正比，即：

m βD(t)=

所以每尾鱼重量的净增长率r(t)为：

S m αβ-r(t)=

因此，建立微分方程为：

dm

S m dt

αβ-= 因为该微分方程涉及多个变量间的数量关系，所以我们暂时无法求解该微分方程。但是要想解决此微分方程还需要更多的信息，例如，每尾鱼表面积与其重量间的关系，一旦此关系确定，便可轻松解出每尾鱼的质量随时间的变化，即m(t)。

（2）

用控制网眼的办法不捕小鱼，假设t=T 时开始捕捞，且单位时间的捕捞率为E ，依题意建立微分方程：

,()dn

kn En t T dt

=--≥ 因此得:

()()0()t E t T n t n e e λλ--+-=

所以单位时间的捕捞鱼的尾数为En(t)，因此从T 时刻开始的总捕捞量为：

()()T

y m t En t dt ∞

=?

问题就转化为求E 和λ的值，使得y 最大，由于条件不足导致m(t)求解不出，因此无法求出y 的具体解释式。

4.P213.2 雨滴的速度v 与空气密度ρ、粘滞系数μ和重力加速度g 有关，其中粘滞系数的定义是：运动物体在空气中受的摩擦力与速度梯度和接触面积的乘积成正比，比例系数为粘滞系数，用量纲分析方法给出速度v 的表达式。 解：

雨滴速度问题中涉及的物理量：雨滴的速度v ，空气密度ρ，粘滞系数μ,重力加速度g ，长度γ。要寻找的关系是：

(,,,)v g ψγρμ=

更一般的将各个物理量之间的关系写作：

0),,,,(=g v f μργ

这里没有因变量与自变量之分，进而设：

3

5

1

2

4

.................(1)y y y y y v g πγρμ=

其量纲表达式为：

0130110002[]LM T []L MT ,[]L MT []LM T []LM T g νρμγ-----=====,,，

其中L ，M ，T 是基本量纲。 因此量纲表达式可以写成:

5120000130110002(LM T )(L MT )(L MT (LM T (LM T )y y y L M T -----=34y y ))

根据量纲原则可写成：

???

?

?=---=+=+--+020035414354321y y y y y y y y y y

量纲矩阵为：

11311()001

1

()10012()()()()

()()L A M T v g γρμ--????=????---??

解得方程的基本解为：

1211(1,,0,0,)22......................(2)31(0,,1,1,)22Y Y ?

=--???

?=---??

将（2）代入（1）可得两个相互独立的无量纲量

???==-----2/11

2/322

/12/11g g v μργ

πγπ 为了得到形如(,,,)v g ?γρμ=的关系，取12()πψπ=，其中

ψ是某个函数，

所以（2）式为：

1/21/23/211/2()v g g γψγρμ-----=

于是：

3/211/21/21/2()()v g g ψγρμγ---=

作业三

5.P248.13 一个岛屿上栖居着食肉爬行动物和哺乳动物，又长着茂盛的植物。爬行动物以哺乳动物为食物，哺乳动物又依赖植物生存。在适当假设下建立三者关系的模型，求其平衡点。 解：

)(1t x 、)(2t x 、)(3t x 分别表示植物、哺乳动物、食肉爬行动物在时刻t 的数量。假设不考虑植物、哺乳动物和食肉爬行动物对自身的阻滞增长作用。

设1r 为植物的固有增长率，而哺乳动物的存在使植物的增长率减少，设减小的程度与捕食者数量成正比，于是建立植物数量的模型：

)()

(21111x r x dt

t dx λ-= 比例系数1λ反映了哺乳动物消耗植物的能力。

哺乳动物离开植物无法生存，设其死亡率为2r ，则哺乳动物独自存在时有：

222)

(x r dt

t dx -= 而植物的存在可以为哺乳动物提供食物，但是食肉爬行动物的存在使哺乳动物数量减少，设减少的程度与食肉爬行动物数量成正比，于是建立哺乳动物数量模型：

)()

(312222x x r x dt

t dx μλ-+-=

其中比例系数2λ反映了植物对哺乳动物的供养能力，μ反映了食肉爬行动物掠取哺乳动物的能力。

食肉爬行动物离开动物无法生存，设其死亡率为3r ，则食肉爬行动物独自存在时有：

333)

(x r dt

t dx -= 而哺乳动物的存在可以为食肉爬行动物提供食物，于是(4)式右端应加上哺乳动物对食肉爬行动物的增长作用，设为3λ，于是建立食肉爬行动物的数量模型： )()

(23333x r x dt

t dx λ+-= 比例系数3λ反映了哺乳动物对食肉爬行动物的供养能力。

综上所述，建立如下微分方程组模型

????

??

???+-=-+-=-=)()

()()

()()

(2333331222221111x r x dt t dx x x r x dt t dx x r x dt t dx λμλλ

求得微分方程组的平衡点为

)0,,(

),0,0,0(1

1

22

21λλr r P P 其中平衡解)0,0,0(1P 对是没有意义的。

6.P43

7.9 一个服务网络由k 个工作站v1，v2，…，vk 依次串接而成，当某种服务请求到达工作站vi 时， vi 能够处理的概率为 pi ，转往下一站vi+1处理的概率为 qi （i=1,2,... ,k-1,设 qk=0），拒绝处理的概率为 ri ，满足pi + qi + ri =1。试构造马氏链模型，确定到达 vi 的请求平均经过多少工作站才能获得接受或拒绝处理的结果，被接受和拒绝的概率各多大。 解：

用随机变量i X 表示第i 站对请求服务的处理方式，i X =1表示接受请求

i X =2表示拒绝请求，()i ,,k =12L ，用()a i 1表示第i 站接受请求的概率，()a i 2表

示第i 站拒绝请求的概率。i q 表示第i 站转移至下一站的转移概率。分析可知，第i +1站处理请求的概率和第i 站处理请求的概率以及转移概率有关，由此可得

()()i i i i i i a i p q q q a i r q q q +-+-+=???

+=??L L 1111

2111

11 其中，()a p =111 ，()a r =211，由（1）可以计算出k 个站各自接受和拒绝服务请

()i k k k a X p p q p q q p q q q --==++++1213211211L L

服务请求被拒绝处理的概率为

()i k k k a X p p q p q q r q q q --==++++121321121

2L L

将服务请求到达工作站i v 记做状态i ，i ,,k =12L ，设v 0表示请求被拒绝，v 00

表示请求被接受，于是转移概率矩阵为:

()()()

()()()

()()()()()()()()k k k k k p r q p r p p r q k p r k k k ---??

??????????

=???

???

-??

??-????

O O M M O

O

O

M L 111221111

00001000102010121000 转移矩阵i p =1的状态i 称为吸收状态，如果马氏链中至少含有一个吸收状态，并且从每一个非吸收状态出发，能以正的概率经有限次转移到达某一个吸收

状态，那么这个马氏链称为吸收链。吸收链可以写成简单的标准形式，若有r 个吸收状态，k r -个非吸收状态，则转移矩阵p 可表示为

r r I p R Q ???

=?

???

0 其中k r -阶子方阵Q 的特征值λ满足λ<1，这要求子阵()k r r R -?中必含有非零元素，以满足任一非吸收状态经有限次转移能够到达某一个吸收状态的条件，这样Q 就不是随机矩阵，它至少存在一个小于1的行和，且如下定理成。由于

()I Q -可逆因此：

()s s M I Q Q ∞

-==-=∑1

记元素全为1的列向量()T

e ,,=111L ,则y Me =的第i 个分量是从第i 个非吸收状态出发，被某个吸收状态吸收的平均转移次数。

所以有

()k k q q q q q q q q M I Q ---???

????

?=-=?

???

???

?

121112121

1101L L

L O

M O

M

到达i v 的请求获得接受或拒绝时，平均经过的工作站数由y Me =可得：

k i q q q q q q -=++++1121211L L