13 一些基本的数学建模示例.docx

1、13 一些基本的数学建模示例1.3 一些基本的数学建模示例1.3.1 椅子的摆放问题 椅子能在不平的地面上放稳吗？下面用数学建模的方法解决此问题。模型准备 仔细分析本问题的实质，发现本问题与椅子腿、地面及椅子腿和地面是否接触有关。如果把椅子腿看成平面上的点，并引入椅子腿和地面距离的函数关系就可以将问题1与平面几何和连续函数联系起来，从而可以用几何知识和连续函数知识来进行数学建模。为讨论问题方便，我们对问题进行简化，先做出如下3个假设：模型假设1、椅子的四条腿一样长，椅子脚与地面接触可以视为一个点，四脚连线是正方形（对椅子的假设）2、地面高度是连续变化的，沿任何方向都不出现间断。（对地面的假设）

2、3、椅子放在地面上至少有三只脚同时着地，（对椅子和地面之间关系的假设）根据上述假设做本问题的模型构成：模型构成用变量表示椅子的位置，引入平面图形及坐标系如图1-1。图中A、B、C、D为椅子的四只脚，坐标系原点选为椅子中心，坐标轴选为椅子的四只脚的对角线。于是由假设2，椅子的移动位置可以由正方形沿坐标原点旋转的角度来唯一表示，而且椅子脚与地面的垂直距离就成为的函数。注意到正方形的中心对称性，可以用椅子的相对两个脚与地面的距离之和来表示这对应两个脚与地面的距离关系，这样，用一个函数就可以描述椅子两个脚是否着地情况。本题引入两个函数即可以描述椅子四 图1-1个脚是否着地情况。记函数f()为椅脚A和C

3、与地面的垂直距离之和。函数g()为椅脚B 和D与地面的垂直距离之和。则显然有f()0、 g()0，且它们都是的连续函数（假设2）。由假设3，对任意的，有f()、 g()至少有一个为0，不妨设当=0时，f(0)0、 g(0)=0，故问题1可以归为证明如下数学命题：数学命题（问题1的数学模型） 已知f()、 g()都是的非负连续函数，对任意的，有f() g()=0，且f(0) 0、 g(0)=0 ，则有存在0，使f(0)= g(0)=0。模型求解证明：将椅子旋转90，对角线AC与BD互换，由f(0)0、 g(0)=0 变为f(/2) =0、 g(/2) 0 构造函数 h()=f() - g(),

4、则有h(0) 0和h(/2) 1）块积木后，再加一块积木的怎样叠放问题。下面我们就n+1（n1）块积木的叠放问题来讨论。假设增加的一块积木插入最底层积木后，我们选择这底层积木的最右端为坐标原点建立如图坐标系（见图1-5）。考虑上面的n块 图1-5积木的重心关系。我们把上面的n块积木分成两部分：1） 从最高层开始的前n-1块积木，记它们的水平重心为x1,总质量为n-12） 与最底层积木相连的第n块积木, 记它的水平重心为x2,质量为1此外，我们也把上面的n块积木看作一个整体，并记它的重心水平坐标，显然n块积木的质量为n。那么，在保证平衡的前提下，上面的n块积木的水平重心应该恰好在最底层积木的右端

5、，即有=0；假设第n块积木超过最底层积木右端的最大前伸距离为z,同样在保证平衡的前提下，从最高层开始的前n-1块积木的总重心的水平坐标为z，即有 x1=z，而第n块积木的水平重心在距第n块积木左端的处，于是在图1-5的坐标系下，有第n块积木的水平重心坐标为x2= 。由重心的关系，有于是有，对三块积木n=2, 第3块积木的右端到第1块积木的右端距离最远可以前伸 对四块积木n=3, 第4块积木的右端到第1块积木的右端距离最远可以前伸设从第n+1块积木的右端到第1块积木的右端最远距离为，则有当时，有。这说明，随着积木数量的无限增加，最顶层的积木可以前伸到无限远的地方。简评: 本题给出的启示是：当问题

