发明名称 一种用于动态稳定控制的多平衡点非线性系统的设计方法
摘要 本发明提供的一种用于动态稳定控制的多平衡点非线性系统的设计方法,包括:1)建立同步机组模型,并确定类梯度的定义;2)判断同步机组模型的类梯度性;若同步机组模型是混沌的,3)对步骤1)同步机组模型进行控制器设计;4)将步骤3)的控制器进行电力系统动态稳定分析与控制,从而保证整个电力系统的动态稳定性。本发明所利用的动态稳定控制方法没有涉及到之前类似方法中的线性化,因此大大降低了该方法中产生的累积误差。并且对以后的实际工程具有指导作用。
申请公布号 CN102799111B 申请公布日期 2016.08.03
申请号 CN201210167908.2 申请日期 2012.05.25
申请人 中国电力科学研究院;北京理工大学 发明人 徐式蕴;汤涌;孙华东
分类号 G05B17/00(2006.01)I 主分类号 G05B17/00(2006.01)I
代理机构 北京安博达知识产权代理有限公司 11271 代理人 徐国文
主权项 一种用于动态稳定控制的多平衡点非线性系统的设计方法,其特征在于,所述方法包括如下步骤:1)建立同步机组模型,并确定类梯度的定义;通过KYP引理将类梯度性的频域条件转化为等价的线性矩阵不等式;所述建立同步机组模型:<img file="FDA0000961236380000011.GIF" wi="1316" he="351" />其中,线性部分的传递函数为m×m维的矩阵K(s):K(s)=C<sup>T</sup>(A‑sI)<sup>‑1</sup>B将其简化为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mover><mi>&sigma;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>&eta;</mi></mrow></mtd></mtr><mtr><mtd><mrow><mover><mi>&eta;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mrow><mo>(</mo><mo>-</mo><msub><mi>&alpha;</mi><mn>1</mn></msub><msub><mi>y</mi><mn>1</mn></msub><mo>-</mo><msub><mi>&alpha;</mi><mn>2</mn></msub><msub><mi>y</mi><mn>2</mn></msub><mo>)</mo></mrow><mi>sin</mi><mi>&sigma;</mi><mo>+</mo><mrow><mo>(</mo><msub><mi>&alpha;</mi><mn>3</mn></msub><msub><mi>y</mi><mn>3</mn></msub><mi>cos</mi><mi>&sigma;</mi><mo>-</mo><msub><mi>&alpha;</mi><mn>4</mn></msub><mi>sin</mi><mi>&sigma;</mi><mi>cos</mi><mi>&sigma;</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mover><msub><mi>y</mi><mn>1</mn></msub><mo>&CenterDot;</mo></mover><mo>=</mo><msub><mi>&alpha;</mi><mn>5</mn></msub><mo>-</mo><msub><mi>&alpha;</mi><mn>6</mn></msub><msub><mi>y</mi><mn>1</mn></msub><mo>+</mo><msub><mi>&alpha;</mi><mn>7</mn></msub><msub><mi>y</mi><mn>2</mn></msub><mo>+</mo><msub><mi>&alpha;</mi><mn>8</mn></msub><mi>cos</mi><mi>&sigma;</mi></mrow></mtd></mtr><mtr><mtd><mrow><mover><msub><mi>y</mi><mn>2</mn></msub><mo>&CenterDot;</mo></mover><mo>=</mo><msub><mi>&alpha;</mi><mn>9</mn></msub><msub><mi>y</mi><mn>1</mn></msub><mo>-</mo><msub><mi>&alpha;</mi><mn>10</mn></msub><msub><mi>y</mi><mn>2</mn></msub><mo>+</mo><msub><mi>&alpha;</mi><mn>11</mn></msub><mi>cos</mi><mi>&sigma;</mi></mrow></mtd></mtr><mtr><mtd><mrow><mover><msub><mi>y</mi><mn>3</mn></msub><mo>&CenterDot;</mo></mover><mo>=</mo><mo>-</mo><msub><mi>&alpha;</mi><mn>12</mn></msub><msub><mi>y</mi><mn>3</mn></msub><mo>+</mo><msub><mi>&alpha;</mi><mn>13</mn></msub><mi>sin</mi><mi>&sigma;</mi></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000961236380000012.GIF" wi="1717" he="486" /></maths>其中<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mo>-</mo><msub><mi>&alpha;</mi><mn>6</mn></msub></mrow></mtd><mtd><msub><mi>&alpha;</mi><mn>7</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>&alpha;</mi><mn>9</mn></msub></mtd><mtd><mrow><mo>-</mo><msub><mi>&alpha;</mi><mn>10</mn></msub></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><msub><mi>&alpha;</mi><mn>12</mn></msub></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mi>B</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>&alpha;</mi><mn>8</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>&alpha;</mi><mn>11</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><msub><mi>&alpha;</mi><mn>13</mn></msub></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000961236380000013.GIF" wi="750" he="200" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>C</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>&alpha;</mi><mn>1</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>&alpha;</mi><mn>2</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><msub><mi>&alpha;</mi><mn>3</mn></msub></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mi>a</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>&alpha;</mi><mn>5</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo><mi>f</mi><mrow><mo>(</mo><mi>&sigma;</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mi>c</mi><mi>o</mi><mi>s</mi><mi>&sigma;</mi></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mi>s</mi><mi>i</mi><mi>n</mi><mi>&sigma;</mi></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000961236380000014.GIF" wi="766" he="199" /></maths><img file="FDA0000961236380000015.GIF" wi="588" he="55" />σ、η、y都对应于同步机组中的物理参数,是随时间变化的向量,其中向量y=[y<sub>1</sub>,y<sub>2</sub>,y<sub>3</sub>]<sup>T</sup>;A、B、C是同步机组模型中的已知参数矩阵,是常数矩阵;f(σ)、<img file="FDA0000961236380000016.GIF" wi="182" he="63" />g(η,σ)是同步机组模型中关于σ的周期非线性函数;α为一般化的同步机组模型中的已知参数;α<sub>1</sub>‑α<sub>13</sub>为方程系数;2)判断同步机组模型的类梯度性;若同步机组模型具有类梯度性,则电力动态系统稳定;所述判断同步机组模型的类梯度性包括:假设μ<sub>1</sub>≥0,存在矩阵P=P<sup>T</sup>使得引理2中的假设条件[1]、[4]、[5]以及下列时域线性矩阵不等式条件<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msup><mi>A</mi><mi>T</mi></msup><mi>P</mi><mo>+</mo><mi>P</mi><mi>A</mi></mrow></mtd><mtd><mrow><mi>P</mi><mi>B</mi><mo>-</mo><msup><mi>A</mi><mi>T</mi></msup><mi>C</mi></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>B</mi><mi>T</mi></msup><mi>P</mi><mo>-</mo><msup><mi>C</mi><mi>T</mi></msup><mi>A</mi></mrow></mtd><mtd><mrow><mo>-</mo><mrow><mo>(</mo><msup><mi>C</mi><mi>T</mi></msup><mi>B</mi><mo>+</mo><msup><mi>B</mi><mi>T</mi></msup><mi>C</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>&lt;</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000961236380000021.GIF" wi="1310" he="143" /></maths>满足,则模型(1)为类梯度的;其中,μ<sub>1</sub>为矩阵参数;引理2:令μ<sub>1</sub>≥0,且下列条件满足:[1]矩阵K(0)是对角矩阵;[2]对于所有ω∈R,都有<img file="FDA0000961236380000022.GIF" wi="359" he="123" />[3]<img file="FDA0000961236380000023.GIF" wi="574" he="119" />[4]函数f(σ)具有连续二阶导数,且在任何区间(θ<sub>1</sub>,θ<sub>2</sub>)内都有f'(σ)≠0;[5]<img file="FDA0000961236380000024.GIF" wi="318" he="90" />那么非线性系统的模型(1)是类梯度的;若同步机组模型是混沌的,则进行步骤3);3)对步骤1)同步机组模型进行控制器设计;4)将步骤3)的控制器进行电力系统动态稳定分析与控制。
地址 100192 北京市海淀区清河小营东路15号