发明名称 含有非线性的飞行器控制器设计方法
摘要 本发明公开了一种含有非线性的飞行器控制器设计方法,用于解决现有的鲁棒控制理论缺乏设计步骤难以直接设计飞行控制器的技术问题。技术方案是给出含有不确定非线性项系统的鲁棒稳定可解条件,直接利用线性系统状态反馈的闭环期望极点选择,并根据所有闭环期望极点的实部全部为负数的特点,给出了限定条件不等式直接设计反馈矩阵。使得本研究领域的工程技术人员可以对风洞或飞行试验得到的含有不确定非线性项的飞行器模型直接设计飞行控制器,解决了当前研究只给出鲁棒稳定性不等式而无法直接设计飞行控制器的技术问题。
申请公布号 CN102929138A 申请公布日期 2013.02.13
申请号 CN201210381072.6 申请日期 2012.10.10
申请人 西北工业大学 发明人 史忠科
分类号 G05B13/04(2006.01)I;G05D1/00(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种含有非线性的飞行器控制器设计方法,其特征在于包括以下步骤:步骤一、在给定高度、马赫数条件下通过风洞或飞行试验得到含有不确定非线性项的飞行器模型为: <mrow> <mover> <mi>x</mi> <mo>&CenterDot;</mo> </mover> <mo>=</mo> <msub> <mi>A</mi> <mn>0</mn> </msub> <mi>x</mi> <mo>+</mo> <msub> <mi>A</mi> <mi>non</mi> </msub> <msub> <mi>f</mi> <mi>&Delta;</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>+</mo> <mo>[</mo> <mi>B</mi> <mo>+</mo> <mi>&Delta;B</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>]</mo> <mi>u</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>式中,x∈Rn,u∈Rm分别为状态和输入向量,A0、B、Anon为已知的系数矩阵,fΔ(x)、ΔB(x)为未知非线性项;按照不同的飞行区域和飞行机动大小,将(1)式右端的状态项表达成: <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>A</mi> <mn>0</mn> </msub> <mi>x</mi> <mo>+</mo> <msub> <mi>A</mi> <mi>non</mi> </msub> <msub> <mi>f</mi> <mi>&Delta;</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>A</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>+</mo> <mi>&Delta;</mi> <msub> <mi>A</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <mi>B</mi> <mo>+</mo> <mi>&Delta;B</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>B</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>+</mo> <mi>&Delta;</mi> <msub> <mi>B</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mi>min</mi> </mrow> </msub> <mo>&le;</mo> <mi>x</mi> <mo>&lt;</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mi>max</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>=</mo> <mn>1,2</mn> <mo>,</mo> <mo>.</mo> <mo>.</mo> <mo>.</mo> <mo>,</mo> <mi>r</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow>式中,A0i、B0i为已知的常数矩阵,ΔA0i ΔB0i为未知矩阵,ximin、ximax分别为第i个飞行区域所对应的状态最小和最大值,r为正整数,i为下标;在飞行区域ximin≤x<ximax内,飞行控制器为:u(t)=‑Kix(t)式中,Ki为常数反馈矩阵;带入(1)式中,有: <mrow> <mover> <mi>x</mi> <mo>&CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>[</mo> <mrow> <mo>(</mo> <msub> <mi>A</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>B</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <msub> <mi>K</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mi>&Delta;</mi> <msub> <mi>A</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>-</mo> <mi>&Delta;</mi> <msub> <mi>B</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <msub> <mi>K</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>]</mo> <mi>x</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow>步骤二、选取(A0i‑B0iK0i)的特征值各不相同且实部为负,设计反馈矩阵Ki使得满足条件: <mrow> <msub> <mi>&Lambda;</mi> <mi>i</mi> </msub> <mo>></mo> <msubsup> <mi>M</mi> <mi>i</mi> <mi>T</mi> </msubsup> <msup> <mrow> <mo>(</mo> <mi>&Delta;</mi> <msub> <mi>A</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>-</mo> <mi>&Delta;</mi> <msub> <mi>B</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <msub> <mi>K</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mi>T</mi> </msup> <msubsup> <mi>M</mi> <mi>i</mi> <mrow> <mo>-</mo> <mi>T</mi> </mrow> </msubsup> <msubsup> <mi>M</mi> <mi>i</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <mi>&Delta;</mi> <msub> <mi>A</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>-</mo> <mi>&Delta;</mi> <msub> <mi>B</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <msub> <mi>K</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>M</mi> <mi>i</mi> </msub> <mo>;</mo> </mrow>该控制器使得 <mrow> <mover> <mi>x</mi> <mo>&CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>[</mo> <mrow> <mo>(</mo> <msub> <mi>A</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>B</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <msub> <mi>K</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mi>&Delta;</mi> <msub> <mi>A</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>-</mo> <mi>&Delta;</mi> <msub> <mi>B</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <msub> <mi>K</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>]</mo> <mi>x</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow>鲁棒稳定;式中,Mi为线性变换矩阵, <mrow> <msubsup> <mi>M</mi> <mi>i</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>A</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>B</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <msub> <mi>K</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>M</mi> <mi>i</mi> </msub> <mo>=</mo> <mi>diag</mi> <mo>[</mo> <msub> <mi>&sigma;</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> <mo>+</mo> <mi>j</mi> <msub> <mi>&omega;</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>&sigma;</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> </msub> <mo>+</mo> <mi>j</mi> <msub> <mi>&omega;</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> </msub> <mo>,</mo> <mo>.</mo> <mo>.</mo> <mo>.</mo> <mo>,</mo> <msub> <mi>&sigma;</mi> <mi>in</mi> </msub> <mo>+</mo> <mi>j</mi> <msub> <mi>&omega;</mi> <mi>in</mi> </msub> <mo>]</mo> <mo>,</mo> </mrow>σik,ωik(k=1,2,…,n)为实数,jωik(k=1,2,…,n)表示虚数,diag为对角矩阵符号, <mrow> <msub> <mi>&Lambda;</mi> <mi>i</mi> </msub> <mo>=</mo> <mi>diag</mi> <mo>[</mo> <msubsup> <mi>&sigma;</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> <mn>2</mn> </msubsup> <mo>,</mo> <msubsup> <mi>&sigma;</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> <mn>2</mn> </msubsup> <mo>,</mo> <mo>.</mo> <mo>.</mo> <mo>.</mo> <mo>,</mo> <msubsup> <mi>&sigma;</mi> <mi>in</mi> <mn>2</mn> </msubsup> <mo>]</mo> <mo>;</mo> </mrow>ΔA0i‑ΔB0iKi通常假设为ΔA0i‑ΔB0iKi=HiFiWi,Hi,Wi均假设为已经矩阵,0<Fi≤I,I=diag[1,1,…,1]为单位阵。
地址 710072 陕西省西安市友谊西路127号