6、涉及到较多对象时，对考虑的进行合理的分类进行解决，往往会使问题变得清晰。此外，一些看似不可能的事情其实并非不可能。1.3.4 四足动物的身长和体重关系问题 四足动物的躯干（不包括头尾）的长度和它的体重有什么关系？这个问题有一定的实际意义。比如，生猪收购站的人员或养猪专业户，如果能从生猪的身长估计它的重量可以给他们带来很大方便。模型准备四足动物的生理构造因种类不同而异，如果陷入生物学对复杂的生理结构的研究，将很难得到什么有价值的模型。为此我们可以在较粗浅的假设的基础上，建立动物的身长和体重的比例关系。本问题与体积和力学有关，搜集与此有关的资料得到弹性力学中两端固定的弹性梁的一个结果：长度为L的圆

7、柱型弹性梁在自身重力f作用下, 弹性梁的最大弯曲v与重力f和梁的长度立方成正比，与梁的截面面积s和梁的直径d平方成反比，即利用这个结果，我们采用类比的方法给出假设。模型假设1.设四足动物的躯干（不包括头尾）为长度为L、断面直径为d的圆柱体，体积为m。2.四足动物的躯干（不包括头尾）重量与其体重相同，记为f。3.四足动物可看做一根支撑在四肢上的弹性梁,其腰部的最大下垂对应弹性梁的最大弯曲，记为v。模型构成 根据弹性理论结果及重量与体积成正比关系， 有：由正比关系的传递性，得上式多一个变量v，为替代变量v，注意到是动物躯干的相对下垂度，从生物进化观点，讨论相对下垂度有： 太大，四肢将无法支撑，此种

8、动物必被淘汰； 太小，四肢的材料和尺寸超过了支撑躯体的需要，无疑是一种浪费，也不符合进化理论。因此从生物学的角度可以确定，对于每一种生存下来的动物，经过长期进化后，相对下垂度已经达到其最合适的数值，应该接近一个常数（当然，不同种类的动物，常数值不同）。于是可以得出 d2 L3, 再由f sL和s d2 得f L4，由此得到四足动物体重与躯干长度的关系 它就是本问题的数学模型。模型应用 如果对于某一种四足动物，比如生猪，可以根据统计数据确定公式中的比例常数k而得到用该类动物的躯体长度估计它体重公式。简评: 发挥想象力，利用类比方法，对问题进行大胆的假设和简化是数学建模的一个重要方法。不过，使用此

9、方法时要注意对所得数学模型进行检验。此外，从一系列的比例关系着手推导模型可以使推导问题大为简化。1.3.5 圆杆堆垛问题 把若干不同半径的圆柱形钢杆水平地堆放在一个长方体箱子里，若已知每根杆的半径和最底层各杆的中心坐标，怎样求出其它杆的中心坐标？模型准备 本问题是一个解析几何问题，利用解析几何的有关结论既可。模型假设1. 箱中最底层的杆接触箱底或紧靠箱壁2. 除最底层之外的箱中每一根圆杆都恰有两根杆支撑3. 箱中的钢杆至少有两层以上模型构成本问题如果把箱中所有钢杆一起考虑会带来不便。现把本问题分解为已知三个圆杆的半径和两根支撑杆的坐标来求另一个被支撑杆坐标的三杆堆垛问题。如果三杆堆垛问题解决了

10、，则我们可以利用它依次求得箱中其它所有圆杆的坐标了。虽说涉及的是空间物体，但可以用其堆垛的横截面图化为平面问题来解决。设三个圆杆中两根支撑杆的半径分别为Rl, Rr，对应坐标为(xl,yl),(xr,yr),被支撑杆的半径和坐标为Rt和（xt,yt）。连接三根圆杆的中心获得一个三角形，用a,b,c表示对应的三条边。另用两个支撑圆杆的中心做一个直角三角形，如图1-6由几何知识和三角公式有：xt =xl+acos(+)=xl+a(coscos-sinsin)yt=yl+asin(+)=yl+a(sincos+cossin)这里计算公式中涉及的数据由如下公式获得：a=Rl+ Rtb=Rr+ Rtd=

11、xr- xle=yr- ylc=(d 2+e 2)1/2cos=d/c ; sin=e/ccos=(a2+c2-b2) /2ac 图 1-6 三杆横截面图及辅助三角形 sin=(1-cos2)1/2 在编程计算支撑钢杆的坐标时，为了能快速求出（xt,yt），可以按如下顺序编程计算求解：a=Rl+ Rtb=Rr+ Rtd=xr- xle=yr- ylc=(d 2+e 2)1/2csb=d/c snb=e/ccsa=(a2+c2-b2) /2acsna=(1-csa2)1/2 xt=xl +a (csacsb-snasnb)yt=yl+a (snacsb+csasnb) 有了如上三杆问题的求解，对

