[发明专利]一种自抗扰控制器的鲁棒保性能设计与参数整定方法

摘要

本发明公开一种自抗扰控制器的鲁棒保性能设计与参数整定方法，包括步骤：1、建立不确定线性系统状态空间模型；2、建立全维扩张状态观测器，其分为：步骤21、建立系统的扩张状态空间模型和步骤22、建立全维扩张状态观测器；3、建立状态反馈控制结构；4、分析系统鲁棒稳定性和鲁棒性能上界；5控制系统参数寻优。本发明适用于最小相位系统与非最小相位系统的自抗扰控制器，以保证控制器系统在不确定性条件下，能够达到给定性能指标下的最优控制性能。

主权项

一种自抗扰控制的鲁棒保性能设计与参数整定方法，其特征在于：包括以下步骤：步骤1、建立不确定线性系统状态空间模型，表示成如下形式：xo(t)=(Ao+ΔA)xo+(Bo+ΔB)uy(t)=Coxo---(1)]]>其中是系统状态向量；Ao，Bo，Co是适当维数的常数矩阵，其描述了系统的名义模型，即忽略模型不确定后的系统模型；ΔA，ΔB是适当维数的结构化不确定矩阵函数；u(t)∈R是控制输入；y(t)∈R是系统输出；假定系统是能控能观的，且所考虑的不确定性是范数有界的，所述不确定矩阵函数具有以下形式:[ΔA(t) ΔB(t)]＝ME(t)[N0 N1]其中M、N0和N1是适当维数的已知矩阵，反映了不确定性的结构信息；E(t)是适当维数的未知矩阵，它可以是时变的，反映了系统模型的参数不确定性，且满足ET(t)E(t)≤I：步骤2、建立全维扩张状态观测器，其包括以下步骤：步骤21、建立系统的扩张状态空间模型：选取系统状态变量组xn+1＝f，其中f是由系统内部状态和输入信号构成等效扰动，并扩张成新状态，则式(1)可表示为以下状态方程组的形式：其中xE＝[x1 x2 … xn xn+1]T是系统状态向量，b0是系统增益，BE=0...010(n+1)×1T,FE=00...01(n+1)×1T,CE=10..01×(n+1);]]>步骤22、建立全维扩张状态观测器，表示成如下形式：z·=AEz+BEb0u+L(y-yz)=(AE-LCE)z+BEb0u+LyyZ=CEz---(2)]]>其中z＝[z1 … zn+1]T是观测器状态向量，L＝[l1 … ln+1]T是观测器增益向量，yz是观测器的输出，观测器增益使得观测器的状态z1,z2,…zn跟踪系统的状态变量x1,x2,…xn,xn+1，zn+1可观测系统的扩张状态，即系统的等效扰动，根据自抗扰控制理论，对扩张状态观测器进行极点配置s＝‑ωo，建立观测器增益L计算公式：|sI‑(AE‑LCE)|＝(s+ωo)n+1   (3)；步骤3、建立状态反馈控制结构在等效扰动估计和补偿作用下，系统近似补偿为串联积分形式，同时引入状态反馈，得到系统的状态反馈控制率：u=1b0(a~n(r-z1)-a~n-1z2-...-a~1zn-zn+1)=1b0(a~nr-Kz)---(4)]]>其中，在系统的等效扰动补偿和状态反馈作用下，闭环系统可近似达到理想的传递函数：Gcl(s)=a~nsn+a~1sn-1+...+a~n-1s+a~n=(ωcs+ωc)n]]>其中ωc可以是闭环系统的带宽参数，且[a~1a~2...a~n-1a~n]=[ωcα1ωc2α2...ωcn-1αn-1ωcnαn],αi=n!i!(n-i)!]]>结合(1)(2)(4)式，推导闭环系统的状态空间模型，由不确定性系统和全维扩张状态观测器构成的闭环系统的状态空间方程描述为：x·o=(Ao+ΔA)xo+(Bo+ΔB)1b0(a~nr-Kz)z·=LCoxo+(AE-LCE-BEK)z+a~nBEry=Coxoyz=CEz---(5)]]>在所述闭环系统中，K,L和b0为未知参数，其中，ωc根据系统响应速度要求进行设定，即K能通过(4)式确定；而L和b0则通过结合系统的鲁棒稳定性条件和控制性能进行设计；步骤4、分析系统鲁棒稳定性和鲁棒性能上界，其包括以下步骤：步骤41、确定系统的鲁棒稳定性条件定义另一个状态向量，令r(t)＝0，此时闭环系统为其中A‾=Ao+ME(t)N0-1b0(Bo+ME(t)N1)KLCoAE-LCE-BEK---(6)]]>闭环鲁棒稳定条件为：对于不确定线性系统(1)，如果存在对称正定矩阵和给定控制参数L和b0，使得对所有允许的不确定性，以下矩阵不等式(7)成立，则闭环系统是鲁棒渐进稳定的：Ψ11+N0TN0Ψ12-1b0N0TN1KP11M*Ψ22+1b02(N1CF)TN1KP12TM**-I<0---(7)]]>其中Ψ12=-1b0P11BoK+CoTLTP22+P12(AE-LCE-BEK),]]>Ψ22=P22(AE-LCE-BEK)+(AE-LCE-BEK)TP22-1b0P12TBoK-1b0KTBoTP12.]]>步骤42：确定系统的鲁棒性能上界对不确定系统(1)，定义以下二次型性能指标：J=∫0∞[xoT(t)Qxo(t)+uT(t)Ru(t)]dt---(8)]]>其中Q和R是给定的正定加权矩阵，给定合适的L和b0，如果存在对称正定矩阵P，使得对所有允许的不确定性，下列矩阵不等式成立：Ψ11+Q+N0TN0Ψ12-1b0N0TN1KP11M*Ψ22+1b02CFT(R+N1TN1)KP12TM**-I<0---(9)]]>则状态反馈控制率u(t)是系统(1)的一个保性能控制律，相应的系统性能指标满足不失一般性，假定初始状态是一个满足的随机变量，结合性能指标的期望值，得到即得到性能指标上界；步骤5、控制系统参数(ωo,b0)寻优根据步骤1‑4，闭环系统中ωo，b0是由满足(9)式约束条件的优化问题确定，描述为：minωo,b0{minPTrace(P)}]]>s.t.(9)      (11)所述优化问题为两层嵌套寻优，其中外层参数寻优采用优化函数实现(ωo,b0)的迭代；内层参数寻优则针对每个迭代点(ωo,b0)，采用线性矩阵不等式工具箱求解(9)式中的最优P矩阵，得到性能指标最优的ωo，b0参数，再通过(3)式得到观测器增益参数L，最终得到具有最优性能指标的鲁棒自抗扰控制器。