发明名称 用于兆瓦级逆变系统的极值优化自整定数字PID控制方法
摘要 本发明公开一种用于兆瓦级逆变系统的极值优化自整定数字PID控制方法。该方法为将兆瓦级逆变系统的模型输入、系统误差信号和可调参数初始值作为数字PID控制器输入;将PID控制器控制参数<i>K<sub>p</sub></i>、<i>K<sub>i</sub></i>和<i>K<sub>d</sub></i>编码成二进制序列,再通过二进制编码极值优化求解器对二进制序列进行迭代优化,直到满足工程设计人员预先设定的终止条件,从而自动整定出最优的二进制序列,通过二进制解码转化为PID控制器最优控制参数,并传输至兆瓦级逆变系统实现实时在线的优化控制。该发明的有益效果为:相比基于经验规则的工程整定方法,本发明可以实时在线整定和优化PID控制器参数;相比遗传算法等其它智能优化方法,本发明需设计的可调参数更少,实施更简单,且优化整定效率更高。
申请公布号 CN103365213B 申请公布日期 2015.10.21
申请号 CN201310300513.X 申请日期 2013.07.15
申请人 温州大学 发明人 曾国强;陈杰;戴瑜兴;张正江;郑崇伟;陆康迪;彭文文;彭子舜
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 周烽
主权项 一种用于兆瓦级逆变系统的极值优化自整定数字PID控制方法,其特征在于,该方法为通过数据采集和初始化设置模块将兆瓦级逆变系统的模型输入、系统误差信号和可调参数初始值作为数字PID控制器输入;通过二进制编码操作将PID控制器控制参数K<sub>p</sub>、K<sub>i</sub>和K<sub>d</sub>编码成二进制序列,再通过二进制编码极值优化求解器对二进制序列进行迭代优化,直到满足工程设计人员预先设定的终止条件,从而自动整定出最优的二进制序列,通过二进制解码转化为PID控制器最优控制参数K<sub>p</sub>、K<sub>i</sub>和K<sub>d</sub>,并将上述最优控制参数K<sub>p</sub>、K<sub>i</sub>和K<sub>d</sub>传输至兆瓦级逆变系统实现实时在线的优化控制;该方法具体包括以下步骤:步骤1:数据采集与初始化设置:采集兆瓦级逆变系统模型输入、控制系统误差信号,并对二进制编码极值优化方法可调参数进行初始化设置;步骤2:二进制编码操作:将PID控制器3个控制参数编码成一串长度为L的二进制序列,其中<img file="FDA0000783183730000011.GIF" wi="217" he="128" />l<sub>i</sub>为第i个控制参数的二进制编码长度,是PID控制器优化参数精度的重要影响因素之一;步骤3:初始解生成模块:随机产生初始二进制序列解S,按照式(1)计算初始适应度C(S),并设置当前最好解S<sub>best</sub>=S和对应的适应度C<sub>best</sub>=C(S<sub>best</sub>);<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>C</mi><mrow><mo>(</mo><mi>S</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>w</mi><mn>1</mn></msub><mo>|</mo><mi>e</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>|</mo><mo>+</mo><msub><mi>w</mi><mn>2</mn></msub><msup><mi>u</mi><mn>2</mn></msup><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>+</mo><msub><mi>w</mi><mn>3</mn></msub><msub><mi>t</mi><mi>u</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi></mrow></mtd><mtd><mrow><mi>e</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>w</mi><mn>1</mn></msub><mo>|</mo><mi>e</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>|</mo><mo>+</mo><msub><mi>w</mi><mn>2</mn></msub><msup><mi>u</mi><mn>2</mn></msup><mo>(</mo><mi>k</mi><mo>)</mo><mo>+</mo><msub><mi>w</mi><mn>4</mn></msub><mo>|</mo><mi>e</mi><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><mo>)</mo><mo>+</mo><msub><mi>w</mi><mn>3</mn></msub><msub><mi>t</mi><mi>u</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi></mrow></mtd><mtd><mrow><mi>e</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000783183730000012.GIF" wi="1673" he="260" /></maths>其中,T<sub>s</sub>为采样周期,n为采样次数且n=INT(t/T<sub>s</sub>),e(k)为第k个采样时刻的误差信号,u(k)为第k个采样时刻的PID控制器输出信号,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>K</mi><mi>p</mi></msub><mi>e</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>K</mi><mi>i</mi></msub><msub><mi>T</mi><mi>s</mi></msub><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mi>k</mi></munderover><mi>e</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><mfrac><msub><mi>K</mi><mi>d</mi></msub><msub><mi>T</mi><mi>s</mi></msub></mfrac><mo>&lsqb;</mo><mi>e</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mi>e</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>&rsqb;</mo><mo>,</mo></mrow>]]></math><img file="FDA0000783183730000013.GIF" wi="1162" he="174" /></maths>t<sub>u</sub>为系统在单位阶跃信号作用下输出响应曲线的上升时间,w<sub>1</sub>,w<sub>2</sub>,w<sub>3</sub>,w<sub>4</sub>均为权重系数,且w<sub>4</sub>&gt;&gt;w<sub>1</sub>;步骤4:随机翻转当前解S中第i个二进制位的数值,保持其它二进制位不变,从而生成解S<sub>i</sub>,并计算相应的适应度C(S<sub>i</sub>);步骤5:局部适应度计算与排序操作:计算第i个二进制位的局部适应度λ<sub>i</sub>=C(S<sub>i</sub>)-C<sub>best</sub>,并对λ<sub>i</sub>数值按照从大到小的顺序进行排序;步骤6:选择操作:按照幂律概率分布P(k<sub>1</sub>)∝k<sub>1</sub><sup>‑τ</sup>,选择λ<sub>i</sub>的排序号,其中1≤k<sub>1</sub>≤L,其中τ为幂律分布可调参数,并将对应的二进制位标记为x<sub>j</sub>;步骤7:变异操作和新解生成机制:翻转二进制位x<sub>j</sub>的数值,保持当前解S中其它位数值不变,从而生成新解S<sub>new</sub>;步骤8:比较C<sub>best</sub>和C(S<sub>new</sub>)的大小:若C<sub>best</sub>≤C(S<sub>new</sub>),则保持S<sub>best</sub>和C<sub>best</sub>不变;否则,S<sub>best</sub>=S<sub>new</sub>,C<sub>best</sub>=C(S<sub>new</sub>);步骤9:无条件地接受S=S<sub>new</sub>;步骤10:判断是否满足终止条件,若是,则转向下一步骤;否则,则返回到步骤4;步骤11:数据转化与输出模块:输出最好解S<sub>best</sub>和对应的适应度C<sub>best</sub>,将二进制序列S<sub>best</sub>解码转化为对应的PID控制参数K<sub>p</sub>、K<sub>i</sub>和K<sub>d</sub>,并将所述PID控制参数K<sub>p</sub>、K<sub>i</sub>和K<sub>d</sub>传输至兆瓦级逆变系统实现实时在线控制。
地址 325035 浙江省温州市茶山高教园区