简介：围绕2003年全国大学生数学建模竞赛D题「抢渡长江」的参考论文，面向备战数学建模竞赛的选手、指导教师及想掌握最优控制建模思路的学习者，重点解决恒定与分段水流条件下如何选择游泳角度θ(t)，使渡江总时间最短的问题。压缩包内共1个doc文档，约124KB，内容为命题人给出的参考答案，便于直接对照题意梳理建模脉络。文档从运动方程出发，论证最优路径为直线，推导斜率z与到达时间T的关系，并由二次方程求根及可行性条件筛选合理解；结合H=1160m、L=1000m、v=1.89m/s等数据给出角度、速度与用时的示例计算。文档还讨论垂直游动所需速度阈值，以及将水流分为三段后动态调整方向、用枚举表逼近最短总时间的做法。已有777人学习，适合用来复盘完整推导、核对数值结果与理解复杂水流下的建模取舍。
1. 抢渡长江这道 D 题，参考论文真正值钱的是哪几页
翻开 2003 年全国赛 D 题的参考论文，多数人直奔公式页，把最优路线那条解析表达式抄进笔记就合上了。这道题问的是渡江选手在水流里怎么选方向、能不能到岸、游得更快是不是成绩就一定更好，本质是带流场的路径优化，数学建模里最典型的“路线 + 时间”型问题。抄完公式你会发现，换一道题还是不会做，因为参考论文里真正值钱的东西不在那几行：前面那张假设表决定了模型能走多远，中间那段流速剖面的处理决定了结论对不对，最后用实测成绩做检验的那一章决定了这篇论文能走到什么位置。下面按这三块拆，把能抄作业的代码、要调的参数和容易翻车的地方摆出来。
2. 把渡江问题翻译成可计算的模型：坐标系、流场与三个决策变量
拿到题目先别急着写公式。这道题的文字描述里塞了大量背景，真正有用的信息只有三类：江面几何（宽度、起终点位置）、水流分布（沿江宽方向的流速剖面数据）、选手能力（相对水的游速）。把这三类信息翻译成符号，剩下的就是标准化操作。
2.1 坐标系先定死：把二维问题压成一维积分
取出发点为原点，x 轴沿江岸指向下游为正，y 轴垂直岸线指向对岸为正，江宽记作 L，终点坐标记作 (X_B, L)。选手位置 (x(t), y(t))，相对水的游速为 v，游动方向与 x 轴正向的夹角为 θ。运动学方程只有两行：
dx/dt = v·cosθ + u(y) dy/dt = v·sinθ第一行里 u(y) 是水流沿 x 方向的速度，它把选手往下游推；第二行是横渡分量，只有它能把人送到对岸。到这里大部分人会直接上二维数值积分，然后卡在“怎么搜路线”上。正确做法是换自变量：因为 u 只随 y 变化，把 y 当自变量，把两式相除，得到
dx/dy = (v·cosθ + u(y)) / (v·sinθ) T = ∫₀ᴸ dy / (v·sinθ(y))二维轨迹问题就此变成一个一维积分问题：给定一条 θ 随 y 变化的策略，横渡时间是一个一维积分，落点是一个一维积分，约束只有一条 x(L) = X_B。这是整篇参考论文里最关键的化简，缺了这一步，后面的数值结果很难解释清楚。
提示：写参考论文时把这段推导单独成节，并在符号表里明确 θ 的零点与正方向。角度定义含混是这道题最常见的失分点，评审看路线图时第一眼就是判断角度是不是从正确的轴量起。
对应的假设表也是在这时候列的。下面这张表是这道题常见的假设及其代价，写论文时可以直接照着改：
| 假设 | 内容 | 放松后要改什么 |
|---|---|---|
| H1 | 选手相对水的速度大小恒定，方向可任意调整 | 引入体力衰减 v(t)，模型变成时变最优控制 |
| H2 | 水流只有沿岸方向分量，无横向分量 | 加入 v_y(y)，横渡时间不再与路线解耦 |
| H3 | 流场沿江宽方向一维分布，沿程不变 | 变成 u(x,y)，需要二维流场数据 |
| H4 | 选手入水点与出水点固定，不允许沿岸陆上移动 | 允许陆上移动后要加入走行速度与切换点优化 |
| H5 | 忽略出发与到岸的加速过程 | 加入加减速段，短距离赛程下影响不可忽略 |
H2 和 H4 是最值得在论文里讨论的两条。真实渡江比赛里，选手往往会先贴着岸边游一段再斜穿，原因就是岸边流速低；如果你的模型不允许沿岸方向“免费”前进，就解释不了这个现象。
2.2 流速剖面怎么拟合：分段线性、抛物线，还是直接插值
题目给的是一张测点数据表：若干个横向位置 y_i 对应的流速 u_i，通常呈现“两岸小、江心大”的形状。处理方式有三种，选择取决于你后面要做什么。
第一种是抛物线拟合 u(y) = 4·u_max·y·(L−y)/L²。好处是解析可积、形式漂亮、参数只有一个 u_max；坏处是它强制岸边流速为零，而实测数据里岸边往往还有一个不小的残余流速，另外真实剖面的峰值位置通常不在正中间，抛物线拟合会把峰值锁死在 L/2。
第二种是分段线性插值。保形、不会插出负流速、能忠实反映岸边低速带的宽度，缺点是导数不连续，做梯度类优化时目标函数会有轻微抖动。实际做这道题，我一般用这一种，因为“岸边低速带有多宽”恰恰是决定最优路线的关键量，不能被一个光滑函数抹平。
第三种是分段三次样条。光滑性最好，但容易在数据稀疏处过冲，出现物理上不可能的负流速或者双峰。要用的话必须做单调性检查。
不管选哪种，落笔前先算两个派生量：峰值流速 u_max 和低速带宽度 δ（定义为 u < 0.3·u_max 的横向范围）。这两个数直接决定后面可行性判据的松紧，也是灵敏度分析里最该扰动的参数。
import numpy as np # 示意数据：第一列是归一化横向位置 y/L，第二列是流速 (m/s) # 实际做题时替换成题目给的测点表 PROFILE = np.array([ [0.00, 0.35], [0.10, 0.90], [0.25, 1.35], [0.50, 1.70], [0.75, 1.45], [0.90, 0.95], [1.00, 0.30], ]) def make_u(L): yn = PROFILE[:, 0] un = PROFILE[:, 1] def u(y): # np.interp 默认做线性插值，超出范围时钳制到端点值 return np.interp(np.asarray(y) / L, yn, un) return u L = 1000.0 # 江宽，示意值 u = make_u(L) print("峰值流速", PROFILE[:, 1].max(), "低速带宽度占比", PROFILE[PROFILE[:, 1] < 0.3 * PROFILE[:, 1].max()][:, 0])这段代码只做一件事：把测点表包装成一个可调用的函数 u(y)，后面所有积分和优化都复用它。注意 np.interp 在区间外做的是端点钳制而不是外推，这正是我们要的行为，避免优化过程中把 y 推到负数时拿到荒谬的流速。
2.3 可行性判据：什么速度下根本到不了终点
在动手写数值优化之前，先用均匀流情形推一个解析判据，这个判据能解释后面很多现象。假设 u 为常数、θ 也恒定（时不变系统的最优控制就是常数控制），从 (0,0) 到 (X_B, L) 要求
(v·cosθ + u)·T = X_B v·sinθ·T = L两式相除，记 k = X_B / L：
cosθ − k·sinθ = −u/v左边是一个振幅为 √(1+k²) 的余弦型表达式，因此有解的必要条件是
v ≥ u / √(1 + k²)这个不等式就是可行性判据。代一组示意数值：u = 1.5 m/s、k = 0.8，临界速度是 1.171 m/s。选手速度低于这个值，无论怎么调整方向，到达对岸时都会落过终点下游，模型无解。
判据还说明两件事。终点越偏下游（k 越大），可行性越宽松；终点偏上游（k 为负）时虽然 √(1+k²) 形式不变，但需要顶着水流游，实际耗时急剧上升。把这条判据写进论文，比单纯给一张“最优路线图”有说服力得多，因为它回答的是“什么情况下这场比赛对某类选手根本不成立”。
3. 用 Python 复现参考论文的计算：均匀流解析解到非均匀流数值解
模型立住之后，剩下的就是算。建议按“解析解 → 常数角打靶 → 逐点最优控制”三级递进，每一级都能和上一级对照，出问题时能立刻定位是哪一层写错了。
3.1 均匀流的解析解：先算一个能手算校验的基准
把 2.3 的推导解出来，得到最优角的闭式表达。设 A = √(1+k²)、φ = arctan(k)，则
θ = arccos(−u / (v·A)) − φ T = L / (v·sinθ)import numpy as np def uniform_route(v, u, L, xb): """均匀流场下最短时间路线的解析解。 v : 选手相对水的速度 (m/s) u : 沿 x 方向的均匀流速 (m/s)，正值为下游 L : 江宽 (m) xb : 终点相对出发点的下游偏移 (m) 返回 (theta_rad, T_sec)；无可行解时返回 (None, None) """ k = xb / L A = np.hypot(1.0, k) if v * A <= abs(u): # 触及可行性边界，无解 return None, None phi = np.arctan(k) theta = np.arccos(-u / (v * A)) - phi T = L / (v * np.sin(theta)) return theta, T if __name__ == "__main__": L, xb = 1000.0, 800.0 # 示意参数，非原题数据 <p> <a href="https://download.csdn.net/download/zichengyixuan/3258104" style="color:#ec7500;font-size:14px;"> 本文还有配套的精品资源，点击获取 </a> <img alt="menu-r.4af5f7ec.gif" src="https://csdnimg.cn/release/wenkucmsfe/public/img/menu-r.4af5f7ec.gif" style="width:16px;margin-left:4px;vertical-align:text-bottom;cursor:text;"> </p>
