《数学建模简明教程》课件第2章.ppt

1、1 1第二章 用初等数学方法建模2.1 比例与函数比例与函数2.2 关于自然数的奇偶性关于自然数的奇偶性2.3 量纲分析法量纲分析法2 2 本章介绍的模型比较简单，只要具有比例、函数、奇偶性、状态转移等数学知识就可以构造和求解模型，我们称其为初等数学方法建模.应当指出，解决实际问题时应尽可能地用简单的数学模型，着重于问题的解决，而不在于采用了多么高深的数学方法.从这个意义上说，培养良好的数学思维能力往往比学习更多更深的知识更为有用.3 3本节给出利用比例和函数建立数学模型的例子.我们将会看到，在日常生活中，到处都会遇到应用数学方法来解决的问题.2.1 比例与函数比例与函数4 42.1.1 四足

2、动物的身长和体重关系问题四足动物躯干(不包括头尾)的长度和它的体重有什么关系？这个问题有一定的实际意义.比如，生猪收购站的人员或养猪专业户如果能从生猪的身长估计它的重量，则可以给他们带来很大的方便.四足动物的生理构造因种类不同而异，如果陷入生物学对复杂的生理结构的研究，将很难得到什么有价值的模型.为此，我们可以在较粗浅的假设的基础上，建立动物的身长和体重的比例关系.本问题与体积和力学有关，搜集与此有关的资料得到弹性力学中两端固定的弹性梁的一个结果：5 5长度为L的圆柱型弹性梁在自身重力f作用下，弹性梁的最大弯曲v与重力f和梁的长度立方成正比，与梁的截面面积s和梁的直径d的平方成反比，即32fL

3、vsd6 6利用这个结果，我们采用类比的方法给出如下假设：(1)四足动物的躯干(不包括头尾)的长度为L，断面直径为d的圆柱体，体积为m；(2)四足动物的躯干(不包括头尾)重量与其体重相同，记为f；(3)四足动物可看做一根支撑在四肢上的弹性梁，其腰部的最大下垂对应弹性梁的最大弯曲，记为v.根据弹性理论结果及重量与体积成正比的关系，有：fm，msL7 7由正比关系的传递性，得(2.1.1)式(2.1.1)中多了一个变量v，为替代变量v，注意到是动物躯干的相对下垂度，从生物进化的观点，讨论相对下垂度有：太大，四肢将无法支撑，此种动物必被淘汰；太小，四肢的材料和尺寸超过了支撑躯体的需要，无疑是一种浪费

4、，也不符合进化理论.443422sLLvLvsddLdvLvLvL8 8因此从生物学的角度可以确定，对于每一种生存下来的动物，经过长期进化后，相对下垂度 已经达到其最合适的数值，应该接近一个常数(当然，不同种类的动物，此常数值不同).于是可以得出d2L3，再由fsL和sd2得fL4，由此得到四足动物体重与躯干长度的关系fkL4(2.1.2)式(2.1.2)就是本问题的数学模型.vL9 9如果对于某一种四足动物，比如生猪，可以根据统计数据确定公式中的比例常数k，那么就可得到用该类动物的躯体长度估计其体重的公式.发挥想象力，利用类比方法，对问题进行大胆的假设和简化是数学建模的一个重要方法.不过，使

5、用此方法时要注意对所得数学模型进行检验.此外，从一系列的比例关系着手推导模型可以使推导过程大为简化.10 102.1.2 公平席位分配问题公平席位分配问题学校学生会为了协调各系的工作，需给各系分配学生会成员名额.如何分配才算合理呢?例如某校甲系100名学生，乙系60名学生，丙系40名学生，学生会成员共设20个名额，显然公平而又简单的分配方案是按学生人数比例分配.那么甲、乙、丙三个系分别应有10、6、4个学生会成员.倘若丙系有6名学生转入其他两系学习，甲系成为103名学生，乙系成为63名学生，丙系成为34名学生，仍按比例分配，就会出现小数，但成员数必须是整数，一个自然的想法就是“四舍五入”，即“

