主权项 |
一种基于量化的广域电力系统输出反馈控制方法,其特征在于:所述控制方法包括以下步骤:1)在广域测量系统当中,信息的采集是通过PMU来完成的,PMU以固定的采样周期对电力系统状态量进行测量,以T为采样周期的闭环电力系统离散化状态空间模型为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = '{' close = ''><mtable><mtr><mtd><mi>x</mi><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo><mo>=</mo><mi>A</mi><mi>x</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>+</mo><mi>B</mi><mi>u</mi><mo>(</mo><mi>k</mi><mo>)</mo></mtd></mtr><mtr><mtd><mi>y</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>=</mo><mi>C</mi><mi>x</mi><mo>(</mo><mi>k</mi><mo>)</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000747948040000011.GIF" wi="529" he="157" /></maths>其中,k表示第k个传输周期kT;A,B表示电机的常数矩阵;C为系统的输出矩阵;x(k+1)表示广域电力系统在第k+1周期的状态变量;u(k)为系统第k周期的输入变量;y(k)表示系统在第k周期的输出变量;2)首先根据输出信号对系统的状态进行观测,然后设计基于观测器的输出反馈控制器对系统进行控制,并在电力系统闭环系统中加入量化器,对PMU的输出量进行量化,采用对数量化器:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open = '' close = ''><mtable><mtr><mtd><mrow><msub><mi>y</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>y</mi><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><msub><mi>u</mi><mi>i</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><mi>δ</mi></mrow></mfrac><msub><mi>u</mi><mi>i</mi></msub><mo><</mo><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>≤</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>-</mo><mi>δ</mi></mrow></mfrac><msub><mi>u</mi><mi>i</mi></msub><mo>;</mo><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>></mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mi>f</mi><mrow><mo>(</mo><mo>-</mo><mi>y</mi><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo><mo>,</mo></mrow></mtd><mtd><mrow><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo><</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000747948040000012.GIF" wi="1124" he="317" /></maths>相应的量化等级定义为:U={±u<sub>(i)</sub>,u<sub>(i)</sub>=ρ<sup>i</sup>u<sub>(0)</sub>,i=±1,±2,…}∪{±u<sub>(0)</sub>}∪{0}0<ρ<1,u<sub>(0)</sub>>0其中,y<sub>q</sub>(k)为量化后的输出变量,u<sub>(i)</sub>为第i级量化水平、u<sub>(0)</sub>为初始量化水平,ρ为量化密度,<img file="FDA0000747948040000021.GIF" wi="229" he="140" />带有量化器的观测器结构为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mi>A</mi><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>B</mi><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>L</mi><mo>(</mo><msub><mi>y</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>y</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow>]]></math><img file="FDA0000747948040000022.GIF" wi="871" he="94" /></maths>其中,<img file="FDA0000747948040000023.GIF" wi="168" he="71" />为k+1时刻系统状态的估计值、L为观测器增益、<img file="FDA0000747948040000024.GIF" wi="111" he="71" />为观测器的输出变量;3)采用Lyapunov稳定性理论及线性矩阵不等式方法可得到保证整个广域电力系统在量化误差存在的情况下保持稳态运行的充分条件:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open = '[' close = ']'><mtable><mtr><mtd><mrow><mo>-</mo><mi>P</mi><mo>+</mo><msup><mi>ϵδ</mi><mn>2</mn></msup><msup><mi>C</mi><mi>T</mi></msup><mi>C</mi></mrow></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mi>Q</mi></mrow></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd></mtr><mtr><mtd><mrow><mi>P</mi><mi>A</mi><mo>+</mo><mi>Y</mi></mrow></mtd><mtd><mrow><mo>-</mo><mi>Y</mi></mrow></mtd><mtd><mrow><mo>-</mo><mi>P</mi></mrow></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mi>Q</mi><mi>A</mi><mo>+</mo><mi>Z</mi><mi>C</mi></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mi>Q</mi></mrow></mtd><mtd><mo>*</mo></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msup><mi>Z</mi><mi>T</mi></msup></mtd><mtd><mrow><mo>-</mo><mi>ϵ</mi><mi>I</mi></mrow></mtd></mtr></mtable></mfenced><mo><</mo><mn>0</mn></mrow>]]></math><img file="FDA0000747948040000025.GIF" wi="1009" he="387" /></maths>其中,P为Lyapunov函数的状态权重矩阵,Q为Lyapunov函数的误差权重矩阵,ε为1维参数变量,矩阵Y为Y=PBK,矩阵Z为Z=QL,通过Matlab中的LMI工具箱可解出矩阵Y和Z,则所求的控制器增益为K=(PB)<sup>‑1</sup>Y,再将得到的控制器加入系统中,使系统保持稳定。 |