发明名称 内含微小电感或电容的复杂电路系统建模与控制方法
摘要 一种内含微小电感或电容的复杂电路系统建模与控制方法,复杂电路系统控制技术领域。该方法基于UDTFSPM,结合谱范数与线性矩阵不等式方法,为被控CCSWSIOCs设计RFSOFC,实现CCSWSIOCs的高精度控制。根据CCSWSIOCs的动力学模型,建立其不确定性连续时间模糊奇异摄动模型,选择适当的采样周期,采用零阶保持器方法,对所获连续时间模糊奇异摄动模型进行离散化,获得CCSWSIOCs的UDTFSPM,在此基础上设计RFSOFC。优点在于,解决现有建模与控制方法无法消除CCSWSIOCs的内部微小电感或电容引起的不稳定或稳态误差大的问题,大幅提高CCSWSIOCs的控制性能。将本发明应用于Van der Pol电路系统的高精度控制,进行仿真验证,表明提出方法的有效性。
申请公布号 CN103645646A 申请公布日期 2014.03.19
申请号 CN201310646527.7 申请日期 2013.12.04
申请人 冶金自动化研究设计院 发明人 陈金香
分类号 G05B17/02(2006.01)I 主分类号 G05B17/02(2006.01)I
代理机构 北京华谊知识产权代理有限公司 11207 代理人 刘月娥
主权项 1.一种内含微小电感或电容的复杂电路系统建模与控制方法,其特征在于:步骤一、根据CCSWSIOCs的动力学方程,建立被控CCSWSIOCs的不确定性连续时间模糊奇异摄动模型;将CCSWSIOCs的变化缓慢或能够直接测量的状态变量看作为慢变量,小参数相关或变化较快的状态变量看作为快变量,采用扇区非线性方法,建立CCSWSIOCs的不确定性连续时间模糊奇异摄动模型:规则i:如果ξ<sub>1</sub>(t)是φ<sub>i1</sub>,...,ξ<sub>g</sub>(t)是φ<sub>ig</sub>,那么<![CDATA[<math><mrow><msub><mi>E</mi><mi>&epsiv;</mi></msub><mover><mi>x</mi><mo>.</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msub><mi>A</mi><mi>ci</mi></msub><mo>+</mo><mi>&Delta;</mi><msub><mi>A</mi><mi>ci</mi></msub><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>B</mi><mi>ci</mi></msub><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>Dw</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></math>]]></maths>y(t)=Hx(t)                              (1)其中,<![CDATA[<math><mrow><msub><mi>E</mi><mi>&epsiv;</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>I</mi><mrow><mi>n</mi><mo>&times;</mo><mi>n</mi></mrow></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>&epsiv;</mi><msub><mi>I</mi><mrow><mi>m</mi><mo>&times;</mo><mi>m</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>x</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mi>H</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>I</mi><mrow><mi>n</mi><mo>&times;</mo><mi>n</mi></mrow></msub></mtd><mtd><msub><mn>0</mn><mrow><mi>n</mi><mo>&times;</mo><mi>m</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>x<sub>s</sub>(t)∈R<sup>n</sup>为慢变量,x<sub>f</sub>(t)∈R<sup>m</sup>为快变量,u(t)∈R<sup>q</sup>为控制输入,y(t)∈R<sup>l</sup>为系统输出,w(t)∈R<sup>q</sup>为外扰,φ<sub>i1</sub>,...,φ<sub>ig</sub>(i=1,2,...,r)均为模糊集合,ξ<sub>1</sub>(t),..,ξ<sub>g</sub>(t)为可测量的系统变量,A<sub>ci</sub>,B<sub>ci</sub>,D为合适维数矩阵,ΔA<sub>ci</sub>为合适维数的不确定性矩阵,ε是奇异摄动参数;步骤二、建立被控CCSWSIOCs的UDTFSPM;控制系统中的传感器和执行器均采用时间驱动方式,且二者采用相同的采样时间T<sub>s</sub>,在零阶保持器的作用下,将以上连续时间模型(1),离散化为UDTFSPM:规则i:如果ξ<sub>1</sub>(k)是φ<sub>i1</sub>,...,ξ<sub>g</sub>(k)是φ<sub>ig</sub>,那么<![CDATA[<math><mrow><mover><mi>x</mi><mo>.</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>E</mi><mi>&epsiv;</mi></msub><mrow><mo>(</mo><msub><mi>A</mi><mi>i</mi></msub><mo>+</mo><mi>&Delta;</mi><msub><mi>A</mi><mi>i</mi></msub><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>E</mi><mi>&epsiv;</mi></msub><msub><mi>B</mi><mi>i</mi></msub><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>E</mi><mi>&epsiv;</mi></msub><mi>Dw</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></math>]]></maths>y(k)=Hx(k)                  (2)其中,ΔA<sub>i</sub>为适当维数的不确定性矩阵,<![CDATA[<math><mrow><msub><mi>A</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>E</mi><mi>&epsiv;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msup><mi>e</mi><mrow><msubsup><mi>E</mi><mi>&epsiv;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>A</mi><mi>ci</mi></msub><msub><mi>T</mi><mi>s</mi></msub></mrow></msup><mo>,</mo><msub><mi>B</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>E</mi><mi>&epsiv;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msubsup><mo>&Integral;</mo><mn>0</mn><mi>h</mi></msubsup><msubsup><mi>E</mi><mi>&epsiv;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msup><mi>e</mi><mrow><msubsup><mi>E</mi><mi>&epsiv;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>A</mi><mi>ci</mi></msub><mi>&tau;</mi></mrow></msup><mi>d&tau;</mi><msub><mi>B</mi><mi>ci</mi></msub></mrow></math>]]></maths>给定[x(k);u(k);w(k)],应用标准模糊推理方法,得到全局UDTFSPM:<![CDATA[<math><mrow><mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>E</mi><mi>&epsiv;</mi></msub></mrow><mrow><mo>(</mo><mi>A</mi><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>+</mo><mi>&Delta;A</mi><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>E</mi><mi>&epsiv;</mi></msub><mi>B</mi><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>E</mi><mi>&epsiv;</mi></msub><mi>Dw</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></math>]]></maths>y(k)=Hx(k)                    (3)其中,隶属度函数<![CDATA[<math><mrow><msub><mi>&mu;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&xi;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&xi;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>r</mi></munderover><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&xi;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&xi;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Pi;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>g</mi></munderover><msub><mi>&phi;</mi><mi>ij</mi></msub><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo><mo>,</mo></mrow></math>]]></maths>φ<sub>ij</sub>(ξ<sub>j</sub>(k))为ξ<sub>j</sub>(k)在φ<sub>ij</sub>中的隶属度,设w<sub>i</sub>(ξ(k))≥0,for i=1,2,…,r,r为规则数,μ<sub>i</sub>(ξ(k))≥0,<img file="FDA0000430077030000017.GIF" wi="218" he="71" />为了便于记录我们令μ<sub>i</sub>=μ<sub>i</sub>(ξ(k)),<![CDATA[<math><mrow><mi>A</mi><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>r</mi></munderover><msub><mi>&mu;</mi><mi>i</mi></msub><msub><mi>A</mi><mi>i</mi></msub><mo>,</mo><mi>B</mi><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>r</mi></munderover><msub><mi>&mu;</mi><mi>i</mi></msub><msub><mi>B</mi><mi>i</mi></msub><mo>,</mo><mi>&Delta;A</mi><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>r</mi></munderover><msub><mi>&mu;</mi><mi>i</mi></msub><mi>&Delta;</mi><msub><mi>A</mi><mi>i</mi></msub></mrow></math>]]></maths>步骤三、基于UDTFSPM(3),对被控对象设计RFSOFC设计如下RFSOFC,其模糊规则前件与系统(3)的模糊规则前件相同,u(k)=G(μ)y(k)                 (4)其中,<img file="FDA0000430077030000021.GIF" wi="259" he="90" />G<sub>i</sub>为控制器增益;步骤四、建立闭环系统模型;针对被控系统模型(3),应用控制率(4),获得闭环系统模型:x(k+1)=E<sub>ε</sub>(A(μ)+B(μ)G(μ)H+ΔA(μ))x(k)+E<sub>ε</sub>Dw(k)              (5)步骤五、采用谱范数方法和线性矩阵不等式方法,推导出RFSOFC存在的充分条件,推导过程不要求知道系统不确定性参数的上确界。下面是求解RFSOFC增益的线性矩阵不等式组:Θii&lt;0          (i=1,2,...r)                   (6)<![CDATA[<math><mrow><msup><mi>&epsiv;</mi><mn>2</mn></msup><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Omega;</mi></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>+</mo><mi>&epsiv;</mi><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Xi;</mi></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><msub><mi>&Theta;</mi><mi>ii</mi></msub><mo>&lt;</mo><mn>0</mn><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mi>r</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>Θ<sub>ij</sub>+Θ<sub>ji</sub>&lt;0            (1≤i&lt;j≤r)                   (8)<![CDATA[<math><mrow><msup><mi>&epsiv;</mi><mn>2</mn></msup><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Omega;</mi></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>+</mo><mi>&epsiv;</mi><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Xi;</mi></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><msub><mi>&Theta;</mi><mi>ij</mi></msub><mo>+</mo><msub><mi>&Theta;</mi><mi>ji</mi></msub><mo>&lt;</mo><mn>0</mn><mrow><mo>(</mo><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&lt;</mo><mi>j</mi><mo>&le;</mo><mi>r</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,<![CDATA[<math><mrow><msub><mi>&Gamma;</mi><mi>ij</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Lambda;</mi><mo>-</mo><mi>&beta;I</mi></mtd><mtd><mo>*</mo></mtd></mtr><mtr><mtd><msub><mi>A</mi><mi>i</mi></msub><mi>Y</mi><mo>+</mo><msub><mi>B</mi><mi>i</mi></msub><msub><mi>W</mi><mi>j</mi></msub></mtd><mtd><mo>-</mo><mi>Z</mi></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>W<sub>j</sub>=[W<sub>ij</sub>  0<sub>q×m</sub>],<![CDATA[<math><mrow><mi>&Lambda;</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>Z</mi><mn>11</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mi>&gamma;Y</mi><mo>-</mo><msup><mi>&gamma;Y</mi><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths><![CDATA[<math><mrow><mi>&Omega;</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>Z</mi><mn>22</mn></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mi>&Xi;</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mo>*</mo></mtd></mtr><mtr><mtd><msub><mi>Z</mi><mn>21</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>γ(0&lt;γ≤1),β为大于零的常数,γ,β的值可由设计者选择(设计者通过选择合适的γ,β值得到最优控制器增益),<![CDATA[<math><mrow><mi>Z</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>Z</mi><mn>11</mn></msub></mtd><mtd><msubsup><mi>Z</mi><mn>21</mn><mi>T</mi></msubsup></mtd></mtr><mtr><mtd><msub><mi>Z</mi><mn>21</mn></msub></mtd><mtd><msub><mi>Z</mi><mn>22</mn></msub></mtd></mtr></mtable></mfenced><mrow><mo>(</mo><msub><mi>Z</mi><mn>11</mn></msub><mo>&Element;</mo><msup><mi>R</mi><mrow><mi>n</mi><mo>&times;</mo><mi>n</mi></mrow></msup><mo>,</mo><msub><mi>Z</mi><mn>22</mn></msub><mo>&Element;</mo><msup><mi>R</mi><mrow><mi>m</mi><mo>&times;</mo><mi>m</mi></mrow></msup><mo>)</mo></mrow></mrow></math>]]></maths>为对称正定矩阵,<![CDATA[<math><mrow><mi>Y</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>Y</mi><mn>11</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>Y</mi><mn>22</mn></msub></mtd></mtr></mtable></mfenced><mrow><mo>(</mo><msub><mi>Y</mi><mn>11</mn></msub><mo>&Element;</mo><msup><mi>R</mi><mrow><mi>n</mi><mo>&times;</mo><mi>n</mi></mrow></msup><mo>,</mo><msub><mi>Y</mi><mn>22</mn></msub><mo>&Element;</mo><msup><mi>R</mi><mrow><mi>m</mi><mo>&times;</mo><mi>m</mi></mrow></msup><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>控制器增益:<img file="FDA0000430077030000029.GIF" wi="209" he="42" />for   i=1,2,…,r.                   (10)步骤六、将所得控制器Matlab代码传化为C语言代码,植入控制器。控制器采用事件驱动方式,当采样数据到达控制器时,控制器立刻进行计算,并将控制信号传给执行器,执行器按照固定的采样周期读取控制信号,生成控制输入,作用于被控CCSWSIOCs,从而实现CCSWSIOCs的高精度控制。
地址 100071 北京市丰台区西四环南路72号