简介本资源面向2025年电工杯数学建模竞赛参赛队伍及建模学习者聚焦B题‘城市垃圾分类运输的路径优化与调度’这一现实痛点问题提供从建模思路、算法实现到成果呈现的一站式解决方案。压缩包共含多类核心文件包括Word格式无水印论文含完整解题逻辑、多目标模型构建、GIS路径仿真与结果分析、Python与MATLAB双版本可运行代码覆盖数据清洗、车辆路径VRP建模、遗传/模拟退火算法求解及动态可视化、结构化结果表格与原始/处理后数据集整体大小为910.11MB。已有255人下载学习适用于需快速掌握复杂物流调度建模流程、复现高分论文结构、理解多约束条件如分类时效性、载重限制、时间窗下路径协同优化机制的中高级建模者。所有代码模块化注释详尽论文内容符合竞赛规范支持直接提交或针对性修改显著降低备赛门槛与试错成本。 拿到2025年电工杯B题的时候我心里其实先松了一口气——城市垃圾分类运输的路径优化与调度这个题目听起来很唬人但本质上就是我们一直在练的车辆路径问题VRP换了一层场景外衣。真正让人头疼的不是建模本身而是怎么把垃圾分类的行业约束翻译进数学模型以及怎么在有限时间内跑出一套能自圆其说、有收敛趋势的算法结果。这道题把三个子问题揉在了一起垃圾产生量的数据处理、收运车辆的多路径规划、以及全局的车辆调度决策。三个子问题环环相扣数据量大、约束杂、评价指标多稍有遗漏就会整体跑偏。这篇文章我不打算讲那种从零开始的建模教程而是直接复盘我拿到题目之后从拆题到建模、从算法到代码、从结果到论文的完整链路。中间会穿插大量我个人踩过的坑和觉得特别值得注意的细节希望能给正在备赛或者打算参加建模竞赛的同学一个可复用的打法。1. 拿到题目后我先把垃圾分类运输翻译成了数学模型1.1 本质是带容量约束和时间窗的VRP变种垃圾分类运输表面看是民生工程但从优化角度说它就是标准的车辆路径问题。赛题通常会给出一批垃圾收运点每个点有坐标、有各类垃圾的产生量给出处理厂或中转站的位置给出车辆的数量、载重、工作时间等参数。要做的事就是安排几辆车、各自去收哪些点、按什么顺序收最后回到处理厂使总成本最低同时满足所有约束。这个描述一出来参加过建模的同学应该已经条件反射了——这就是经典的带容量约束的车辆路径问题CVRP。再加上我推测赛题里会包含收运时间窗口比如厨余垃圾需要日产日清、有些小区只允许夜间收运所以它又升级成了带时间窗的车辆路径问题VRPTW。如果再算上垃圾分成四五个品类不同品类不能用同一辆车混装那就是多车型、多品类的异质车队VRP。别看题目包装得花里胡哨把问题归到VRP家族之后整个解题思路就清晰了一大半。后面所有的工作——数据预处理、模型构建、算法设计都围绕车辆-路径-时间窗-品类这四个关键词展开。1.2 子问题解耦预测、规划、调度三个层面分开处理这类题目千万不要一上来就想建一个全能模型把所有东西同时优化。因为以竞赛的时间体量同时优化所有层级的决策会让模型复杂度爆炸算法根本收敛不动。我的做法是把问题拆成三个层面第一层垃圾量数据的处理与预测。如果赛题给了历史垃圾量时间序列可能需要预测未来某一时段各收运点的垃圾产生量。如果只给了一个静态数值那这层就简化为数据清洗和统计分析。第二层收运路径规划。已知每个点的垃圾量或预测值、车辆容量和各类约束确定每辆车的访问顺序。这是整个题目的核心也是算法最能拉开差距的地方。第三层车辆调度决策。在路径规划的基础上决定每辆车几点出发、是否执行多趟任务以及如何处理突发的新增需求。这一层在竞赛中可以适当简化重点放在能用约束说清楚。把问题解耦之后每一层的目标就非常明确了。我在论文里的写法也是按这个逻辑递进的先交代数据怎么处理再建立路径优化模型最后把调度规则作为一个扩展场景来讨论。评阅人看了也能一眼看出你的思路是清晰的。1.3 赛题里最容易被忽略的隐藏约束这类题目有个很坑的地方很多约束不会直接写出来而是藏在数据的编排方式里。比如不同垃圾品类是分开统计的每个收运点某一类垃圾的产生量为0很正常说明这个点不需要该类别的车去收处理厂可能有多个而且不同处理厂只接收特定品类车辆工作时段是有限制的比如早高峰时段不能上路有些收运点之间可能存在单向道路导致距离矩阵不对称。这些都是我在读题和翻数据时特别注意的地方。如果不做这一步直接套标准CVRP模板大概率会在细节上翻车。建议所有做这道题的同学拿到数据后先花半天时间做探索性数据分析把每一个字段的含义、取值范围、缺失情况都摸一遍再做建模。2. 数据坑位盘点这题90%的坑都在数据理解上2.1 坐标转换和距离计算矩阵赛题给的位置信息有的是经纬度有的是平面坐标有的甚至只给一个编号对应到地图上的点。这个处理直接决定整个算法的输入质量。如果你拿到的是经纬度坐标记得用Haversine公式计算球面距离而不是简单套欧氏距离。尤其在城市尺度下经纬度一度对应的地面距离在不同纬度差异很大直接用欧氏距离算出来的路径完全不可信。我一般会写一个距离矩阵工具函数输入点坐标列表输出任意两点间的距离矩阵和行驶时间矩阵。如果赛题给了道路网络信息那就要用最短路算法比如Floyd或Dijkstra先算出实际路网距离再用这个距离矩阵作为后续算法的输入。import numpy as np import math def haversine_dist(lon1, lat1, lon2, lat2): R 6371.0 # 地球半径单位km dlon math.radians(lon2 - lon1) dlat math.radians(lat2 - lat1) a math.sin(dlat / 2) ** 2 math.cos(math.radians(lat1)) * \ math.cos(math.radians(lat2)) * math.sin(dlon / 2) ** 2 return 2 * R * math.asin(math.sqrt(a)) def build_distance_matrix(points): n len(points) dist_mat np.zeros((n, n)) for i in range(n): for j in range(n): if i ! j: dist_mat[i][j] haversine_dist( points[i][lon], points[i][lat], points[j][lon], points[j][lat] ) return dist_mat这里有个小经验车辆行驶时间不只是距离除以速度那么简单城市工况下红灯、堵车会导致平均速度远低于道路限速。我在模型里给时间估算留了一个安全系数比如平均车速取20-25km/h并且设置一个固定的装卸服务时间参数比如每点10分钟这样算出来的总时间会更符合实际。2.2 垃圾量数据的周期特征如果赛题给了历史垃圾量数据那处理的时候要充分考虑周期性。城市垃圾产量通常有明显的周内波动和节假日效应——周末厨余垃圾和生活垃圾往往比工作日多节假日前后的可回收物也会冲高。这些规律在做预测和后续调度时非常关键。我拿到这类数据后第一件事不是急着建模而是画时间序列图按星期几分组看箱线图看有没有明显异常值。如果有异常点需要判断是记录错误还是真实的突发事件比如大型活动后的垃圾量激增。如果是记录错误可以直接剔除或修正如果是真实事件把它作为特殊需求单独处理。在此基础上我会用一些相对简单但稳健的方法做预测——移动平均、指数平滑、或者小规模的时间序列回归。除非数据量特别大且有明显复杂的非线性模式否则不建议在预测环节用太复杂的深度学习模型因为竞赛时间有限而且后面算法环节才是拿分重点。预测步骤只要能给出合理、可解释、误差可接受的数值就可以了。2.3 数据的量级决定了算法的选择赛题点数不同解题策略完全不同。我遇到过几十个点的小规模数据和几百个点的大规模数据两种情况。点少的时候可以用精确算法或者穷举思路做验证算法性能不是瓶颈点多的时候必须上启发式算法否则计算时间完全没法接受。一个很实际的自测标准如果收运点在50个以内精确算法比如分支定界或者Gurobi直接求解MILP模型完全可以得到一个全局最优或近似最优解如果超过100个点直接上Gurobi很容易卡死这时候就要把重心放在设计高效的元启发式算法上。我在做题前会根据数据量快速判断该走哪条路避免在一开始就选错算法路线。3. 模型构建目标函数和约束条件的取舍策略3.1 目标函数三层级设计垃圾分类运输的优化目标表面上是总成本最低但细拆起来其实是多层级的。我的做法是设置一个分层目标函数不同层级的优先级不同第一优先级使用的车辆数量最少。在车辆路径问题里减少一辆车节省的固定成本购车、司机、维护通常远大于缩短行驶距离带来的可变成本节省。所以第一目标一定是尽量减少出车的台数。第二优先级总行驶距离最短。在车辆数固定的前提下尽可能让所有车辆的总行驶里程最小对应的是油耗、磨损、司机工时等可变成本。第三优先级各个车辆的负载均衡和时间均衡。这一点评阅人很看重因为如果算法算出某一辆车干重活、另一辆车只收两个点就回来了在实际场景中是不可用的方案。我会在目标函数里加一个均衡性的惩罚项或者在算法的每一代评估中加一个负载均衡的指标。在论文里三个优先级可以写成加权求和的形式也可以写成分层优化的形式。如果求加权和注意权重系数要设置合理否则会出现车辆数多但距离短的方案比车辆数少但距离长的方案得分更高的反直觉结果。我采用的是先按车辆数排序再按总距离排序的字典序优化思路更稳妥。3.2 约束条件的合理取舍模型里的约束条件不是越多越好因为每加一个约束求解难度都会上一个台阶。我在建模时把约束分成了三类第一类是硬性约束违反就不可能执行。比如容量约束——任何一辆车装载的垃圾量不能超过车辆载重上限时间窗约束——必须在规定的时间窗口内到达收运点品类匹配约束——某类垃圾只能用对应的专用车辆运输车辆路径闭合约束——车辆从处理厂出发最终必须回到处理厂。第二类是可以软化的约束用一个惩罚项代替。比如时间窗可以设置一个软时间窗——允许迟到但是每迟到一分钟有一个惩罚成本这样算法在搜索时更有弹性避免因为个别点的约束太紧导致整个解不可行。实际项目中这是非常常用的手段。第三类是可选的博弈性约束。比如每辆车每天最多跑两趟这类限制如果写进模型会大大增加调度层面的复杂度。我在初版模型里一般不写死而是在路径层解完之后再做可达性验证如果确实需要车辆复用再单独扩展。3.3 决策变量的表示方法VRP类模型的决策变量通常是0-1变量x_ijk表示车辆k是否从点i行驶到点j。这个变量在数学上很标准但如果你用的是遗传算法等启发式方法实际上并不需要显式定义这个变量——你只需要一个染色体编码来表示每辆车的服务顺序解码的时候就能自动得到对应的路径集合和总成本。我在代码里主要用的是基于路径的编码方式用一串整数表示车辆访问的收运点顺序用分隔符比如0区分不同车辆的任务段。这种编码方式直观、容易实现交叉变异解码速度也快。4. 主算法选型为什么是节约算法构造解 遗传算法加局部搜索4.1 精确算法的天花板很多第一次做竞赛的同学会想着用Gurobi或者CPLEX直接求解MILP模型但VRP是NP-hard问题点一多就完全跑不动。即便只有50个点完整的VRPTW模型的整数变量数量也会达到几千甚至上万精确求解器可能要跑几小时甚至数天。竞赛只有三天时间不可能靠精确算法拿结果。我的建议是可以用Gurobi求解小规模算例来做结果验证和对比但主算法一定得是启发式。4.2 构造初始解改进型节约算法C-W算法好的启发式算法需要一个靠谱的初始解。如果初始解质量太差后续的进化优化会花大量时间在修复而不是改进上。VRP里最经典的初始解构造方法是节约算法Clarke-Wright Savings Algorithm核心思想很直观先给每个收运点单独派一辆车去服务然后计算任意两个点合并到同一条路线时能节省的距离按节约值从大到小排序逐步合并路线直到容量或时间窗约束不允许为止。def clarke_wright_savings(dist_mat, demands, capacity): # 初始化每个点单独一条路径 routes [[i] for i in range(1, len(dist_mat))] # 计算节约值 savings [] for i in range(1, len(dist_mat)): for j in range(i 1, len(dist_mat)): saving dist_mat[0][i] dist_mat[0][j] - dist_mat[i][j] savings.append((saving, i, j)) savings.sort(reverseTrue) # 按节约值从大到小合并路径注意检查容量约束 ...在实际实现时我不会完全按照经典C-W算法走而是做两个改进一是在合并时加入时间窗检查如果合并后到达时间不满足任何一端的点的时间窗要求就放弃这次合并二是合并时尝试多个插入位置不仅限于把两条路径的首尾相接这样可以进一步压缩距离。4.3 为什么用遗传算法而不是模拟退火或粒子群VRP的搜索空间是离散且组合爆炸的适合用群体智能算法。遗传算法、模拟退火、粒子群、蚁群都可以做。我个人更推荐遗传算法GA原因是它的并行性天然适合VRP这种多路径结构种群中的每个个体就是一套完整的路径方案交叉变异操作可以直接在编码层面操作路线片段非常契合VRP的解结构。当然纯遗传算法也容易早熟收敛所以我的方案是遗传算法 局部搜索的混合框架也就是常说的Memetic Algorithm。每一代在所有个体完成交叉变异之后对最优的几个个体做局部搜索强化——用2-opt交换路径内两条边的连接方式、Or-opt移动路径中的一小段到其他位置、Swap交换两条路径中的两个点这几个邻域算子反复搜索直到没有改进为止。4.4 参数设置的具体经验值我这次用的参数范围可以分享给大家参考种群规模100-200交叉概率0.8-0.9变异概率0.1-0.2迭代代数300-500精英保留数5-10局部搜索频率每代对最优的3-5个个体执行这些参数不是拍脑袋定的。种群太小容易早熟太大会拖慢每一代的计算速度迭代代数太少还没收敛太多则后期几乎没有改进纯浪费计算资源。建议做一组小规模的参数敏感性测试画出收敛曲线观察一下再最终确定。5. 代码实现数据结构、算法流程与关键代码5.1 数据结构设计这道题的数据结构设计直接影响后续算法的开发效率。我用类的方式封装节点、车辆和路径让代码更清晰也方便后期扩展。class Node: def __init__(self, node_id, category, amount, x, y, tw_min0, tw_max240): self.id node_id self.category category self.amount amount self.x x self.y y self.tw_min tw_min # 最早服务时间 self.tw_max tw_max # 最晚服务时间 class Vehicle: def __init__(self, vehicle_id, capacity, category0, fixed_cost1.0): self.id vehicle_id self.capacity capacity self.category category # 0表示通用1/2/3/4表示对应品类 self.fixed_cost fixed_cost class Route: def __init__(self, vehicle): self.vehicle vehicle self.nodes [0] # 从depot出发 self.load 0 self.cost 0 def add_node(self, node): self.nodes.append(node.id) self.load node.amount self.cost ... # 根据距离矩阵累加把每个点的时间窗作为Node的一部分有一个好处后续做解码和时间窗校验时直接遍历route.nodes查对应Node的属性即可不需要额外维护复杂的映射表。这个设计在后面大量调试中帮我省了很多时间。5.2 遗传算法的整体流程完整算法流程我用一个框架表示读取并预处理数据计算距离矩阵与行驶时间矩阵。对每个垃圾品类分别构造初始路径集合因为不同品类不能混装。将初始路径集合编码成染色体初始化种群。循环迭代计算每个个体的适应度总成本 惩罚项执行锦标赛选择执行交叉、变异执行精英保留对精英个体做局部搜索。迭代结束后取最优个体解码为路径方案。后处理计算各项指标总车辆数、总距离、平均装载率、时间窗违反率生成结果表。5.3 遗传算子设计细节编码方式我采用的是基于整数序列的自然数编码。比如有10个收运点、3辆车一个染色体可以表示为[2, 5, 0, 1, 3, 6, 0, 4, 7, 8, 9]其中0是分隔符代表车辆切换。但这里有个细节不同品类的收运点混合编码会导致解码时出现某辆车拉错品类的错误。所以我给染色体加了一个品类段标记确保每条路径段只包含同品类节点。交叉算子我用的是顺序交叉Order Crossover和部分匹配交叉PMX两种交替使用。PMX适合处理有约束的路径段交换顺序交叉则更擅长保留节点的相对顺序。变异算子则包括交换变异随机交换两个基因位、插入变异把一个点插入到另一个位置和逆转变异反转一段子路径。每次变异后要立即做容量和时间窗可行性校验如果不可行就放弃这次变异重新选一段尝试。下面给一个交叉算子的简化示例展示PMX的核心逻辑def pmx_crossover(p1, p2): size len(p1) # 随机选择两个交叉点 start, end sorted(random.sample(range(size), 2)) child [-1] * size # 复制父本1的交叉区间 child[start:end1] p1[start:end1] # 填充其余位置先找父本2中不在child区间内的基因 for i in range(size): if child[i] -1: candidate p2[i] while candidate in child[start:end1]: candidate p2[p1.index(candidate)] child[i] candidate return child这个实现的逻辑比较简洁但要注意在真实代码里编码里的0分隔符会干扰交叉操作所以交叉前需要把0去掉、按路径段拼接或者对分隔符做特殊处理。我处理的办法是分层编码染色体只编码每个车辆段的节点序列车辆分配另用一个数组表示这样交叉时车辆分配和节点顺序就能分开进化。5.4 局部搜索的实现要点局部搜索是提升解质量的关键。我用三个算子2-opt在一条路径内部选择两条不相交的边断开后反向连接。这个算子在路径优化里极其经典几乎适用于所有VRP变种。具体实现时把路径看成一个环形结构从depot出发又回到depot所以需要对中间段做反转而不是简单交换两个点。Or-opt把一条路径中的连续1到3个点取出来插入到这条路径或另一条路径的其他位置。这个操作对改善局部拥堵非常有效实现起来也不复杂。Swap交换两条路径中的各一个点或者一条路径内的两个点。适用于车辆间的负载均衡调整。局部搜索的停止条件是连续k次迭代没有改进k我一般设为100。每次局部搜索后重新计算路径的容量和总成本如果有改进就更新。def two_opt(route, dist_mat): improved True best_route route[:] best_cost compute_cost(route, dist_mat) while improved: improved False for i in range(1, len(route) - 2): for j in range(i 2, len(route)): new_route route[:i] route[i:j][::-1] route[j:] new_cost compute_cost(new_route, dist_mat) if new_cost best_cost - 1e-9: best_route new_route best_cost new_cost improved True route best_route return best_route, best_cost在实际代码里还有个加速优化不需要每次重新计算整条路径的成本只需要计算断点附近的那几段边的变化量即可。但对竞赛规模的数据来说全量重算是可以接受的代码简洁优先。6. 结果怎么呈现才能对得起你的努力6.1 收敛曲线与分析算法跑完第一件事是画收敛曲线和路径图。收敛曲线要能证明你的算法是正常迭代的而不是一开始就乱跳到末尾。具体来说画出每一代的最优适应度和平均适应度两条曲线最优适应度应该是单调下降或先快速下降后趋于平稳平均适应度围绕最优值波动但整体也呈下降趋势。如果最优曲线在早期就完全不动了说明算法大概率陷入了局部最优就要考虑调整变异概率或者增加局部搜索力度。路径图方面我会画一个带坐标的地图散点图用不同颜色区分不同车辆的路线并标注depot和各收运点。这张图一定要画得清晰因为评阅人大概率第一眼就看结果图路径规划效果好不好一目了然。如果路径交叉严重即便总距离指标不错印象分也会打折扣。我的经验是画图时把od的标记放大、路线线条加粗、图例注明车辆编号和车辆载重利用率让每辆车跑的全貌一看就懂。6.2 对比实验你至少要做三组竞赛评阅里对比实验是拉开档次的关键。光说自己算法好是不够的得在同一套数据上对比才行。我做了三组对比第一组无算法调度当前人工方案vs 改进后优化方案。如果题目数据里给了类似现有方案的参考值这组对比最直观直接说明你的方法能省多少钱、少用多少车、少跑多少公里。如果题目没有给可以用一个简单的按距离最近邻贪心构造的基线方案做对比体现算法优化的必要性和效果。第二组单一算法 vs 你的改进算法。比如纯遗传算法和GA局部搜索在相同迭代次数下的对比重点展示混合算法在收敛速度和解的质量上的优势。第三组小规模精确解 vs 启发式解。如果收运点数较少用Gurobi求出精确最优解再和你的启发式结果对比说明启发式算法取得的解和最优解的差距在可接受范围内比如5%以内。这组实验最有说服力直接给评审一个量化感知你的算法离最优解还有多远。6.3 灵敏度分析怎么设计灵敏度分析是很多同学容易忽略但特别加分的部分。我做了三组灵敏度分析车辆载重变动-20%、-10%、10%、20%对总车辆数和总距离的影响时间窗松紧程度扩大或缩小20%对可行性的影响垃圾量预测误差±10%对方案鲁棒性的影响。这些分析一方面展示了对模型的深入理解另一方面也回应了现实场景中参数波动的风险。如果题目给了多个场景的数据比如不同片区或不同季节的数据还可以做多场景的对比分析说明算法在不同规模、不同分布下的泛化能力。这是拿创新点分的好地方。6.4 论文每部分的写作要点摘要部分要四段式第一段用两三句话说清楚背景和问题是什么第二段写你建立了什么模型、用了什么方法第三段给出核心结果用数据说话第四句轻轻一句点出方案的实践价值和可推广性。摘要里的数一定不能拍脑袋要和正文结果表完全对应评阅人会专门抽查。问题分析和模型建立部分重点写清楚为什么这样建模和哪些假设是不合理的。不要只罗列数学公式要交代每一条约束的现实意义——比如车辆装载量不能超过额定载重对应的是车辆行驶安全。符号说明要全公式编号要有序这是细节分。算法设计部分是评阅人重点看的段落。建议给每一个算法模块配上流程图或伪代码伪代码的每一行都要和后面代码核心逻辑对应这样审稿人读起来会非常顺。参数表要单独列出来包括参数名称、含义、取值、调节依据。结果部分要注意图文并茂。一张总结果对比表、一张路径可视化图、一张收敛曲线图三个基础图必须有。好的论文在结果部分会很克制每个数据都有明确的作用不堆砌。7. 给下一届选手的几条实在建议代码规范这块从第一天就按工程标准写。变量名别用a、b、c这种文件夹按data、src、results、figures分好结果输出统一用CSV保存。比赛的第三天晚上你可能要改十几次参数重跑几十组实验如果代码一团乱光debug就能耗掉你半天时间。时间分配上我的建议是第一天上午快速读题和数据下午完成数据预处理和基线模型第二天全力做算法和实验晚上开始写论文框架第三天上午补全实验和分析下午集中写完整论文、统一格式、校对摘要。千万不要花两天时间搭模型最后只剩半天写论文——那样你再好的算法也白搭。最后说一个我个人特别受益的习惯做任何模型改动前先把当前版本的结果完整备份好并记录下来改了什么、效果如何。这样在论文讨论部分写我们尝试了xxx方法结果yyy的时候每一步都有真实依据而不是靠回忆拼凑。建模竞赛比的从来不只是算法它比的是谁在有限时间里更系统地解决问题。希望这篇复盘能帮你少走一些弯路祝你比赛顺利。本文还有配套的精品资源点击获取 垃圾分类收运路径优化全解析：从VRP建模到遗传算法求解实战
简介本资源面向2025年电工杯数学建模竞赛参赛队伍及建模学习者聚焦B题‘城市垃圾分类运输的路径优化与调度’这一现实痛点问题提供从建模思路、算法实现到成果呈现的一站式解决方案。压缩包共含多类核心文件包括Word格式无水印论文含完整解题逻辑、多目标模型构建、GIS路径仿真与结果分析、Python与MATLAB双版本可运行代码覆盖数据清洗、车辆路径VRP建模、遗传/模拟退火算法求解及动态可视化、结构化结果表格与原始/处理后数据集整体大小为910.11MB。已有255人下载学习适用于需快速掌握复杂物流调度建模流程、复现高分论文结构、理解多约束条件如分类时效性、载重限制、时间窗下路径协同优化机制的中高级建模者。所有代码模块化注释详尽论文内容符合竞赛规范支持直接提交或针对性修改显著降低备赛门槛与试错成本。 拿到2025年电工杯B题的时候我心里其实先松了一口气——城市垃圾分类运输的路径优化与调度这个题目听起来很唬人但本质上就是我们一直在练的车辆路径问题VRP换了一层场景外衣。真正让人头疼的不是建模本身而是怎么把垃圾分类的行业约束翻译进数学模型以及怎么在有限时间内跑出一套能自圆其说、有收敛趋势的算法结果。这道题把三个子问题揉在了一起垃圾产生量的数据处理、收运车辆的多路径规划、以及全局的车辆调度决策。三个子问题环环相扣数据量大、约束杂、评价指标多稍有遗漏就会整体跑偏。这篇文章我不打算讲那种从零开始的建模教程而是直接复盘我拿到题目之后从拆题到建模、从算法到代码、从结果到论文的完整链路。中间会穿插大量我个人踩过的坑和觉得特别值得注意的细节希望能给正在备赛或者打算参加建模竞赛的同学一个可复用的打法。1. 拿到题目后我先把垃圾分类运输翻译成了数学模型1.1 本质是带容量约束和时间窗的VRP变种垃圾分类运输表面看是民生工程但从优化角度说它就是标准的车辆路径问题。赛题通常会给出一批垃圾收运点每个点有坐标、有各类垃圾的产生量给出处理厂或中转站的位置给出车辆的数量、载重、工作时间等参数。要做的事就是安排几辆车、各自去收哪些点、按什么顺序收最后回到处理厂使总成本最低同时满足所有约束。这个描述一出来参加过建模的同学应该已经条件反射了——这就是经典的带容量约束的车辆路径问题CVRP。再加上我推测赛题里会包含收运时间窗口比如厨余垃圾需要日产日清、有些小区只允许夜间收运所以它又升级成了带时间窗的车辆路径问题VRPTW。如果再算上垃圾分成四五个品类不同品类不能用同一辆车混装那就是多车型、多品类的异质车队VRP。别看题目包装得花里胡哨把问题归到VRP家族之后整个解题思路就清晰了一大半。后面所有的工作——数据预处理、模型构建、算法设计都围绕车辆-路径-时间窗-品类这四个关键词展开。1.2 子问题解耦预测、规划、调度三个层面分开处理这类题目千万不要一上来就想建一个全能模型把所有东西同时优化。因为以竞赛的时间体量同时优化所有层级的决策会让模型复杂度爆炸算法根本收敛不动。我的做法是把问题拆成三个层面第一层垃圾量数据的处理与预测。如果赛题给了历史垃圾量时间序列可能需要预测未来某一时段各收运点的垃圾产生量。如果只给了一个静态数值那这层就简化为数据清洗和统计分析。第二层收运路径规划。已知每个点的垃圾量或预测值、车辆容量和各类约束确定每辆车的访问顺序。这是整个题目的核心也是算法最能拉开差距的地方。第三层车辆调度决策。在路径规划的基础上决定每辆车几点出发、是否执行多趟任务以及如何处理突发的新增需求。这一层在竞赛中可以适当简化重点放在能用约束说清楚。把问题解耦之后每一层的目标就非常明确了。我在论文里的写法也是按这个逻辑递进的先交代数据怎么处理再建立路径优化模型最后把调度规则作为一个扩展场景来讨论。评阅人看了也能一眼看出你的思路是清晰的。1.3 赛题里最容易被忽略的隐藏约束这类题目有个很坑的地方很多约束不会直接写出来而是藏在数据的编排方式里。比如不同垃圾品类是分开统计的每个收运点某一类垃圾的产生量为0很正常说明这个点不需要该类别的车去收处理厂可能有多个而且不同处理厂只接收特定品类车辆工作时段是有限制的比如早高峰时段不能上路有些收运点之间可能存在单向道路导致距离矩阵不对称。这些都是我在读题和翻数据时特别注意的地方。如果不做这一步直接套标准CVRP模板大概率会在细节上翻车。建议所有做这道题的同学拿到数据后先花半天时间做探索性数据分析把每一个字段的含义、取值范围、缺失情况都摸一遍再做建模。2. 数据坑位盘点这题90%的坑都在数据理解上2.1 坐标转换和距离计算矩阵赛题给的位置信息有的是经纬度有的是平面坐标有的甚至只给一个编号对应到地图上的点。这个处理直接决定整个算法的输入质量。如果你拿到的是经纬度坐标记得用Haversine公式计算球面距离而不是简单套欧氏距离。尤其在城市尺度下经纬度一度对应的地面距离在不同纬度差异很大直接用欧氏距离算出来的路径完全不可信。我一般会写一个距离矩阵工具函数输入点坐标列表输出任意两点间的距离矩阵和行驶时间矩阵。如果赛题给了道路网络信息那就要用最短路算法比如Floyd或Dijkstra先算出实际路网距离再用这个距离矩阵作为后续算法的输入。import numpy as np import math def haversine_dist(lon1, lat1, lon2, lat2): R 6371.0 # 地球半径单位km dlon math.radians(lon2 - lon1) dlat math.radians(lat2 - lat1) a math.sin(dlat / 2) ** 2 math.cos(math.radians(lat1)) * \ math.cos(math.radians(lat2)) * math.sin(dlon / 2) ** 2 return 2 * R * math.asin(math.sqrt(a)) def build_distance_matrix(points): n len(points) dist_mat np.zeros((n, n)) for i in range(n): for j in range(n): if i ! j: dist_mat[i][j] haversine_dist( points[i][lon], points[i][lat], points[j][lon], points[j][lat] ) return dist_mat这里有个小经验车辆行驶时间不只是距离除以速度那么简单城市工况下红灯、堵车会导致平均速度远低于道路限速。我在模型里给时间估算留了一个安全系数比如平均车速取20-25km/h并且设置一个固定的装卸服务时间参数比如每点10分钟这样算出来的总时间会更符合实际。2.2 垃圾量数据的周期特征如果赛题给了历史垃圾量数据那处理的时候要充分考虑周期性。城市垃圾产量通常有明显的周内波动和节假日效应——周末厨余垃圾和生活垃圾往往比工作日多节假日前后的可回收物也会冲高。这些规律在做预测和后续调度时非常关键。我拿到这类数据后第一件事不是急着建模而是画时间序列图按星期几分组看箱线图看有没有明显异常值。如果有异常点需要判断是记录错误还是真实的突发事件比如大型活动后的垃圾量激增。如果是记录错误可以直接剔除或修正如果是真实事件把它作为特殊需求单独处理。在此基础上我会用一些相对简单但稳健的方法做预测——移动平均、指数平滑、或者小规模的时间序列回归。除非数据量特别大且有明显复杂的非线性模式否则不建议在预测环节用太复杂的深度学习模型因为竞赛时间有限而且后面算法环节才是拿分重点。预测步骤只要能给出合理、可解释、误差可接受的数值就可以了。2.3 数据的量级决定了算法的选择赛题点数不同解题策略完全不同。我遇到过几十个点的小规模数据和几百个点的大规模数据两种情况。点少的时候可以用精确算法或者穷举思路做验证算法性能不是瓶颈点多的时候必须上启发式算法否则计算时间完全没法接受。一个很实际的自测标准如果收运点在50个以内精确算法比如分支定界或者Gurobi直接求解MILP模型完全可以得到一个全局最优或近似最优解如果超过100个点直接上Gurobi很容易卡死这时候就要把重心放在设计高效的元启发式算法上。我在做题前会根据数据量快速判断该走哪条路避免在一开始就选错算法路线。3. 模型构建目标函数和约束条件的取舍策略3.1 目标函数三层级设计垃圾分类运输的优化目标表面上是总成本最低但细拆起来其实是多层级的。我的做法是设置一个分层目标函数不同层级的优先级不同第一优先级使用的车辆数量最少。在车辆路径问题里减少一辆车节省的固定成本购车、司机、维护通常远大于缩短行驶距离带来的可变成本节省。所以第一目标一定是尽量减少出车的台数。第二优先级总行驶距离最短。在车辆数固定的前提下尽可能让所有车辆的总行驶里程最小对应的是油耗、磨损、司机工时等可变成本。第三优先级各个车辆的负载均衡和时间均衡。这一点评阅人很看重因为如果算法算出某一辆车干重活、另一辆车只收两个点就回来了在实际场景中是不可用的方案。我会在目标函数里加一个均衡性的惩罚项或者在算法的每一代评估中加一个负载均衡的指标。在论文里三个优先级可以写成加权求和的形式也可以写成分层优化的形式。如果求加权和注意权重系数要设置合理否则会出现车辆数多但距离短的方案比车辆数少但距离长的方案得分更高的反直觉结果。我采用的是先按车辆数排序再按总距离排序的字典序优化思路更稳妥。3.2 约束条件的合理取舍模型里的约束条件不是越多越好因为每加一个约束求解难度都会上一个台阶。我在建模时把约束分成了三类第一类是硬性约束违反就不可能执行。比如容量约束——任何一辆车装载的垃圾量不能超过车辆载重上限时间窗约束——必须在规定的时间窗口内到达收运点品类匹配约束——某类垃圾只能用对应的专用车辆运输车辆路径闭合约束——车辆从处理厂出发最终必须回到处理厂。第二类是可以软化的约束用一个惩罚项代替。比如时间窗可以设置一个软时间窗——允许迟到但是每迟到一分钟有一个惩罚成本这样算法在搜索时更有弹性避免因为个别点的约束太紧导致整个解不可行。实际项目中这是非常常用的手段。第三类是可选的博弈性约束。比如每辆车每天最多跑两趟这类限制如果写进模型会大大增加调度层面的复杂度。我在初版模型里一般不写死而是在路径层解完之后再做可达性验证如果确实需要车辆复用再单独扩展。3.3 决策变量的表示方法VRP类模型的决策变量通常是0-1变量x_ijk表示车辆k是否从点i行驶到点j。这个变量在数学上很标准但如果你用的是遗传算法等启发式方法实际上并不需要显式定义这个变量——你只需要一个染色体编码来表示每辆车的服务顺序解码的时候就能自动得到对应的路径集合和总成本。我在代码里主要用的是基于路径的编码方式用一串整数表示车辆访问的收运点顺序用分隔符比如0区分不同车辆的任务段。这种编码方式直观、容易实现交叉变异解码速度也快。4. 主算法选型为什么是节约算法构造解 遗传算法加局部搜索4.1 精确算法的天花板很多第一次做竞赛的同学会想着用Gurobi或者CPLEX直接求解MILP模型但VRP是NP-hard问题点一多就完全跑不动。即便只有50个点完整的VRPTW模型的整数变量数量也会达到几千甚至上万精确求解器可能要跑几小时甚至数天。竞赛只有三天时间不可能靠精确算法拿结果。我的建议是可以用Gurobi求解小规模算例来做结果验证和对比但主算法一定得是启发式。4.2 构造初始解改进型节约算法C-W算法好的启发式算法需要一个靠谱的初始解。如果初始解质量太差后续的进化优化会花大量时间在修复而不是改进上。VRP里最经典的初始解构造方法是节约算法Clarke-Wright Savings Algorithm核心思想很直观先给每个收运点单独派一辆车去服务然后计算任意两个点合并到同一条路线时能节省的距离按节约值从大到小排序逐步合并路线直到容量或时间窗约束不允许为止。def clarke_wright_savings(dist_mat, demands, capacity): # 初始化每个点单独一条路径 routes [[i] for i in range(1, len(dist_mat))] # 计算节约值 savings [] for i in range(1, len(dist_mat)): for j in range(i 1, len(dist_mat)): saving dist_mat[0][i] dist_mat[0][j] - dist_mat[i][j] savings.append((saving, i, j)) savings.sort(reverseTrue) # 按节约值从大到小合并路径注意检查容量约束 ...在实际实现时我不会完全按照经典C-W算法走而是做两个改进一是在合并时加入时间窗检查如果合并后到达时间不满足任何一端的点的时间窗要求就放弃这次合并二是合并时尝试多个插入位置不仅限于把两条路径的首尾相接这样可以进一步压缩距离。4.3 为什么用遗传算法而不是模拟退火或粒子群VRP的搜索空间是离散且组合爆炸的适合用群体智能算法。遗传算法、模拟退火、粒子群、蚁群都可以做。我个人更推荐遗传算法GA原因是它的并行性天然适合VRP这种多路径结构种群中的每个个体就是一套完整的路径方案交叉变异操作可以直接在编码层面操作路线片段非常契合VRP的解结构。当然纯遗传算法也容易早熟收敛所以我的方案是遗传算法 局部搜索的混合框架也就是常说的Memetic Algorithm。每一代在所有个体完成交叉变异之后对最优的几个个体做局部搜索强化——用2-opt交换路径内两条边的连接方式、Or-opt移动路径中的一小段到其他位置、Swap交换两条路径中的两个点这几个邻域算子反复搜索直到没有改进为止。4.4 参数设置的具体经验值我这次用的参数范围可以分享给大家参考种群规模100-200交叉概率0.8-0.9变异概率0.1-0.2迭代代数300-500精英保留数5-10局部搜索频率每代对最优的3-5个个体执行这些参数不是拍脑袋定的。种群太小容易早熟太大会拖慢每一代的计算速度迭代代数太少还没收敛太多则后期几乎没有改进纯浪费计算资源。建议做一组小规模的参数敏感性测试画出收敛曲线观察一下再最终确定。5. 代码实现数据结构、算法流程与关键代码5.1 数据结构设计这道题的数据结构设计直接影响后续算法的开发效率。我用类的方式封装节点、车辆和路径让代码更清晰也方便后期扩展。class Node: def __init__(self, node_id, category, amount, x, y, tw_min0, tw_max240): self.id node_id self.category category self.amount amount self.x x self.y y self.tw_min tw_min # 最早服务时间 self.tw_max tw_max # 最晚服务时间 class Vehicle: def __init__(self, vehicle_id, capacity, category0, fixed_cost1.0): self.id vehicle_id self.capacity capacity self.category category # 0表示通用1/2/3/4表示对应品类 self.fixed_cost fixed_cost class Route: def __init__(self, vehicle): self.vehicle vehicle self.nodes [0] # 从depot出发 self.load 0 self.cost 0 def add_node(self, node): self.nodes.append(node.id) self.load node.amount self.cost ... # 根据距离矩阵累加把每个点的时间窗作为Node的一部分有一个好处后续做解码和时间窗校验时直接遍历route.nodes查对应Node的属性即可不需要额外维护复杂的映射表。这个设计在后面大量调试中帮我省了很多时间。5.2 遗传算法的整体流程完整算法流程我用一个框架表示读取并预处理数据计算距离矩阵与行驶时间矩阵。对每个垃圾品类分别构造初始路径集合因为不同品类不能混装。将初始路径集合编码成染色体初始化种群。循环迭代计算每个个体的适应度总成本 惩罚项执行锦标赛选择执行交叉、变异执行精英保留对精英个体做局部搜索。迭代结束后取最优个体解码为路径方案。后处理计算各项指标总车辆数、总距离、平均装载率、时间窗违反率生成结果表。5.3 遗传算子设计细节编码方式我采用的是基于整数序列的自然数编码。比如有10个收运点、3辆车一个染色体可以表示为[2, 5, 0, 1, 3, 6, 0, 4, 7, 8, 9]其中0是分隔符代表车辆切换。但这里有个细节不同品类的收运点混合编码会导致解码时出现某辆车拉错品类的错误。所以我给染色体加了一个品类段标记确保每条路径段只包含同品类节点。交叉算子我用的是顺序交叉Order Crossover和部分匹配交叉PMX两种交替使用。PMX适合处理有约束的路径段交换顺序交叉则更擅长保留节点的相对顺序。变异算子则包括交换变异随机交换两个基因位、插入变异把一个点插入到另一个位置和逆转变异反转一段子路径。每次变异后要立即做容量和时间窗可行性校验如果不可行就放弃这次变异重新选一段尝试。下面给一个交叉算子的简化示例展示PMX的核心逻辑def pmx_crossover(p1, p2): size len(p1) # 随机选择两个交叉点 start, end sorted(random.sample(range(size), 2)) child [-1] * size # 复制父本1的交叉区间 child[start:end1] p1[start:end1] # 填充其余位置先找父本2中不在child区间内的基因 for i in range(size): if child[i] -1: candidate p2[i] while candidate in child[start:end1]: candidate p2[p1.index(candidate)] child[i] candidate return child这个实现的逻辑比较简洁但要注意在真实代码里编码里的0分隔符会干扰交叉操作所以交叉前需要把0去掉、按路径段拼接或者对分隔符做特殊处理。我处理的办法是分层编码染色体只编码每个车辆段的节点序列车辆分配另用一个数组表示这样交叉时车辆分配和节点顺序就能分开进化。5.4 局部搜索的实现要点局部搜索是提升解质量的关键。我用三个算子2-opt在一条路径内部选择两条不相交的边断开后反向连接。这个算子在路径优化里极其经典几乎适用于所有VRP变种。具体实现时把路径看成一个环形结构从depot出发又回到depot所以需要对中间段做反转而不是简单交换两个点。Or-opt把一条路径中的连续1到3个点取出来插入到这条路径或另一条路径的其他位置。这个操作对改善局部拥堵非常有效实现起来也不复杂。Swap交换两条路径中的各一个点或者一条路径内的两个点。适用于车辆间的负载均衡调整。局部搜索的停止条件是连续k次迭代没有改进k我一般设为100。每次局部搜索后重新计算路径的容量和总成本如果有改进就更新。def two_opt(route, dist_mat): improved True best_route route[:] best_cost compute_cost(route, dist_mat) while improved: improved False for i in range(1, len(route) - 2): for j in range(i 2, len(route)): new_route route[:i] route[i:j][::-1] route[j:] new_cost compute_cost(new_route, dist_mat) if new_cost best_cost - 1e-9: best_route new_route best_cost new_cost improved True route best_route return best_route, best_cost在实际代码里还有个加速优化不需要每次重新计算整条路径的成本只需要计算断点附近的那几段边的变化量即可。但对竞赛规模的数据来说全量重算是可以接受的代码简洁优先。6. 结果怎么呈现才能对得起你的努力6.1 收敛曲线与分析算法跑完第一件事是画收敛曲线和路径图。收敛曲线要能证明你的算法是正常迭代的而不是一开始就乱跳到末尾。具体来说画出每一代的最优适应度和平均适应度两条曲线最优适应度应该是单调下降或先快速下降后趋于平稳平均适应度围绕最优值波动但整体也呈下降趋势。如果最优曲线在早期就完全不动了说明算法大概率陷入了局部最优就要考虑调整变异概率或者增加局部搜索力度。路径图方面我会画一个带坐标的地图散点图用不同颜色区分不同车辆的路线并标注depot和各收运点。这张图一定要画得清晰因为评阅人大概率第一眼就看结果图路径规划效果好不好一目了然。如果路径交叉严重即便总距离指标不错印象分也会打折扣。我的经验是画图时把od的标记放大、路线线条加粗、图例注明车辆编号和车辆载重利用率让每辆车跑的全貌一看就懂。6.2 对比实验你至少要做三组竞赛评阅里对比实验是拉开档次的关键。光说自己算法好是不够的得在同一套数据上对比才行。我做了三组对比第一组无算法调度当前人工方案vs 改进后优化方案。如果题目数据里给了类似现有方案的参考值这组对比最直观直接说明你的方法能省多少钱、少用多少车、少跑多少公里。如果题目没有给可以用一个简单的按距离最近邻贪心构造的基线方案做对比体现算法优化的必要性和效果。第二组单一算法 vs 你的改进算法。比如纯遗传算法和GA局部搜索在相同迭代次数下的对比重点展示混合算法在收敛速度和解的质量上的优势。第三组小规模精确解 vs 启发式解。如果收运点数较少用Gurobi求出精确最优解再和你的启发式结果对比说明启发式算法取得的解和最优解的差距在可接受范围内比如5%以内。这组实验最有说服力直接给评审一个量化感知你的算法离最优解还有多远。6.3 灵敏度分析怎么设计灵敏度分析是很多同学容易忽略但特别加分的部分。我做了三组灵敏度分析车辆载重变动-20%、-10%、10%、20%对总车辆数和总距离的影响时间窗松紧程度扩大或缩小20%对可行性的影响垃圾量预测误差±10%对方案鲁棒性的影响。这些分析一方面展示了对模型的深入理解另一方面也回应了现实场景中参数波动的风险。如果题目给了多个场景的数据比如不同片区或不同季节的数据还可以做多场景的对比分析说明算法在不同规模、不同分布下的泛化能力。这是拿创新点分的好地方。6.4 论文每部分的写作要点摘要部分要四段式第一段用两三句话说清楚背景和问题是什么第二段写你建立了什么模型、用了什么方法第三段给出核心结果用数据说话第四句轻轻一句点出方案的实践价值和可推广性。摘要里的数一定不能拍脑袋要和正文结果表完全对应评阅人会专门抽查。问题分析和模型建立部分重点写清楚为什么这样建模和哪些假设是不合理的。不要只罗列数学公式要交代每一条约束的现实意义——比如车辆装载量不能超过额定载重对应的是车辆行驶安全。符号说明要全公式编号要有序这是细节分。算法设计部分是评阅人重点看的段落。建议给每一个算法模块配上流程图或伪代码伪代码的每一行都要和后面代码核心逻辑对应这样审稿人读起来会非常顺。参数表要单独列出来包括参数名称、含义、取值、调节依据。结果部分要注意图文并茂。一张总结果对比表、一张路径可视化图、一张收敛曲线图三个基础图必须有。好的论文在结果部分会很克制每个数据都有明确的作用不堆砌。7. 给下一届选手的几条实在建议代码规范这块从第一天就按工程标准写。变量名别用a、b、c这种文件夹按data、src、results、figures分好结果输出统一用CSV保存。比赛的第三天晚上你可能要改十几次参数重跑几十组实验如果代码一团乱光debug就能耗掉你半天时间。时间分配上我的建议是第一天上午快速读题和数据下午完成数据预处理和基线模型第二天全力做算法和实验晚上开始写论文框架第三天上午补全实验和分析下午集中写完整论文、统一格式、校对摘要。千万不要花两天时间搭模型最后只剩半天写论文——那样你再好的算法也白搭。最后说一个我个人特别受益的习惯做任何模型改动前先把当前版本的结果完整备份好并记录下来改了什么、效果如何。这样在论文讨论部分写我们尝试了xxx方法结果yyy的时候每一步都有真实依据而不是靠回忆拼凑。建模竞赛比的从来不只是算法它比的是谁在有限时间里更系统地解决问题。希望这篇复盘能帮你少走一些弯路祝你比赛顺利。本文还有配套的精品资源点击获取