1. 赛题回顾与核心挑战解析每年九月的全国大学生数学建模竞赛对很多理工科学生来说都是一场“硬仗”。2023年的B题题目是“多波束测线问题”初看之下很多同学可能会觉得这又是一个典型的优化问题无非是设计测线、覆盖区域、计算效率。但真正上手后你会发现这道题远不止于此它更像是一个披着优化外衣的“系统工程”问题考察的是从物理模型理解、到数学抽象、再到算法实现、最后到结果分析的完整链条。我指导过几支队伍也复盘了网上一些优秀论文发现这道题的“坑”和“亮点”都非常集中。今天我就以一个过来人的视角拆解一下这道题的核心脉络、常见误区以及那些能让论文脱颖而出的关键点。这道题背景是海洋测绘使用多波束声呐对一片矩形海域进行测量。声呐会向海底发射一个扇形的波束这个波束在垂直于航迹的方向上横向有一个覆盖宽度而沿着航迹方向纵向的覆盖则由相邻测点决定。题目给出了波束开角、海水深度、测线间距等关键参数。核心任务有两个第一在给定重叠率约束下设计测线布局使得总测线长度最短即测量效率最高第二讨论测线方向与海底坡面的夹角对测量覆盖率的影响。这听起来很清晰但难点在于如何将“重叠率”这个物理约束精准地转化为数学模型以及如何处理海底地形起伏带来的复杂性。2. 物理模型到数学模型的精准转化重叠率是核心这是整道题的第一个分水岭。很多队伍在这里就出现了理解偏差导致后续模型全盘皆错。题目中提到的“重叠率”不是指两条相邻测线在海底投影区域的面积重叠比例而是指相邻波束在海底脚印边缘的横向重叠程度。这是一个一维线性的概念而不是二维面积的概念。2.1 波束脚印与覆盖宽度的计算多波束声呐的单个波束在海底会形成一个近似椭圆形的“脚印”。其横向宽度垂直于航迹方向是建模的关键。在平坦海底假设下这个宽度 ( W ) 可以通过简单的几何关系求得 [ W 2 \times D \times \tan(\theta / 2) ] 其中( D ) 是海水深度( \theta ) 是波束的开角。这里要注意单位统一深度 ( D ) 和宽度 ( W ) 通常以米为单位。然而题目中海底是有坡度的。当测船航迹与海底等深线即坡度方向存在夹角时实际的有效深度和覆盖宽度会发生变化。这是第一个容易忽略的细节。你需要建立一个坐标系将海底坡度、测线方向、波束平面三者之间的关系用向量和三角函数表达出来。实际计算覆盖宽度时深度 ( D ) 应取测线中心点处的海水深度并考虑坡度引起的沿波束横向的水深变化。一个简化的处理是使用中心点深度进行计算但要在模型假设中明确指出这一点并讨论其带来的误差。2.2 重叠率约束的数学表达设相邻两条测线的间距为 ( d )。为了保证测量数据无缝拼接且不过度冗余要求相邻波束脚印的边缘重叠一部分。设要求的最小重叠率为 ( \eta )例如10%。 那么数学模型应该是 [ \frac{W - d}{W} \ge \eta ] 即 [ d \le (1 - \eta) \times W ] 这个不等式意味着为了满足最小重叠率测线间距 ( d ) 不能大于覆盖宽度 ( W ) 的 ( (1-\eta) ) 倍。这里一个巨大的坑是有些同学错误地理解为 ( d \ge \eta \times W )这完全颠倒了。重叠率是“重叠部分占单条覆盖宽度的比例”所以间距越小重叠部分越大重叠率越高。我们的目标是在满足最小重叠率的前提下让间距 ( d ) 尽可能大因为间距越大所需测线越少总长度越短但不能超过 ( (1-\eta)W ) 这个上限。因此在平坦海底且测线平行布置的情况下最优的测线间距就是取等号 [ d_{opt} (1 - \eta) \times W ] 这样既能刚好满足重叠率要求又能使测线数量最少。3. 问题一矩形区域最优测线布局设计问题一是在一个给定的矩形海域内设计一组平行的测线使得在满足重叠率约束下所有测线的总长度最短。这本质上是一个一维条带填充Strip Packing问题。3.1 模型建立与求解确定单条测线长度矩形区域长 ( L )宽 ( W_{area} )。如果测线方向平行于长边则单条测线长度 ( l_{track} L )如果平行于宽边则 ( l_{track} W_{area} )。题目通常要求测线平行于矩形的一边我们需要比较两种方向的总长度。计算所需测线条数在确定了测线方向假设平行于长边和最优间距 ( d_{opt} ) 后我们需要用测线覆盖整个区域的宽度 ( W_{area} )。第一条测线通常从区域边界开始布放。 所需测线数量 ( N ) 需要向上取整 [ N \lceil \frac{W_{area}}{d_{opt}} \rceil 1 ] 这里“1”和取整函数ceil的运用是关键。因为第一条测线放在边界最后一条测线必须覆盖另一条边界即使最后剩余宽度小于 ( d_{opt} )也需要一条测线去覆盖。更严谨的计算是( N \lceil (W_{area} - W/2) / d_{opt} \rceil 1 )其中 ( W/2 ) 是考虑波束半宽覆盖从边界开始。但很多优秀论文采用了更简洁的ceil(W_area / d_opt) 1模型并在假设中说明已将边界考虑在内。计算总长度总长度 ( L_{total} N \times l_{track} )。方向对比分别计算测线平行于长边和平行于宽边两种情况下的总长度选择总长度更短的那个方案作为最优布局。一个重要的洞察是由于 ( d_{opt} ) 是固定的总长度主要取决于 ( N \times l_{track} )。而 ( N ) 与区域宽度垂直于测线的方向成正比( l_{track} ) 就是区域长度平行于测线的方向。因此最优策略是让测线平行于区域较长的那一边。这样可以减少测线条数 ( N )虽然单条线变长了但通常总长度会更优。这可以通过简单的代数比较证明应在论文中清晰展示。3.2 编程实现与结果验证这部分需要编程计算。以Python为例核心代码结构如下import math # 参数输入 L 4.0 # 区域长 (海里需换算) W_area 3.0 # 区域宽 (海里) D 120.0 # 中心深度 (米) theta math.radians(120) # 波束开角转弧度 eta 0.1 # 最小重叠率 # 单位换算1海里1852米 L_m L * 1852 W_area_m W_area * 1852 # 计算覆盖宽度 W_beam 2 * D * math.tan(theta / 2) # 计算最优间距 d_opt (1 - eta) * W_beam # 情况1测线平行于长边L track_length_1 L_m num_tracks_1 math.ceil(W_area_m / d_opt) 1 # 简化模型 total_length_1 track_length_1 * num_tracks_1 # 情况2测线平行于宽边W_area track_length_2 W_area_m num_tracks_2 math.ceil(L_m / d_opt) 1 total_length_2 track_length_2 * num_tracks_2 # 比较结果 if total_length_1 total_length_2: print(f最优方案测线平行于长边。总长度{total_length_1:.2f}米) else: print(f最优方案测线平行于宽边。总长度{total_length_2:.2f}米)注意在实际论文中必须将结果换算回海里单位并给出清晰的布局示意图。示意图可以用Matplotlib绘制显示矩形区域和等间距的平行测线。4. 问题二坡度影响与测线方向优化问题二引入了海底坡度并问测线方向与坡面方向即等深线法线方向的夹角 (\alpha) 如何影响测量。这是本题的深化和难点考察对模型物理意义的理解。4.1 坡度如何影响覆盖宽度当存在坡度为 (\beta) 的斜面且测线方向与坡面方向夹角为 (\alpha) 时声呐波束平面与海底斜面的交线几何变得复杂。此时波束在海底的横向覆盖宽度不再简单地是 (2D\tan(\theta/2))。建立三维几何模型需要建立空间直角坐标系。设海底斜面为平面其法向量由坡度 (\beta) 和走向决定。测线方向是水平面上的一个方向。波束平面是垂直于测线方向、且包含测线的一个垂直平面。计算有效深度与投影波束平面与海底斜面的交线决定了实际照射到的海底区域。这条交线的横向跨度就是有效覆盖宽度。它依赖于中心点处的海水深度 (D)。坡度角 (\beta)。夹角 (\alpha)。波束开角 (\theta)。推导覆盖宽度公式通过空间解析几何可以推导出覆盖宽度 (W(\alpha)) 的表达式。一个相对简洁的近似公式适用于小坡度是 [ W(\alpha) \approx \frac{2D\tan(\theta/2)}{\cos\beta - \sin\beta \cdot \tan(\theta/2) \cdot \sin\alpha} ] 这个公式显示覆盖宽度 (W) 是夹角 (\alpha) 的函数。当 (\alpha 0^\circ)测线平行于等深线即垂直于坡度方向时(\sin\alpha0)宽度最大当 (\alpha 90^\circ)测线垂直于等深线即沿坡度方向时(\sin\alpha1)宽度最小。4.2 影响分析与优化方向基于上述公式我们可以分析夹角 (\alpha) 的影响对单条测线覆盖宽度的影响如公式所示(W(\alpha)) 随 (\alpha) 增大而减小。这意味着沿着坡度方向测量(\alpha) 大时单次扫描覆盖的海底条带更窄。对测线总长度的影响由于最优间距 (d_{opt} (1-\eta)W(\alpha))当 (W(\alpha)) 变小时(d_{opt}) 也同比减小。为了覆盖同样宽度的区域所需的测线条数 (N) 会增加。同时单条测线长度不变假设区域矩形且测线平行于某边。因此总测量长度 (L_{total}) 会随着 (\alpha) 的增大而显著增加。也就是说沿着坡度方向测量效率最低。对测量质量覆盖率均匀性的影响这是更高阶的考量。当 (\alpha) 不为0时波束脚印在斜坡不同位置的水深不同导致脚印形状和大小发生变化可能引起边缘区域的照射强度不均、信噪比下降从而影响数据质量。在论文中可以定性讨论这一点或者定义一个简单的“质量系数”如基于水深变化率来量化。优化建议从纯效率总长度最短角度最优的测线方向是平行于海底等深线即 (\alpha 0^\circ)。这样覆盖宽度最大测线间距最大所需测线数量最少。在实际海洋测绘中这也确实是常用策略因为垂直于等深线航行可能导致船位控制困难和水深剧烈变化。4.3 模型求解与可视化对于问题二需要编程计算不同 (\alpha) 角下的 (W(\alpha)) 和对应的总长度 (L_{total}(\alpha))并绘制关系曲线。import numpy as np import matplotlib.pyplot as plt # 参数 D 120.0 theta np.radians(120) beta np.radians(1.5) # 坡度1.5度转弧度 eta 0.1 L_m 4 * 1852 W_area_m 3 * 1852 # 定义覆盖宽度函数近似公式 def coverage_width(alpha_deg): alpha np.radians(alpha_deg) denominator np.cos(beta) - np.sin(beta) * np.tan(theta/2) * np.sin(alpha) # 避免除零理论上cos(beta)足够大 W 2 * D * np.tan(theta/2) / denominator return W # 定义总长度函数假设测线平行于长边 def total_length(alpha_deg): W_beam coverage_width(alpha_deg) d_opt (1 - eta) * W_beam N math.ceil(W_area_m / d_opt) 1 return N * L_m # 计算不同角度下的值 alphas np.arange(0, 91, 10) # 0到90度间隔10度 widths [coverage_width(a) for a in alphas] lengths [total_length(a) for a in alphas] # 绘制结果 fig, (ax1, ax2) plt.subplots(1, 2, figsize(12, 4)) ax1.plot(alphas, widths, bo-) ax1.set_xlabel(测线与坡面夹角 α (度)) ax1.set_ylabel(覆盖宽度 W(α) (米)) ax1.set_title(覆盖宽度随夹角变化) ax1.grid(True) ax2.plot(alphas, lengths, rs-) ax2.set_xlabel(测线与坡面夹角 α (度)) ax2.set_ylabel(总测线长度 (米)) ax2.set_title(总测线长度随夹角变化) ax2.grid(True) plt.tight_layout() plt.show()通过这样的图表可以直观展示 (\alpha) 对效率和覆盖宽度的影响有力支撑“应使测线平行于等深线”的结论。5. 论文写作提分关键与常见误区数学建模竞赛模型和算法只占一部分论文表述同样至关重要。针对B题以下几个点是区分普通论文和优秀论文的关键。5.1 模型假设的清晰与合理必须明确列出所有假设并说明其合理性。对于本题关键假设包括海水声速均匀恒定影响可忽略。海底坡度均匀且为平面简化模型便于分析。波束脚印形状近似为矩形或椭圆用于简化重叠率计算。测量船沿直线匀速航行且定位无偏差。忽略地球曲率对于小范围海域合理。 这些假设不仅要说出来最好能在论文开头的“问题重述与假设”部分用列表形式清晰呈现。5.2 符号说明的规范性使用规范的数学符号。例如(D): 海水深度 (米)(\theta): 波束开角 (度或弧度)(\eta): 最小重叠率 (无量纲)(d): 相邻测线间距 (米)(W): 波束覆盖宽度 (米)(\alpha): 测线方向与坡面方向的夹角 (度)(\beta): 海底坡度 (度) 所有符号在首次出现时注明并可在附录或正文前集中列出。5.3 对“最优解”的深入讨论不要仅仅给出一个数字。要讨论你的解为什么是最优的。对于问题一论证为什么测线平行于长边是最优的或证明两种方向对比后哪个更优。可以画出示意图展示测线布局。对于问题二不仅要给出 (\alpha0) 最优的结论还要解释其物理意义为什么平行于等深线时效率最高从几何和测量质量两个角度阐述。灵敏度分析这是重要的加分项。讨论如果关键参数如深度 (D)、坡度 (\beta)、重叠率要求 (\eta)发生微小变化你的最优方案总长度、最优夹角会如何变化这体现了模型的稳健性。例如可以计算深度增加10%对总长度的影响。5.4 常见误区与扣分点重叠率模型错误如前所述这是最致命的错误。务必理解其物理意义和一维线性本质。忽略单位换算题目中海域尺寸是海里深度是米。计算时必须统一单位通常全换算为米最后结果再根据需要转回海里。单位混乱会直接导致结果错误。对坡度影响处理过于简单或完全忽略在问题二中如果仍然使用平坦海底的宽度公式会被认为没有理解题目深意。即使推导不出精确的 (W(\alpha)) 公式也应通过文字和示意图定性分析坡度带来的影响并提出合理的简化计算思路。缺乏示意图一道几何性这么强的题目论文中必须有清晰的示意图。包括波束几何示意图、测线布局俯视图、坡度影响下的三维关系图。一图胜千言。模型评价与推广部分空洞不要只说“模型优点很多缺点很少”。要具体。优点可以是模型简洁、物理意义清晰、计算效率高。缺点可以是假设海底为平面与实际复杂地形不符忽略了海流对船迹的影响等。推广可以提到模型可应用于机载激光雷达测绘、无人机航拍路径规划等类似条带覆盖问题。6. 从解题到竞赛给参赛同学的建议回顾2023年B题它是一道非常“正”的建模题没有偏难怪的算法但非常考验基本功阅读理解、物理建模、数学推导、编程实现、论文写作。针对这类题型我的建议是组队与分工要明确且互补。团队里最好有一个人擅长快速理解物理背景并建立数学模型通常由数学或物理专业的同学担任一个人擅长编程实现和数值计算计算机、信工专业一个人擅长写作和绘制图表任何专业但需细心、文笔好。在三天时间里定期同步进度防止有人“掉队”或方向跑偏。第一天上午必须吃透题目。不要急着敲代码。花2-3个小时所有人一起反复读题在白板上画图讨论每一个概念如重叠率、波束开角、坡度方向到底指什么直到达成一致理解。这是避免方向性错误的最重要环节。论文写作要贯穿始终。不要等到最后一天晚上才开始写。从第一天确定模型开始就同步撰写“问题分析”、“模型假设”、“符号说明”。第二天编程出结果后立即将结果填入论文并撰写“模型求解”部分。最后一天主要用于打磨摘要、进行灵敏度分析、撰写模型评价与推广以及整体的格式、图表美化。摘要一定要最后写精炼地概括问题、方法、模型、结果和结论。结果的可视化至关重要。对于本题至少要有的图波束几何示意图、矩形海域测线布局图、覆盖宽度/总长度随夹角α变化的关系曲线图。使用Matplotlib、Visio或PPT如果擅长绘制清晰、专业的图表并配上详细的图注。保持良好心态。竞赛期间肯定会遇到困难比如公式推不下去、程序跑不出结果。这时不要慌张团队内部快速讨论如果一条路走不通及时调整简化模型。记住一个完整且自洽的简化模型远胜过一个残缺不全的复杂模型。完成比完美更重要。这道“多波束测线问题”很好地体现了数学建模的核心用数学工具解决实际问题。它不需要你知道多么高深的算法但需要你耐心、细致、逻辑清晰并能将解决方案有效地传达给他人。希望这篇分析能帮助未来参赛的同学不仅看懂这道题更能掌握处理这类问题的方法论。