6、去掉尾数取整”.而这样的话，常常导致名额多余或不够分配，更严重的是，这种似乎公平的分配方法有时会出现不公平的结果.表2-1和表2-2分别是学生会成员为20个名额和21个名额时的分配表.11 11表 2-112 12 表 2-213 13从上述两表明显可以看到，当学生会成员名额增加1个时，丙系反而减少了1个名额，这个结果对丙系来讲太不公平了.这说明，传统分配方法存在严重缺陷，要解决这一问题，必须寻找公平分配的衡量指标，建立新的分配办法.14 14我们先就A，B两方公平分配席位情形加以说明.设A，B两方人数为p1，p2，占有席位分别为n1，n2，则表示两方每个席位所代表的人数.显然，当且仅当时，名

7、额分配才是公平的.但是一般说来，它们不会相等，这表明席位分配不公平，直观的想法是用数值来表示双方的不公平程度，称为绝对不公平值，但绝对不公平值往往无法区分不公平程度.1212ppnn和1212ppnn1212ppnn15 15例如当p1120，p2100，n1n210时，有：当双方人数增加为p11020，p21000，n1，n2仍为10时，有：121212 102ppnn1212102 1002ppnn16 16 常识告诉我们，尽管它们有着相同的绝对不公平值，但前一种不公平程度明显大于后一种不公平程度，所以绝对不公平值不是好的衡量标准.由此我们想到要用相对标准，仍设p1，p2为A，B两方的人数

8、；n1，n2为双方分配的席位，我们来定义“相对不公平”的概念.若，则称：1212ppnn12121222(,)Appnnr n npn17 17为对A的相对不公平值.若，则称：为对B的相对不公平值.1212ppnn21211211(,)Bppnnr n npn18 18现在我们用确立的衡量分配不公平程度的数量指标rA，rB来制定席位的分配方案.当席位总名额增加1个时，应该给A还是B呢?不失一般性，设，这时对A不公平，当再增加1个席位名额时，有下面三种情形：(1)表明即使A再增加1个名额，仍对A不公平，当然这一名额应给A；1212ppnn12121ppnn19 19(2)表明A增加1个名额就对B

9、不公平，这时对B的相对不公平值为：(3)，即给B增加1个名额时将对A不公平.这时对A的相对不公平值为：12121ppnn211212(1)(1,)1Bp nr nnp n1212(1)ppnn+121221(1)(,1)1Ap nr n np n+=-2020 公平分配席位的原则是使相对不公平值尽可能地小，所以若rB(n11，n2)rA(n1，n21)则增加的名额应给A.反之，则应该给B.注意到rB(n11，n2)rA(n1，n21)等价于：(2.1.3)22212211(1)(1)ppn nn npm，这时由供给曲线S上的点A1决定商品数量为q1.由于q1qm，按需求曲线D，消费者认为价格下

10、至p1才能购买.一旦价格降到p1，生产者按供给曲线S将产量减少成q2，以减少损失.由于q2pm又刺激生产者大量生产，使产量上升为q3，这一过程如此下去，即A1A2A3A4M，达到市场平衡.如果供给曲线S和需求曲线D如图2-3所示，类似分析可知，市场供求关系将按A1A2A3A4M的方向远离平衡点M.3030为什么会出现供求关系有时趋于平衡点，有时偏离平衡点，这是由什么因素来决定的呢？从图2-2和图2-3可以看出，当需求曲线D的斜率绝对值小于供给曲线S的斜率值时，市场趋于平衡，M点是稳定平衡点；反之，当需求曲线D的斜率绝对值大于供给曲线S的斜率值时，市场不稳定，M点是不稳定平衡点.现进一步分析不稳

11、定的原因.两条曲线的切线斜率K实际上表示商品价格随商品数量变化而变化的程度.|KD|KS，表明消费者对商品价格的敏感度比生产者要高，商品稍少一点，人们便去抢购，导致商品价格很快发生大的变化，因此KD较大，容易引起供求的不平稳.31 31图 2-23232图 2-33333如何解决这一问题，即怎样才能总有|KD|KS成立呢？一种办法是控制物价.比如价格有微小的改变，则KD0；另一种办法是控制市场上商品的数量，不能缺少或积压商品，这样，KS.总之，上述两种办法都可以使|KD|KS，从而保证了市场的供需平稳.3434数的奇偶性即数是奇数还是偶数的性质.自然数的奇偶性在数学建模中有着广泛的应用.下面我