12、多于三杆的问题可以按支撑关系和先后顺序依次求出所有其它杆的坐标。例如，如果长方体箱子中有6根圆杆，已知1，2，3号的圆杆在箱底，4号杆由1，2号杆支撑，5号杆由2，3号杆支撑，6号杆由4，5号杆支撑，则可以调用如上三杆问题的算法先由1，2号杆算出4号杆坐标，接着再用2，3号杆算出5号杆坐标，最后用4，5号杆算出6号杆坐标。简评：本题建立计算模型的关键是把原题分解为一组等价的子问题从而是问题得以简化，研究者可以通过讨论子问题的求解获得原问题的解决。这种处理问题的方法可以使复杂问题变得简单有效，它是处理一些有规律复杂问题的常用方法。1.3.6 公平的席位分配问题 席位分配在社会活动中经常遇到，如：

13、人大代表或职工学生代表的名额分配和其他物质资料的分配等。通常分配结果的公平与否以每个代表席位所代表的人数相等或接近来衡量。目前沿用的惯例分配方法为按比例分配方法，即： 某单位席位分配数 = 某单位总人数比例总席位 如果按上述公式参与分配的一些单位席位分配数出现小数,则先按席位分配数的整数分配席位,余下席位按所有参与席位分配单位中小数的大小依次分配之。这种分配方法公平吗？下面来看一个学院在分配学生代表席位中遇到的问题：某学院按有甲乙丙三个系并设20个学生代表席位。它的最初学生人数及学生代表席位为 系名 甲 乙 丙 总数 学生数 100 60 40 200 学生人数比例 100/200 60/20

14、0 40/200 席位分配 10 6 4 20后来由于一些原因，出现学生转系情况，各系学生人数及学生代表席位变为 系名 甲 乙 丙 总数 学生数 103 63 34 200 学生人数比例 103/200 63/200 34/200 按比例分配席位 10.3 6.3 3.4 20 按惯例席位分配 10 6 4 20 由于总代表席位为偶数，使得在解决问题的表决中有时出现表决平局现象而达不成一致意见。为改变这一情况，学院决定再增加一个代表席位，总代表席位变为21个。重新按惯例分配席位,有 系名 甲 乙 丙 总数 学生数 103 63 34 200 学生人数比例 103/200 63/200 34/2

15、00 按比例分配席位 10.815 6.615 3.57 21 按惯例席位分配 11 7 3 21这个分配结果出现增加一席后，丙系比增加席位前少一席的情况，这使人觉得席位分配明显不公平。这个结果也说明按惯例分配席位的方法有缺陷，请尝试建立更合理的分配席位方法解决上面代表席位分配中出现的不公平问题。模型构成 先讨论由两个单位公平分配席位的情况，设 单位 人数 席位数 每席代表人数单位A p1 n1 单位B p2 n2 要公平，应该有=， 但这一般不成立。注意到等式不成立时有 若 ，则说明单位A 吃亏(即对单位A不公平 ) 若，定义 为对A的相对不公平值；若 ，即对单位A不公平，再分配一个席位时，

16、关于,的关系可能有 1. ,说明此一席给A后,对A还不公平;2. ,说明此一席给B后,对A不公平,不公平值为4. ,不可能 上面的分配方法在第1和第3种情况可以确定新席位的分配，但在第2种情况时不好确定新席位的分配。用不公平值的公式来决定席位的分配，对于新的席位分配，若有rB(n1+1,n2)rA (n1,n2+1)则增加的一席应给A ，反之应给B。对不等式 rB(n1+1,n2)rA (n1,n2+1)进行简单处理，可以得出对应不等式引入公式于是知道增加的席位分配可以由Qk的最大值决定，且它可以推广到多个组的一般情况。用Qk的最大值决定席位分配的方法称为Q值法。对多个组（m个组）的席位分配Q

17、值法可以描述为： 1先计算每个组的Q值： Qk , k=1,2,m 2求出其中最大的Q值Qi（若有多个最大值任选其中一个即可） 3将席位分配给最大Q值Qi对应的第i组。这种分配方法很容易编程处理。模型求解 先按应分配的整数部分分配，余下的部分按Q值分配。 本问题的整数名额共分配了席，具体为： 甲 10.815 n1 =10 乙 6.615 n2 =6 丙 3.570 n3 =3对第席的分配，计算Q值Q1=1032/(1011) = 96.45 ; Q2=632/(67)= 94.5; Q3 =342/(34)=96.33因为Q1最大，因此第席应该给甲系; 对第席的分配，计算Q值Q1=1032/

