1、第三章 初等数学模型所谓初等数学模型主要是指建立模型所用的数学知识和方法主要是初等的，而不是高等的。在解决实际问题的过程中，往往主要是是看解决问题的效果和应用的结果如何，而不在于用了初等的方法还是高等的方法，对于数学建模也是这样。本章介绍了量纲分析法、比例与函数建模法，并给出了相应的一些模型。第一节 量纲分析法量纲分析提出于 20 世纪初，是物理学中常用的一种定性分析方法，也是在物理领域中建立数学模型的一个有力工具。它是在经验和实验的基础上, 利用物理定律的量纲齐次原则，确定各物理量之间的关系。1.1 量纲齐次原则许多物理量是有量纲的，有些物理量的量纲是基本的，另一些物理量的量纲则可以由基本量
2、纲根据其定义或某些物理定律推导出来。例如在动力学中，把长度 ， 质量 和时lm间 的量纲作为基本量纲，记为t； TtMmLl,而速度 的量纲可表示为 . fv， 力 21,MLTfLv在国际单位制中，有 7 个基本量：长度、质量、时间、电流、温度、光强度和物质的量，它们的量纲分别为 L、M、T、I、 、J、和 N，称为基本量纲。任一个物理量 q 的量纲都可以表成基本量纲的幂次之积， q量纲齐次性原则 用数学公式表示一个物理定律时，等式两端必须保持量纲一致。量纲分析就是在保证量纲一致的原则下，分析和探求物理量之间关系。先看一个具体的例子，再给出。1.2 量纲分析的一般方法例 1 （单摆运动）质量
3、为 的小球系在长度为 的线的一端，线的另一端固定，小ml球偏离平衡位置后，在重力 作用下做往复摆动，忽略阻力，求摆动周期 的表达式。g t解：在这个问题中有关的物理量有 设它们之间有关系式lt,-（1.1）321lt其中 为待定常数，入为无量纲的比例系数，取（1.1）式的量纲表达式有,整理得： -（1.2）321glmt 321TLM由量纲齐次原则应有 -1203（1.3）解得： 代入（31）得 -,2,1,0321 glt（1.4）（1.4）式与单摆的周期公式是一致的1.3 Buckingham Pi 定理下面我们给出用于量纲分析建模的 Buckingham Pi 定理，定理（Bucking
4、ham Pi 定理） 设 n 个物理量 之间存在一个函数关系nxx,21-（1.5）0,21xf为基本量纲， 。 的量纲可表示为 mx1 1xnixijji ,21 矩阵 称为量纲矩阵，若 的秩 则（1.5）式与下式等价，mnijA)(A,rRak0)21rnF其中 F 为一个未定的函数关系， 为无量纲量，且 可表示为ss-(1.6)(1sinisx而 为线性齐次方程组 的基本解向量. )()(2) sns 0TA利用 Pi 定理建模，关键是确定与该问题相关的几个基本量纲的无量纲量。rn21第二节 航船的阻力长 吃水深度 的船以速度 航行，若不考虑风的影响，航船受到的阻力 除依赖于lhv f船
5、的诸变量 以外，还与水的参数密度 P，粘性系数 , 以及重力加速度 g 有v、 关。我们利用 pi 定理分析 f 和上述物理之间的关系1 航船问题中涉及的物理量及其量纲为2132LTgMvLhlTf我们要寻求的关系式为 -（2.1）0),(gvhlf这些物理量中涉及到的基本量纲为 L、M、T2 写出量纲距阵TLAT 2102133rank3 解齐次线性方程组 , 可得 个基本解向量 TA4ArankTT）、（ ）、（ 、 、 01210)()4(3)21由（1.6）式，可给出 4 个无量量纲-（2.2）124321vflglh由 Pi 定理， （2.1）等价于下列方程 -（2.3）0),(43
6、这里 是未定的函数由（2.2） ， （2.2）可得阻力 f 的显式表达式，-（2.4）),(3212vlf其中 是由（2.3）可得到的未知的函数关系,在力学上, 称为 Froude 数，记为 Fr ;lgv称为 Reynold 数，记为 Re , 因此（2.4）又可写为lv-（2.5）Re),(2Frvlfhl4. 下面我们利用物理模拟进一步确定航船在水中的阻力。 设: 分别表示模型和原型中gvlfghlf 、和、 的各物理量，由(2.5) 有 ),l(2vvlf),gl,(2hlf当无量纲量 -（2.6） vlvll成立时 , 可得 -（2.7）2lvf则此时由模型船的阻力 ，及其它的 ；可
7、确定原型船的阻力f 、 vll。f下面,我们讨论一下(2.6)成立的条件,如果在实验中采用跟实际同样的水质 ,则 又;pg故可得 ： -（2.8）lvlvhl 要使得(2.8)成立 , 必有 ；也即模型船与原型船一样大，这显然排除了hl,物理模拟的可行性。若考虑选用不同的水质， 仍设 则(2.6) 式,ts化为-（2.7）lvlvhl由(2.7)可得 ，23l若按 1：20 的比例， ，显然无法找到如此小的粘性系数的液体。01实际上的一种近似处理方法是，在一定条件下 Re 数的影响很小，这样可近似得到，)lg,(2vhvlf类似地分析，只要 即有 lvl-（2.8）3lf由（2.8）式就容易确
8、定原型船的阻力 f第三节 抛射问题下面我们通过一个例子，介绍如何使用无量纲化方法简化模型。抛射问题：在某星球表面以初速度 竖直向上发射火箭，记星球半径为 ，星球表面vr重力加速度为 ，忽略阻力，讨论发射高度 随时间 T 的变化规律。gx模型建立：设 轴竖直向上， 时 ，火箭和星球质量分别记为 和 ，x0t 1m2由牛顿第二定律和万有引力定律可得: -（3.1）211)(rxmk以 代入（2.17）可得 gx)0(,)( 2grk故得如下初值问题-（3.2）vxyr)0(2(3.2)式的解可以表为 也即发射高度是以 为参数的 的)(grt、 gvr、 t函数,下面我们采用无量纲化方法化简方程（3
9、.2） ， 显然抛射问题中的基本量纲为 ， 而 ， ， ，LTLxTtLr，1LTv2LTg所谓无量纲化是指，对(3.2)式中的 和 分别构造且有相同的参数组合 和 ，使得xt cxt新变量 00tx为无量纲量，其中 ， 称为特征尺度或参考尺度；把方程（3.2）化为 对 cxt x的微分方程，即可简化模型，如何寻找特征尺度，这里我们以 为例，首先写出参数t ct的量纲距阵 gvr、 A210t 的量纲向量为 记为 : T),(0求解线性方程组 得通解: 0ATTk)1,2()0,1(任取 ，即得到一种特征尺度，例如 得 ， 得 k krvtc1k得 同理可得 的几种特征尺度 等21;1vgtc
10、 1rgtc x2,g以下，我们利用不同的 和 化简（3.2） x1. 令 ; 则1,rvtxcc 1,rvtrx由 , tdv2tdr代入(3.2)可得: -（3.3）1)0(,22xrgv(3.3)式的解可表为 ,含一个独立参数且 为无量纲量 .)(t 2. 令 , 类似地可将(3.2)化为 :,rxc1gtc-（3.4）rgvx2,)0(,)3. 令 可将(3.3)化为 ,120gvx;10vt-（3.4）1)0(,22xrgv按照现有科技能力, ,秒米8rgv在（3.4）中令 ，则有 ,10-（3.5）1)0(x(2.22) 的解为: ，tt2)(代回原变量 得: , -（3.6）vt
11、gtx1(3.6)式恰为假定火箭运动过程中所受星球引力不变的运动方程。无量纲化是用数学工具研究物理问题时常用的方法，恰当地选择特征尺度不仅可以减少参数的个数，而且可以帮助人们决定舍弃哪些次要因素。虽然量纲分析法可以得到一些重要有用的结论，但它也有很大的局限性，它要用于对实际问题的定性分析。思考题： 在（3.3）和(3.4)式中,令 会得到什么结果 ?0第四节 比例与函数建模本节介绍的几个模型,都是利用基本的比例关系与函数建立起的数学模型。4.1 动物体型问题问题： 某生猪收购站，需要研究如何根据生猪的体长（不包括头尾）估计其体重？模型假设：1 将四足动物的躯干（不含头尾）视为质量为 m 的圆柱
12、体，长度为 ，截面面积l，直径为 d，图 4-1s图 4-12 把圆柱体的躯干看作一根支撑在四肢上的弹性梁，动物在体重 f 作用下的最大下垂为 ，即梁的最大弯曲，根据弹性力学弯曲度理论有： -（4.1） 23sdfl3 以生物进化学的角度，可认为动物的相对下垂度 已达到一个最合适的数值，l也即 为常数。l模型建立： ； - -（4.2） 42dsldf42由（1）式 可令 为比例常数231sflk1由（2）式 -2224)4(lflds（4.3） flk5144151llf令 由假设 3， 为常数 lk4k-（4.4）lf因此生猪的体重与体长的四次方成正比，在实际工作中，工作人员可由实际经验及
13、统计数据找出常数 K，则可近似地由生猪的体长估计它的体重。4.2 双重玻璃的功效问题： 房间居室的窗户有的是双层的，即在窗户上装两层玻璃，且中间留有一定的空隙，试比较双层玻璃窗与单层玻璃窗的热量流失？模型假设：1。设双层玻璃窗的两玻璃的厚度都为 d，两玻璃的间距为 L；单层玻璃窗的玻璃厚度为 2d，所用玻璃材料相同，如图 4-22假设窗户的封闭性能很好，两层玻璃之间的空气不流动，即忽略热量的对流，只考虑热量的传导.3. 室内温度 和室外温度 保持不变，热传导过程处于稳定状态，即单1T2位时间通过单位面积的热量为常数4. 玻璃材料均匀，热传导系数为常数.图 4-2模型建立：对于厚度为 d 的均匀
14、介质，两侧温度差为 ，则单位时间由温度高的一T侧向温度低的一侧通过单位面积的热量 Q 满足 k 为热传导系数 TkQ设玻璃的热传导系数为 ，空气的热传导系数为 1 2k（1） 先考虑单层玻璃的单位时间，单位面积的热量传导-（4.5）dTk211（2） 考虑双层玻璃情形此时热量先通过厚度为 d 的玻璃传导到两层玻璃的夹层空气中，再通过空气传导，再通过厚度为 d 的玻璃传导；设内层玻璃的外侧温度为 ，外层玻aT璃的内侧温度为 ;则有:bT-（4.6）dTklkdQba21212 由（4.6）式可得 -（4.7） )(221TdlkTbba记 ls1则 -)(222TsTbb（4.8） - -)()
15、( 2212bb（4.9） -)()(212TSTb（4.10）-)()2(121TdkSQ（4.11） 考虑两者之比 -s21（4.12）显然 , 也即双层玻璃的热量损失较小. 12Q模型分析与应用： 常用玻璃的热传导系数 ，csmJk3311084而不流通,干燥空气的热传导系数 ，c425.若取 , 则 ， hdlhS316故 -（4.13）Q812若取 ， 则 又此可见双层玻璃的保暖效果是相当可观的。4h312我国北方寒冷地区的建筑物，通常采用双层玻璃；由（4.13）式 h=4 时, , h 再大，热量传递的减少就不明显了，再考虑到墙体的厚度；所以123Q建筑规范通常要求 .44.3 席
16、位分配模型问题：某校有 200 名学生，甲系 100 名，乙系 60 名，丙系 40 名；若学生会中学生代表有 20 个席位，则公平又简单的分法应各有 10，6，4 个席位。若丙系有 6 名学生分别转入甲、乙两系各 3 人，此时各系的人数为 103，63，34；按比例席位分配应为 10.3，6.3 和3.4，出现了小数，19 个整数席位分配完后，最后一席留给小数部分最大的丙系，分别为10，6，4。为方便提案表决，现增加 1 席共 21 席，按比例计算甲、乙、丙三系分别占有10.815，6.615，3.570；按上面的分法应分别为 11，7，3；这样虽然增加了一个席位，但丙系的席位反而减少一席，
17、因此这种分法显然是不合理的，请给出一个比较公平的席位分配方案？问题分析： 席位分配问题，当出现小数时，无论如何分配都是不完全公平的。那么一个比较公平的分法是应该找到一个不公平程度最低的方法，因此首先要给出不公平程度的数量化，然后考虑使之最小的分配方案。模型建立：一、讨论不公平程度的数量化设 A，B 两方人数分别为 ；分别占有 和 个席位，则两方每个席位所代21,p1n表的人数分别为 和 。1np2我们称 为绝对不公平值。例： ， ， ；21np120p210,21n则 又 ， ， 则 ；2110221n21p由上例可知，用绝对不公平程度作为衡量不公平的标准，并不合理，下面我们给出相对不公平值。
18、若 ，则称 为对 A 的相对不公平值， 记为 ；21np1221np ),(21nrA若 则称 为对 B 的相对不公平值 ，记为 ；21np12112np ),(21nrB上例中，相对不公平值分别为：0.2 和 0.02，可见相对不公平值较合理。二、 下面我们用相对不公平值建立模型，设，A，B 两方人数分别为 ；分别占有 和 个席位现在增加一个席位，21,p1n应该给 A 还是 B ？不妨设 ，此时对 A 不公平，下面分二种情形21n（1） ，这说明即使 A 增加 1 席，仍对 A 不公平, 故这一席应给 A。21np（2） ， 说明 A 方增加 1 席时，将对 B 不公平，此时计算对 B 的
19、相对21不公平值)(),(2121npnrB若这一席给 B，则对 A 的相对不公平值为)(),(1221rA本着使得相对不公平值尽量小的原则，若 -（4.14）)1,(),1(22nrnrAB则增加的 1 席给 A 方，若 -（4.15）),(),(212B则增加的 1 席给 B 方由（4.14）式可得 ： )()(122npn由（4.15）式可得 ： )()(122记 ： 则增加的 1 席，应给 Q 值大的一方 )1(iinpQ第一种情形，显然也符合该原则，现在将上述方法推广到 方分配席位的情况 方人数为 已占有 席 miAipini,2,1计算 则将增加的 1 席分配应给 值最大的一方)1
20、(2iinpQ下面考虑原问题：前 19 席的分配没有争议，甲系得 10 席，乙系得 6 席，丙系得 3 席第 20 席的分配4.96)10(321Q5.)(23.96)1(43Q故第 20 席分配给甲系第 21 席的分配4.80)1(3213.96,5.942QQ故第 21 席分配给乙系甲、乙、丙三系各分得 11，6，4 席，这样丙系保住它险些丧失的 1 席。思考题：1 若单层玻璃窗的玻璃厚度也是 ，结果将如何？d2 怎样讨论三层玻璃的功效？3 怎样讨论双层玻璃的隔音效果？4 比利时 分配方案：)(Hondt将甲、乙、丙三系的人数都用 去除，将商从大到小排列，取前 21 个最,321大的，这 21 个中各系占有几个，就分给几个席位，你认为这种方法合理吗？5 学校共有 1000 名学生，235 人住在 A 楼，333 人住在 B 楼，432 人住在 C 楼。学生们要组成一个 10 人委员会，使用 值方法及 方法给出分配方案。如QHondt果委员会为 15 人，分配方案是什么？
