数模竞赛中部分几何物理问题解析

数模竞赛中部分几何物理问题解析谭劲英14/07/181.CUMCM-1995A：一个飞行管理问题2.CUMCM-2000D：空洞探测3.CUMCM-2010A：储油罐的变位识别与罐容表标定CUMCM-1995A：一个飞行管理问题在约10000m高空的某边长160km的正方形区域内，经常有若干架飞机作水平飞行，区域内每架飞机的位置和速度向量均由计算机记录其数据，以便进行飞行管理。当一架欲进入该区域的飞机到达边界区域边缘时，记录其数据后，要立即计算并判断是否会与其区域内的飞机发生碰撞。如果会碰撞，则应计算如何调整各架(包括新进入的)飞机飞行的方向角，以避免碰撞。现假设条件如下：1)不碰撞的标准为任意两架飞机的距离大于8km；2)飞机飞行方向角调整的幅度不应超过30度；3)所有飞机飞行速度均为每小时为800km；4)进入该区域的飞机在到达区域边缘时，与区域内飞机的距离应在60km以上；5)最多考虑6架飞机；6)不必考虑飞机离开此区域后的状况。请你对这个避免碰撞的飞行管理问题建立数学模型。列出计算步骤，对以下数据进行计算(方向角误差不超过0.01度)，要求飞机飞行方向角调整的幅度尽量小。设该区域4个顶点坐标为(0,0)，(160,0)，(160,160)，(0,160)。记录数据为：飞机编号横坐标x纵坐标y方向角(度)1150140243285852363150155220.54145501595130150230新进入0052注：方向角指飞行方向与x轴正向的夹角。初始位置：时刻t飞机的位置：两架飞机的距离（平方）：两架飞机不碰撞的条件222)()()(tjttjtijyyxxtrii,sin,cos00iitiitivtyyvtxxi064)()(2trtfijij（0≤t≤Tij）Ti为第i架飞机飞出区域的时刻：不碰撞条件：000),,(iiiyxiii0不必考虑在区域外的碰撞！两架飞机都在区域中的时间：),min(jiijTTT00000000000000000000tan,223tan,23,sin,tan,23tan,2,cos,tan,2tan,20,sin,tan,223tan,20,cosiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiixDyorxyifvyxyorxyDifvxxyDorxDyDifvyDxDyorxDyDifvxDT具体来看，第i架飞机在区域内的时间：飞机飞出区域的时刻0i.)(2ijijijijijczbztf整理：fij(t)的最小值(-bij2/4+cij)；此时其中：.2sin4*jiijijvbt不碰撞条件的等价表述最后，优化模型为0*ijt若fij(t)大于等于０肯定成立ijijTt*若fij(t)大于等于０等价于0)(ijijTfijijTt*0若fij(t)大于等于０等价于0)(*ijtfij,042ijijcb000),,(iiiyxiii0.61iiMin其他目标：调整后的方向角：总的调整量最小：||max6,...,1iiMin最大调整量最小：初始位置与方向角：山体隧道坝体等的某些内部结构可用弹性波测量来确定。简化问题可叙述为，一块均匀介质构成的矩形平板内有一些充满空气的空洞。在平板的两个邻边分别等距地设置若干波源，在他们的对边对等地安放同样多的接收器，记录弹性波由每个波源到达对边上每个接收器的时间。根据弹性波在介质和在空气中不同的传播速度来确定板内空洞的位置。CUMCM-2000D：空洞探测具体问题：一块240(米)×240(米)的平板ABCD：在AB边等距地设置7个波源Pi(i=1,…,7)，在CD边等距地设置7个接收器Qj(j=1,…,7)，记录由Pi发出的弹性波到达Qj的时间tij(秒)；在AD边等距地设置7个波源Ri(i=1,…,7)，在BC边等距地设置7个接收器Sj(j=1,…,7)，记录由Ri发出的弹性波到达Sj的时间τij(秒)。已知弹性波在介质和空气中的传播速度分别为2880(米/秒)和320(米/秒)，且弹性波沿板边缘的传播速度与在介质中的传播速度相同。CBADP2Q4R3S6TP=(tij)tijQ1Q2Q3Q4Q5Q6Q7P1.0611.0895.1996.2032.4181.4923.5646P2.0989.0592.4413.4318.4770.5242.3805P3.3052.4131.0598.4153.4156.3563.1919P4.3221.4453.4040.0738.1789.0740.2122P5.3490.4529.2263.1917.0839.1768.1810P6.3807.3177.2364.3064.2217.0939.1031P7.4311.3397.3566.1954.0760.0688.1042TR=(τij)τijS1S2S3S4S5S6S7R1.0645.0602.0813.3516.3867.4314.5721R2.0753.0700.2852.4341.3491.4800.4980R3.3456.3205.0974.4093.4240.4540.3112R4.3655.3289.4247.1007.3249.2134.1017R5.3165.2509.3214.3256.0904.1874.2130R6.2749.3891.5895.3016.2058.0841.0706R7.4434.4919.3904.0786.0709.0914.0583要求：(1)确定该平面内空洞的位置。(2)只根据Pi发出的弹性波到达Qj的时间tij能确定空洞的位置吗？讨论在同样能够确定空洞位置的前提下，减少波源和接收器的方法。分析：弹性波沿平板边缘的理论传播时间：t=240/2880=0.0833(秒)弹性波沿平板边缘的实际传播时间：t11=.0611,t77=.1042,τ11=.0645,τ77=.0583题目中已假设“弹性波沿板边缘的传播速度与在介质中的传播速度相同”。观测数据的最大绝对误差为d=0.025秒。可以认为，0.025*320=8(米）以下的空洞是探测不出的。假设1.观测数据有测量误差。观测数据除测量误差外是可靠的。2.波在传播过程中沿直线单向传播，且不考虑波的反射、折射以及干涉等现象。3.空气密度和介质密度都均匀。4.“弹性波”在传播过程中没有能量损失。其波速仅与介质有关，且在同一均匀介质中波速不变。弹性波沿板边缘的传播速度与在介质中的传播速度相同。5.假设平板可划分化为网格，空洞定位于每个网格单元内，空洞大小大致相同。波线与网格交线长度的计算(k,l);52120,;711,40kjikjikjikjikbijkl，且或如果或如果记波源Pi与接收器Qj决定的波线与每个单元（k，l）的交线长度为bijkli=j时，123456654321PiQj决定的直线方程：(j-i）y=6(x-40(i-1))i=j以外的情况单元(k,l)左边缘直线方程x=40(k-1)波线与单元(k,l)左边缘对应交点的y坐标为y1ijkl=240(k-i)/(j-i)，其中l-1≤6(k-i)/(j-i)≤l(k,l)波线与网格交线长度的计算PiQj决定的直线方程：(j-i）y=6(x-40(i-1))i=j以外的情况单元(k,l)右边缘直线方程x=40k波线与单元(k,l)右边缘对应交点的y坐标为y2ijkl=240(k+1-i)/(j-i)，其中l-1≤6(k+1-i)/(j-i)≤l(k,l)波线与网格交线长度的计算PiQj决定的直线方程：(j-i）y=6(x-40(i-1))i=j以外的情况单元(k,l)下边缘直线方程y=40(l-1)波线与单元(k,l)下边缘对应交点的y坐标为y3ijkl=40(l-1)，其中0≤6(i-k)-(i-j)(l-1)≤6(k,l)波线与网格交线长度的计算PiQj决定的直线方程：(j-i）y=6(x-40(i-1))i=j以外的情况单元(k,l)上边缘直线方程y=40l波线与单元(k,l)上边缘对应交点的y坐标为y4ijkl=40l，其中0≤6(i-k)-(i-j)l≤6(k,l)波线与网格交线长度的计算交线在y轴的投影长度(交点条件最多只有2个成立)i=j以外的情况dyijkl=max(y1ijkl,y2ijkl,y3ijkl,y4ijkl)-min(y1ijkl,y2ijkl,y3ijkl,y4ijkl)由相似三角形关系QjABCDPiRiSjPjdyijkl(k,l)bijklEFGbijkl=aijdyijkl/240)(624024022jiaiji=j也成立波线与网格交线长度的计算由对称性，RiSj与单元(k,l)的交线长度ci,j,k,l=bj,i,l,7-k波线与网格交线长度的计算参量、变量：xkl：单元(k,l)是否为空洞(1：是；0：否)aij：波源Pi与接收器Qj，或Ri与Sj之间的距离Pij：经过介质的长度,qij经过空气的长度tij(同样ij):传播时间观测值优化模型（拟合／回归）71k,lklijklijijijijxbaqap71k,lklijklijxbq若没有误差：tij=pij/v1+qij/v2271171//vxbvxbatk,lklijklk,lklijklijij同理：271171//vxcvxcak,lklijklk,lklijklijij模型：22711712271171////minvxcvxcavxbvxbatk,lklijklk,lklijklijijk,lklijklk,lklijklijij优化模型（拟合／回归）计算结果空洞X(P2,Q2)1X(P2,Q3)1X(P2,Q5)1X(P3,Q2)1X(P3,Q3)1X(P3,Q4)1X(P4,Q4)1X(P5,Q3)1通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。CUMCM-2010A：储油罐的变位识别与罐容表标定图1是一种典型的储油罐尺寸及形状示意图，其主体为圆柱体，两端为球冠体。图2是其罐体纵向倾斜变位的示意图。图3是罐体横向偏转变位的截面示意图。请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。（1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为的纵向变位两种情况做了实验，实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。4.1（2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度和横向偏转角度）之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。附件