前段时间有个刚入门的同学问我教科书上把节点电价LMP公式抄了一遍又一遍但真让他自己写一个能算节点出清电价的小程序他连从哪下手都不知道。这个问题很典型——电力市场理论学习与实际出清程序之间缺的往往不是公式而是一份注释清楚、能直接改参数运行的代码。这篇就聊一聊我用一个三节点算例程序理解节点出清电价LMP的过程把模型、代码和踩过的坑都摆出来适合刚接触电力市场、想让公式落地的读者。程序备注写得很细几乎每行都能对应到教材里的某个概念你拿到手之后可以直接改负荷、改线路容量、改机组报价观察电价怎么变。1. 三节点算例程序解决什么问题从“背公式”到“看价格生成”1.1 节点电价不是“平均成本”很多初学者最容易误解的一件事是把实时电价理解成“系统里所有发电报价的平均值”。你问他什么是LMP他能背出“边际机组的报价加上阻塞分量”但你让他解释为什么同样一小时、同一个系统节点1的电价是25元节点3的电价要涨到40元他就开始含糊了。问题出在大多数人只看了结论没看过价格在优化模型里是怎么“长出来”的。节点出清电价不是财务人员事后摊出来的成本它是在安全约束经济调度结果里由约束条件的边际信息直接生成的数字。什么意思就是优化模型求完最优解之后每个约束背后都带了一个影子价格节点电价就是这些影子价格组合出来的结果。这份程序最好的地方就在于它把优化模型的每一个约束都写在明面上哪一行对应功率平衡哪一行对应机组上下限哪一行对应线路潮流约束清清楚楚。你跑完一遍再把求解器返还的对偶变量打印出来电价是怎么来的一目了然。1.2 为什么非要自己跑一遍程序看书和跑程序是两种完全不同的学习状态。看书的时候你是在接受别人的推导公式跳一步你就跟丢了跑程序的时候所有变量、约束、参数都是你自己定义的求解器任何一个警告你都得正面回答。我刚接触电力市场时也犯过懒总觉得自己知道原理就行。后面发现不行因为对话里全是“阻塞盈余”“影子价格”“对偶变量”这些东西你要是不亲手算出来别人问细节你马上就露馅。跑通这个三节点程序之后我最大的收获不是学会了一个工具而是把脑子里那些漂浮的概念全部锚定了下来。而且这个程序刻意控制得很小三个节点、两条发电机、一条负荷几分钟就能跑完。小算例有个优势每个数字你都能用手大概验算一遍程序结果对不对、物理现象合不合理你能迅速判断。换成几百个节点的大系统出了错你连错在哪都定位不了。1.3 这套程序适合谁如果你满足下面任何一个条件这套程序就值得你花半天时间跑一跑刚接触电力市场学完LMP定义但没写过出清模型在教材里看到过安全约束经济调度但不知道约束条件怎么影响价格想读《机组运行约束对节点电价的影响》这类文献但被里面的模型公式劝退准备做电力市场相关的课程设计、毕业设计需要一个基础框架不需要你有很强的编程基础。程序用的Python建模部分用CVXPY这类优化库你只要能看懂列表、字典、循环能区分“变量”和“参数”就够了。关键是每行代码都有注释注释写的是“这一行对应教材里的哪条约束”你对着书看代码就变成了公式的翻译器。2. 出清模型的最小骨架目标函数、三组约束与对偶价格2.1 为什么要用直流潮流近似电力系统交流潮流是非线性的含电压、无功、功角、网损一大堆变量。初学者如果直接面对交流最优潮流等于刚学会加减法就去做微积分。出清价学习的起步阶段完全没必要把交流模型拉进来。这里用的是直流潮流DC-OPF它做了三个关键近似假设所有节点电压幅值都等于1、线路电阻远小于电抗所以忽略网损、功角差足够小所以正弦值近似等于角度本身。做完这些简化之后线路有功潮流变成了节点相角的线性函数线路功率等于两端相角差除以电抗。这三条假设在物理上不够精确但在理解市场机制时完全够用。因为节点电价之所以有空间差异核心来自线路潮流约束——也就是阻塞。直流模型恰好把这种阻塞效应完整保留了下来而把无功、电压这些与电价形成机理关系不大的因素暂时放到一边。你先把直流模型的LMP搞清楚再往AC模型走就会发现网损分量只是多出来的一项无功价格只是在告诉你“电力市场不只有有功一个维度”。地基没打好就直接上交流模型那才是真正的灾难。2.2 目标函数和三类约束逐条拆解这个三节点程序里出清模型的经济含义是在保证安全的前提下找到总购电成本最低的发电组合。目标函数写出来就是min Σ(c_i * P_i)其中c_i是机组i的报价单位是元/MWhP_i是机组i的中标出力单位是MW。求和就是所有机组的发电成本之和。注意这里用的是“报价”而不是“成本”市场出清按报价最小化来排这是电力市场的规则设定。光有目标函数还不够还要有约束否则模型会给出不切实际的解。程序里有三组基本约束第一组是机组出力上下限也就是每台发电机的技术边界。一台机组不能低于最小技术出力运行也不能超过额定容量写成约束就是P_min P_i P_max第二组是系统功率平衡。任意时刻所有机组出力之和必须等于系统总负荷这里不考虑网损所以就是Σ(P_i) 总负荷第三组是线路潮流约束。每条输电线路都有热稳定极限直流模型下线路功率是相角的线性函数必须落在正反两个方向的最大容量内-F_max F_l F_max这三组约束就是出清模型的骨架。很多教材里的“安全约束经济调度”SCED本质上就是在这个骨架上增加各种细节爬坡约束、备用约束、网络安全约束、机组启停状态等。2.3 拉格朗日乘子如何变成节点电价节点电价不是从目标函数里直接读出来的它是优化问题里约束条件的“影子价格”。学过最优化的人都知道带约束的极值问题可以写成拉格朗日函数每个约束对应一个拉格朗日乘子。功率平衡约束的乘子物理含义是“在系统里多负荷1MWh总购电成本会增加多少”。这个数字是所有节点电价的公共基准通常叫能量价格。机组出力约束的乘子含义是“机组容量多1MW能带来的价值”这在真实市场里对应容量价值。线路潮流约束的乘子含义是“这条线路多1MW传输容量能节省多少系统成本”这就是阻塞价格。有了这些乘子节点电价的计算公式就顺理成章了LMP_i 能量价格 Σ(线路l对节点i的功率转移因子 * 线路l的阻塞价格)功率转移因子PTDF描述的是节点i增加1MW注入会在线路l上产生多大的功率变化。它是由网络拓扑和线路阻抗决定的可以预先算好。这样每个节点的电价都由两部分构成全系统共享的能量价格加上该节点独特的阻塞价格。2.4 为什么说“约束越紧价格信号越强”理解影子价格最好的方式就是想清楚如果某个约束从“紧”变成“松”价格会发生什么变化。功率平衡约束永远都是紧的因为负荷必须满足所以能量价格永远存在。机组出力上限约束在机组满发时才是紧的此时机组不再是边际机组它的容量价值通过约束乘子体现。线路容量约束在满载时才是紧的一旦某条线路满送它就会在所有相关节点之间制造电价差。这就是标题里那句“机组运行约束对节点电价的影响”的由来。任何一条机组约束只要它成为瓶颈就会以乘子的形式进入节点电价。程序里你只要把一台机组的出力上限设得很低它很快就会被拉满然后你会在运行结果里看到它对电价的影响。3. 程序怎么写成数据、建模、取对偶变量的关键代码3.1 先规划好数据再写约束写程序之前先想清楚需要哪些数据。一个最简出清算例至少要包含节点、机组、负荷、线路阻抗和线路容量五样东西。我的三节点算例数据是这么设定的节点1和节点2各挂一台机组节点3带100MW负荷。G1在节点1报价25元/MWh容量100MWG2在节点2报价40元/MWh容量100MW。线路有三条1-3线路容量40MW这是最关键的一条它会造成阻塞2-3线路容量80MW1-2线路容量60MW但这回把它的电抗设得很大让功率几乎不走这条线方便理解。数据定义用Python来写大概是这个风格# 机组参数所在节点、报价、最大出力、最小出力 gens [ {bus: 0, cost: 25.0, Pmax: 100.0, Pmin: 0.0}, # G1节点1 {bus: 1, cost: 40.0, Pmax: 100.0, Pmin: 0.0}, # G2节点2 ] # 节点负荷节点3负荷100MW索引号2 loads {2: 100.0} # 线路参数[首端节点、末端节点、电抗(标幺值)、容量] lines [ (0, 2, 0.1, 40.0), (1, 2, 0.1, 80.0), (0, 1, 10.0, 60.0), ]这种用字典和列表组织数据的方式对初学者最友好。后面约束怎么写、结果怎么解读全部可以回到这三行数据结构上找答案。3.2 建模核心变量、目标、约束的写法接下来用CVXPY把优化模型写出来。这里面最重要的是搞清楚哪些是决策变量哪些是已知参数决策变量是各机组出力P和节点相角theta其余都是已知参数。目标函数那一行就是在求和所有机组的费用。注意这里我用线性报价也就是每度电价格恒定没有分段报价。分段报价会让模型稍微复杂一点但原理是一样的等线性版本跑通之后再改不迟。约束部分建议先在草稿纸上写好数学形式再一行行翻译成代码。功率平衡是一个等式约束代码写一行机组上下限是边界约束可以直接写进变量的定义里线路潮流约束是这堆代码里最需要理解的要先根据线路阻抗形成导纳矩阵再求出功率转移矩阵。核心建模代码大致长这样import cvxpy as cp import numpy as np n_bus 3 n_gen len(gens) # 决策变量机组出力、节点相角 P cp.Variable(n_gen) theta cp.Variable(n_bus) # 目标函数总购电成本最小化 cost 0 for g, data in zip(gens): cost data[cost] * P[g] objective cp.Minimize(cost) constraints [] # 机组出力上下限 for g, data in zip(gens): constraints.append(P[g] data[Pmin]) constraints.append(P[g] data[Pmax]) # 节点功率平衡注入-流出负荷 # 这里需要按节点聚合机组出力再减去线路流出功率 # 简化写法先做全系统功率平衡 total_load sum(loads.values()) constraints.append(cp.sum(P) total_load) # 参考节点相角为0 constraints.append(theta[0] 0) # 线路潮流约束Bf theta 的计算略结果即各线路功率 # flow Bf theta # constraints.append(flow line_capacity) # constraints.append(-flow line_capacity) prob cp.Problem(objective, constraints) prob.solve()这里Bf矩阵是怎么来的需要先形成直流潮流的导纳矩阵B然后对B求逆得到相角与注入功率的关系。我在程序里已经把这些计算封装成函数并逐行注释了你不需要死抠矩阵推导但要知道它的作用把节点注入功率映射成线路潮流。有了这个映射阻塞约束才能写进模型。3.3 从求解结果里取对偶变量模型求完解之后出力P和相角theta可以直接打印。但要拿到节点电价必须从约束里取对偶变量也就是影子价格。CVXPY里每个约束都带有dual_value属性求解成功之后就能访问。关键要分清取哪一条功率平衡约束的dual_value就是能量价格线路潮流上限约束的dual_value就是阻塞价格的正向部分线路潮流下限约束的dual_value就是阻塞价格的反向部分程序里取对偶变量的代码大概是这样if prob.status cp.OPTIMAL: print(最优出清结果) for g, data in zip(gens): print(fG{g1} 出力: {P[g].value:.2f} MW) print(功率平衡约束影子价格能量价格:) print(flambda {constraints[bal_idx].dual_value:.4f} 元/MWh) # 线路阻塞价格 上限约束乘子 - 下限约束乘子 mu np.array([constraints[ub_idx].dual_value - constraints[lb_idx].dual_value for ...])因为线路约束写成了正反两个方向所以要相减得到一个净阻塞价格。然后节点电价用PTDF矩阵乘以阻塞价格向量加上能量价格就得到每个节点的LMP。完整的代码里这一段的注释写得非常详细因为这是整个程序最精华的部分。3.4 注释该怎么写才对初学者有用写这份程序时我给自己定了一个规矩每行代码都要能回答两个问题——这行在做什么以及为什么要这么做。比如说机组出力上限约束那一行注释不能只写“设置机组上限”而应该写“对应机组运行约束中的最大技术出力满发时其对偶变量会成为容量影子价格”。这样的注释才有信息量。再比如功率平衡约束注释要写明“等式约束乘子即为系统能量价格是LMP的共同基准”。这种注释等于把教材重点重新标注了一遍读代码就是在复习概念。反过来最没有用的注释是“这行定义变量”“这行设置参数”——废话。初学者看了跟没看一样。真正好的注释应该能让人在脱离教材的情况下只靠读代码就把出清机制理顺。这也是为什么我在程序开头专门写了一段README性质的大注释把三节点网络拓扑、报价表、预期现象全部讲清楚。3.5 跑通后的典型输出长什么样以我设的算例程序跑完后输出的核心结果大概是这个风格节点LMP元/MWh节点125.0节点240.0节点340.0机组中标出力MWG140.0G260.0线路潮流MW传输极限MW状态1-340.040.0满载2-360.080.0未满载1-2接近060.0未满载看到这个结果你应该马上产生两个问题为什么G1报价25元却只出了40MW为什么节点1电价25元而节点3电价40元这两个问题之间存在因果关系——正因为节点1到节点3的线路卡在40MW便宜机组不能随心所欲地把电送到负荷中心系统不得不从节点2调用报价40元的高价机组节点3的电价因此被抬到边际成本水平而节点1因为“便宜电送不出去”电价被压在25元。这个现象如果只是看教科书你可能知道“阻塞导致节点电价升高”但它很难在你脑子里形成一个生动的画面。程序输出出来的数字直接把这个画面钉死了。4. 跑通算例后做两个实验容量变化如何改写节点电价4.1 实验一线路容量够大时LMP全网相等跑通三节点程序之后我第一个建议你做的实验特别简单把1-3线路的容量从40MW改成200MW。重新运行一遍你会发现所有节点的LMP几乎全部变成25元/MWh。原因也很直观线路容量不再是瓶颈G1这台25元的便宜机组可以全力发电100MW全部送出还绰绰有余完全没有必要调用40元的G2。此时系统边际机组就是G1能量价格25元阻塞价格全部归零所有节点共享同一个电价。这个实验最大的价值在于帮助你形成“节点价差从哪来”的判断标准当网络没有约束时LMP处处相等当网络出现约束时价差出现。所以下次你看到哪个节点电价特别高第一反应应该是去查这个节点和便宜电源之间的输电通道是不是满了而不是疑惑“是不是哪台机组在乱报价”。4.2 实验二把线路容量收紧看价格分离第二个实验反过来把1-3线路容量从200MW再一步一步往回调150、100、80、50、40。每调一次记录一次LMP你会发现规律是这样的线路容量足够宽时节点电价全部25元随着容量收窄到一定程度G1不能全出力G2开始介入节点3的价格开始攀升等到线路容量低于40后价格分离格局彻底形成。这个渐进实验比直接给出40MW的最终结果更能说明问题。它演示的是“约束从松到紧、影子价格从零到非零”的连续过程。你亲自调参数跑一遍就会理解为什么说约束的松紧程度才是电价背后的真正推手。我在程序里把这组实验做成了一个循环函数一行行注释说明“每次只改线路容量其他参数不变”。如果想更直观还可以让程序在每次运行后打印一条消息标出当前线路满载状态和节点3电价。整个过程我建议你自己敲一遍不要直接复制粘贴敲代码的过程会强迫你理解变量之间的引用关系。4.3 实验三改变机组报价看谁决定节点电价第三个实验更有意思把G2的报价从40元改成60元其他不动重新求解。你会看到节点3的LMP不再是40元而是变成了60元。这个实验说明在1-3线路阻塞的格局下节点3的边际供应者已经变成G2G2的报价直接决定了节点3的电价。再改一次把G1报价从25元改成35元结果节点1的LMP又会跟着变。这个现象引出一个重要结论节点电价由“该节点边际供应机组的报价”加“阻塞分量”决定不是由所有机组报价的平均值决定。在某些节点电价可能长期由一台高价机组锁定即便全系统有很多便宜机组。这也就是为什么电力市场里要关注“市场力”——在阻塞情况下个别机组对局部电价几乎有定价权。4.4 用LMP分解解释“为什么这台机组满发还带低了整个电价”程序的输出里我专门放了一段LMP分解的汇总表把每个节点电价拆成能量价格和阻塞价格两部分。拿三节点算例来说能量价格由阻塞网络里的边际机组报价决定。当线路满载时G2是边际机组所以能量价格为40元/MWh。节点1的阻塞价格为负代表它受网络约束无法把电送到节点3因此节点电价低于能量价格节点3的阻塞价格接近零或略正所以电价等于或略高于能量价格。这个分解表建议你打印出来对照着看。很多教材里讲LMP分解会用数学公式推一遍但你真正理解它的方式是看到同一个数字在程序里被拆成两个低维分量。拆完之后你会发现所谓节点电价本质上就是一组约束条件的边际成本叠加。5. 程序调试中的五类典型坑对偶符号、参考节点、单位与求解器5.1 对偶变量符号反了这是所有初学者第一个踩进去的坑我自己也不例外。用CVXPY求对偶变量约束写法不同dual_value的正负号也不同你要是直接拿去当阻塞价格算LMP算出来的节点电价可能完全离谱。判断符号对不对有一个笨但可靠的方法手动构造一个只有一台机组、一条线路的极简系统用笔算出理论价格值再对比程序输出。只要符号反了程序价格一定会出现负值或者价差方向颠倒。还有个特征是当某条线路满载时阻塞乘子应该为正如果某条满载线路的乘子输出成负数大概率是符号问题。CVXPY文档里对偶变量的定义方向需要仔细看不同版本的求解器行为也略有差异。我在程序里刻意把所有含对偶变量的代码都写了类型注释并且用一个断言语句检查满载线路的乘子是否非负这样程序出错时会直接报警而不是让你默默拿到一个错误的电价。5.2 机组上下限约束的对偶变量怎么解释很多人只盯线路阻塞乘子忽略了机组出力上下限约束的乘子。实际上机组约束的乘子才是阅读《机组运行约束对节点电价的影响》这类文献时必须掌握的概念。举一个最直观的情况一台机组被出力上限卡住、满发运行。此时它的出力上限约束是紧的对应的乘子通常为正值物理含义是“如果这台机组能再多发1MW系统总成本能下降多少”。这个值越大说明这台机组在当前状态下越宝贵它的容量价值越高。出力下限约束的乘子则相反。机组被强压着不能低于最小技术出力运行时下限约束也会产生乘子但它往往意味着“这台机组本来不想发这么多但因为技术原因必须维持出力”其乘子有时为负值对应的是市场中的“最小出力约束补偿”问题。程序输出里我把机组约束乘子也列了出来并用注释标明“该值不等于机组报价它是约束的边际价值”。初学者往往把“报价”和“电价”混为一谈搞清机组约束乘子之后这个概念自然就分开了。5.3 线路功率方向带来的“哪条线满载”误解直流潮流的线路功率有方向程序里如果把两条线路潮流求出来一条为正一条为负初学者经常弄不清到底哪条线满了。以我三节点程序里的1-3线和2-3线为例它们采用正方向都是从低编号节点流向高编号节点。如果某条线路的实际功率是负的比如从节点3流向节点2那么“接近正值容量上限”和“接近负值容量上限”两个约束都得检查。很多时候一条线路看起来“没满”是因为它在反方向输送而程序判断的是绝对值是否超过容量。这个坑在查看程序输出时会很明显。我在代码注释里写了一个检查函数线路满载判定不是看潮流正数而是看潮流绝对值是否接近容量上限。千万别只看一个数字就下结论要先确认方向。5.4 单位与基准值标幺值换算错了价格全错电力系统计算常用标幺值节点电价却要落到元/MWh这种有名值。如果程序内部功率用标幺值报价却用元/MWh求出来的影子价格单位会乱套。我这个三节点程序里功率统一用MW有名值阻抗用标幺值但只影响转移因子不影响成本单位这样对初学者最不容易出错。你如果在其他材料里看到标幺值模型想拿到有名值价格一定要找到系统的基准功率然后把所有涉及功率的单位统一换算。这个坑极其隐蔽因为模型能正常求解出力的相对关系也对唯独电价数值不对。一个简单的自检方法把功率平衡约束的影子价格和边际机组报价对比在无阻塞情况下两者应该相等。如果不相等先查单位再查符号别急着怀疑模型结构。5.5 求解器报错时先看哪里程序跑不通时第一反应不要是逐行检查代码而是先看求解器给出的状态信息。CVXPY常见的状态有optimal、infeasible、unbounded三种分别对应的排查方向完全不同。optimal代表求出了最优解这时看看对偶变量是否合理。infeasible代表约束条件互相矛盾比如总负荷超过所有机组上限之和、线路容量低到系统无解优先检查数据是否有物理矛盾。unbounded代表模型缺少必要约束比如没有给机组设出力上限模型会无限增大便宜机组出力来凑目标函数这时功率平衡约束甚至可能无法成立。我在程序里加了一行注释专门提醒初学者先看prob.status再看约束的dual_value最后才看P的值。很多人一上来就打印P程序报错还在找P为什么不对方向就搞反了。6. 从节点电价程序走向机组运行约束与研究文献6.1 加入爬坡约束从单时段走向多时段三节点程序是单时段模型只能算某一时刻的出清结果。真实现货市场是连续运行的多时段模型机组上一小时出力这一小时能不能突然大幅变化不能于是有了爬坡约束。爬坡约束写成数学形式就是相邻时段出力之差不能超过规定的爬坡速率。加入这个约束之后模型的时间耦合就出现了之前那个时段的出力会限制下个时段的出力影子价格也会跨时段传递。对应到电价上你会看到某些时段即使负荷没变化电价也会因为机组爬坡受限而上涨。给三节点程序加爬坡约束是很自然的扩展路径只需要把变量从一维变成两维机组、时段把约束方程改成相邻时段差的形式。我建议你在这个版本的代码里先加一个24时段、每小时负荷变化的数据看看系统峰谷时段的LMP曲线长什么样这比直接啃多时段最优潮流的数学公式直观得多。6.2 加入机组启停约束最小运行时间与最小停机时间机组运行约束不只包括上下限和爬坡还有启停状态约束。机组不能刚开机就马上停机也不能刚停机就立刻开机这就是最小运行时间和最小停机时间。这些约束之所以重要是因为它们引入了状态变量模型从线性规划变成混合整数规划求解难度会上升一个台阶。但它们的影子价格同样会进入出清电价体系影响系统边际成本。标题里提到的《机组运行约束对节点电价的影响》讲的正是这类约束如何通过乘子改变节点价格。如果你在程序里加入启停变量建议从最简单的两台机组开始观察一台机组有最小运行时间约束另一台没有然后改变负荷序列看哪些时段电价会出现跳变。这种“加一个约束就观察一次电价变化”的思路比通读十篇文献都管用。6.3 参考文献应该怎么读提到这部分说说标题里那篇《机组运行约束对节点电价的影响》类的文献。很多初学者拿到这种论文习惯从头开始啃结果被第一节的模型符号劝退。我的建议是第一步先看摘要和结论搞清楚文章里“运行约束”指的是哪几类约束是出力上下限、爬坡、启停还是备用容量。第二步找到文章里有没有算例部分通常会有一个小规模系统比如三节点或五节点给出了机组参数和结果。第三步把你手头这份三节点程序改造成论文里的算例把参数填进去对比程序算出的节点电价和论文列出的结果。这样做的好处是你把一篇文献变成了自己程序的验证数据。你不需要完全复现它全部理论只需要抓住它的算例就能理解论文里那些公式在实际运行中对应的是什么。我当年就是靠这个方法把一堆关于机组组合和节点电价关系的文献逐个读下来的。6.4 从DC到AC网损、无功与更多真实细节直流模型把网损近似成零节点电价里没有网损分量。真实系统中电能从电源送到负荷中心必然产生损耗网损分量会让送端的节点电价更低、受端的节点电价更高。如果你继续往AC最优潮流方向走需要补充的东西是节点电压幅值变量、无功功率平衡方程、交流线路潮流方程、变压器分接头等。模型复杂度会大幅上升但LMP分解从三分量变成四分量多了网损分量价格形成机理终于完整。使用这个程序的一个好处是它为你搭好了骨架将来无论是加爬坡、加启停还是切换成AC模型底层的数据结构、结果读取方式都不用大改。很多真实系统里的电价分析工具走的也是类似的路子。6.5 一条更稳的学习路线总结一下这几步的顺序这是我试过之后觉得最平滑的路径第一步先跑通三节点DC-OPF能把LMP打印出来。第二步做容量和报价的敏感性实验理解LMP分解。第三步把单时段改成多时段加入爬坡约束。第四步用论文算例验证自己的程序。第五步再考虑交流模型或更复杂的市场机制。每一步都有一个可运行的代码版本兜底不会出现“概念没懂还要硬啃下一篇”的情况。你越往后面走越会发现自己对前面基础模型的理解在自动加深这就是程序驱动学习的正循环。最后再分享一点我个人的体会刚接触电力市场的人很容易高估理论的难度、低估实践的必要性。公式里的每个符号到底在算什么程序里一跑就露馅。所以别怕代码写得不够漂亮先让它跑起来再逐步把注释补充整齐。等你什么时候能拿着这套程序给别人从头到尾讲清楚每个约束、每个乘子为什么存在你对节点出清电价的理解就真正过关了。