12、们介绍这方面的几个例子2.2 关于自然数的奇偶性关于自然数的奇偶性35352.2.1 铺瓷砖问题铺瓷砖问题要用40块正方形瓷砖铺设如图2-4所示图形的地面，但商店里只有长方形瓷砖，每块大小等于正方形的两块.一人买了20块长方形瓷砖，结果无论怎样摆弄，都无法完整铺好.问题在于用20块长方形瓷砖正好铺成图2-4所示地面的可能性是否存在，只有可能性存在才谈得上用什么方法铺设的问题.为此，在图2-4上黑、白相间染色，我们发现共有19个白格和21个黑格，一块长方形瓷砖可盖住一白一黑两格，所以铺上19块长方形瓷砖后，不管用什么方式总要剩下2个黑格没有铺上，而一块长方形瓷砖是无法盖住如图所示的两个黑格的，唯

13、一的方法是把最后一块长方形瓷砖一分为二.3636图 2-43737解决铺瓷砖问题中所用的方法在数学上称为奇偶校验，即如果两个数都是奇数或偶数，则称其具有相同的奇偶性.如果一个数是奇数，另一个数是偶数，则称其具有相反的奇偶性.在铺瓷砖问题中，同色的两个格子具有相同的奇偶性，异色的两个格子具有相反的奇偶性，长方形瓷砖显然只能覆盖具有相反奇偶性的一对方格.因此，把19块长方形瓷砖在地面上铺好后，只有在剩下的两个方格具有相反的奇偶性时，才有可能把最后一块长方形瓷砖铺上.现在由于剩下的两个方格具有相同的奇偶性，所以无法铺上最后一块长方形瓷砖.这就从理论上证明了用20块长方形瓷砖铺好如图2-4所示图形是不

14、可能的，任何改变铺设方法的努力都是徒劳的.3838奇偶校验法巧妙简单，富有创造力在估计事情不可能成立时，可考虑使用奇偶性这一方法来论证.39392.2.2 菱形十二面体上的菱形十二面体上的H路径问题路径问题沿一菱形十二面体各棱行走，要寻找一条这样的路径：经过各顶点恰好一次，这个问题被称为Hamilton路径问题.如图2-5所示.我们利用奇偶校验法证明：在菱形十二面体上没有Hamilton路径.事实上，菱形十二面体每个顶点的度或者是3，或者是4.(所谓顶点的度，是指通过这一顶点的棱数.)另外，仔细观察，每个3度顶点被4度顶点所包围，即3度顶点的三条棱中的每一条棱都与一个4度顶点关联；反之，每个4

15、度顶点也都被3度顶点包围.因此，一条Hamilton路径要通过所有的顶点必须依次轮流通过3度顶点和4度顶点.4040图 2-541 41如果这个多面体的顶的个数是奇数，则3度顶点和4度顶点必须只能差一个，若是偶数，则必须恰好相等.但菱形12面体的顶点是14个，4度顶点只有6个，所以Hamilton路径是不可能存在的.42422.2.3 自然数的因子个数与狱吏问题自然数的因子个数与狱吏问题将自然数n分解为两个因子的乘积，其因子个数用d(n)表示，则d(n)有的为奇数，有的为偶数，例如d(3)2，因为3有两个因子1，3；d(4)3，因为4有三个因子1，2，4.这当中有无规律可循呢?考察一下表2-3

16、再来分析.表 2-34343从表2-3中我们发现，当自然数是完全平方数，即n1，4，9，16时，d(n)为奇数.这暗示有如下定理：定理定理2.1 正整数M的因子个数d(n)为奇数的充要条件是M是一个完全平方数.定理的证明很简单，只需注意M的因子通常是成双成对出现的，即若abM，则a和b都是M的因子，所以d(M)为偶数，并且仅当M是完全平方数时，为整数，且因此产生一个额外因子，使d(M)成为奇数.MMMM=4444下面用这一结果研究一个有趣的问题狱吏问题.某王国按大赦条款，让一狱吏n次通过一排锁着的M间牢房，每通过一次按照所定规则转动门锁.牢房的门锁的结构是这样的：每转动一次，原来锁着的锁被打开

