三问递进数学建模：分段弹道轨迹与误差分析-开发者社区

简介：《第八届MathorCup数学建模竞赛优秀论文C3289》是一份面向数学建模竞赛参赛者、指导教师及军事运筹方向学习者的获奖论文PDF，围绕陆基导弹打击航母的赛题展开，重点解决分段轨道建模、移动目标跟踪和命中率评估等问题。资源包仅含1个PDF文件，大小约511KB，篇幅集中、结构完整，便于直接阅读与打印。论文按赛题三问组织：问题1建立发射段、中段与末段的分段静态轨道模型，并给出MATLAB轨道曲线；问题2引入三维坐标，用Excel散点与SPSS三维散点图刻画动态轨迹；问题3讨论系统误差、随机误差及指令延迟对命中率的影响。目前已有109人学习，适合用来揣摩优秀论文的建模假设、公式表达、图表呈现与论文写作逻辑，也可作为同类赛题的分段建模和误差分析参考。

1. 从一道三问递进的数学建模题说起

数学建模挑战赛里有一类题特别容易把人卡住：题目给一张受力表、几个时间节点，要求你先把运动轨迹切成几段分别建模，再把二维的结论升到三维，最后回头算误差和命中概率。这三问看着是三件事，其实是一条链——第一问定几何骨架，第二问加一个维度，第三问检验前两问的假设到底能撑住多大的不确定度。某届赛事里一份拿奖的论文处理的就是这类题：陆基发射平台对一个海面移动平台做分段弹道建模，发射段抛物线、中段两段变推力、末段平抛，全篇没有一句玄学，靠的全是把每个分段的受力写清楚、把每个衔接点的状态量对上。这套拆法对做毕设、做课设、打竞赛的人都通用，下面把它拆开讲透。

2. 赛题拆解与总体建模路线：三问递进、六条假设和它们埋下的误差

2.1 已知条件怎么变成可计算的量

题面给的是一堆"人话"条件：发射点在某经纬度、目标初始位置在某经纬度、目标以 32 节向正南航行、发射初速竖直向上 500 m/s、末段不超过 20 s、末段最高速度 1000 m/s、目标平台长度 335 m。这些东西必须先统一成米、秒、弧度，否则后面每一个方程都会带着量纲错误往下传。

经纬度到平面距离的换算有个坑：不能直接用大圆距离当斜距，因为后面的三维模型需要分别知道"沿经线方向"和"沿纬线方向"两个分量。我一般用球面近似做局部分解，取题目给定的地球半径 6300 km，纬线方向乘一个平均纬度的余弦修正：

import math R = 6_300_000.0 # 题目给定的地球半径，单位 m phi1, lam1 = 27.5, 120.5 # 发射点纬度/经度（度） phi2, lam2 = 25.65, 123.75 # 目标初始位置纬度/经度（度） dphi = math.radians(abs(phi1 - phi2)) # 纬差 dlam = math.radians(abs(lam1 - lam2)) # 经差 phi_m = math.radians((phi1 + phi2) / 2) # 平均纬度，修正纬线长度 a = R * dphi # 经向距离分量 b = R * dlam * math.cos(phi_m) # 纬向距离分量 L = math.hypot(a, b) # 两点斜距 theta = math.atan2(b, a) # 连线与正南方向的夹角 print(f"经向 {a/1000:.2f} km, 纬向 {b/1000:.2f} km") print(f"斜距 {L/1000:.2f} km, theta = {math.degrees(theta):.2f} 度")

跑出来经向约 203 km、纬向约 320 km、斜距约 378.8 km、夹角约 57.5°。这个 378.8 km 不是凑出来的数——它正好等于后面弹道模型算出来的总水平射程，也就是说题目给的初速、段时长、末段速度这几个参数是配套设计过的。做这类题，第一步先把这类"闭合关系"找出来，后面一旦算出来的射程和它差出几十公里，你立刻知道哪一段的参数抄错了。

速度单位同样要换干净。32 节 = 32 × 1.852 km/h = 59.264 km/h = 16.4622 m/s。这个数会在三维模型里反复出现，写错一位小数，几百秒的飞行时间累下来就是十几公里的偏差。

2.2 三问不是三道题，是一条递进的链

第一问假设目标静止，本质是在求一条"从发射点到固定点"的曲线。它要交付的东西只有三样：每一段的时长、每一段结束时的位置和速度、以及拼起来的整条轨迹方程。这三样东西一旦定死，第二问就变成了"把这条曲线整体抬升到三维空间，再沿目标运动方向加一个匀速分量"，不需要重新推受力。