18、(1112)=80.37 ; Q2 =632/(67)=94.5; Q3 =342/(34)=96.33因为Q3最大，因此第1席应该给丙系最后的席位分配为： 甲席乙席 丙席注：若一开始就用Q值分配，以n1=n2=n3 =1逐次增加一席，也可以得到同样的结果。简评：本题给出的启示是对涉及较多对象的问题，可以先通过研究两个对象来找出所考虑问题的一般的规律，这也是科学研究的常用方法。请对一般情况编程。1.3.7 中国人重姓名问题由于中国人口的增加和中国姓名结构的局限性，中国人姓名相重的现象日渐增多，特别是单名的出现，使重姓名问题更加突出，其引出的矛盾也日益突出。例如有时在一个60人的集体里，重姓名竟

19、有两对！若在更大范围里调查，重姓名的严重情况将会更多。重姓名现象引起的误会与带来的弊端是众所周知的，伴随着我国经济文化的高速发展和对外交往的扩大，重姓名引起的问题将更加突出，可以说有效地克服重姓名问题也即中国人姓名改革是迫在眉睫！因此，用合理的方法对中国人姓名进行改革是非常必要的。请尝试提出一个合理且可以有效解决此问题的中国人取名方案。模型准备先来看一下中国姓名的结构和取名习惯，中国的姓名是由姓和名来组成的。姓在前名在后，目前姓大约有5730个但常用姓只有2077个左右。名通常由两个字组成，较早时名的第一个字体现辈份，由一个家族的族谱决定，后一个字可任选。近来随着家族概念的淡化，名子中已经无辈

20、份之意，人们为图好记，只有一个字的单字名在增多。组合数学中的乘法原理和鸽笼原理可以非常简单地解释拥有十几亿人口的中国重名现象。姓名是由汉字排列而成，构成姓名的汉字多，则姓名总数就多。要想有效地克服重姓名问题，就该增加姓名的汉字数。因此，本问题可以用排列组合理论来解决。模型假设1 中国的所有姓名共有N个,其中姓有S个；2 取名的方法和习惯不改变，即姓名中应该父亲姓氏在姓名首位模型构成 靠机械地增加名字的个数解决重姓名问题或完全改变现有的姓名是不明智和不可取的。应该采用兼顾现有姓名习惯来做这件事。为扩大姓名集合并考虑到中国姓名的特色和兼顾原有取名习惯，利用排列组合的理论，这里提出如下体现父母姓的复

21、姓名方式来解决重姓名问题。引入的中国姓名取名方法称为FM取名方法, 这里 F和M分别是中文父(f)母(m)二字的拼音首字母，同时也是英语父(father)母(mother)二词的首字母，它表示父母之意，“FM姓名”是与父母有关的姓名。一个FM姓名的结构为：主姓名. 辅姓名这里主姓名就是现在的人们所用的姓名，而辅姓名可以只是母亲的姓，也可以是用母亲姓起的另一个姓名，不过这个姓名要名在前姓在后以区别于主姓名，中间的 . 是间隔号，如果用？和分别表示父母姓，* 和 表示对应的两个名字，则FM姓名可表示为？* . 或 ？* . FM姓名的结构的取名方法为： 给一个人取一个FM名是很简单的，只要按以前的习惯用父母姓各取一个姓名，然后按FM姓名结构的要求就得到FM姓名。例如：父亲姓王，母亲姓孙，给孩子取的名字是中华和靖，则孩子的FM姓名为：王中华靖孙如果只取一个名字，则孩子的FM姓名为：王中华孙显见，这种“FM姓名”对原姓名改动较小，也无须重新翻字典去找新的名字，易于在人口普查时对全国公民的姓名做统一改动。在一般场合或不易引起混淆的情况下，直接使用主姓名或原来的姓名即可，但是在正式场合，如译写著作，论文及发明创造，申请专利等署名是应写完整的“FM姓名”，这有利于中国姓名朝正规化、合理化的多字姓名过渡。模型分析 由假设1，按排序原则，在“FM姓名体系”下