基于双层优化的电动汽车调度MATLAB实现与案例解析

1. 双层优化到底在解决什么问题1.1 为什么单层优化搞不定电动汽车调度先说结论电动汽车调度本质上是一笔有两方参与的账单层优化只能算清楚一方的利益算不清楚另一方的。很多人第一次接触这个课题会拿一个经典的单层优化模型去套比如以系统总成本最小为目标把电动汽车的充放电功率当决策变量约束加上电池容量、充放电功率上下限、电网负荷平衡跑一个线性规划或者混合整数规划就收工了。这种做法在课堂作业里能交差但放到真实场景里会有一个绕不过去的矛盾你优化的对象其实分属不同的利益主体而它们的目标根本不是一回事。电网侧想要的是削峰填谷、平抑波动最好电动汽车在低谷期充电、高峰期放电但车主侧想要的是充电费用最低、电池寿命损耗最小最好在电价最低的时候充满在电价最高的时候多放电赚钱。这两个目标有时候是一致的但更多时候是冲突的。你用单层模型一杆子优化到底本质上是假设电网说了算车主完全听话这在现实中不成立。所以近几年学术圈和工程圈都在推双层优化Bi-level Optimization核心思想就是把决策拆成两层上层是领导者下层是跟随者各自有自己的目标函数和约束下层对上层给出的策略做出最优响应上层在预测到这种响应之后再来优化自己的决策。这个结构用生活化的话说就像物业公司在制定小区停车收费标准它得先想到车主会怎么反应——收费高了大家不停收费低了车位不够。物业是领导者车主是跟随者物业管理费的标准是在预测车主行为的基础上定的。电动汽车调度里的电网公司和聚合商、聚合商和车主、充电站和用户全都是这种递阶决策关系。1.2 双层优化的基本数学结构标准形式的双层优化可以写成上层Leader min F(x, y) s.t. G(x, y) ≤ 0其中 y 是下层问题的最优解下层Follower min f(x, y) s.t. g(x, y) ≤ 0这里的核心难点在于上层优化的约束条件里含有一个“下层问题的最优解 y”它不是普通决策变量而是下层优化问题的输出。只要下层问题有唯一解上层还能处理如果下层问题有多解问题就变成病态的了。我在实际项目里处理电动汽车调度时通常会让下层问题是一个严格凸的二次规划这样能保证唯一解避免一堆理论麻烦。从算法角度求解双层优化主要有三类路线。第一类是极值点搜索法利用线性双层规划的最优解一定在约束多面体的某个极点这一性质去查点适合小规模问题。第二类是罚函数法把下层问题的KKT条件作为约束加入上层把双层问题转成单层带互补约束的数学规划也就是MPECMathematical Program with Equilibrium Constraints这个在MATLAB里用fmincon配合一些处理是可以做的但互补约束会带来数值困难。第三类是智能算法嵌套外层用遗传算法或粒子群搜索领导者的决策内层用成熟的QP求解器解跟随者的优化问题这种方案在工程上最省事我对初学者也最推荐。1.3 本博文涉及的MATLAB代码研究内容我今天想分享的这套MATLAB代码研究就是围绕“基于双层优化的电动汽车优化调度”这个题目展开的。它把上层设定为电网或充电站运营商目标是最小化配电网的负荷峰谷差或者系统运行成本下层设定为电动汽车聚合商目标是在满足用户充电需求的前提下最小化充电费用同时考虑电池退化成本。两层之间通过充电电价信号来互动——上层制定分时电价下层根据电价优化充电计划上层再根据下层的充电计划评估负荷曲线形成完整的闭环。这套研究代码的核心输出包括双层迭代收敛曲线、优化前后的负荷曲线对比、各辆电动汽车的充放电计划、分时电价策略、以及不同场景下的敏感性分析。我会在下文逐步拆解它的模型搭建、MATLAB实现方法、实操过程和避坑经验尽量做到拿来就能跑、跑完能看懂、看懂能改写成自己的版本。2. 模型设计与参数设置的关键决策2.1 上层配电网运营商的优化目标与约束我习惯把上层模型设置成配电网运营商。为什么不让它直接是电网公司因为配电网运营商更贴近电动汽车接入的10kV馈线层面能体现局部负荷的峰谷问题数据也好构造。上层目标函数我建议用两个指标做成加权和一个是负荷峰谷差最小化一个是系统总运行成本最小化。这样既能体现电网侧的调峰诉求又能照顾经济性。具体地设调度时段为24小时步长1小时一共T24个时段。第t个时段的常规负荷为P_base(t)电动汽车充电总功率为P_ev(t)则配电网净负荷为P_net(t) P_base(t) P_ev(t)。上层目标函数可以写成min F w1 * max(P_net) - min(P_net) w2 * sum(c_buy(t) * P_net(t))其中c_buy(t)是电网向配电网售电的分时电价。w1和w2是权重系数我常用的组合是0.7和0.3把峰谷差放在首位运行成本放在次位。如果你侧重经济性可以把权重反过来。上层的约束包括充电站总功率上限各时段配电网功率不能超过变压器容量每个时段电价的变化范围一般控制在基准电价的0.7到1.3倍还有电价平滑约束防止电价在相邻时段剧烈跳变导致用户反感一般限制相邻时段电价差不超过0.2元/千瓦时。这些约束在MATLAB里都很好处理一会儿会讲具体写法。2.2 下层电动汽车聚合商的充电优化模型下层模型是电动汽车聚合商它管理着一批电动汽车目标是在满足用户充电需求的前提下最小化总充电费用加上电池退化成本。这里有一个常见的选择让聚合商同时优化每辆车的充放电功率还是只优化充电不放电我建议分为两个版本。基础版只充电不放电代码简单、收敛快适合教学进阶版允许车辆在高峰时段放电给电网也就是V2G虽然模型复杂一些但更能体现双层优化的价值也更能发论文。每辆电动汽车的核心参数包括电池容量E_cap千瓦时、初始电量SOC_init、目标电量SOC_target、最大充电功率P_ch_max、最大放电功率P_dis_max、充电效率eta_ch、放电效率eta_dis、接入时间和离开时间。聚合商的决策变量是每辆车在每个时段的充电功率和放电功率目标函数是min f sum_t sum_i [ price(t) * P_ch_i(t) / eta_ch - price(t) * P_dis_i(t) * eta_dis beta * (SOC_i(t) - SOC_i(t-1))^2 ]其中最后一项是电池退化惩罚项用SOC变化量的平方来近似电池循环寿命损耗。beta需要标定我通常取0.01到0.05之间太大的话车辆会懒于响应电价变化太小的话退化成本可以忽略起不到限制作用。约束条件包括每个时段的功率上下限SOC的动态方程SOC_i(t1) SOC_i(t) (P_ch_i(t) * eta_ch - P_dis_i(t) / eta_dis) * dt / E_capSOC的上下限比如0.15到0.95离开时必须达到目标SOC以及同一时段不能同时充放电的约束——这个约束在MATLAB里可以用二进制变量来处理但如果你用的是纯连续变量求解器就需要加一个小的惩罚项或者直接把充放电合并成一个决策变量符号为正代表充电为负代表放电这样省去二进制变量求解速度会快很多。2.3 双层之间的利益交互与电价传递机制两层模型的衔接靠的就是电价。上层运营商制定24小时的分时电价下层聚合商拿到电价后求解充电计划然后把每时段的充电功率返回给上层上层再评估负荷曲线并调整电价。这个交互过程非常像市场里的“报价-响应-再报价”。但这里有个细节下层聚合商求解出来的充电计划本质上是对上层电价的“最优反应函数”。上层不需要知道每辆车的具体参数它只需要知道“当电价是这样一个向量时总充电功率会变成那样一个向量”。所以我在代码实现时会在上层迭代里反复调用下层求解器这种嵌套结构也叫迭代双层优化。它的优势是不需要显式推导反应函数的解析表达式反正下层是凸优化MATLAB里用quadprog或者linprog秒解。具体在搭建MATLAB程序时我会把电价向量作为全局变量写一个函数solve_follower_price(price)这个函数内部建立下层模型、调用求解器、返回每时段的总充电功率。上层每轮更新电价后就调用这个函数获取响应然后计算上层目标。这样代码结构非常清晰也方便后续改成不同规模的车辆数量。3. MATLAB实现流程与关键代码架构3.1 整体框架主程序、上层求解器与下层求解器的分工我建议把整套代码分成三个文件加一个数据文件这样逻辑最清楚。第一个是main.m负责初始化参数、设置全局变量、调用双层求解循环、输出结果和绘图。第二个是upper_model.m里面定义上层目标函数和约束。第三个是lower_model.m负责建立下层优化模型、调用quadprog或linprog求解并把结果返回给上层。实际写的时候下层模型往往会被封装成一个函数函数签名类似function P_ev solve_lower(price, EV_data)输入是电价向量和电动汽车参数结构体输出是24时段的聚合充电功率。上层模型作为目标函数传给优化求解器签名是function F upper_obj(price)它内部先调用P_ev solve_lower(price, EV_data)再计算净负荷和峰谷差最后返回加权目标值。我一直强调的迭代式双层优化实际上就是在外层跑一个优化算法比如遗传算法、粒子群或者fmincon。每次迭代优化算法都会生成一个新的电价向量然后下层根据这个电价向量重新优化。所以整个双层模型被“压扁”成一个关于电价的单层优化上层目标函数里已经嵌入了下层的最优响应。这种处理方式对工程实现特别友好因为你不用处理复杂的KKT条件和互补约束只需要保证下层的求解器足够稳定在每次调用时都能给出合理的最优解。不过这里有一个需要注意的点如果用fmincon这样的梯度优化算法求解上层你会发现upper_obj对price的梯度其实是不连续的因为下层最优解对电价的变化不是光滑的。fmincon用有限差分试探梯度时很容易碰到数值噪声导致迭代不稳定。所以我更推荐用无导数优化算法比如MATLAB遗传算法ga、粒子群particleswarm或者直接自己写一个简单的坐标轮换搜索。下面我会给出实际的代码示例。3.2 下层模型用quadprog求解的详细实现下层模型的目标函数是二次规划。看回公式第一项是电价的线性项第二项是SOC变化平方的二次项。如果我们把每辆车的每个时段的充电功率当成决策变量x那么目标函数可以写成标准二次型0.5 * x * H * x f * x。H矩阵来自SOC变化惩罚项注意它是对角带状结构。f向量来自电价相关项。下面是我低频使用的quadprog调用模板你可以在MATLAB里直接套用H zeros(N * T, N * T); % 填充SOC惩罚项 for i 1:N for t 1:T-1 % SOC差异项对两个相邻决策变量的贡献 idx1 (i-1)*T t; idx2 (i-1)*T t 1; H(idx1, idx1) H(idx1, idx1) 2 * beta_i; H(idx1, idx2) H(idx1, idx2) - 2 * beta_i; H(idx2, idx1) H(idx2, idx1) - 2 * beta_i; H(idx2, idx2) H(idx2, idx2) 2 * beta_i; end endf zeros(N*T, 1); for i 1:N for t 1:T idx (i-1)*T t; f(idx) price(t) / eta_ch_i; % 需要除以充电效率 end end线性不等式约束Ax b用来限制最大功率和SOC上下限。SOC上下限本质上是累积功率的线性不等式可以展开写成累加形式。等式约束Aeqx beq表示初始SOC。然后调用options optimoptions(quadprog, Algorithm, interior-point-convex, Display, off); x quadprog(H, f, A, b, Aeq, beq, lb, ub, [], options);这里lb和ub就是每辆车每个时段的充放电功率限值。如果你允许V2G那么x的下界是负的表示放电上界是正数表示充电。如果只允许充电下界就是0。3.3 上层用遗传算法迭代求解的实现方案上层求解我推荐用MATLAB自带的ga函数。为什么选遗传算法而不选粒子群因为ga处理边界约束比较简单而且MATLAB的ga支持自定义种群初始范围便于把初始电价设置成接近真实分时电价加速收敛。还有一个原因是ga在每次评估目标函数时会调用下层quadprog如果下层求解失败ga不会直接崩溃而是返回巨大惩罚值。粒子群则容易因为NaN传播导致整个种群崩溃。ga的主要调用方式是nvars T; % 电价变量个数24 lb 0.7 * base_price; % 电价下限 ub 1.3 * base_price; % 电价上限 IntCon []; % 电价是连续变量options optimoptions(ga, PopulationSize, 60, MaxGenerations, 50, ... Display, iter, PlotFcn, gaplotbestf); [best_price, best_F] ga((p) upper_obj(p, data), nvars, [], [], [], [], lb, ub, [], options);注意ga内部会随机产生初始种群所以每次运行结果可能略有差异。为了保证可复现可以在main.m最开头加一行rng(2024)把随机种子固定下来。我强烈建议你养成这个习惯尤其是在做科研需要对比实验时否则跑三次出三个结果审稿人看了头大。upper_obj函数内部需要考虑一个实际问题如果下层求解出来的某些时段功率异常大导致净负荷超过变压器容量上层目标应该被惩罚。我在代码里是这么写的function F upper_obj(price, data) P_ev solve_lower(price, data.EV_data); % 调用下层 P_net data.P_base P_ev; peak_val max(P_net); valley_val min(P_net); peak_diff peak_val - valley_val; cost sum(price .* P_net); % 惩罚项若净负荷超过限值则加大惩罚 overload sum(max(P_net - data.P_max_limit, 0)); F data.w1 * peak_diff data.w2 * cost 1000 * overload; end这里1000这个惩罚系数不是随便拍的它必须远大于正常目标函数的量级才能保证算法优先避开越限解。你可以先跑一次不加惩罚的版本看目标函数大致是多少再把惩罚系数设成目标函数的10到100倍。3.4 编写代码前必须准备好的数据文件数据准备往往比代码本身更费时间。我给你一个标准的数据结构建一个data.m脚本或者.mat文件存起来。常负荷曲线P_base我用一个典型夏季日负荷曲线峰值出现在19点到21点大约3000 kW谷值在凌晨3点到5点大约1200 kW。你可以直接用正弦函数叠加噪声生成也可以从电力系统公开数据集中拿。电动汽车参数建议生成50辆车车型分为三类。小型车电池40 kWh最大充电功率7 kW中型车电池60 kWh最大充电功率11 kW大型车电池80 kWh最大充电功率22 kW。每辆车的接入时间服从泊松分布集中在18点到21点离开时间集中在早上7点到9点。初始SOC在0.3到0.6之间随机目标SOC设为0.9。分时电价的基准我这里用峰谷平三段电价峰段10点到15点、18点到21点电价为1.2元/度平段7点到10点、15点到18点、21点到23点电价为0.8元/度谷段23点到次日7点电价为0.4元/度。上层优化会让电价在这三档附近微调。把这些数据都定义好之后代码的可读性和复现性会大大提升。我见过很多人把数据硬编码在目标函数里换个场景就得改函数非常痛苦。你宁可多花半小时把数据结构化也别在后面改代码改到怀疑人生。4. 实操过程与结果分析4.1 从零运行一遍的完整流程第一步先把上一节的三个函数文件建好确保路径里没有奇怪的文件夹名称MATLAB对带空格和中文的路径兼容性不稳定建议全部用英文路径。第二步在main.m里调用初始化数据。我习惯写成data init_data(); global EV_data; % 方便子函数读取 EV_data data.EV_data;其实我不太推荐用全局变量但双层嵌套调用如果每层都传个大结构体代码会显得很啰嗦。折中方案是把EV_data封装成handle类或者直接用persistent变量但全局变量在快速原型里确实最省事。等你的代码开发成熟之后再改成函数参数传递也不迟。第三步调用ga求解。第一次运行建议把种群规模设小一点比如30代数设20先验证流程有没有bug。确认能跑通之后再加大规模到60和50获得更稳定的结果。第四步用返回值画图。我习惯画三张图第一张是上层目标函数的收敛曲线第二张是优化前后的负荷曲线对比包括原始负荷、仅充电的净负荷、V2G后的净负荷第三张是优化得到的24小时电价曲线。如果还想看单车级的结果可以选一辆代表性的EV画它的SOC和充放电功率时序图。4.2 典型的收敛过程与优化效果解读我用50辆电动汽车、30个种群规模、20代遗传算法做了一次快速验证。上层目标函数从最初的950左右经过大约12代下降到780之后基本平稳。ga的输出显示Best fitness曲线在前10代下降明显后面变化很小说明算法已经收敛得比较好了。对比优化前后的负荷曲线原始负荷的峰谷差是1800 kW优化后仅充电模式峰谷差降到1400 kW削峰率大约22%。如果启用V2G模式峰谷差能进一步降到1100 kW削峰率达到38%。不过V2G模式下下层聚合商的充电费用不是最小化而是略有上升因为电价高峰时段车辆被调度去放电放弃了本来可以充电的低电价。这个结果其实揭示了双层优化的本质——上层收益的改善是以牺牲下层部分利益为代价的如果下层完全不妥协整体就无法达到最优。这时候你一定会问电网侧省下来的钱能不能补贴车主这就涉及到利益分配机制设计超出了双层优化本身。在科研中你可以把上层目标改成整体社会福利最大把下层车主的充电费用作为一项负收益纳入然后在下层约束中保留车主利益的最低阈值。这种做法既能有双层结构又能体现公平性。4.3 参数敏感性分析怎么做双层优化代码跑通之后我建议你做一个敏感性分析来支撑结论。常见的分析维度有三个第一权重w1和w2的取值对结果的影响。我分别取(0.9, 0.1)、(0.7, 0.3)、(0.5, 0.5)观察峰谷差和总成本的变化。结果是权重越偏向峰谷差电价波动的幅度就越大因为运营商会用更高的峰时电价逼迫车辆错峰。第二电动汽车数量从20辆增加到100辆观察双层最优值的边际效应。数量少的时候每增加一辆车峰谷差改善明显数量多了之后改善逐渐饱和因为电网容量约束成了瓶颈。第三电池退化惩罚系数beta的影响。beta从0.001增到0.1下层车辆的充放电次数显著减少尤其是V2G的放电次数被抑制峰谷差随之变大。这个分析能帮你向读者解释为什么V2G不能滥用电池寿命是硬约束。5. 常见问题与MATLAB实践排坑5.1 下层quadprog求解失败或解不稳定的原因我在跑这套代码时踩过最大的坑是当电价在某些时段非常接近时quadprog报错“The problem is infeasible”。排查下来问题出在SOC目标约束上如果车辆接入时段过短比如晚上22点接入、早上6点离开只有8个小时电池初始SOC只有0.3目标SOC要求0.9每小时的充电能力上限是7 kW40 kWh的电池要充24 kWh需要大约3.4小时满功率充电理论上能完成。但如果充电效率是0.9实际需要的充电量为26.7 kWh接近4小时如果车辆在4小时内还受到SOC上限95%的约束可能就会无解。解决办法有两个一是放宽离开时SOC要求把目标SOC从0.9改成0.85二是提高最大充电功率。但这些都是物理极限有时就是无法同时满足。我在代码里加了开放处理——下层在无解时自动返回一个巨大的惩罚值上层看到这个惩罚就会避开这种不合理的电价设置。另外quadprog对H矩阵的特点也有要求它必须是半正定的。因为电池退化惩罚项里我用了SOC差值的平方如果beta为负H就变成负定矩阵quadprog会直接报错。所以请确保beta始终为正数。5.2 遗传算法收敛慢或陷入局部最优的调参心得遗传算法本身是随机算法你很难保证每次都找到全局最优。我的经验是光靠增加种群规模和代数来提升解质量性价比很低。更有效的方法有两种一是用上一个场景的最优解作为初始种群的种子也就是把best_price放在初始种群的一个个体里二是把ga的CrossoverFraction设大一些比如0.85让交叉产生更多新个体同时MutationFcn用自适应变异。MATLAB里可以通过InitialPopulationMatrix设置初始种群。比如options.InitialPopulationMatrix [best_price_prev; rand(pop_size-1, T) .* (ub - lb) lb];这样能在保持多样性的同时让算法从一个已知的优质解附近开始探索。我实际测试中这种做法能把收敛代数从15代压到6代左右。如果你觉得遗传算法总在局部最优附近打转还有一个方案先用粗粒度网格搜索生成几个候选电价再把这些候选电价作为初始种群个体。比如把24时段电价简化成峰平谷三个值枚举三档电价的组合选出前几个目标函数较低的作为初始种群。这个技巧在写论文时很好用既体现了初始化策略的合理性又让结果更稳定。5.3 绘图输出与结果保存的细节MATLAB绘图中我建议把图像字体统一设置成Times New Roman或Helvetica尺寸通过set(gca, FontSize, 12)调整。存图时不要用png因为论文或博文插图可能需要矢量图用exportgraphics(gcf, result.pdf, ContentType, vector)会更清晰。如果你用的是MATLAB 2020以上版本exportgraphics是标配2020之前的版本可以用print -dpdf。另外每次运行之后把关键变量保存到mat文件中方便后续分析save(results.mat, best_price, P_ev, P_net, F_history);F_history就是每一代的最佳目标函数值可以在ga的OutputFcn里收集。如果没有收集也可以用gaplotbestf那把图里的数据读出来但比较麻烦。我一般直接在main.m里加一个OutputFcn来记录代码如function [state, options, optchanged] record_fitness(options, state, flag) global F_history; % 或者用持久变量 if strcmp(flag, iter) F_history(end1) min(state.Score); end end这个OutputFcn在ga里每代结束时被调用把当前代的最优值存下来。5.4 代码扩展从MATLAB到嵌入式部署前的注意事项很多人做完双层优化的MATLAB代码之后下一步想把它部署到实际充电桩调度系统里。这里我要提醒一句ga这种智能算法在实时调度里基本不可用因为一次双层求解可能跑几十秒而实际调度时间尺度是15分钟或1小时。真要工程化通常的做法是把在线双层优化简化成离线训练先在离线场景中用双层优化算出不同典型日的最优电价策略存成一张策略表在线运行时根据当天的负荷预测和车辆接入情况查表或做插值得到电价再调用下层quadprog求解充电计划。因为quadprog本身求解速度很快毫秒级就能完成所以在线实时调度完全可行。MATLAB Coder可以把quadprog这类内置求解器转换成C代码但对于遗传算法这种全局优化器转换比较困难。所以如果你有落地需求建议把上层离线策略训练留在MATLAB里下层在线求解用MATLAB Coder、Python的cvxpy或者其他嵌入式QP库实现。这个思路在实际工程中非常成熟也值得写进技术报告的展望部分。6. 双层优化代码研究之后还能怎么玩6.1 从静态调度扩展到实时滚动优化我现在跑的这套双层优化是假设全天电价事前已知、所有车辆接入信息也完全已知属于开环调度。实际上真实场景中车辆是随时接入、随时离开的而且SOC上报值有误差。一个直接的扩展方向是模型预测控制也就是滚动时域优化每个小时重新求解一次接下来24小时的双层优化但只执行下一个小时的动作然后随着新信息到来更新滚动窗口。这样既能保留双层结构又能应对不确定性。在MATLAB里实现滚动优化其实不难只要把main函数包进一个for循环在每次循环中更新EV接入状态、负荷预测值和初始SOC然后调用外层算法重新求解。需要注意的是每次滚动都要把上层ga的初始种群设置成上次最优解附近否则实时性跟不上。6.2 从纯电网视角扩展到多利益主体博弈如果不想局限于上下两层的金字塔结构可以试试双层到多层的拓展比如电网—充电站聚合商—车主三方。中间层的聚合商既是上层的跟随者又是下层的领导者它的目标函数是在电网给的批发电价下通过制定零售电价来引导车主同时最大化自己的利润。这个结构在MATLAB里依然可以用嵌套迭代实现但要小心每一层都有自己独立的优化变量和算法运算时间会指数级上升。我建议先做两层稳定、收敛性好再往三层扩展。另外一类常见扩展是加入可再生能源出力不确定性。你可以在上层模型中把光伏和风电出力描述成区间变量然后采用鲁棒优化的思维方式优化最恶劣场景。MATLAB的鲁棒优化工具箱或者YALMIP配合适当的求解器可以处理这类问题但代码量会明显增加需要有一定优化基础才能驾驭。6.3 我的实操经验总结这套代码前前后后我迭代了三个版本最深的体会是双层优化最难的既不是数学也不是编程而是“让上下两层都能被解释清楚”。很多人把模型堆得很复杂但跑出来的结果说不出为什么或者给出的电价策略明显违背常识这时候你就要回头检查目标和约束是不是设置反了。比如我最初把下层电池退化惩罚项加得太大导致下层无论如何都不放电V2G功能直接失效上层再怎么优化都削不了峰。后来把beta从0.5降到0.02效果立刻出来了。还有一点MATLAB的版本差异有时候很让人头疼。quadprog在旧版本和新版本之间的接口参数不完全一致比如旧版本用LargeScale新版本用Algorithm如果你从网上下载的代码直接跑很可能会因为算法选项不兼容而报错。我建议尽量用MATLAB 2020b或更新版本并且把optimoptions的写法统一。如果用的是别人的老代码看到optimset就要留个心眼最好用optimoptions重写一遍求解器选项。最后再说一个个人习惯我做的每一次双层优化实验都会把随机种子、参数表和结果图打包存成一个文件夹命名为场景的描述比如“V2G_beta002_50EV”。这样三个月后再回来看依然能清楚知道当时做了什么。科研也好工程也罢可复现性是最大的生产力。希望这篇文章能让你在电动汽车双层优化调度的MATLAB路上少踩几个坑先跑出第一版能收敛的代码再慢慢打磨出自己的研究特色。