数学模型第三版课后习题集答案解析

《数学模型》作业解答

第二章(1)（2008年9月16日）

1． 学校共1000名学生，235人住在A 宿舍，333人住在B 宿舍，432人住在C 宿舍.学生们

要组织一个10人的委员会，试用下列办法分配各宿舍的委员数：

（1）. 按比例分配取整数的名额后,剩下的名额按惯例分给小数部分较大者;

（2）. §1中的Q 值方法；

（3）.d ’Hondt 方法：将A 、B 、C 各宿舍的人数用正整数n=1,2,3,……相除，其商数如下表：

将所得商数从大到小取前10个（10为席位数），在数字下标以横线，表中A 、B 、C 行有横线的数分别为2，3，5，这就是3个宿舍分配的席位.你能解释这种方法的道理吗？

如果委员会从10个人增至15人，用以上3种方法再分配名额，将3种方法两次分

配的结果列表比较.

解：先考虑N=10的分配方案，

,432 ,333 ,235321===p p p

∑==31.1000i i p 方法一（按比例分配） ,35.23111==∑=i i p

N

p q ,33.33122==∑=i i p N p q 32.43133==∑=i i

p N p q 分配结果为： 4 ,3 ,3321===n n n

方法二（Q值方法）

9个席位的分配结果（可用按比例分配）为：

4 ,3 ,2321===n n n

第10个席位：计算Q 值为

,17.92043223521=?=Q ,75.92404333322=?=Q 2.93315

44322

3=?=Q 3Q 最大，第10个席位应给C.分配结果为 5 ,3 ,2321===n n n

方法三（d ’Hondt 方法）

此方法的分配结果为：5 ,3 ,2321===n n n

此方法的道理是：记i p 和i n 为各宿舍的人数和席位（i=1,2,3代表A 、B 、C 宿舍）.i i n p 是每席位代表的人数，取,,2,1 =i n 从而得到的

i i n p 中选较大者，可使对所有的,i i

i n p 尽量接近. 再考虑15=N 的分配方案，类似地可得名额分配结果.现将3种方法两次分配的结果列表如下：

2． 试用微积分方法，建立录像带记数器读数n 与转过时间的数学模型.

解： 设录像带记数器读数为n 时，录像带转过时间为t.其模型的假设见课本.

考虑t 到t t ?+时间内录像带缠绕在右轮盘上的长度，可得,2)(kdn wkn r vdt π+=两边积分，得 ??+=n

t dn wkn r k vdt 00)(2π )22 2

n wk k(r n πvt +=∴ .2 22

n v

k w n v rk t ππ+=∴ 第二章(2)（2008年10月9日）

15．速度为v 的风吹在迎风面积为s 的风车上，空气密度是ρ ，用量纲分析方法确定风车获得的功率P 与v 、S 、ρ的关系.

解: 设P 、v 、S 、ρ的关系为0),,,(=ρs v P f ， 其量纲表达式为:

[P]=32-T ML , [v ]=1-LT ,[s ]=2L ,[ρ]=3-ML ,这里T M L ,,是基本量纲.

量纲矩阵为：

A=)

??????????---ρ()()

()()()()(001

310013212s v P T M L

齐次线性方程组为： ??

???=--=+=-++030

032221414321y y y y y y y y 它的基本解为)1,1,3,1(-=y

由量纲i P 定理得 1131ρπs v P -=， 113ρλs v P =∴ ， 其中λ是无量纲常数.

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

解：设v ,

ρ,μ,g 的关系为(f v ,ρ,μ,g )=0.其量纲表达式为[v ]=LM 0T -1，[ρ]=L -3MT 0，[μ]=MLT -2（LT -1L -1）-1L -2=MLL -2T -2T=L -1MT -1，[g ]=LM 0T -2,其中L ，M ，T 是基本量纲.

量纲矩阵为

A=)

()()()()

()()(210101101131g v T M L μρ??????????----- 齐次线性方程组Ay=0 ，即

?????==+=+02y -y - y -0

y y 0y y -3y -y 431

324321 的基本解为y=(-3 ,-1 ,1 ,1)

由量纲i P 定理 得 g v μρπ13--=. 3ρ

μλg v =∴，其中λ是无量纲常数.

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

解：设v ,ρ,μ,γ，g 的关系为0),,,,(=g v f μργ.其量纲表达式为

[v ]=LM 0T -1，[ρ]=L -3MT 0，[μ]=MLT -2（LT -1L -1）-1L -2=MLL -2T -2T=L -1MT -1，[γ]=LM 0T 0 ，[g ]=LM 0T -2

其中L ，M ，T 是基本量纲.

量纲矩阵为

A=)()()()()

()()()(21001

110011311g v T M L μργ??????????----- 齐次线性方程组Ay=0 即 ??

???=---=+=+--+020035414354321y y y y y y y y y y 的基本解为

??

???---=--=)21,1,1,23,0()21,0,0,21,1(21y y 得到两个相互独立的无量纲量

???==-----2/112/32

2

/12/11g g v μργπγπ 即 1212/12/31,--==πμργπγg g v . 由0),(21=Φππ , 得 )(12

1-=π?π ∴ )(12/12/3-=

μργ?γυg g , 其中?是未定函数.

20.考察阻尼摆的周期，即在单摆运动中考虑阻力，并设阻力与摆的速度成正比.给出周期的表达式，然后讨论物理模拟的比例模型，即怎样由模型摆的周期计算原型摆的周期. 解：设阻尼摆周期t ,摆长l , 质量m ,重力加速度g ，阻力系数k 的关系为

0),,,,(=k g m l t f

其量纲表达式为：

1

12120000000)(]][[][,][,][,][,][-----======LT MLT v f k T LM g MT L m T LM l T M L t

10-=MT L ， 其中L ，M ，T 是基本量纲.

量纲矩阵为

A=)

()()()()()()()(120011010001010k g m l t T M L ????

??????-- 齐次线性方程组

?????=--=+=+0200541

5342y y y y y y y 的基本解为

??

???--=-=)1,21,1,21,0()0,21,0,21,1(21Y Y 得到两个相互独立的无量纲量

∴g l t =1π, )(21π?π=, 2/12

/12mg

kl =π ∴)(2/12

/1mg kl g l t ?=

，其中?是未定函数 . 考虑物理模拟的比例模型，设g 和k 不变，记模型和原型摆的周期、摆长、质量分别为

t ,'t ；l ,'l ;m ,'m . 又)(2

/12

/1g m l k g l t '''='? 当无量纲量l l m m '='时， 就有 l l l g g l t

t '=?'='. 《数学模型》作业解答

第三章1（2008年10月14日）

???==---2

2/112/112/12/1ππk g m l g tl

1.在3.1节存贮模型的总费用中增加购买货物本身的费用,重新确定最优订货周期和订货

批量．证明在不允许缺货模型中结果与原来的一样，而在允许缺货模型中最优订货周期和订货批量都比原来结果减少．

解：设购买单位重量货物的费用为k ,其它假设及符号约定同课本．

01 对于不允许缺货模型，每天平均费用为：

kr rT c T c T C ++=2

)(21 2221

r c T

c dT dC +-= 令0=dT

dC ， 解得 r c c T 21*2= 由rT Q = ， 得212c r c rT Q =

=** 与不 …… 此处隐藏：4916字，全部文档内容请下载后查看。喜欢就下载吧 ……