发明名称 基于优化迭代算法的时滞广域电力系统控制器
摘要 本发明属于电力系统广域时滞控制领域,具体涉及一种基于优化迭代算法的时滞广域电力系统控制器,在广域环境下控制信号在传输和处理过程中产生的时滞现象会对电力系统的稳定运行产生难以忽略的负面影响,本发明基于Lyapunov稳定性理论,研究在时滞影响下,基于状态反馈的广域电力系统控制器的设计问题,通过构造Lyapunov-Krasovskii泛函并对泛函进行求导,得到了非线性矩阵不等式稳定性判据;然后将不等式中的非线性项做线性化处理,使其转化为隶属于线性矩阵不等式的锥补问题,在迭代求解时,利用优化算法对迭代次数进行了优化,平衡了迭代时间与时滞上界的关系,最后通过仿真算例验证了所得控制器不仅具有较低的保守性,而且具有较快的响应速度。
申请公布号 CN104932256A 申请公布日期 2015.09.23
申请号 CN201510247402.6 申请日期 2015.05.15
申请人 河南理工大学 发明人 钱伟;高超;李赓;李冰锋;党建军
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 郑州科维专利代理有限公司 41102 代理人 亢志民;张欣棠
主权项 一种基于优化迭代算法的变时滞广域电力系统控制器,其特征在于:该控制器是基于以下步骤实现的:步骤1、建立广域电力系统变时滞状态空间模型:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mover><mi>x</mi><mo>.</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mi>A</mi><mo>+</mo><mi>&Delta;A</mi><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><msub><mi>A</mi><mn>1</mn></msub><mo>+</mo><msub><mi>&Delta;A</mi><mn>1</mn></msub><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>d</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><mi>Bu</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mfenced open='' close=''><mtable><mtr><mtd><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>&phi;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd><mtd><mi>t</mi><mo>&Element;</mo><mo>[</mo><mo>-</mo><mi>h</mi><mo>,</mo><mn>0</mn><mo>]</mo></mtd></mtr></mtable></mfenced></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000717579090000011.GIF" wi="935" he="143" /></maths>其中:x(t)为电力系统状态变量,x(t‑d(t))为经过时间延时后的状态变量;控制输入量u为附加的励磁输入ΔV<sub>s</sub>;φ(t)为[‑h,0]上连续的初始相量函数,ΔA、ΔA<sub>1</sub>为系统的扰动项目;[ΔA,ΔA<sub>1</sub>]=DF(t)[E,E<sub>1</sub>],其中F(t)满足条件:F<sup>T</sup>(t)F(t)≤I;步骤2、设置状态反馈控制器u(t)=Kx(t),其中K为状态反馈增益,以保证闭环系统<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>x</mi><mo>.</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mrow><mo>(</mo><mi>A</mi><mo>+</mo><mi>&Delta;A</mi><mo>+</mo><mi>BK</mi><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><msub><mi>A</mi><mn>1</mn></msub><mo>+</mo><msub><mi>&Delta;A</mi><mn>1</mn></msub><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>d</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000717579090000015.GIF" wi="852" he="87" /></maths>为渐近稳定;步骤3、给出时滞电力系统稳定性定理;定理1:针对任意的时滞d(t)满足0≤d(t)≤h,<img file="FDA0000717579090000016.GIF" wi="197" he="84" />如果存在矩阵L=L<sup>T</sup>>0,W=W<sup>T</sup>≥0,R=R<sup>T</sup>≥0,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>Y</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>Y</mi><mn>11</mn></msub></mtd><mtd><msub><mi>Y</mi><mn>12</mn></msub></mtd></mtr><mtr><mtd><mo>*</mo></mtd><mtd><msub><mi>Y</mi><mn>22</mn></msub></mtd></mtr></mtable></mfenced><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo></mrow>]]></math><img file="FDA0000717579090000012.GIF" wi="358" he="147" /></maths>以及任意合适维数矩阵M<sub>1</sub>,M<sub>2</sub>,V和一标量λ>0,使得下列矩阵不等式成立,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>&Pi;</mi><mn>11</mn></msub></mtd><mtd><msub><mi>&Pi;</mi><mn>12</mn></msub></mtd><mtd><mi>h</mi><mrow><mo>(</mo><msup><mi>LA</mi><mi>T</mi></msup><mo>+</mo><msup><mi>V</mi><mi>T</mi></msup><msup><mi>B</mi><mi>T</mi></msup><mo>+</mo><msup><mi>&lambda;DD</mi><mi>T</mi></msup><mo>)</mo></mrow></mtd><mtd><msup><mi>LE</mi><mi>T</mi></msup></mtd></mtr><mtr><mtd><mo>*</mo></mtd><mtd><msub><mi>&Pi;</mi><mn>22</mn></msub></mtd><mtd><msup><msub><mi>hLA</mi><mn>1</mn></msub><mi>T</mi></msup></mtd><mtd><msubsup><mi>LE</mi><mn>1</mn><mi>T</mi></msubsup></mtd></mtr><mtr><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>-</mo><mi>hR</mi><mo>+</mo><msup><mi>&lambda;h</mi><mn>2</mn></msup><msup><mi>DD</mi><mi>T</mi></msup></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>-</mo><mi>&lambda;I</mi></mtd></mtr></mtable></mfenced><mo>&lt;</mo><mn>0</mn><mo>,</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>Y</mi><mn>11</mn></msub></mtd><mtd><msub><mi>Y</mi><mn>12</mn></msub></mtd><mtd><msub><mi>M</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mo>*</mo></mtd><mtd><msub><mi>Y</mi><mn>22</mn></msub></mtd><mtd><msub><mi>M</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&mu;</mi><mo>)</mo></mrow><msup><mi>LR</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>L</mi></mtd></mtr></mtable></mfenced><mo>&GreaterEqual;</mo><mn>0</mn></mrow>]]></math><img file="FDA0000717579090000013.GIF" wi="1483" he="314" /></maths>式中;Π<sub>11</sub>=LA<sup>T</sup>+AL+BV+V<sup>T</sup>B<sup>T</sup>+M<sub>1</sub>+M<sub>1</sub><sup>T</sup>+W+hY<sub>11</sub>+λDD<sup>T</sup>Π<sub>12</sub>=A<sub>1</sub>L‑M<sub>1</sub>+M<sub>2</sub><sup>T</sup>+hY<sub>12</sub>Π<sub>22</sub>=‑M<sub>2</sub>‑M<sub>2</sub><sup>T</sup>‑(1‑μ)W+hY<sub>22</sub>则闭环系统<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mover><mi>x</mi><mo>.</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mrow><mo>(</mo><mi>A</mi><mo>+</mo><mi>&Delta;A</mi><mo>+</mo><mi>BK</mi><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><msub><mi>A</mi><mn>1</mn></msub><mo>+</mo><msub><mi>&Delta;A</mi><mn>1</mn></msub><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>d</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000717579090000017.GIF" wi="847" he="78" /></maths>是渐近稳定的;步骤4、定理1中的LR<sup>‑1</sup>L为非线性项,采用锥补线性化算对其进行线性化处理,并利用优化算法<img file="FDA0000717579090000014.GIF" wi="449" he="126" />对迭代次数进行优化;步骤5、根据步骤3、4的给出的算法,利用MATLAB软件求得保证广域电力系统的最大时滞上界h和状态反馈控制器K。
地址 454000 河南省焦作市高新区世纪路2001号河南理工大学