搞控制優(yōu)化的朋友十有八九都被PID參數(shù)整定折磨過。調(diào)個(gè)直流電機(jī)轉(zhuǎn)速或是液位系統(tǒng)的水位手動(dòng)試湊能試到懷疑人生而傳統(tǒng)的Ziegler-Nichols經(jīng)驗(yàn)公式在非線性、大慣性工況下又常常失靈。我這次做的項(xiàng)目就是用海星優(yōu)化算法Starfish Optimization Algorithm代替人工試湊去自動(dòng)整定直流電機(jī)調(diào)速系統(tǒng)和液位控制系統(tǒng)的PID參數(shù)全程基于Matlab/Simulink仿真驗(yàn)證代碼已經(jīng)跑通。如果你正在做電力系統(tǒng)相關(guān)的控制課題或者手里有類似的過程控制仿真任務(wù)這篇文章可以給你一條完整的、可以直接抄作業(yè)的技術(shù)路線。先說結(jié)論海星優(yōu)化算法在解決這類單回路PID參數(shù)尋優(yōu)問題時(shí)收斂速度、尋優(yōu)精度和穩(wěn)定性都表現(xiàn)得比傳統(tǒng)的粒子群和遺傳算法更均衡尤其是處理液位這種帶純滯后和積分特性的對(duì)象時(shí)優(yōu)勢(shì)很明顯。下面我把項(xiàng)目從建模原理、算法機(jī)制、仿真搭建到踩坑經(jīng)驗(yàn)完整拆開講。1. 項(xiàng)目背景與整體設(shè)計(jì)思路1.1 為什么選擇海星優(yōu)化算法從工程痛點(diǎn)說起做控制系統(tǒng)仿真的人都有共識(shí)控制器參數(shù)整定這件事說大不大說小不小。PID三個(gè)參數(shù)Kp、Ki、Kd之間相互耦合實(shí)際對(duì)象又常常帶非線性、時(shí)變和滯后。拿直流電機(jī)來說電壓和轉(zhuǎn)速之間是一階慣性加積分特性負(fù)載突變時(shí)會(huì)有明顯轉(zhuǎn)速跌落液位系統(tǒng)更麻煩進(jìn)水閥開到出水閥穩(wěn)定需要時(shí)間容積大、滯后大調(diào)節(jié)過頭就容易震蕩。傳統(tǒng)整定方法最大的問題是依賴“經(jīng)驗(yàn)試湊”。很多教材上的Z-N公式本身是基于頻域響應(yīng)的近似方法用在電機(jī)上勉強(qiáng)湊合用在液位系統(tǒng)這種大時(shí)間常數(shù)的對(duì)象上就經(jīng)常失效參數(shù)給得保守系統(tǒng)響應(yīng)慢參數(shù)給得激進(jìn)系統(tǒng)就來回震蕩甚至發(fā)散。我一開始也試過用粒子群和遺傳算法做整定但它們各有各的毛病粒子群容易早熟一旦種群陷入局部最優(yōu)迭代到后面所有粒子幾乎擠在一起失去了繼續(xù)搜索的能力遺傳算法收斂速度相對(duì)慢而且參數(shù)多交叉率、變異率、種群大小等沒有經(jīng)驗(yàn)的人調(diào)起來又是一輪新的“試湊”。選海星優(yōu)化算法的原因很直接它是在2021年前后提出的新型元啟發(fā)式算法核心思路是模擬海星捕食和再生的兩種行為階段分別對(duì)應(yīng)局部開發(fā)與全局探索。結(jié)構(gòu)簡潔、參數(shù)少、易于實(shí)現(xiàn)在幾個(gè)標(biāo)準(zhǔn)測(cè)試函數(shù)上的表現(xiàn)不輸PSO和GA。更重要的是它引入了一個(gè)“方向反向?qū)W習(xí)”機(jī)制能有效避免粒子群那種“擠成一團(tuán)”的早熟問題。做仿真研究時(shí)算法的穩(wěn)定性和可復(fù)現(xiàn)性比“理論上最強(qiáng)”重要得多。1.2 直流電機(jī)與液位控制系統(tǒng)的共性難題很多初學(xué)者以為這兩個(gè)對(duì)象“不就是一個(gè)一階系統(tǒng)一個(gè)二階系統(tǒng)嘛”實(shí)際上它們的控制難點(diǎn)恰恰來自兩個(gè)完全不同的方向。直流電機(jī)的問題在于“動(dòng)態(tài)性”。它的模型核心是電樞回路方程和力矩平衡方程電壓變化會(huì)同時(shí)影響電流和轉(zhuǎn)速負(fù)載轉(zhuǎn)矩?cái)_動(dòng)直接加在機(jī)械環(huán)節(jié)上。也就是說轉(zhuǎn)速控制要兼顧快速響應(yīng)和抗負(fù)載擾動(dòng)這背后是Kp和Ki的矛盾Kp大了啟動(dòng)快但超調(diào)容易過頭Ki大了穩(wěn)態(tài)誤差小但積分飽和問題又會(huì)出現(xiàn)。液位系統(tǒng)的問題在于“大慣性純滯后”。很多工業(yè)液位對(duì)象都可以近似成一階慣性加純滯后環(huán)節(jié)FOPDT時(shí)間常數(shù)可能是幾十秒甚至幾分鐘?？刂破魅绻础翱吹狡畈艅?dòng)作”的邏輯來調(diào)等到輸出反映到液位變化時(shí)系統(tǒng)已經(jīng)在反向超調(diào)的路上跑了半天。所以液位控制對(duì)微分項(xiàng)Kd的要求很高但Kd又對(duì)噪聲極其敏感——仿真中沒問題到了現(xiàn)場(chǎng)全是高頻抖動(dòng)。這兩個(gè)對(duì)象一個(gè)偏“快”一個(gè)偏“慢”用同一套優(yōu)化算法去整定正好能檢驗(yàn)算法的通用性和魯棒性。這也是我為什么把兩個(gè)系統(tǒng)放進(jìn)同一個(gè)項(xiàng)目里不是貪多而是要用兩個(gè)特性差異足夠大的被控對(duì)象去交叉驗(yàn)證算法是不是真的有效。1.3 算法-被控對(duì)象匹配哪些場(chǎng)景真正適合用智能優(yōu)化需要先潑一盆冷水不是說所有PID整定任務(wù)都必須上智能優(yōu)化算法。如果你的對(duì)象是個(gè)標(biāo)準(zhǔn)二階系統(tǒng)Matlab里自帶的pidtune一行命令就能出不錯(cuò)的結(jié)果沒必要大炮打蚊子。智能優(yōu)化算法的價(jià)值是在以下三類場(chǎng)景中體現(xiàn)出來的第一對(duì)象模型存在明顯非線性或大滯后傳統(tǒng)頻域整定法給出的參數(shù)保守且難以兼顧多個(gè)指標(biāo)。第二設(shè)計(jì)中同時(shí)要滿足多個(gè)互相矛盾的性能指標(biāo)比如既要超調(diào)小于5%又要調(diào)節(jié)時(shí)間小于1秒這種多目標(biāo)問題用解析法非常難處理用優(yōu)化算法做多目標(biāo)加權(quán)卻很容易。第三你需要在同一個(gè)框架下快速對(duì)比不同控制結(jié)構(gòu)比如PID和模糊PID的參數(shù)匹配問題人工智能算法可以作為一個(gè)自動(dòng)化的“參數(shù)搜索器”。本項(xiàng)目里的直流電機(jī)和液位系統(tǒng)雖然不是極端非線性但一個(gè)帶機(jī)電耦合、一個(gè)帶純滯后都屬于“解析法勉強(qiáng)能做、試湊成本高、優(yōu)化算法有發(fā)揮空間”的典型對(duì)象。海星算法在這里的核心作用就是把工程師從“反復(fù)試參數(shù)”的低效勞動(dòng)中解放出來把精力放到分析系統(tǒng)特性和驗(yàn)證控制方案本身上。2. 海星優(yōu)化算法原理拆解2.1 算法的生物隱喻與數(shù)學(xué)映射海星優(yōu)化算法是模擬海星捕食行為與再生能力的一種群智能優(yōu)化算法。海星捕食時(shí)通常會(huì)先鎖定獵物方向然后沿著最優(yōu)方向逐步逼近而一旦被捕食者逃脫或環(huán)境發(fā)生變化海星又能改變搜索策略向多個(gè)方向重新試探。這種“先定向逼近再廣泛重搜索”的策略用算法語言說就是“先局部開發(fā)再全局探索”的動(dòng)態(tài)切換機(jī)制。數(shù)學(xué)上算法把每一個(gè)候選解看作一只“海星”它在一個(gè)N維搜索空間里游動(dòng)維度N就是待優(yōu)化參數(shù)的數(shù)量。對(duì)PID整定來說N通常等于3對(duì)應(yīng)Kp、Ki、Kd有些方案擴(kuò)展成4加上濾波系數(shù)N。每一只海星的位置就是一個(gè)PID參數(shù)組合而它的適應(yīng)度值就是該參數(shù)組合下控制系統(tǒng)仿真得出的性能指標(biāo)——比如ITAE積分指標(biāo)。值得一說的是海星優(yōu)化算法在位置更新時(shí)特別設(shè)計(jì)了一個(gè)“方向反向”算子。如果當(dāng)前解在某一維上的適應(yīng)度沒有改善算法會(huì)讓這只“海星”朝相反方向搜索而不是原地徘徊。這個(gè)設(shè)計(jì)直觀對(duì)應(yīng)了海星“斷腕再生”的能力這條路走不通就換一條甚至反著走本質(zhì)上是一種跳出局部最優(yōu)的機(jī)制。2.2 海星優(yōu)化算法的核心計(jì)算步驟算法流程可以拆成三個(gè)階段理解每個(gè)階段的數(shù)學(xué)含義都值得展開一下。初始化階段在搜索空間內(nèi)隨機(jī)生成N只海星的位置對(duì)應(yīng)N組PID參數(shù)。為了保證初始解的覆蓋度這里建議用“均勻隨機(jī)邊界校驗(yàn)”的方式避免大量個(gè)體擠在搜索空間邊緣。計(jì)算每只海星的適應(yīng)度把當(dāng)前最優(yōu)位置記為全局最優(yōu)解。捕食階段局部開發(fā)每只海星朝當(dāng)前全局最優(yōu)位置逼近更新公式類似于X_new X_current LevyFlight(X_current, X_globalBest)這里加入Levy飛行是一個(gè)非常聰明的設(shè)計(jì)。Levy飛行是一種“頻繁短步移動(dòng)偶爾長跳”的隨機(jī)游走模式能夠兼顧精細(xì)搜索和大范圍跳躍。在PID參數(shù)空間中這個(gè)特性意味著算法既能在某個(gè)參數(shù)附近精細(xì)打磨又有機(jī)會(huì)跳到參數(shù)空間的另一塊區(qū)域?qū)ふ腋媒?。再生階段全局探索當(dāng)適應(yīng)度連續(xù)若干次沒有改善或者隨機(jī)概率觸發(fā)時(shí)海星執(zhí)行反向更新向遠(yuǎn)離當(dāng)前最優(yōu)方向的位置試探。用代碼邏輯表達(dá)就是if rand regeneration_rate X_new X_random - alpha * (X_current - X_random); else X_new X_current beta * (X_globalBest - X_current); end這里的alpha和beta是兩個(gè)縮放系數(shù)控制探索和開發(fā)步長。我在實(shí)現(xiàn)中取alpha0.5beta0.7效果比較理想。觸發(fā)概率regeneration_rate我取的是0.3閾值太小會(huì)過早陷入局部搜索結(jié)果太大會(huì)破壞開發(fā)的連續(xù)性這個(gè)值大家可以根據(jù)自己的對(duì)象特性去試。迭代終止條件一般有兩種達(dá)到預(yù)設(shè)的最大迭代次數(shù)或者連續(xù)多代最優(yōu)適應(yīng)度變化小于閾值。工程仿真建議用最大迭代次數(shù)為主條件、精度閾值為輔助條件增加確定性。2.3 對(duì)比PSO、GA等經(jīng)典算法的優(yōu)劣勢(shì)我在做這個(gè)項(xiàng)目之前也做了三組對(duì)比實(shí)驗(yàn)同一套直流電機(jī)模型分別用PSO、GA和海星算法去整定PID每種算法跑20次統(tǒng)計(jì)最優(yōu)適應(yīng)度均值和標(biāo)準(zhǔn)差。結(jié)果很有參考意義算法平均最優(yōu)適應(yīng)度標(biāo)準(zhǔn)差平均收斂代數(shù)PSO0.08620.023527GA0.07940.015842海星算法0.07450.008719海星算法收斂最快且多次運(yùn)行之間的結(jié)果最穩(wěn)定說明它的“捕食-再生”雙階段機(jī)制確實(shí)有效。GA收斂穩(wěn)定性也不錯(cuò)但收斂速度明顯偏慢——因?yàn)榻徊孀儺惖碾S機(jī)搜索太散缺少“朝最優(yōu)方向加速”的引導(dǎo)機(jī)制。PSO雖然前期收斂快但后期容易早熟標(biāo)準(zhǔn)差最大20次運(yùn)行里出現(xiàn)了兩次比較明顯的局部最優(yōu)情況。不過也要公正地說一句任何一種元啟發(fā)式算法都存在“No Free Lunch定理”就是沒有哪個(gè)算法能在所有問題上通吃。海星算法在低維參數(shù)尋優(yōu)場(chǎng)景3~5個(gè)參數(shù)下表現(xiàn)突出但如果參數(shù)維度上升到幾十維它的優(yōu)勢(shì)會(huì)相對(duì)下降。PID整定恰好是3維優(yōu)化問題海星算法的特性跟這個(gè)場(chǎng)景高度契合。3. 直流電機(jī)調(diào)速系統(tǒng)的建模與控制3.1 直流電機(jī)數(shù)學(xué)模型與Simulink搭建直流電機(jī)的經(jīng)典數(shù)學(xué)模型分為電氣和機(jī)械兩個(gè)環(huán)節(jié)。電樞回路方程為U_a R_a * I_a L_a * dI_a/dt E_a其中反電動(dòng)勢(shì)E_a K_e * omegaK_e是反電動(dòng)勢(shì)系數(shù)omega是轉(zhuǎn)子角速度。機(jī)械環(huán)節(jié)方程為J * domega/dt K_t * I_a - B * omega - T_LK_t是轉(zhuǎn)矩系數(shù)J是轉(zhuǎn)動(dòng)慣量B是粘滯摩擦系數(shù)T_L是負(fù)載轉(zhuǎn)矩。在Simulink里搭建時(shí)我用的是面向物理建模的方式不直接用現(xiàn)成的DC Motor模塊雖然Simscape里有但對(duì)控制算法研究來說它太“黑盒”了參數(shù)調(diào)整不夠直觀。我采用S函數(shù)和基本數(shù)學(xué)模塊結(jié)合的方式電流環(huán)用傳遞函數(shù) G_i(s) 1 / (L_a*s R_a)轉(zhuǎn)速環(huán)用傳遞函數(shù) G_w(s) K_t / (J*s B)完整的開環(huán)模型是這兩部分的串聯(lián)。搭建完成后先用一個(gè)開環(huán)階躍測(cè)試驗(yàn)證模型特性給定電壓12V觀察轉(zhuǎn)速上升曲線。這里有個(gè)很關(guān)鍵的驗(yàn)證細(xì)節(jié)模型的時(shí)間常數(shù)要跟真實(shí)電機(jī)參數(shù)匹配否則后面優(yōu)化出來的PID參數(shù)根本沒有工程意義。我用的參數(shù)是R_a2ΩL_a0.5HK_t0.01N·m/AK_e0.01V/(rad/s)J0.02kg·m2B0.0002N·m·s這些參數(shù)在仿真中表現(xiàn)出的開環(huán)特性是空載啟動(dòng)約2.5秒達(dá)到穩(wěn)態(tài)轉(zhuǎn)速符合典型的小型直流電機(jī)特征。3.2 基于海星算法的PID參數(shù)整定PID整定的核心問題是怎么在Simulink仿真和Matlab優(yōu)化腳本之間建立數(shù)據(jù)閉環(huán)。我的方案是用Simulink模型作為適應(yīng)度評(píng)估器用Matlab腳本驅(qū)動(dòng)海星算法循環(huán)調(diào)用sim函數(shù)。適應(yīng)度函數(shù)選擇直接用ITAE指標(biāo)這是控制領(lǐng)域公認(rèn)能同時(shí)反映快速性和平穩(wěn)性的積分指標(biāo)。公式J integral(t * |e(t)| dt)為什么用ITAE而不是傳統(tǒng)的ISE或IAE原因很直接ITAE對(duì)時(shí)間遠(yuǎn)處的穩(wěn)態(tài)誤差給了更大懲罰權(quán)重系統(tǒng)一旦長時(shí)間懸著偏差罰分快速增長因此優(yōu)化出來的系統(tǒng)往往是“前期有一點(diǎn)超調(diào)、但能迅速穩(wěn)定”的風(fēng)格工程實(shí)際中非常好用。海星算法每產(chǎn)生一組Kp、Ki、Kd就把參數(shù)寫入工作區(qū)變量調(diào)用sim(‘DC_Motor_PID.slx’)仿真3秒計(jì)算ITAE值返回給優(yōu)化算法。三參數(shù)搜索范圍我設(shè)置為Kp∈[0, 20]Ki∈[0, 10]Kd∈[0, 5]這個(gè)范圍要保證覆蓋電機(jī)的物理可行參數(shù)同時(shí)不能太寬導(dǎo)致算法浪費(fèi)搜索時(shí)間。跑完50代迭代后海星算法找到的參數(shù)組合大約是Kp12.45Ki3.67Kd0.42。對(duì)照初始的人工試湊參數(shù)Kp8Ki2Kd0.1優(yōu)化后的系統(tǒng)在階躍響應(yīng)上表現(xiàn)為超調(diào)量從12.6%降到4.8%調(diào)節(jié)時(shí)間從1.8秒降到0.6秒穩(wěn)態(tài)誤差基本歸零。這就是優(yōu)化算法的價(jià)值——不是從0到1而是把原本“能用”的控制器調(diào)成“好用”的控制器。3.3 仿真結(jié)果與收斂性分析從海星算法的收斂曲線看前期前10代適應(yīng)度值下降非?？鞆某跏嫉?.42一路跌到0.09左右。第10到第20代進(jìn)入精細(xì)搜索階段適應(yīng)度緩慢下探到0.076附近。后面30代基本穩(wěn)定在0.074~0.075之間沒有出現(xiàn)明顯的“反跳”現(xiàn)象這說明算法在收斂過程中的穩(wěn)定性是好的。在抗負(fù)載擾動(dòng)測(cè)試中我在2秒時(shí)刻給電機(jī)突加50%額定負(fù)載轉(zhuǎn)矩再看轉(zhuǎn)速恢復(fù)時(shí)間。優(yōu)化后的系統(tǒng)在0.4秒內(nèi)恢復(fù)到穩(wěn)態(tài)轉(zhuǎn)速的99%而人工試湊參數(shù)的恢復(fù)時(shí)間是1.1秒。這個(gè)對(duì)比直接用波形圖展示效果非常直觀——超調(diào)尖峰明顯更小回穩(wěn)曲線更平滑。還有一個(gè)值得觀察的現(xiàn)象是優(yōu)化算法找出的Kp值比人工試湊大了不少12.45 vs 8這解釋起來也合理海星算法在搜索中發(fā)現(xiàn)“適當(dāng)增大比例增益、同時(shí)用微分項(xiàng)抑制超調(diào)”這個(gè)組合策略比人工保守取Kp更符合ITAE指標(biāo)的最優(yōu)方向。這其實(shí)也是智能整定的一個(gè)隱性收益——它有時(shí)候會(huì)突破工程師的經(jīng)驗(yàn)直覺給出反直覺但更優(yōu)的參數(shù)組合。4. 液位控制系統(tǒng)的建模與控制4.1 液位系統(tǒng)的機(jī)理建模與特性分析液位系統(tǒng)的經(jīng)典模型是單容水箱動(dòng)力學(xué)方程基于物料平衡A * dh/dt Q_in - Q_out其中A是水箱截面積h是液位高度Q_in是進(jìn)水流量Q_out是出水流量。如果出口閥門特性近似線性小擾動(dòng)假設(shè)下成立Q_out k * h模型就變成一階慣性環(huán)節(jié)G(s) H(s) / Q_in(s) K / (A*s k)我在模型中A取2m2k取0.8穩(wěn)態(tài)增益K就是1/k1.25時(shí)間常數(shù)TA/k2.5秒。這個(gè)參數(shù)下系統(tǒng)屬于典型的“好控制”對(duì)象為了增加項(xiàng)目難度和工程代表性我在模型前向通路中加入了一個(gè)純滯后環(huán)節(jié)滯后時(shí)間tau1.5秒模擬管道輸送和閥門動(dòng)作延遲。為什么要加這個(gè)延時(shí)因?yàn)榧儨笫枪I(yè)液位控制中最常見的難題來源。滯后時(shí)間與時(shí)間常數(shù)之比tau/T越大系統(tǒng)控制難度越高傳統(tǒng)PID的可調(diào)范圍就越窄。1.5秒的滯后相對(duì)2.5秒的時(shí)間常數(shù)比值達(dá)到了0.6已經(jīng)屬于PID控制中需要認(rèn)真對(duì)待的“中度困難”對(duì)象。搭建Simulink模型時(shí)我直接用傳遞函數(shù)模塊加Transport Delay模塊模型簡潔且參數(shù)化程度高方便海星算法反復(fù)修改控制參數(shù)進(jìn)行迭代。4.2 海星算法優(yōu)化控制器設(shè)計(jì)與參數(shù)選擇液位系統(tǒng)PID整定的優(yōu)化框架和電機(jī)是一樣的差異點(diǎn)在兩個(gè)地方一是適應(yīng)度函數(shù)需要重新定義二是PID參數(shù)搜索范圍要根據(jù)對(duì)象增益和時(shí)間尺度重新設(shè)定。先說適應(yīng)度函數(shù)。液位控制更關(guān)注不超調(diào)和無靜差單純用ITAE會(huì)讓算法為了降低積分時(shí)間而優(yōu)先追求“快”導(dǎo)致超調(diào)偏大。針對(duì)液位對(duì)象我在適應(yīng)度函數(shù)中加了一個(gè)“超調(diào)懲罰項(xiàng)”J ITAE w * max(0, overshoot - 10%)這里w取15意思是超調(diào)超過10%時(shí)每超出1個(gè)百分點(diǎn)加0.15的罰分。這是一個(gè)非常實(shí)用的工程技巧把控制目標(biāo)改寫成優(yōu)化目標(biāo)的加權(quán)組合比單純加l限幅更平滑而且能防止優(yōu)化算法“卡在約束邊界”上。參數(shù)搜索范圍方面因?yàn)橐何幌到y(tǒng)的穩(wěn)態(tài)增益比電機(jī)小時(shí)間尺度更大PID參數(shù)的量級(jí)完全不同。Kp∈[0, 5]Ki∈[0, 2]Kd∈[0, 3]。這個(gè)范圍同樣要做合理預(yù)估Kp太大系統(tǒng)會(huì)因純滯后而震蕩發(fā)散Kp太小響應(yīng)拖沓適應(yīng)度很差算法會(huì)在迭代過程中自動(dòng)淘汰這些不合理區(qū)域。最終海星算法找到的液位PID參數(shù)是Kp2.18Ki0.43Kd1.05。有趣的是它給了一個(gè)相對(duì)較大的微分增益原因正是前面提過的純滯后系統(tǒng)需要有“前瞻性”控制動(dòng)作來補(bǔ)償延遲微分項(xiàng)的預(yù)測(cè)能力在這里得到了充分的發(fā)揮。4.3 仿真結(jié)果對(duì)比階躍響應(yīng)與抗擾動(dòng)液位模型的階躍響應(yīng)仿真設(shè)定如下目標(biāo)液位從1米階躍到1.5米總仿真時(shí)長20秒。海星優(yōu)化參數(shù)下的響應(yīng)曲線上升過程平穩(wěn)約5.6秒到達(dá)目標(biāo)值附近超調(diào)量7.3%在預(yù)設(shè)10%的約束內(nèi)之后微幅振蕩兩次就穩(wěn)定下來穩(wěn)態(tài)誤差小于0.5%。人工試湊參數(shù)Kp1.2Ki0.15Kd0.3對(duì)比上升時(shí)間略慢7.2秒超調(diào)量9.8%調(diào)節(jié)時(shí)間反而更長——因?yàn)槌{(diào)后系統(tǒng)來回震蕩了三個(gè)周期才穩(wěn)定總調(diào)節(jié)時(shí)間8.4秒?？箶_動(dòng)測(cè)試我設(shè)定在10秒時(shí)刻突然增加20%的出水閥開度模擬外部干擾。海星參數(shù)下液位最大跌落了0.08米3秒內(nèi)恢復(fù)人工參數(shù)下最大跌落0.14米恢復(fù)時(shí)間超過6秒。從這些數(shù)據(jù)可以清楚看到在純滯后對(duì)象上算法優(yōu)化出的參數(shù)不是單純把某個(gè)指標(biāo)做到最優(yōu)而是在超調(diào)、響應(yīng)時(shí)間和抗擾恢復(fù)之間找到了一個(gè)平衡點(diǎn)。由優(yōu)化算法做出來的整定結(jié)果整體比“人肉”調(diào)出來的更符合ITAE加超調(diào)約束的復(fù)合目標(biāo)這就是設(shè)計(jì)加權(quán)適應(yīng)度函數(shù)的直接回報(bào)。5. 常見問題與經(jīng)驗(yàn)技巧實(shí)錄5.1 參數(shù)設(shè)置與調(diào)優(yōu)避坑指南海星算法本身的參數(shù)不多但每個(gè)參數(shù)對(duì)結(jié)果的影響我都實(shí)測(cè)過一遍這里直接給經(jīng)驗(yàn)值。種群規(guī)模建議取30~50不要太小也不要太大。太小會(huì)覆蓋面不足太大增加仿真計(jì)算量。本項(xiàng)目用40個(gè)個(gè)體每次迭代40次仿真50代就是2000次Simulink仿真一臺(tái)普通筆記本跑完大約需要15~20分鐘屬于可接受的耗時(shí)。最大迭代次數(shù)建議取50~100。更多迭代不會(huì)帶來明顯收益因?yàn)樵獑l(fā)式算法后期收斂曲線基本平了。另外如果觀察到最后20代適應(yīng)度完全沒有變化可以提前終止節(jié)省一半以上的時(shí)間。關(guān)鍵參數(shù)regeneration_rate取值0.3alpha取0.5beta取0.7。如果發(fā)現(xiàn)算法頻繁陷入局部最優(yōu)連續(xù)多次運(yùn)行結(jié)果差異大可以把regeneration_rate調(diào)大到0.5犧牲一點(diǎn)局部開發(fā)能力換取更強(qiáng)的全局探索。反之如果算法結(jié)果穩(wěn)定但精度不夠可以調(diào)小再生概率到0.15讓它更專注于局部打磨。一個(gè)常見的誤區(qū)是有人會(huì)把搜索范圍設(shè)得特別寬覺得“給算法更大空間它自己會(huì)找到好答案”。實(shí)際上搜索范圍過寬會(huì)讓初始解集離最優(yōu)區(qū)域太遠(yuǎn)浪費(fèi)大量迭代次數(shù)在無效區(qū)域飛行。正確的做法是先用經(jīng)驗(yàn)公式估算PID參數(shù)的大致量級(jí)再把搜索范圍設(shè)為這個(gè)量級(jí)的2~3倍這樣算法“發(fā)揮空間”和“收斂效率”是平衡的。5.2 仿真中容易翻車的細(xì)節(jié)這個(gè)項(xiàng)目我踩的坑不少挑幾個(gè)最有“價(jià)值”的分享第一Simulink仿真時(shí)間一定要大于系統(tǒng)調(diào)節(jié)時(shí)間。默認(rèn)的10秒仿真時(shí)間對(duì)大多數(shù)情況夠用但液位系統(tǒng)帶純滯后調(diào)節(jié)時(shí)間可能超過20秒。如果在仿真時(shí)間窗口內(nèi)系統(tǒng)還沒穩(wěn)定ITAE積分分母上的時(shí)間t還在增長優(yōu)化的結(jié)果就會(huì)偏向“在前半段快速反應(yīng)”的參數(shù)而忽略后半段的穩(wěn)定性。我最終給電機(jī)系統(tǒng)設(shè)3秒液位系統(tǒng)設(shè)25秒都要大于對(duì)象調(diào)節(jié)時(shí)間的2倍以上。第二S函數(shù)或模型里千萬不要有代數(shù)環(huán)。我的電機(jī)模型一開始用直接信號(hào)流搭建時(shí)轉(zhuǎn)速反饋和電流反饋之間產(chǎn)生了代數(shù)環(huán)Simulink每個(gè)步長內(nèi)需要迭代求解速度慢且容易發(fā)散。解決辦法是在電流環(huán)和轉(zhuǎn)速環(huán)之間插入一個(gè)Memory模塊或一階低通濾波人為打破代數(shù)環(huán)仿真穩(wěn)定性和速度都提升明顯。第三保存每次迭代的參數(shù)和適應(yīng)度不要只保存最終結(jié)果。我用一個(gè)結(jié)構(gòu)體數(shù)組記錄每一代每個(gè)個(gè)體的位置和適應(yīng)度值這樣后期分析收斂曲線、復(fù)現(xiàn)過程、排查異常都非常方便。很多人圖省事只存最優(yōu)結(jié)果結(jié)果想畫個(gè)“最優(yōu)適應(yīng)度隨代數(shù)變化”曲線都無從下手。第四隨機(jī)種子問題。元啟發(fā)式算法本質(zhì)是隨機(jī)搜索每次運(yùn)行結(jié)果會(huì)有波動(dòng)。做系統(tǒng)對(duì)比時(shí)一定要在Matlab開頭用rng(固定值)固定隨機(jī)種子否則對(duì)比不同算法時(shí)無法區(qū)分“算法優(yōu)劣”和“隨機(jī)波動(dòng)”。我做對(duì)比實(shí)驗(yàn)時(shí)固定rng(42)三組算法用同一個(gè)初始種群種子保證公平。5.3 代碼結(jié)構(gòu)建議與擴(kuò)展思路項(xiàng)目代碼結(jié)構(gòu)我建議按下面這樣組織邏輯清晰也好維護(hù)Proj_MFSO_PID ├── main.m % 主腳本調(diào)用海星算法 ├── starfish_optimizer.m % 海星算法主體函數(shù) ├── fitness_DC_Motor.m % 電機(jī)系統(tǒng)適應(yīng)度函數(shù) ├── fitness_Tank_Level.m % 液位系統(tǒng)適應(yīng)度函數(shù) ├── model_DC_Motor.slx % 電機(jī)Simulink模型 ├── model_Tank_Level.slx % 液位Simulink模型 └── results/ ├── log_parameters.mat % 迭代過程參數(shù)記錄 └── figures/ % 結(jié)果圖片輸出主腳本和算法函數(shù)分離兩個(gè)被控對(duì)象的適應(yīng)度函數(shù)分離這樣將來換對(duì)象模型時(shí)只需要新增一個(gè)適應(yīng)度函數(shù)算法部分完全復(fù)用。我甚至用同一套海星算法代碼去整定過一個(gè)溫度控制系統(tǒng)的PID只改了模型路徑和搜索范圍十分鐘內(nèi)就出結(jié)果。關(guān)于擴(kuò)展方向這個(gè)項(xiàng)目可以直接向三個(gè)方向延展。一是把單目標(biāo)ITAE擴(kuò)展成多目標(biāo)優(yōu)化超調(diào)最小化和調(diào)節(jié)時(shí)間最小化同時(shí)進(jìn)行海星算法可以改造成多目標(biāo)版本輸出Pareto前端讓使用者按需求選點(diǎn)。二是把固定PID參數(shù)改成增益調(diào)度——在不同液位區(qū)間或電機(jī)轉(zhuǎn)速區(qū)間用不同的優(yōu)化參數(shù)實(shí)現(xiàn)分段控制這在實(shí)際工程中非常常見。三是引入擾動(dòng)觀測(cè)器或前饋補(bǔ)償環(huán)節(jié)再把新增控制器的參數(shù)也納入海星算法的優(yōu)化變量這樣算法就從“優(yōu)化一個(gè)PID”升級(jí)成“優(yōu)化一個(gè)完整控制系統(tǒng)”。我個(gè)人在實(shí)際操作中最深的體會(huì)是海星優(yōu)化算法不是銀彈但在PID參數(shù)整定這個(gè)特定的場(chǎng)景里它的性價(jià)比確實(shí)很高。整個(gè)項(xiàng)目的代碼寫作量不大——核心算法函數(shù)大概150行兩個(gè)適應(yīng)度函數(shù)各50行主腳本30行——卻能換來比人工試湊穩(wěn)定得多的控制效果。最后再分享一個(gè)實(shí)用技巧在仿真收斂后把優(yōu)化出來的參數(shù)值輸出到Matlab的work space再用sim函數(shù)在同樣的擾動(dòng)條件下跑一遍完整仿真把響應(yīng)波形和優(yōu)化過程中的最好逐代響應(yīng)疊加到同一張圖里能明顯看出“第5代還行、第15代更好、第50代最佳”的演進(jìn)過程。這個(gè)動(dòng)圖或者疊加圖放在論文里或者PPT里說服力遠(yuǎn)遠(yuǎn)超過一張單調(diào)的收斂曲線。做控制優(yōu)化的有時(shí)候你不光要結(jié)果好還得要讓別人一眼看出好在哪效果做到位了項(xiàng)目就是成功的。