发明名称 飞行器多个时间滞后模型逼近及控制器设计方法
摘要 本发明公开了一种飞行器多个时间滞后模型逼近及控制器设计方法,用于解决现有的鲁棒控制理论缺乏设计步骤难以直接设计飞行控制器的技术问题。技术方案是给出时变系统分段鲁棒稳定可解条件,直接利用线性系统状态反馈的闭环期望极点选择,并根据所有闭环期望极点的实部全部为负数的特点,给出了限定条件不等式直接设计反馈矩阵。使得本研究领域的工程技术人员对风洞或飞行试验得到的含有多个时间滞后不确定性飞行器模型直接设计飞行控制器,解决了当前研究只给出鲁棒稳定性不等式而无法直接设计飞行控制器的技术问题。
申请公布号 CN102929136A 申请公布日期 2013.02.13
申请号 CN201210380995.X 申请日期 2012.10.10
申请人 西北工业大学 发明人 史忠科
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 1.一种飞行器多个时间滞后模型逼近及控制器设计方法,其特征在于包括以下步骤:步骤一、在给定高度、马赫数条件下通过风洞或飞行试验得到含有r个时间滞后不确定性的飞行器模型为:<maths num="0001"><![CDATA[<math><mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mrow><mo>(</mo><msub><mi>A</mi><mn>0</mn></msub><mo>+</mo><mi>&Delta;</mi><msub><mi>A</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>x</mi><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>r</mi></munderover><msub><mi>A</mi><mi>&tau;i</mi></msub><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>&tau;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>+</mo><mi>Bu</mi></mrow></math>]]></maths>式中,x∈R<sup>n</sup>,u∈R<sup>m</sup>分别为状态和输入向量,A<sub>0</sub>,A<sub>τi</sub>(i=1,2,…,r),B为已知常系数矩阵,τ<sub>i</sub>(i=1,2,…,r)为未知延迟时间,ΔA<sub>0</sub>为系数矩阵未知部分;全文符号相同;<img file="FDA00002234789800012.GIF" wi="1760" he="275" />迭代求解N,式中:x<sub>real</sub>(t)为实际系统状态响应;将模型表达成:<maths num="0002"><![CDATA[<math><mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mrow><mo>(</mo><mi>A</mi><mo>+</mo><mi>&Delta;A</mi><mo>)</mo></mrow><mi>x</mi><mo>+</mo><mi>Bu</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,x∈R<sup>n</sup>,u∈R<sup>m</sup>分别为状态和输入向量,A=(A<sub>0</sub>+A<sub>τ</sub>N),B为已知的系数矩阵,ΔA为系数矩阵未知部分;选择飞行控制器为:u=-Kx式中,K为反馈矩阵;带入(1)式中,有:<maths num="0003"><![CDATA[<math><mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mrow><mo>(</mo><mi>A</mi><mo>-</mo><mi>BK</mi><mo>+</mo><mi>&Delta;A</mi><mo>)</mo></mrow><mi>x</mi></mrow></math>]]></maths>步骤二、选取(A-B K)的特征值各不相同且实部为负,设计反馈矩阵K使得满足条件:Λ&gt;M <sup>T</sup>(ΔA)<sup>T</sup> M<sup>-T</sup> M<sup>-1</sup>ΔA M;该控制器使得<maths num="0004"><![CDATA[<math><mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mrow><mo>(</mo><mi>A</mi><mo>-</mo><mi>BK</mi><mo>+</mo><mi>&Delta;A</mi><mo>)</mo></mrow><mi>x</mi></mrow></math>]]></maths>鲁棒稳定;式中,M为线性变换矩阵,M<sup>-1</sup>(A-B K)M=diag[σ<sub>1</sub>+jω<sub>1</sub>,σ<sub>2</sub>+jω<sub>2</sub>,…,σ<sub>n</sub>+jω<sub>n</sub>],σ<sub>i</sub>,ω<sub>i</sub>(i=1,2,…,n)为实数,jω<sub>i</sub>(i=1,2,…,n)表示虚数,diag为对角矩阵符号,<maths num="0005"><![CDATA[<math><mrow><mi>&Lambda;</mi><mo>=</mo><mi>diag</mi><mo>[</mo><msubsup><mi>&sigma;</mi><mn>1</mn><mn>2</mn></msubsup><mo>,</mo><msubsup><mi>&sigma;</mi><mn>2</mn><mn>2</mn></msubsup><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup><mo>]</mo><mo>;</mo></mrow></math>]]></maths>ΔA-ΔBK通常假设为ΔA-ΔBK=HFW,H,W均假设为已经矩阵,0&lt;F≤I,I=diag[1,1,…,1]为单位阵。
地址 710072 陕西省西安市友谊西路127号