如何评价2022年第十四届电工杯数学建模竞赛？
A、B题思路和部分代码已经更新！！！大家可以点赞收藏。主页自行领取gzh：UST数模社
<hr/>A题：高比例风电电力系统储能运行及配置分析 A题是一个风力系统储能问题，需要相关的专业知识，本题考量的是风力发电接入容量、供电成本，储能容量等；需要考虑建立多目标优化模型来求解，后面会更新相对应的详细思路
B 题 5G 网络环境下应急物资配送问题
B题是一个典型的路径优化问题，考察的是“配送车辆+无人机”的配送模式下的路径规划问题，我们可以用到模拟退火算法、遗传算法、粒子群算法来解决该问题。
我们首先可以将其配送路径设置为如下图所示：
VRP模型
下面放一种可供参考的VRP模型来解决问题：
首先我们对模型参数和决策变量进行设定
模型参数：
参数含义N所有点的集合M可以使用无人机配送的点的集合O起止点wi点i的配送重量W配送车辆的载重上限dij点i到点j的距离c无人机配送的距离成本系数B一个充分大的数决策变量：
变量类型含义xi,j0-1卡车是否经过从点 i 到点 j 的路径yi,j0-1无人机是否经过从点 i 到点 j 的路径si0-1点 i 是否由卡车配送zi,jR0+从点 i 到点 j 的路径上的卡车载重，可用于消除子回路we问题1 图 1给出 14 个地点，其中实线代表各地点之间的路线情况。若目前所有应急物资集中在第 9 个地点，配送车辆的最大载重量为 1000 千克，采取配送车辆（无人机不参与）的配送模式。请结合附件 1，建立完成一次整体配送的数学模型，并给出最优方案。
针对问题一，我们的约束条件为
(1) 流入/流出
(1a) 只能使用车辆配送的点，最大载重量为1000kg
物资都集合在了第 9 个点，车辆配送最大载重为 1000kg，那么我们可以设计程序步骤：
Step1：设置优化算法基本参数，设置最大路线数N，设置惩罚系数（该系数后面问通用）
Step2：生成初始个体，每个个体长度为n（地点数），通过randi()得到k条路线，通过randi(,1,n)随机赋予每个点所在的路线编号，并检验是否存在没有经过点的路线（检验时排除第9个点），否则重新randi(,1,n)，然后通过randperm(n)随机赋予各点一个顺序，这样就可以得到不同路线的不同路径
Step3：计算初始个体的目标函数，有了前面的距离矩阵，可以直接调用计算出总时间（肯定要算完整的路线从出发点到终点），同时计算每条路线的运载量（需求量之和），超载部分乘以惩罚系数作为额外的成本
Step4：迭代，这里遗传算法的交叉变异函数需要自己写，交叉过程就是rand<交叉率时，随便选择两个点的顺序交换或者路线编号交换，变异过程就建议重新用randperm(n)生成下路径顺序，或者用randi(,1,n)重新给予路线编号都可以，但交叉和变异过程都得进行检验
Step5：最后输出最优路线及路径
（最短路径还有模拟退火、禁忌搜索、蚁群等算法都可以实现）
问题2 图 2 中实线代表车辆和无人机都可以走的路线，虚线代表只有无人机可以走的路线。应急物资仍然集中在第 9 个地点，配送车辆的最大载重量为 1000 千克，采取“配送车辆+无人机”的配送模式。请结合附件 2，建立完成一次整体配送的数学模型，并给出最优方案。
第二问就在第一问程序基础上增加条件就可以了，逻辑给同学们理一理哈，无人机的速度是优先于车辆的，所以只要满足无人机配送条件的都可以直接通过无人机配送，程序添加第一问Step2最后，针对每条路线中，判断无人机可以连续配送多少个点，配送完了后就到达下一个点等待车辆或者车辆先到达等待无人机，车辆就直接前往无人机配送完后将抵达的下一个点，这里一定要考虑无人机的航程和续航时间，必须到达点上才行，都这只配送一个点或者不配送。
建议同学们增加一个0-1矩阵，维度就是所有点数，用来记录那些点是无人机配送的就行，第二问主要是加个0-1矩阵，然后重新换算下配送时间，第一问整体程序框架未做大改
(1b) 可以使用无人机配送的点
(2) 去除自循环
(3)是否由车辆配送
(4) 子回路去除
(5) 无人机配送的流入点和流出点，符合车辆路径的先后关系
(6) 每辆车只有一架无人机
目标函数
最小化，加权距离成本
VPR部分代码：
from datetime import datetime
import random
import math
from utils.logger import log_vrp as log
dts = datetime.now()
# nodes
NUM_NODE = 50
NUM_NODE_NAV = 10
LIST_NODE_UAV = list(range(NUM_NODE - NUM_NODE_NAV + 1, NUM_NODE + 1))
# coordinates
random.seed(2021)
RANGE_COORDINATE = (0, 100)
LIST_COORDINATE = [(random.randint(RANGE_COORDINATE, RANGE_COORDINATE),
random.randint(RANGE_COORDINATE, RANGE_COORDINATE)) for _ in range(NUM_NODE)]
# distance
ORIGIN = (round((RANGE_COORDINATE - RANGE_COORDINATE) / 2),
round((RANGE_COORDINATE - RANGE_COORDINATE) / 2))
LIST_COORDINATE_ = + LIST_COORDINATE
MAT_DISTANCE = [ for _ in range(0, NUM_NODE + 1)]
range_wave = (0.8, 1.2)
for i in range(0, NUM_NODE + 1):
for j in range(0, NUM_NODE + 1):
distance = round(math.sqrt((LIST_COORDINATE_ - LIST_COORDINATE_) ** 2
+ (LIST_COORDINATE_ - LIST_COORDINATE_) ** 2)
* (range_wave + random.random() * (range_wave - range_wave)), 4)
MAT_DISTANCE = distance
# cost
COST_UAV = 0.1
# weight, load
upper_weight_item = 20
LIST_WEIGHT =
UPPER_LOAD = 100
log.info(msg=&#34;total weight:{}&#34;.format(sum(LIST_WEIGHT)))
dte = datetime.now()
tm = round((dte - dts).seconds + (dte - dts).microseconds / (10 ** 6), 3)
log.info(msg=&#34;random data generating time:{} s&#34;.format(tm))
更多思路：
会提供参考资料和思路模型
页:
[1]