[发明专利]一种基于非光滑控制理论的励磁系统调节器计算方法

摘要

本发明公开一种基于非光滑控制理论的励磁系统调节器计算方法，该方法以同步发电机功率角作为“平坦输出”，基于非光滑控制理论，给出同步发电机自动电压调节器的参数化形式，提高了同步发电机电压调节计算效率；利用和声搜索优化算法，优化自动电压调节器的参数，提高了同步发电机自动电压调节器设计的自动化程度，避免了采用传统PID控制方法设计自动电压调节器时需要反复手动调整PID参数、效率低下的问题。

主权项

1.一种基于非光滑控制理论的励磁系统调节器计算方法，其特征是：该方法包括如下步骤：1)在计算机中编写程序，建立同步发电机系统模型为Pa(δ,E′q)＝Pe(δ,E′q)‑Pm  (1‑d)Ef＝kcu  (1‑g)其中，δ是发电机的功率角，D是阻尼系数，ω是发电机角速度，ωs是同步电动机角速度，v＝ω‑ωs是发电机的相对速度，Pm是机械输入功率，Pe是发电机产生的电功率，Pa＝Pe‑Pm是加速功率，H是单位惯性系数，E′q是瞬时感应电动势，Eq是感应电动势，Ef是励磁绕组等效感应电动势，T′do是瞬时开环电路时间常数，Vs是电网的电压，kc是励磁放大器的增益，u是发电机励磁放大器的输入，xT是变压器的阻抗，xd是d轴的阻抗，x′d是d轴的瞬时阻抗，xL是传输线的阻抗；2)选取δ作为同步发电机系统模型(1‑a)～(1‑g)的“平坦输出变量”F，所谓某个系统的“平坦输出变量”是指该系统的状态变量和控制输入可以表达为“平坦输出变量”及其各阶导数的非线性函数；将同步发电机系统模型的状态变量和控制输入表示为x1＝F  (2‑a)其中，3)记将同步发电机系统模型表示为如下的线性系统模型4)针对线性系统模型(3)，设计如下形式的控制器其中，函数sigα(β)定义为ki>0(i＝1,2,3)可使多项式s3+k3s2+k2s+k1＝0胡尔维茨(Hurwitz)稳定；且有5)将控制器式(4)代入同步电压系统模型式(2～a)到(2～d)，得到同步发电机系统的自动电压调节器公式u的表达式如下其中，k1，k2，k3为待优化的自动电压调节器参数；6)选取自动电压调节器参数优化的性能指标函数f(x)，其中x为待优化参数向量，也即x＝[k1，k2，k3]T；7)采用和声搜索算法，进行自动电压调节器式(5)的参数自动寻优，该算法通过类比音乐和最优化问题的相似性，将求解优化问题的过程，模拟为乐师创作音乐的过程，乐师优化算法通过记忆反复调整乐队中各乐器的音调，更新和声记忆库，即待优化变量的解集，最终得到美妙的和声即最优解，为应用该算法，首先需要选定算法的主要参数，包括和声库规模HMS、和声库取值概率HMCR、音调调节概率PAR和音调调节带宽bw、迭代次数M；8)根据待优化的变量k1、k2和k3的取值范围，记其为xi,令N表示待寻优变量的个数，初始化和声记忆库，采用随机数生成算法，得到HMS组和声，并记为(x1，x2，…，xHMS)，将相应的和声记忆库表示为9)令迭代计数器j＝1，计算和声记忆库中所有和声对应的目标函数值，记为f(x1),f(x2),…,f(xN)；10)采用以下的方式生成一个新的和声：(a)在和声记忆库中选取后，按照一定的概率进行音调微调；(b)在和声库以外，通过随机生成算法在取值范围内，进行随机选择；首先，采用随机数生成算法产生一个随机概率r₁，如果r₁<HMCR，则从和声记忆库式(6)中的选择一个作为x₁′，为采用随机数算法产生的随机概率,对其进行变量微调，微调的幅度为rand×bw,调整算法为其中，PAR是微调概率，bw是微调带宽，rand是微调概率系数，rand1为采用随机数算法产生的随机概率；如果r1>HMCR，则采用随机生成算法，在x1允许的范围内随机产生一个作为x1′，其它待寻优变量x1～xN的处理方式与x1相同，不再赘述；通过这一步骤，得到产生的新和声x′＝[x′1,x′2,...,x′N]T；11)针对x′计算目标函数值，得到f(x′)，记fmin为f(x1),f(x2),…,f(xN)中的最小值，与其对应的和声为xi，如果f(x′)>fmin，则用x′替换xi，得到新的和声记忆库；12)采用计算机程序，判断和声记忆库中各个和声所对应的目标函数值差异是否满足误差范围，若满足误差范围则转步骤14)，否则转入下一步，判断目标函数值差异的函数为13)采用计算机程序，判断j是否大于最大迭代次数M，如果大于，则转入下一步，否则，令迭代次数j＝j+1，转入步骤10)；14)在计算机程序中，停止进行优化计算，输出和声记忆库中的目标函数最大的和声作为参数寻优问题的最优解决方案。