发明名称 基于未知目标轨迹的非线性系统鲁棒自适应跟踪控制方法
摘要 本发明公开了一种基于未知目标轨迹的非线性系统鲁棒自适应跟踪控制方法,包括步骤:步骤一、建立含有执行器故障的多输入多输出非线性系统的数学模型;步骤二、建立估计未知目标理想轨迹的模型,并利用此模型得出估计的目标轨迹逼近理想的目标轨迹;步骤三、设计鲁棒自适应容错控制器。本发明能够得出估计的未知目标轨迹,使其逼近理想的目标轨迹;设计的控制器通过巧妙地对“虚拟控制增益”进行矩阵分解,使得其他控制方法中对增益矩阵的条件限制大大降低;最终,在系统同时存在执行器故障、参数不确定性、外部干扰及追踪目标轨迹未知的情况下仍然可以获得渐近跟踪的稳态性能,使闭环控制系统对不确定模型和未知干扰具有鲁棒自适应作用。
申请公布号 CN106154839A 申请公布日期 2016.11.23
申请号 CN201610831242.4 申请日期 2016.09.19
申请人 重庆大学 发明人 宋永端;张东
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 重庆信航知识产权代理有限公司 50218 代理人 吴彬
主权项 基于未知目标轨迹的非线性系统鲁棒自适应跟踪控制方法,其特征在于:包括以下步骤:步骤一、建立含有执行器故障的多输入多输出非线性系统的数学模型;所述含执行器故障的多输入多输出非线性系统具有如下状态空间形式:<maths num="0001"><math><![CDATA[<mrow><msup><mi>y</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msup><mo>=</mo><mi>F</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>+</mo><mi>G</mi><mrow><mo>(</mo><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><msub><mi>u</mi><mi>a</mi></msub><mo>+</mo><mi>d</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001116547970000011.GIF" wi="843" he="73" /></maths>式中x=[x<sub>1</sub><sup>T</sup>,…,x<sub>q</sub><sup>T</sup>]<sup>T</sup>∈R<sup>n</sup>是系统的整个状态向量,其中<img file="FDA0001116547970000012.GIF" wi="778" he="99" />i=1,2,…,q,并且n<sub>1</sub>+n<sub>2</sub>+…+n<sub>q</sub>=n;<img file="FDA0001116547970000013.GIF" wi="781" he="79" />是系统的输出;F(x)∈R<sup>q</sup>是未知的函数向量;<img file="FDA0001116547970000014.GIF" wi="262" he="66" />是连续可微的未知函数,其中<img file="FDA0001116547970000015.GIF" wi="645" he="101" /><img file="FDA0001116547970000016.GIF" wi="819" he="105" />d(x,t)∈R<sup>q</sup>为系统的不确定性非线性,u<sub>a</sub>∈R<sup>q</sup>为控制输入向量;考虑执行器故障,实际控制输入u<sub>a</sub>与理想控制输入u的关系为:u<sub>a</sub>=ρ(t)u+ε(t),式中ρ=diag{ρ<sub>i</sub>}∈R<sup>q×q</sup>是对角矩阵,ρ<sub>i</sub>为执行器效率因子,并满足0<<u>h</u><sub>i</sub>≤ρ<sub>i</sub>≤1,<u>h</u><sub>i</sub>为ρ<sub>i</sub>的最小值;ε(t)代表控制行为中完全失控的部分并假设为有界;步骤二、建立估计未知目标理想轨迹的模型,并利用此模型得出估计的目标轨迹逼近理想的目标轨迹;对于未知的目标轨迹,采用基于拓展卡尔曼滤波的数学模型去估计,并使其逼近理想的目标轨迹;y<sub>d</sub>(t)=y<sub>EKF</sub>(t)+y<sup>*</sup><sub>guess</sub>(t)<maths num="0002"><math><![CDATA[<mrow><msup><msub><mi>y</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><msup><mi>y</mi><mi>j</mi></msup><mrow><mi>E</mi><mi>K</mi><mi>F</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><msup><mi>y</mi><mrow><mo>*</mo><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></mrow></msup><mrow><mi>g</mi><mi>u</mi><mi>e</mi><mi>s</mi><mi>s</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001116547970000017.GIF" wi="862" he="115" /></maths>式中,y<sub>d</sub>(t)∈R<sup>q</sup>为未知目标理想轨迹的估计值;y<sub>EKF</sub>(t)∈R<sup>q</sup>为用拓展卡尔曼滤波技术预测的理想目标轨迹;y<sup>*</sup><sub>guess</sub>(t)∈R<sup>q</sup>为基于某种已知条件得出的理想轨迹的粗略估计值,如果没有可用的已知条件,此值可以为0;y<sub>d</sub><sup>(j)</sup>(t)是理想轨迹的j阶导数<img file="FDA0001116547970000021.GIF" wi="187" he="90" />的估计值,y<sup>j</sup><sub>EKF</sub>(t)和<img file="FDA0001116547970000022.GIF" wi="299" he="109" />分别为用拓展卡尔曼滤波技术及已知条件预测的对应理想j阶导数轨迹;通过拓展卡尔曼滤波得出的最优估计值对应于模型中的y<sub>EKF</sub>(t)以及y<sup>j</sup><sub>EKF</sub>(t);步骤三、设计鲁棒自适应容错控制器;1)利用估计的未知目标轨迹与系统输出得到跟踪误差E<sub>m</sub>,通过滑模滤波器后得到新的状态变量s<sub>m</sub>;2)控制增益G和执行器效率因子ρ作为整体视作虚拟控制增益,对虚拟控制增益进行矩阵分解,得到已知的矩阵D(x)、U(x)和未知的矩阵S(x);其中矩阵S(x)作为不确定项通过核心函数产生器,同时系统不确定模型和外部干扰不确定项也通过核心函数发生器处理,得到未知的虚拟参数a和可计算的核心函数<img file="FDA0001116547970000023.GIF" wi="124" he="66" />3)核心函数<img file="FDA0001116547970000024.GIF" wi="113" he="65" />的平方通过任意正比例c<sub>1</sub>放大后与之前得到的状态变量s<sub>m</sub>的乘积组成未知的虚拟参数的一部分,再减去估计的虚拟参数<img file="FDA0001116547970000025.GIF" wi="40" he="62" />的c<sub>2</sub>倍,c<sub>2</sub>为任意正常数,得到虚拟参数<img file="FDA0001116547970000026.GIF" wi="43" he="55" />的导数值,最后进行积分运算得到未知的虚拟参数a的估计值<img file="FDA0001116547970000027.GIF" wi="63" he="63" />4)利用得到的状态变量s<sub>m</sub>与核心函数<img file="FDA0001116547970000028.GIF" wi="112" he="69" />的积再乘以未知的虚拟参数a的估计值的‑c<sub>1</sub>倍,得到最后的控制器u;5)控制器u将计算出的控制指令发送给非线性系统的执行器,实现系统输出跟踪理想的目标轨迹。
地址 400030 重庆市沙坪坝区正街174号