17、，而原来开着的锁被锁上.如果通过M次，门锁是打开的，则牢房中的犯人将被赦出；而若门锁仍是锁着的，则牢房中的犯人不能赦出.狱吏转动门锁的规则是：第一次通过牢房时，要转动每一把门锁，即把全部锁打开；第二次通过牢房时，要从第二间开始转动，然后每隔一间转动一次；第三次通过牢房时，从第三间开始，然后每隔两间转动一次；第k次通过牢房时，从第k间开始，每隔k1间转动一次；问题就是如此下去，通过n次后，哪些牢房的锁仍然是打开的?4545观察头10多间牢房的结果，很容易猜到答案.现作一般证明如下：考虑第m间牢房，开始它是锁着的，第一次通过时被转动一次；如果m被2除尽，第二次通过时又被转动一次；如果m被3除尽，第

18、三次通过时还被转动一次；等等.即m每有一个因子，锁就要被转动一次，要使门锁最终是打开的，就必须被转动奇数次，所以充分必要条件是：m有奇数个因子，即m是一个完全平方数.实际上，狱吏转动牢房门锁的规则恰与我们前面提到的自然数因子个数问题相同，因为按锁的转动规则的m因子有几个，门锁就要被转动几次，当然只有是完全平方数的牢房，门锁才能被转动奇数次.数的奇偶性隐藏着很多学问，值得我们细心研究.4646量纲分析法是在物理领域运用实验和经验，根据物理定律中的量纲一致原则，确定各物理量之间的关系，实现数学建模的一种方法.本节先介绍量纲一致原则和Buckingham定理，然后给出量纲分析的应用.2.3 量纲分析

19、法量纲分析法47472.3.1 量纲一致原则量纲一致原则人们在研究事物时，需要对对象进行定性或定量分析，这必然涉及许多物理量，例如长度、质量、密度、速度等，这些表示对象的不同物理特征的量，就构成了不同的量纲，通常记为.如时间的量纲为T，长度的量纲为L等.在这些物理量中，有些物理量的量纲是基本的，而另外的物理量的量纲是由基本量纲推导出来的.例如在力学中，常取质量m、长度l、时间t的量纲为基本量纲，分别记为M、L、T，则速度的量纲可表为：vLT14848加速度的量纲可表为：vLT2而力的量纲可根据牛顿第二定律表示为：fMLT2有些物理常数也有量纲，如万有引力定律 中的引力常数k的量纲是kMLT2L

20、2M2M1L3T2无量纲的物理量，称为量纲为1的量.122m mfkr=4949在国际单位制中，有7个基本量：长度、质量、时间、电流、温度、光强度和物质的量，它们的量纲分别是L，M，T，I，J，N称为基本量纲.将一个物理量q表示成基本量纲的幂次之积qLMTINJ(2.3.1)称为该物理量的量纲，幂指数，和称为量纲指数.在用数学公式描述任一物理规律时，等式两端必须保持量纲一致，即在任一物理方程中，所有的项必须有相同的量纲，这就是量纲一致的原则，也称量纲齐次原则.例如自由落体距离计算公式：5050其中v0是初速度，g是重力加速度，t是时间.在上述物理方程中每一项的量纲都是L，它们具有量纲一致性.量

21、纲分析就是在保证量纲一致的原则上，分析和探求物理量之间的关系.这是一个过程，也是一个分析工具，其基本思路是设想物理量能按有意义的方式进行组合，以减少从实验中为得出有关数据需要的试验次数和测量次数.所以说，量纲分析是简化实验和设计实验时，做相似模拟的基础，2012sv tgt=+51 51同时还是指导对实际问题的讨论及研究各种物理量对问题的影响的工具.由于这些原因，量纲分析就成为对某些问题建立数学模型的实用方法.52522.3.2 量纲分析的应用量纲分析的应用例例1 自由落体运动.设有质量为m的小球、从高度为h的位置落下忽略阻力，求自由落体的速度.在这一问题中出现的物理量有v，m，h，g，设它们

