发明名称 含有分布参数项的飞行器控制器设计方法
摘要 本发明公开了一种含有分布参数项的飞行器控制器设计方法,用于解决现有的鲁棒控制理论缺乏设计步骤难以直接设计飞行控制器的技术问题。技术方案是给出含有分布参数项系统的鲁棒稳定可解条件,直接利用线性系统状态反馈的闭环期望极点选择,并根据所有闭环期望极点的实部全部为负数的特点,给出了限定条件不等式直接设计反馈矩阵。使得本研究领域的工程技术人员对风洞或飞行试验得到含有分布参数项的飞行器模型直接设计飞行控制器,解决了当前研究只给出鲁棒稳定性不等式而无法直接设计飞行控制器的技术问题。
申请公布号 CN102929137B 申请公布日期 2015.03.11
申请号 CN201210381071.1 申请日期 2012.10.10
申请人 西北工业大学 发明人 史忠科
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种含有分布参数项的飞行器控制器设计方法,其特征在于包括以下步骤:步骤一、在给定高度、马赫数条件下通过风洞或飞行试验得到含有结构振动不确定性的飞行器模型为:<maths num="0001" id="cmaths0001"><math><![CDATA[<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><mfrac><mo>&PartialD;</mo><mrow><mo>&PartialD;</mo><mi>&eta;</mi></mrow></mfrac><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>+</mo><mo>{</mo><mi>B</mi><mo>+</mo><mfrac><mo>&PartialD;</mo><mrow><mo>&PartialD;</mo><mi>&eta;</mi></mrow></mfrac><mo>[</mo><mi>&Delta;B</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>]</mo><mo>}</mo><mi>u</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000640206350000011.GIF" wi="1747" he="148" /></maths>式中,x∈R<sup>n</sup>,u∈R<sup>m</sup>分别为状态和输入向量,η为结构振动参数向量,A<sub>0</sub>、B、A<sub>non</sub>为已知的系数矩阵,<img file="FDA0000640206350000012.GIF" wi="539" he="143" />为未知结构振动项;按照不同的飞行区域和结构振动大小,将(1)式右端的状态项表达成:<maths num="0002" id="cmaths0002"><math><![CDATA[<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><mfrac><mo>&PartialD;</mo><mrow><mo>&PartialD;</mo><mi>&eta;</mi></mrow></mfrac><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>&eta;</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><mfrac><mo>&PartialD;</mo><mrow><mo>&PartialD;</mo><mi>&eta;</mi></mrow></mfrac><mo>[</mo><mi>&Delta;B</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>]</mo><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><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><mo>,</mo><msub><mi>&eta;</mi><mrow><mi>i</mi><mi>min</mi></mrow></msub><mo>&le;</mo><mi>&eta;</mi><mo>&lt;</mo><msub><mi>&eta;</mi><mrow><mi>i</mi><mi>max</mi></mrow></msub><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>]]></math><img file="FDA0000640206350000013.GIF" wi="1742" he="296" /></maths>式中,A<sub>0i</sub>、B<sub>0i</sub>为已知的常数矩阵,ΔA<sub>0i</sub>ΔB<sub>0i</sub>为未知矩阵,x<sub>i min</sub>、x<sub>i max</sub>分别为第i个飞行区域所对应的状态最小和最大值,η<sub>i min</sub>、η<sub>i max</sub>分别为第i个飞行区域所对应结构振动参数的最小和最大值,r为正整数,i为下标;在飞行区域和结构振动参数x<sub>i min</sub>≤x&lt;x<sub>i max</sub>,η<sub>i min</sub>≤η&lt;η<sub>i max</sub>范围内,飞行控制器为:u(t)=‑K<sub>i</sub>x(t),式中,K<sub>i</sub>为常数反馈矩阵;带入(1)式中,有:<maths num="0003" id="cmaths0003"><math><![CDATA[<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>]]></math><img file="FDA0000640206350000014.GIF" wi="1206" he="77" /></maths>步骤二、选取(A<sub>0i</sub>‑B<sub>0i</sub>K<sub>i</sub>)的特征值各不相同且实部为负,设计反馈矩阵K<sub>i</sub>使得满足条件:<maths num="0004" id="cmaths0004"><math><![CDATA[<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>]]></math><img file="FDA0000640206350000015.GIF" wi="1514" he="86" /></maths>该控制器使得<maths num="0005" id="cmaths0005"><math><![CDATA[<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>]]></math><img file="FDA0000640206350000016.GIF" wi="1200" he="78" /></maths>鲁棒稳定;式中,M<sub>i</sub>为线性变换矩阵,<maths num="0006" id="cmaths0006"><math><![CDATA[<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><mfenced open='[' close=']'><mtable><mtr><mtd><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></mtd><mtd><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></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo></mtd><mtd><msub><mi>&sigma;</mi><mi>in</mi></msub><mo>+</mo><mi>j</mi><msub><mi>&omega;</mi><mi>in</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000640206350000017.GIF" wi="1932" he="86" /></maths>σ<sub>ik</sub>,ω<sub>ik</sub>,k=1,2,…,n为实数,jω<sub>ik</sub>,k=1,2,…,n表示虚数,diag为对角矩阵符号,<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>&Lambda;</mi><mi>i</mi></msub><mo>=</mo><mi>diag</mi><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>&sigma;</mi><mrow><mi>i</mi><mn>1</mn></mrow><mn>2</mn></msubsup><mo>,</mo></mtd><mtd><msubsup><mi>&sigma;</mi><mrow><mi>i</mi><mn>2</mn></mrow><mn>2</mn></msubsup><mo>,</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo></mtd><mtd><msubsup><mi>&sigma;</mi><mi>in</mi><mn>2</mn></msubsup></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000640206350000018.GIF" wi="928" he="97" /></maths>ΔA<sub>0i</sub>‑ΔB<sub>0i</sub>K<sub>i</sub>通常假设为ΔA<sub>0i</sub>‑ΔB<sub>0i</sub>K<sub>i</sub>=H<sub>i</sub>F<sub>i</sub>W<sub>i</sub>,H<sub>i</sub>,W<sub>i</sub>均假设为已知矩阵,0&lt;F<sub>i</sub>≤I,I=diag[1,1,…,1]为单位阵。
地址 710072 陕西省西安市友谊西路127号