[發(fā)明專利]一種基于自適應(yīng)和聲搜索算法的PI參數(shù)優(yōu)化方法

在審

• 申請?zhí)枺?/td>: 201910147429.6
• 申請日: 2019-02-27
• 公開（公告）號: CN110554599A
• 公開（公告）日: 2019-12-10
• 發(fā)明（設(shè)計）人: 李永麗;趙學(xué)明;孫廣宇;張云柯;曾亮
• 申請（專利權(quán)）人: 天津大學(xué);國網(wǎng)上海市電力公司
• 主分類號: G05B11/42
• 分類號: G05B11/42
• 代理公司: 12201 天津市北洋有限責(zé)任專利代理事務(wù)所
• 代理人: 程毓英

1.一種基于自適應(yīng)和聲搜索算法的PI參數(shù)優(yōu)化方法，設(shè)HMS為和聲庫大小；HM代表和聲記憶庫，是一個HMS*3的矩陣；HMCR代表記憶庫取值概率；PAR代表微調(diào)概率；BW代表微調(diào)帶寬，T_max代表最大迭代次數(shù)。包含以下步驟：

(1)確立變量X₁和X₂與目標(biāo)函數(shù)J_ITAE，其中，X₁代表定電流控制器的比例系數(shù)K_p，X₂代表積分時間常數(shù)K_i；

(2)在變量X₁和X₂的可行域內(nèi)，隨機生成HMS組X₁和X₂，分別放入矩陣HM中的第一列和第二列；

(3)把每一組X₁和X₂分別代入PI控制器進行仿真，記錄系統(tǒng)故障后到恢復(fù)額定狀態(tài)所需時間t_s以及送端母線暫態(tài)電壓的最大值U_1max，送端母線暫態(tài)電壓的最大值U_1max與穩(wěn)態(tài)測量值U₁的差即為e(t)；

(4)計算每一組X₁和X₂的目標(biāo)函數(shù)值J_ITAE，并將目標(biāo)函數(shù)值放入矩陣HM中的第三列中；

(5)X₁’和X₂’，隨機生成一個0-1之間的隨機數(shù)r，若r小于記憶庫取值概率HMCR，則從已經(jīng)生成的HMS組變量中選擇一組新的變量X₁’和X₂’，否則，隨機生成一組新的變量X₁’和X₂’；當(dāng)新的變量X₁’和X₂’從HMS組變量中選擇時，進行下一步的判斷：新生成另一個0-1之間的隨機數(shù)r₁，若隨機數(shù)r₁小于微調(diào)概率PAR，則對新的變量X₁’和X₂’進行微調(diào)操作，即，x₁′＝x₁′+r₁*BW，x₂′＝x₂′+r₁*BW；

(6)把新生成的變量X₁’和X₂’代入PI控制器進行仿真，記錄系統(tǒng)故障后到恢復(fù)額定狀態(tài)所需時間t_s’以及送端母線暫態(tài)電壓的最大值U_1max’，送端母線暫態(tài)電壓的最大值U_1max’與穩(wěn)態(tài)測量值U₁的差即為e(t)；

(7)計算新變量X₁’和X₂’所對應(yīng)的目標(biāo)函數(shù)值J_ITAE’，并與矩陣HM中的J_ITAEi進行比較，若J_ITAE’小于J_ITAEi，則將第i組變量X₁和X₂替換成新的變量X₁’和X₂’；

(8)檢查當(dāng)前迭代次數(shù)gn是否達到最大迭代次數(shù)T_max，若gn<T_max，重復(fù)步驟(5)—(8)；

(9)當(dāng)gn等于最大迭代次數(shù)T_max時，令X₁為當(dāng)前矩陣HM第一列變量的平均值，X₂為當(dāng)前矩陣HM第二列變量的平均值，輸出X₁和X₂的值。