主权项 |
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>ϵ</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>Δ</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>ϵ</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>I</mi><mrow><mi>n</mi><mo>×</mo><mi>n</mi></mrow></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>ϵ</mi><msub><mi>I</mi><mrow><mi>m</mi><mo>×</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>×</mo><mi>n</mi></mrow></msub></mtd><mtd><msub><mn>0</mn><mrow><mi>n</mi><mo>×</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>ϵ</mi></msub><mrow><mo>(</mo><msub><mi>A</mi><mi>i</mi></msub><mo>+</mo><mi>Δ</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>ϵ</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>ϵ</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>ϵ</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msup><mi>e</mi><mrow><msubsup><mi>E</mi><mi>ϵ</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>ϵ</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msubsup><mo>∫</mo><mn>0</mn><mi>h</mi></msubsup><msubsup><mi>E</mi><mi>ϵ</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msup><mi>e</mi><mrow><msubsup><mi>E</mi><mi>ϵ</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>A</mi><mi>ci</mi></msub><mi>τ</mi></mrow></msup><mi>dτ</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>·</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>E</mi><mi>ϵ</mi></msub></mrow><mrow><mo>(</mo><mi>A</mi><mrow><mo>(</mo><mi>μ</mi><mo>)</mo></mrow><mo>+</mo><mi>ΔA</mi><mrow><mo>(</mo><mi>μ</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>ϵ</mi></msub><mi>B</mi><mrow><mo>(</mo><mi>μ</mi><mo>)</mo></mrow><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>E</mi><mi>ϵ</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>μ</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>ξ</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>ξ</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow><mrow><munderover><mi>Σ</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>ξ</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>ξ</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><munderover><mi>Π</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>g</mi></munderover><msub><mi>φ</mi><mi>ij</mi></msub><mrow><mo>(</mo><msub><mi>ξ</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>μ</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>r</mi></munderover><msub><mi>μ</mi><mi>i</mi></msub><msub><mi>A</mi><mi>i</mi></msub><mo>,</mo><mi>B</mi><mrow><mo>(</mo><mi>μ</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>r</mi></munderover><msub><mi>μ</mi><mi>i</mi></msub><msub><mi>B</mi><mi>i</mi></msub><mo>,</mo><mi>ΔA</mi><mrow><mo>(</mo><mi>μ</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>r</mi></munderover><msub><mi>μ</mi><mi>i</mi></msub><mi>Δ</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<0 (i=1,2,...r) (6)<![CDATA[<math><mrow><msup><mi>ϵ</mi><mn>2</mn></msup><mfenced open='[' close=']'><mtable><mtr><mtd><mi>Ω</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>ϵ</mi><mfenced open='[' close=']'><mtable><mtr><mtd><mi>Ξ</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>Θ</mi><mi>ii</mi></msub><mo><</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><0 (1≤i<j≤r) (8)<![CDATA[<math><mrow><msup><mi>ϵ</mi><mn>2</mn></msup><mfenced open='[' close=']'><mtable><mtr><mtd><mi>Ω</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>ϵ</mi><mfenced open='[' close=']'><mtable><mtr><mtd><mi>Ξ</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>Θ</mi><mi>ij</mi></msub><mo>+</mo><msub><mi>Θ</mi><mi>ji</mi></msub><mo><</mo><mn>0</mn><mrow><mo>(</mo><mn>1</mn><mo>≤</mo><mi>i</mi><mo><</mo><mi>j</mi><mo>≤</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>Γ</mi><mi>ij</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>Λ</mi><mo>-</mo><mi>β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>Λ</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>γY</mi><mo>-</mo><msup><mi>γY</mi><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths><![CDATA[<math><mrow><mi>Ω</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>Ξ</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<γ≤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>∈</mo><msup><mi>R</mi><mrow><mi>n</mi><mo>×</mo><mi>n</mi></mrow></msup><mo>,</mo><msub><mi>Z</mi><mn>22</mn></msub><mo>∈</mo><msup><mi>R</mi><mrow><mi>m</mi><mo>×</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>∈</mo><msup><mi>R</mi><mrow><mi>n</mi><mo>×</mo><mi>n</mi></mrow></msup><mo>,</mo><msub><mi>Y</mi><mn>22</mn></msub><mo>∈</mo><msup><mi>R</mi><mrow><mi>m</mi><mo>×</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的高精度控制。 |