我前前后后幫好幾個(gè)課題小組審過(guò)這類題目凡是帶“主從博弈”和“粒子群”兩個(gè)詞的能源系統(tǒng)優(yōu)化十有八九都是圍繞同一個(gè)核心問(wèn)題多個(gè)決策主體各自打自己的小算盤但彼此之間又有上下級(jí)的制約關(guān)系怎么在算法上把這層關(guān)系解出來(lái)。這篇文章我就直接拿“基于粒子群優(yōu)化算法的三方三層主從博弈能源系統(tǒng)優(yōu)化模型”當(dāng)例子把從模型搭建到Matlab代碼實(shí)現(xiàn)的過(guò)程完整拆開講一遍包括我踩過(guò)的坑和調(diào)參經(jīng)驗(yàn)?？赐曛竽阒辽倌芨闱宄碌谝蝗饺龑又鲝牟┺牡臄?shù)學(xué)結(jié)構(gòu)到底怎么落到代碼里第二粒子群算法在這個(gè)模型里承擔(dān)什么角色為什么要把它放在內(nèi)層而不是外層第三整個(gè)Matlab工程該怎么組織才能既跑得通、又方便改。文章適合正在做能源系統(tǒng)優(yōu)化、綜合能源調(diào)度、電力市場(chǎng)策略方向的研究生也適合想用智能算法做博弈問(wèn)題但還沒(méi)理清思路的工程師。1. 先把問(wèn)題的骨架摸清三方三層到底在優(yōu)化什么1.1 主從博弈的本質(zhì)與能源系統(tǒng)場(chǎng)景主從博弈也叫Stackelberg博弈核心思想就是“l(fā)eader先出牌follower看牌后再行動(dòng)”。這個(gè)邏輯放到能源系統(tǒng)里非常貼合現(xiàn)實(shí)能源服務(wù)商售電公司、綜合能源運(yùn)營(yíng)商先定電價(jià)或者熱價(jià)用戶看到價(jià)格之后再調(diào)整自己的用能計(jì)劃運(yùn)營(yíng)商反過(guò)來(lái)又根據(jù)用戶的用能行為去修正定價(jià)策略。兩層決策層層嵌套每一層都是獨(dú)立的優(yōu)化問(wèn)題但目標(biāo)函數(shù)彼此耦合。比傳統(tǒng)單層優(yōu)化麻煩的地方在于你不能直接把所有決策變量丟進(jìn)一個(gè)優(yōu)化器里求解因?yàn)槊總€(gè)主體的目標(biāo)不一樣。運(yùn)營(yíng)商想多賺錢用戶想少花錢兩邊同時(shí)優(yōu)化得不到一個(gè)全局“最優(yōu)解”只能找一個(gè)博弈均衡點(diǎn)——在均衡點(diǎn)上任何一方單獨(dú)改變策略都不能讓自己變得更好。這在數(shù)學(xué)上對(duì)應(yīng)一類特殊的優(yōu)化問(wèn)題上層是帶下層反應(yīng)函數(shù)或KTT條件的約束優(yōu)化。處理方式要么是求下層問(wèn)題的KKT條件塞進(jìn)上層約束里MPEC/EPEC做法要么就是迭代逼近——先給上層一個(gè)策略下層求解后反饋每次根據(jù)反饋調(diào)整上層策略反復(fù)迭代直到收斂。后者對(duì)場(chǎng)景規(guī)模比較大、下層問(wèn)題不好解析求導(dǎo)的情況尤其友好也是粒子群這類元啟發(fā)式算法的用武之地。1.2 三層結(jié)構(gòu)如何對(duì)應(yīng)三方利益“三方”指的是博弈里的三個(gè)決策主體比如綜合能源服務(wù)商上級(jí)、儲(chǔ)能運(yùn)營(yíng)商或微網(wǎng)運(yùn)營(yíng)商中間層、終端用戶底層?！叭龑印敝傅氖菦Q策嵌套的層級(jí)關(guān)系上層定價(jià)中層做容量或調(diào)度決策底層做用能決策。用我自己做過(guò)的項(xiàng)目舉例一個(gè)典型的設(shè)定是這樣上層能源服務(wù)商決策變量是向終端用戶售電的電價(jià)、售熱的熱價(jià)目標(biāo)是最小化自身運(yùn)行成本或最大化利潤(rùn)。它要考慮從上級(jí)電網(wǎng)購(gòu)電的成本、自己運(yùn)營(yíng)分布式能源的成本并預(yù)估用戶對(duì)價(jià)格的響應(yīng)。中層儲(chǔ)能/微網(wǎng)運(yùn)營(yíng)商決策變量是儲(chǔ)能的充放電策略、購(gòu)售電功率目標(biāo)是在服務(wù)商給的結(jié)算價(jià)格和用戶負(fù)荷需求之間做一個(gè)最優(yōu)調(diào)度賺取峰谷價(jià)差或者減少自身購(gòu)電成本。底層終端用戶決策變量是可轉(zhuǎn)移負(fù)荷、可削減負(fù)荷的調(diào)整量目標(biāo)是在電價(jià)信號(hào)下最小化自己的用能成本同時(shí)保證舒適度約束不能打破。三層之間的信息流就是“上層定價(jià)→中層調(diào)度→下層用能”的順次傳遞。用戶調(diào)整完用能計(jì)劃后反饋給中間層一個(gè)負(fù)荷曲線中間層再反饋給上層一個(gè)總購(gòu)電需求上層根據(jù)總需求調(diào)整價(jià)格。這個(gè)過(guò)程循環(huán)往復(fù)每一次循環(huán)都是一次完整的“三層求解”。如果你把三層拆開分別建模會(huì)發(fā)現(xiàn)每一層單獨(dú)拿出來(lái)都不難關(guān)鍵是層與層之間的接口怎么設(shè)計(jì)。接口設(shè)計(jì)得好迭代收斂就快接口設(shè)計(jì)得抽象代碼改起來(lái)也省事。1.3 為什么偏偏選粒子群能源系統(tǒng)主從博弈的下層問(wèn)題還好說(shuō)很多時(shí)候是線性規(guī)劃或二次規(guī)劃用linprog、quadprog都能解。麻煩的是上層。上層的目標(biāo)是利潤(rùn)最大化同時(shí)牽涉到價(jià)格變量、需求響應(yīng)約束、設(shè)備運(yùn)行約束再把下層問(wèn)題的反應(yīng)函數(shù)疊加上來(lái)整個(gè)可行域嚴(yán)重非凸甚至不連續(xù)。傳統(tǒng)基于梯度的優(yōu)化工具在這種問(wèn)題上容易陷進(jìn)局部最優(yōu)或者根本不收斂。粒子群優(yōu)化算法PSO做這種事有三個(gè)天然優(yōu)勢(shì)第一它不要求目標(biāo)函數(shù)可導(dǎo)也不要求約束光滑跑出什么“臟”目標(biāo)值都能用第二全局搜索能力強(qiáng)在處理高維非凸問(wèn)題上比梯度法穩(wěn)得多第三Matlab實(shí)現(xiàn)非常方便三五十行代碼就能跑起一個(gè)基礎(chǔ)的粒子群核心。當(dāng)然粒子群也有短板比如后期容易早熟收斂、速度參數(shù)不好調(diào)但這些完全可以通過(guò)改進(jìn)策略來(lái)規(guī)避。后面的實(shí)操部分我會(huì)講怎么加自適應(yīng)慣性權(quán)重、怎么處理越界粒子、怎么在博弈迭代里嵌PSO才不會(huì)把時(shí)間耗死。2. 模型設(shè)計(jì)與參數(shù)標(biāo)定從數(shù)學(xué)問(wèn)題到可計(jì)算問(wèn)題2.1 數(shù)學(xué)表達(dá)式的設(shè)計(jì)與互補(bǔ)約束處理我見過(guò)不少同學(xué)在模型設(shè)計(jì)階段就把自己困住了。一上來(lái)就試圖把所有主體的目標(biāo)函數(shù)和約束都寫進(jìn)一個(gè)大模型里結(jié)果非線性約束、整數(shù)變量、互補(bǔ)條件堆在一起MATLAB根本跑不動(dòng)然后跑來(lái)問(wèn)是不是算法不行。實(shí)際上思路應(yīng)該反過(guò)來(lái)——先是把“誰(shuí)先動(dòng)、誰(shuí)后動(dòng)、層間傳什么”想清楚再拆分模塊化建模。一個(gè)我常用的切入方式是先把三層主體的目標(biāo)函數(shù)分別寫出來(lái)。上層服務(wù)商的目標(biāo)函數(shù)大致是f_up sum(price_e .* P_user_e price_h .* P_user_h) - C_purchase - C_ope;其中price_e和price_h是上層要定的電、熱價(jià)格向量P_user_e和P_user_h是用戶響應(yīng)后的負(fù)荷向量C_purchase是向電網(wǎng)購(gòu)電的成本C_ope是自己的設(shè)備運(yùn)行維護(hù)成本。中層的調(diào)度目標(biāo)大致是f_mid sum(price_grid .* P_grid) - sum(price_e .* P_user_e) C_bess;這里C_bess是儲(chǔ)能充放電的折舊成本通過(guò)把充放電功率折算成統(tǒng)一的成本系數(shù)實(shí)現(xiàn)。底層用戶的目標(biāo)比較簡(jiǎn)單f_low sum(price_e .* P_user_e price_h .* P_user_h) C_discomfort;C_discomfort是用戶調(diào)整用能計(jì)劃引起的舒適度損失可以當(dāng)成懲罰項(xiàng)。約束條件方面最核心的幾類是功率平衡約束系統(tǒng)內(nèi)發(fā)電/購(gòu)電功率 用戶負(fù)荷 儲(chǔ)能充放電 網(wǎng)損儲(chǔ)能約束荷電狀態(tài)SOC上下限、充放電功率上下限、充放電效率、日始末SOC相等用戶舒適度約束可轉(zhuǎn)移負(fù)荷的轉(zhuǎn)移量范圍、室內(nèi)溫度允許波動(dòng)范圍價(jià)格約束價(jià)格要落在政府指導(dǎo)價(jià)或市場(chǎng)限價(jià)區(qū)間內(nèi)。還有一個(gè)常見的坑是互補(bǔ)約束。比如儲(chǔ)能不能同時(shí)充電和放電如果直接寫P_ch 0和P_dis 0再加二值變量在非商業(yè)求解器里會(huì)拖垮速度。工程上常用三種處理辦法第一種是引入0-1變量需要求解器支持混合整數(shù)第二種是加小規(guī)模罰函數(shù)把“同時(shí)充放”的量做成懲罰項(xiàng)加進(jìn)目標(biāo)函數(shù)第三種是干脆用P_bess P_dis - P_ch單變量建模用效率分段表示避免同時(shí)性。我自己的習(xí)慣是如果儲(chǔ)能數(shù)量不多直接引入一個(gè)0-1變量布爾求解在Matlab里用intlinprog就能做如果儲(chǔ)能多到幾百個(gè)節(jié)點(diǎn)規(guī)模那就用懲罰項(xiàng)近似。這個(gè)取舍對(duì)收斂速度影響很大。2.2 粒子群參數(shù)設(shè)置粒子群算法要設(shè)的參數(shù)不多但每一個(gè)都直接影響收斂行為。標(biāo)準(zhǔn)粒子群的核心公式是速度更新和位置更新v w * v c1 * r1 * (pbest - x) c2 * r2 * (gbest - x); x x v;其中w是慣性權(quán)重c1和c2是學(xué)習(xí)因子r1和r2是[0,1]之間的隨機(jī)數(shù)pbest是粒子個(gè)體歷史最優(yōu)gbest是全局最優(yōu)。我調(diào)試過(guò)幾次之后發(fā)現(xiàn)對(duì)能源系統(tǒng)這樣的高維連續(xù)優(yōu)化問(wèn)題一套比較穩(wěn)的參數(shù)配置是參數(shù)推薦取值說(shuō)明種群規(guī)模30~60狀態(tài)變量多時(shí)往高取但要平衡耗時(shí)迭代次數(shù)100~300主從博弈外層迭代里內(nèi)層要重復(fù)調(diào)用適當(dāng)縮小慣性權(quán)重w0.9→0.4線性遞減前期全局搜索后期局部收斂c1個(gè)體學(xué)習(xí)因子1.5~2.0越大越依賴個(gè)體歷史最優(yōu)c2社會(huì)學(xué)習(xí)因子1.5~2.0越大越依賴群體最優(yōu)速度上限變量范圍的10%~20%防止粒子飛得太離譜一個(gè)很容易被忽略的細(xì)節(jié)是速度初始化。很多初學(xué)者直接把v設(shè)成全零矩陣這樣前幾次迭代粒子幾乎不移動(dòng)白白浪費(fèi)代數(shù)。更好的做法是用變量范圍的比例隨機(jī)初始化速度比如取變量可行域?qū)挾鹊?0%。另一個(gè)細(xì)節(jié)是約束處理。PSO本身不直接處理約束要么用罰函數(shù)把違反約束的量折算進(jìn)目標(biāo)值要么直接對(duì)越界的粒子做“位置修正速度重置”。我自己更推薦速度重置的做法當(dāng)粒子某維越界時(shí)不僅把位置拉回邊界還把該維速度置零。這樣粒子不會(huì)繼續(xù)朝邊界外沖迭代穩(wěn)定性高很多。2.3 博弈迭代與PSO內(nèi)層求解的耦合方式博弈迭代和粒子群不是簡(jiǎn)單的“一個(gè)套一個(gè)”關(guān)系具體耦合方式要看你把哪個(gè)問(wèn)題交給PSO。我常用的做法是上層服務(wù)商的策略價(jià)格用PSO來(lái)尋優(yōu)中層和下層用Matlab的優(yōu)化工具箱linprog或quadprog精確求解。每次PSO迭代產(chǎn)生一組價(jià)格方案就調(diào)用一次中下層求解模塊得到用戶和中間商的最優(yōu)反應(yīng)然后把反應(yīng)結(jié)果帶回上層目標(biāo)函數(shù)算出適應(yīng)度。這個(gè)方案的好處是PSO只負(fù)責(zé)最不光滑的上層問(wèn)題中下層仍然用確定性求解器速度和穩(wěn)定性都有保障。還有一種做法是三層全用PSO形成“外層PSO套中層PSO再套內(nèi)層PSO”的俄羅斯套娃結(jié)構(gòu)。說(shuō)實(shí)話這種做法在論文里寫起來(lái)漂亮但實(shí)踐起來(lái)很難收斂因?yàn)閮?nèi)層的隨機(jī)擾動(dòng)會(huì)不斷傳遞放大最后外層看到的適應(yīng)度噪聲非常大粒子根本分不清哪個(gè)位置是真正的好位置。如果非要用全PSO結(jié)構(gòu)內(nèi)層的迭代次數(shù)至少要壓到50代以內(nèi)而且要固定隨機(jī)種子保證每次評(píng)價(jià)同一方案時(shí)結(jié)果一致不然優(yōu)化完全失去意義。我個(gè)人的經(jīng)驗(yàn)是能用精確求解器的地方不要輕易用元啟發(fā)式。博弈迭代本身已經(jīng)很耗時(shí)了再層層嵌套隨機(jī)算法調(diào)試體驗(yàn)會(huì)非常糟糕。先跑通簡(jiǎn)單方案再逐步考慮增復(fù)雜度比一次性搭大而全的框架靠譜得多。3. Matlab實(shí)現(xiàn)代碼框架與核心環(huán)節(jié)逐步實(shí)現(xiàn)3.1 主循環(huán)架構(gòu)一個(gè)實(shí)際能跑的Matlab工程我建議按“數(shù)據(jù)輸入——參數(shù)初始化——主從迭代——結(jié)果輸出”四個(gè)模塊來(lái)組織。不要把所有代碼堆在一個(gè)文件里否則到后面你會(huì)被自己剛寫的代碼勸退。我一般這樣組織目錄|-- main.m % 主程序入口 |-- data/ | |-- load_data.m % 加載負(fù)荷、價(jià)格、設(shè)備參數(shù) |-- model/ | |-- upper_objective.m % 上層目標(biāo)函數(shù) | |-- middle_solve.m % 中層調(diào)度求解 | |-- lower_solve.m % 底層需求響應(yīng)求解 |-- algorithm/ | |-- pso_main.m % 粒子群主程序 | |-- pso_update.m % 粒子速度位置更新 |-- result/ | |-- plot_result.m % 繪圖與結(jié)果分析主程序的循環(huán)邏輯用一個(gè)while來(lái)實(shí)現(xiàn)博弈迭代。%% main.m 主循環(huán)框架 % 初始化 price_iter price0; % 初始價(jià)格方案上層決策 max_iter 20; % 主從博弈最大迭代次數(shù) tol 1e-3; % 收斂精度 for k 1:max_iter % 1. 給定價(jià)格求解中層調(diào)度問(wèn)題 [P_mid, f_mid] middle_solve(price_iter, system_data); % 2. 給定價(jià)格和中層調(diào)度結(jié)果求解底層用戶響應(yīng)用量 [P_low, f_low] lower_solve(price_iter, system_data); % 3. 把中下層各主體的響應(yīng)結(jié)果返回上層代入上層目標(biāo)函數(shù)的適應(yīng)度計(jì)算 fitness upper_objective(price_iter, P_mid, P_low, system_data); % 4. 用PSO更新價(jià)格策略 [price_new, swarm_state] pso_main(swarm_state, system_data); % 5. 收斂判斷 if norm(price_new - price_iter) tol break; end price_iter price_new; end這里主從博弈疊了個(gè)PSO兩個(gè)循環(huán)清清楚楚所以不容易把邏輯搞混。max_iter設(shè)20次左右的用意是一般而言主從博弈的定價(jià)策略經(jīng)過(guò)十幾次迭代就能趨于穩(wěn)定如果超過(guò)30次還不收斂大概率是模型參數(shù)或者迭代步長(zhǎng)設(shè)置有問(wèn)題跑再多輪也白搭。3.2 PSO核心代碼實(shí)戰(zhàn)粒子群主程序我通常寫成一個(gè)獨(dú)立的m函數(shù)這樣換場(chǎng)景時(shí)可以直接復(fù)制調(diào)用。function [gbest, gbest_fit, history] pso_main(fitness_func, dim, lb, ub, opts) % 粒子群優(yōu)化主函數(shù) % fitness_func: 函數(shù)句柄輸入位置向量輸出適應(yīng)度值 % dim: 決策變量維度 % lb, ub: 變量上下界向量 % opts: 結(jié)構(gòu)體包含 nP, maxIter, wMax, wMin, c1, c2 nP opts.nP; % 種群規(guī)模 maxIter opts.maxIter; % 最大迭代次數(shù) wMax opts.wMax; wMin opts.wMin; c1 opts.c1; c2 opts.c2; % 初始化位置和速度 X repmat(lb, nP, 1) rand(nP, dim) .* repmat(ub - lb, nP, 1); V repmat(lb - ub, nP, 1) .* (0.1 * rand(nP, dim)); % 速度用可行域?qū)挾鹊?0%初始化 % 初始化個(gè)體最優(yōu)和全局最優(yōu) pbest X; pbest_fit arrayfun((i) fitness_func(X(i, :)), 1:nP); gbest pbest(1, :); gbest_fit pbest_fit(1); for i 2:nP if pbest_fit(i) gbest_fit gbest pbest(i, :); gbest_fit pbest_fit(i); end end % 迭代 history zeros(maxIter, 1); for iter 1:maxIter w wMax - (wMax - wMin) * iter / maxIter; % 慣性權(quán)重線性遞減 for i 1:nP r1 rand(1, dim); r2 rand(1, dim); V(i, :) w * V(i, :) c1 * r1 .* (pbest(i, :) - X(i, :)) c2 * r2 .* (gbest - X(i, :)); % 速度限幅 Vmax 0.2 * (ub - lb); V(i, :) max(min(V(i, :), Vmax), -Vmax); % 位置更新 X(i, :) X(i, :) V(i, :); % 越界處理拉回邊界速度置零 for d 1:dim if X(i, d) lb(d) || X(i, d) ub(d) X(i, d) min(max(X(i, d), lb(d)), ub(d)); V(i, d) 0; end end % 適應(yīng)度評(píng)價(jià) fit fitness_func(X(i, :)); if fit pbest_fit(i) pbest_fit(i) fit; pbest(i, :) X(i, :); if fit gbest_fit gbest_fit fit; gbest X(i, :); end end end history(iter) gbest_fit; end end這段代碼有幾個(gè)值得拎出來(lái)細(xì)講的點(diǎn)一是V的初始化不要用全零。全零初速度會(huì)讓前幾輪粒子只在原地小幅震蕩搜索效率極低。用可行域?qū)挾鹊?0%作為初始速度區(qū)間能讓粒子一開始就有足夠的探索動(dòng)能。二是“越界拉回速度置零”。很多新手抄PSO代碼時(shí)直接讓越界粒子的位置等于邊界值但速度不處理結(jié)果粒子下一輪更新時(shí)又沖出邊界在邊界來(lái)回彈跳。把越界維度的速度置零后粒子被“摁”在邊界上需要靠其他維度的社會(huì)/個(gè)體學(xué)習(xí)項(xiàng)把它帶回來(lái)穩(wěn)定性好很多。三是適應(yīng)度函數(shù)統(tǒng)一寫成最小化形式。上層服務(wù)商的利潤(rùn)最大化和目標(biāo)函數(shù)要加負(fù)號(hào)或者整體取倒數(shù)否則粒子群會(huì)朝適應(yīng)度最大的方向飛和代碼邏輯沖突。3.3 中下層求解模塊與目標(biāo)函數(shù)銜接中層的儲(chǔ)能調(diào)度問(wèn)題如果規(guī)模不大直接用linprog就行。以典型的24小時(shí)調(diào)度問(wèn)題為例決策變量是每個(gè)小時(shí)的儲(chǔ)能充放電功率維度48如果充電、放電分開兩個(gè)變量或者24用凈功率單變量建模。function [P_bess, f_mid, exitflag] middle_solve(price_iter, data) % 中層儲(chǔ)能調(diào)度min C_grid * P_grid - revenue C_bess % 采用凈功率建模P_bess 0 表示放電P_bess 0 表示充電 T data.T; % 24小時(shí) P_load data.P_load; % 用戶預(yù)測(cè)負(fù)荷 % 決策變量: [P_bess(1:T), SOC(1:T)] nVar 2 * T; % 目標(biāo)函數(shù)系數(shù)向量 f [zeros(1, T), zeros(1, T)]; % 具體系數(shù)根據(jù)場(chǎng)景調(diào)整 % 等式約束功率平衡 Aeq [eye(T), zeros(T, T)]; beq P_load - price_iter .* 0; % 等等這里要根據(jù)場(chǎng)景修改占位方便說(shuō)明 % 邊界條件 lb [-data.P_bess_max * ones(1, T), 0.2 * ones(1, T)]; ub [data.P_bess_max * ones(1, T), 0.9 * ones(1, T)]; options optimoptions(linprog, Display, off); [x, f_mid, exitflag] linprog(f, A, b, Aeq, beq, lb, ub, options); if exitflag 0 error(中層求解失敗檢查約束是否矛盾); end P_bess x(1:T); end底層用戶需求響應(yīng)部分同理根據(jù)可轉(zhuǎn)移負(fù)荷和可削減負(fù)荷的比例構(gòu)建線性規(guī)劃。如果引入舒適度溫度約束可能需要用到二次規(guī)劃quadprog因?yàn)橛脩舻氖孢m度損失函數(shù)寫成二次形式更自然。中下層寫完以后最要緊的一步是用測(cè)試數(shù)據(jù)單測(cè)。我見過(guò)太多人直接把三層代碼串起來(lái)一跑就是幾個(gè)小時(shí)結(jié)果最后發(fā)現(xiàn)中層功率平衡約束寫反了符號(hào)全部結(jié)果都是錯(cuò)的。單層測(cè)試時(shí)固定另外兩層傳入的數(shù)據(jù)檢查這層求解出來(lái)的變量是否物理上合理花不了多少時(shí)間但能省下大量調(diào)試時(shí)間。4. 常見問(wèn)題與排查實(shí)錄4.1 粒子群早熟收斂結(jié)果明顯不是最優(yōu)解這個(gè)大概是使用PSO過(guò)程里最常碰到的事。能源系統(tǒng)優(yōu)化問(wèn)題往往決策變量多如果所有維度的搜索策略都一樣很快就全體朝某個(gè)局部最優(yōu)靠攏粒子多樣性喪失。我自己的處理思路是先把線性遞減慣性權(quán)重的范圍拉開比如從0.95降到0.35同時(shí)把c1和c2設(shè)成非對(duì)稱比如c11.8、c21.2讓粒子前期更偏向自我搜索不要過(guò)早被群體最優(yōu)帶走。要是跑了好幾遍都一樣那就是初始種群沒(méi)覆蓋好可以嘗試用Sobol序列或者拉丁超立方抽樣生成初始位置避免隨機(jī)抽樣導(dǎo)致種群在可行域里擠成一團(tuán)。還有一個(gè)小技巧是每隔一定迭代代數(shù)對(duì)適應(yīng)度最差的那批粒子做“重新初始化”把它們的速度和位置隨機(jī)重置相當(dāng)于給種群注入新鮮血液。4.2 罰函數(shù)系數(shù)取值敏感輕微調(diào)整結(jié)果差異巨大罰函數(shù)是處理約束時(shí)最簡(jiǎn)單的方式但系數(shù)設(shè)小了約束根本滿足不了設(shè)大了目標(biāo)函數(shù)的真實(shí)梯度又會(huì)被淹沒(méi)。我的建議是不要在全模型用同一個(gè)罰函數(shù)系數(shù)。按約束類型分類懲罰功率平衡約束的懲罰系數(shù)取量級(jí)比較大的數(shù)比如1e4設(shè)備容量約束的懲罰系數(shù)可以小一些比如1e2到1e3。原因很簡(jiǎn)單功率平衡偏差的物理量綱本身就比設(shè)備容量超限的偏差大得多懲罰系數(shù)也必須跟著匹配。如果非要用罰函數(shù)建議把違反約束的量做歸一化處理。比如儲(chǔ)能SOC越界量除以其允許波動(dòng)范圍這樣不同尺度的約束在懲罰項(xiàng)里權(quán)重才一致。4.3 博弈迭代不收斂上層價(jià)格在幾個(gè)值之間來(lái)回震蕩這是主從博弈嵌套優(yōu)化里非常經(jīng)典的現(xiàn)象。價(jià)格迭代時(shí)會(huì)出現(xiàn)在兩個(gè)方向同時(shí)偏移的“鋸齒形”震蕩根源往往是迭代步長(zhǎng)太大或者用戶響應(yīng)函數(shù)對(duì)價(jià)格過(guò)于敏感。兩個(gè)可行的解決辦法一是給價(jià)格迭代加阻尼每次更新時(shí)不完全跳到新值而是取新舊值的加權(quán)平均price_new alpha * price_pso (1 - alpha) * price_old;alpha取0.6左右比較穩(wěn)。二是檢查用戶需求響應(yīng)模型是否合理如果需求彈性系數(shù)過(guò)大價(jià)格稍微一波動(dòng)用戶負(fù)荷就劇烈變化上層看到的目標(biāo)函數(shù)當(dāng)然也會(huì)劇烈變化。很多時(shí)候是模型參數(shù)本身沒(méi)標(biāo)定好不是算法的問(wèn)題。4.4 算法整體運(yùn)行太慢主從迭代加PSO雙重循環(huán)跑不動(dòng)這種雙重循環(huán)結(jié)構(gòu)跑起來(lái)確實(shí)慢尤其是當(dāng)?shù)讓佑胠inprog、中層用quadprog每個(gè)粒子都要調(diào)一次完整中下層求解的時(shí)候。24小時(shí)場(chǎng)景、36個(gè)小時(shí)的場(chǎng)景60個(gè)粒子30次博弈迭代每次迭代里每個(gè)粒子都要調(diào)度一次中下層累計(jì)求解次數(shù)驚人Matlab單線程下可能要跑上幾個(gè)小時(shí)。我優(yōu)化這個(gè)瓶頸的方法主要有三個(gè)減少不必要的粒子評(píng)價(jià)。PSO迭代前期粒子位置變化很大但很多位置的適應(yīng)度差得很明顯根本不影響全局最優(yōu)就沒(méi)必要每次都調(diào)用完整求解器。思路是先把粒子位置做粗糙的可行域檢查明顯不合理的位置直接賦個(gè)很大的懲罰值不進(jìn)求解器。中下層求解器選型上能用線性規(guī)劃就不要上非線性規(guī)劃。很多儲(chǔ)能調(diào)度問(wèn)題在目標(biāo)函數(shù)和約束都是線性的前提下直接用linprog比用fmincon快一個(gè)數(shù)量級(jí)。如果單次循環(huán)時(shí)間實(shí)在降不下去改用并行計(jì)算。Matlab的parfor可以把每個(gè)粒子的適應(yīng)度評(píng)價(jià)分散到多核通信開銷可接受提速基本能到核心數(shù)級(jí)別。還有一個(gè)隱藏很深的坑是主從迭代過(guò)程中每層都在重復(fù)加載數(shù)據(jù)或重新構(gòu)造稀疏矩陣這些“固定開銷”往往被忽略。把不隨迭代變化的數(shù)據(jù)提前算好存入結(jié)構(gòu)體每次只傳遞句柄能省下不少時(shí)間。4.5 常見問(wèn)題速查現(xiàn)象可能原因排查順序與對(duì)策PSO結(jié)果每次跑都不一樣且差異大隨機(jī)性太強(qiáng)種群初始化覆蓋差固定隨機(jī)種子復(fù)現(xiàn)改用拉丁超立方初始化加大種群規(guī)模結(jié)果符合約束但經(jīng)濟(jì)性差罰函數(shù)系數(shù)把目標(biāo)函數(shù)壓沒(méi)了降低懲罰系數(shù)分類型設(shè)置檢查目標(biāo)函數(shù)是否無(wú)意中包含了懲罰項(xiàng)中下層求解器報(bào)錯(cuò)“無(wú)可行解”約束過(guò)強(qiáng)或傳入的接口數(shù)據(jù)異常單獨(dú)測(cè)試中下層輸入輸出先放松一兩個(gè)約束定位矛盾源價(jià)格迭代劇烈震蕩迭代步長(zhǎng)過(guò)大或需求彈性過(guò)大加阻尼減小價(jià)格更新幅度檢查價(jià)格上下限是否過(guò)寬整體耗時(shí)太長(zhǎng)雙重循環(huán)中每個(gè)粒子都調(diào)用完整求解器并行計(jì)算預(yù)檢查減少無(wú)效評(píng)價(jià)盡量把中下層寫成線性規(guī)劃上層目標(biāo)函數(shù)值出現(xiàn)NaN下層求解失敗或罰函數(shù)溢出用exitflag捕獲下層求解狀態(tài)對(duì)NaN直接賦值巨大懲罰值5. 實(shí)操心得與后續(xù)擴(kuò)展調(diào)這個(gè)模型的過(guò)程中我最大的體會(huì)是主從博弈加粒子群最大的難點(diǎn)不在算法本身而在“接口設(shè)計(jì)”。三層模型里每一層的數(shù)據(jù)格式、變量順序、約束矩陣構(gòu)造方式必須從一開始就定好規(guī)范否則代碼寫到后面層與層之間互相傳數(shù)組都對(duì)不上排錯(cuò)排到頭大。一個(gè)小技巧是所有層之間的數(shù)據(jù)傳遞都用結(jié)構(gòu)體統(tǒng)一格式例如layers.up.price、layers.mid.P_bess、layers.low.P_shift這樣在斷點(diǎn)調(diào)試時(shí)一目了然也能避免把行向量傳成列向量的經(jīng)典錯(cuò)誤。另外為了讓中間過(guò)程可解釋建議每一輪主從迭代把三個(gè)層的目標(biāo)函數(shù)值記錄到數(shù)組里最后畫在一張圖中。這樣你能直觀看到三個(gè)主體的利益演化軌跡——哪一方在博弈中占據(jù)優(yōu)勢(shì)、哪一方在讓利一眼便知。這張圖放在論文里也很有說(shuō)服力。關(guān)于后續(xù)擴(kuò)展你還可以往這幾個(gè)方向深挖把粒子群改成多目標(biāo)版本MOPSO同時(shí)考慮經(jīng)濟(jì)性和碳排放兩個(gè)目標(biāo)畫Pareto前沿上層引入多種類型售能主體電、熱、氣形成“多領(lǐng)導(dǎo)者—多跟隨者”的復(fù)雜博弈結(jié)構(gòu)把用戶側(cè)真實(shí)需求響應(yīng)歷史數(shù)據(jù)接進(jìn)來(lái)替代數(shù)學(xué)響應(yīng)模型將PSO換成改進(jìn)的量子粒子群或混合灰狼PSO提升高維下的搜索能力。不過(guò)我的建議是先把基礎(chǔ)的單一PSO版本跑通、跑穩(wěn)把所有模塊驗(yàn)證到位再逐漸往上加復(fù)雜度。我見過(guò)太多人一上來(lái)就想一步到位做多目標(biāo)多主體結(jié)果卡在調(diào)試上兩個(gè)月心態(tài)直接崩掉。一步一步來(lái)先把一個(gè)能收斂的版本攥在手里再考慮升級(jí)的事。最后再分享一個(gè)小技巧無(wú)論你用什么智能算法做博弈優(yōu)化一定要在代碼里固定隨機(jī)種子rng(2024)這種否則你改一次參數(shù)、重跑一遍結(jié)果可能面目全非根本沒(méi)法判斷到底是改動(dòng)生效了還是純隨機(jī)波動(dòng)。固定種子之后算法之間的細(xì)微差別才真正可比較、可復(fù)現(xiàn)。這是我踩過(guò)無(wú)數(shù)次坑之后最想提醒你的一件事。