第二问真正的难点不在物理，在坐标系。第一问是一维水平加一维竖直，X 轴直接指向目标就行；第二问目标会动，如果你还沿着"发射点到初始位置"这条线当 X 轴，目标的运动会斜着穿过坐标面，参数方程会写得很别扭。论文里重建了坐标系：原点仍是发射点，Y 轴取发射点到目标初始位置的水平连线，X 轴取水平面内从 Y 轴顺时针转 90°，Z 轴垂直向上。这样一来，目标沿正南航行这件事，就分解成了 X、Y 两个方向上的匀速分量，每个分量都能单独写进参数方程。

第三问看起来是新增内容，其实是回头审账。你在前面至少做了六条假设：忽略空气阻力、把力简化成牵引力和重力、飞行过程质量不变、控制系统理想无延迟、推力大小恒定方向可调、忽略地球自转。第三问要干的事，就是逐条问"这条假设如果不成立，落点会偏多少"。前两问建得越干净，第三问越好写；反过来，如果前面为了好算而胡乱假设，第三问就会出现"误差大到模型没意义"的尴尬。

2.3 每条假设都在给误差分析埋线

把假设和误差来源做成一张对照表，是这类论文最省事也最显功力的做法：

模型假设	被忽略的物理量	出现在哪一段	影响到什么
忽略空气阻力	与速度平方相关的阻力	全程，末段最严重	射程偏大、轨迹不再是严格抛物线
忽略地球自转	科里奥利惯性力	全程，末段累积明显	横向偏移，量级随飞行时间平方增长
控制系统理想	指令到执行的延迟	中段第二过程	竖直速度归零时刻后移，落点顺航向后偏
推力恒定可调	推力爬升与衰减过程	发射段、中段	分段点位置漂移
质量不变	燃料消耗	全程	加速度被低估
力简化为两个	气动侧力、姿态耦合	全程	高阶小量，通常并入随机误差

这张表的意义在于，它把"误差分析"从一段空话变成了六个可以分别估算的条目。答辩时被问到"你为什么忽略空气阻力"，你指着表说"因为题目明确要求发射段是抛物线，而阻力存在时轨迹不可能是抛物线，这是题目条件倒逼的假设，它的代价在末段体现为射程偏大，量级我在第 6 节估过"——这个回答比"因为阻力复杂"强太多。

提示：这张表最好在建模阶段就填，而不是写完模型再补。填不出来的行，说明你对那条假设的物理含义还没想清楚。

3. 静态分段模型：四段受力表怎么变成一条连续轨迹方程

3.1 发射段：类斜抛与 51 秒的由来

发射段的设定是：初速度竖直向上 500 m/s，同时受到一个大小恒定、方向水平指向目标的牵引力。竖直方向只受重力，水平方向只受牵引力，两个方向独立，这就是一个标准的两维分解。

竖直方向：0 = v0 - g·t1，所以t1 = 500 / 9.8 = 51.02 s，上升到最高点时竖直速度恰好为 0。这个"最高点即发射段结束"不是随便定的，题目说发射段轨迹是抛物线、发射段结束时进入中段，而抛物线顶点正好是竖直速度为零的位置，两者吻合，所以段末时刻被唯一确定。最高点高度h1 = v0²/(2g) = 12755 m。

水平方向：牵引力产生的加速度记作a_F = F/m，这是一个未知量，第一问里它不能自己定，必须由后面的中段约束反解。水平位移L1 = ½·a_F·t1²。

消去时间参数，得到发射段的轨迹方程：

import numpy as np v0, g = 500.0, 9.8 t1 = v0 / g # 51.02 s h1 = v0**2 / (2 * g) # 12755 m def y_stage1(x, aF): """发射段：以水平位移 x 表示的高度 y""" return v0 * np.sqrt(2 * x / aF) - (g / aF) * x # 当 aF = 9.96 时，系数为 224.06 与 0.9839 print(y_stage1(np.array([0.0, 13003.0]), 9.96))

这里有个容易被忽略的自洽性：轨迹方程里√x的系数是v0·√(2/a_F)，一次项系数是g/a_F，两者之比等于v0·√(2a_F)/g。把 a_F = 9.96 代进去，比值是 227.7，与论文给出的 224.06 / 0.9839 完全一致。这种"系数之间能互相验证"的检查，做一次只要两分钟，但能挡住后面 80% 的抄错。