22、之间有关系式：(是无量纲的比例系数)取量纲表达式312vm h gaaal=312 vmhgaaa=5353再将vLT1，mM，hL，gLT2代入得按量纲一致原则应有233112LTMLTaaaa+-=1233012-1aaaa=+=5454其解为10，所以得到再经实验测定常数，于是得到自由落体的速度.从中还可看出，自由落体的速度与质量大小无关.2312aa=vghl=2l=2vgh=5555现给出著名的Buckingham定理:定理2.2 设n个物理量x1，x2，x3，xn之间存在一个函数关系f(x1，x2，x3，xn)0(2.3.2)X1，X2，Xm是基本量纲，其中mn.则xi的量纲可表示

23、为:1(1,2,)ijmiijxXina=5656 如果矩阵A(ij)mn的秩为r，则式(2.3.2)与式F(1，2，nr)0等价，式中F是一个未定函数关系，s是无量纲量，且s可表为：其中是线性齐次方程组A0的基本解，(1，2，n)T.()1(1,2,)simsijxsnrbp=-()()()()12(,)ssssnbbbb=5757根据定理，应用量纲分析构造无量纲量的一般方法概括如下：(1)确定与问题有关的物理量(变量与常数)，记作x1，x2，x3，xn，以及这个问题的基本量纲，记作X1，X2，Xm，并把每个物理量纲用基本量纲表示为:式中ij由已知的物理定律确定1(2.3.3)ijmiijx

24、Xa=5858(2)设(2.3.4)则式(2.3.4)的量纲表达式为：再应用量纲一致性原则，得1imiixbp=1011 nijiimmiijjXXa bp=照10,(1,2,)nijiijma b=5959 (3)解齐次线性方程组：如果系数矩阵(ij)mn的秩为r，则方程组有nr个基本解，记为于是得到x1，x2，x3，xn之间的nr个关系式式中s是无量纲量，1，2，nr即为所求的nr个无量纲量.10,(1,2,)nijiijma b=()()()()12(,)(1,2,)ssssnjnrbbbb=-()1(1,2,)simsijxsnrbp=-6060例例2 航船的阻力问题.长l、吃水深度h

25、的船以速度v航行，若不考虑风的影响，那么航船受到的阻力f除依赖船的诸变量l、h、v以外，还与水的参数密度、粘性系数以及重力加速度g有关.下面用量纲分析法确定阻力f和这些物理量之间的关系.(1)航船问题涉及的物理量有：阻力f、船长l、吃水深度h、水的粘性系数、水的密度、船速v、重力加速度g，要寻求的关系记作j(f，l，h，v，g)0(2.3.5)61 61(2)这是一个力学问题，基本量纲选为L、M、T，上述各物理量的量纲表示为223112 (2.3.6)fLMTlLhLvLTL ML MTgLTrm-=6262其中的量纲由基本关系得到.这里p是压强(单位面积受的力)，即：pLMT2L2L1MT2

26、 v是流速，x是长度，有所以111 vLTLTx-=vpxm=12111 /vpL MTTL MTxm-轾 犏=犏 臌6363 (3)式(2.3.6)的量纲矩阵为并可求出R(A)3.3 71111311 1000110 2001012 ()()()()()()()LMTflhvgrm骣-=-桫A6464(4)解齐次方程A0，可得nr734，即有4个基本解，可取为(2.3.7)T1T2T3T4(0110000)(0102001)(0101110)(1202100)=-=-=-=-6565 (5)式(2.3.7)给出4个无量纲量:(2.3.8)从而得到与式(2.3.7)等价的方程112213221

27、4lhlvglvflvpprmprp-=6666j(1，2，3，4)0(2.3.9)式(2.3.8)与式(2.3.9)表达了航船问题中各物理量间的全部关系，这里j是未定的函数.6767(6)为得到阻力f的显式表达式，由式(2.3.9)及式(2.3.8)中的4可写出fl2v2j(1，2，3)其中j表示一个未定函数.在流体力学中称为Froude数，3称为Reynold数，分别计为(2.3.10)vglFr,Revlvglrm=6868则阻力f又表示为(2.3.11)上式就是应用量纲分析确定的航船阻力与各物理量之间的关系.式中函数的形式虽无从知道，但后面会看到这个表达式在物理模拟中的用途.从本例我们看到，对实际问题应用量纲分析时，其前提就是要找出所有与问题有关的物理量，包括变量与常量，既不能遗漏，也不能多余，因为包含了无关变量或丢掉了必需变量，都会使构造的无量纲量出现错误和矛盾.22(,Fr,Re)lfl vhrj=