发明名称 基于有限状态机的自动编队飞行控制方法
摘要 本发明公开了一种基于有限状态机的自动编队飞行控制方法,包括如下步骤:步骤一:基于有限状态自动机的原理,针对编队飞行训练任务要求,对编队飞行中可能遇到的情况进行划分、归纳,确定出僚机可能出现的飞行模式。步骤二:根据步骤一中飞行模式,确定各飞行模式之间的转换条件。步骤三:结合步骤一中归纳总结出的飞行模式和步骤二中确定的各个飞行模式之间的转换条件,列出状态转换表。步骤四,为步骤一中所述的各飞行模式设计控制方法,为编队队形保持设计控制率。本发明提出了僚机在编队飞行过程中在不同飞行模式间进行自动转换的方法,有效模拟僚机飞行员的智能;设计了编队队形保持飞行模式下的控制率,能够满足编队队形保持的需求。
申请公布号 CN101464692A 申请公布日期 2009.06.24
申请号 CN200910077120.0 申请日期 2009.01.16
申请人 北京航空航天大学 发明人 龚光红;孟光磊;马耀飞
分类号 G05D1/00(2006.01)I 主分类号 G05D1/00(2006.01)I
代理机构 北京永创新实专利事务所 代理人 周长琪
主权项 1、一种基于有限状态机的自动编队飞行控制方法,其特征在于,包含以下步骤:步骤一:基于有限状态自动机的原理,针对编队飞行训练任务要求,对编队飞行中遇到的情况进行划分、归纳,确定出僚机出现的飞行模式;步骤二:根据步骤一中飞行模式,确定各飞行模式之间的转换条件;步骤三:结合步骤一中归纳总结出的飞行模式和步骤二中确定的各个飞行模式之间的转换条件,列出状态转换表;步骤四,为步骤一中所述的各飞行模式设计控制方法,为编队队形保持设计控制率;飞机模型控制量的初始值为δ<sub>T</sub>=δ<sub>e</sub>=δ<sub>a</sub>=δ<sub>r</sub>=0,landGearSwitch=0,peedBreakSwitch=1;δ<sub>T</sub>:油门杆输入量,<maths num="0001"><![CDATA[<math><mrow><msub><mi>&delta;</mi><mi>T</mi></msub><mo>&Element;</mo><mrow><mo>[</mo><msubsup><mi>&delta;</mi><mi>min</mi><mi>T</mi></msubsup><mo>,</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>T</mi></msubsup><mo>]</mo></mrow><mo>,</mo></mrow></math>]]></maths><img file="A200910077120C00022.GIF" wi="65" he="52" />和<img file="A200910077120C00023.GIF" wi="69" he="52" />分别是δ<sub>T</sub>的最大值和最小值;δ<sub>e</sub>:俯仰舵偏角,<maths num="0002"><![CDATA[<math><mrow><msub><mi>&delta;</mi><mi>e</mi></msub><mo>&Element;</mo><mrow><mo>[</mo><msubsup><mi>&delta;</mi><mi>min</mi><mi>e</mi></msubsup><mo>,</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>e</mi></msubsup><mo>]</mo></mrow><mo>,</mo></mrow></math>]]></maths><img file="A200910077120C00025.GIF" wi="63" he="47" />和<img file="A200910077120C00026.GIF" wi="68" he="47" />分别是δ<sub>e</sub>的最大值和最小值;δ<sub>a</sub>:滚转舵偏角,<maths num="0003"><![CDATA[<math><mrow><msub><mi>&delta;</mi><mi>a</mi></msub><mo>&Element;</mo><mrow><mo>[</mo><msubsup><mi>&delta;</mi><mi>min</mi><mi>a</mi></msubsup><mo>,</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>a</mi></msubsup><mo>]</mo></mrow><mo>,</mo></mrow></math>]]></maths><img file="A200910077120C00028.GIF" wi="65" he="46" />和<img file="A200910077120C00029.GIF" wi="69" he="46" />分别是δ<sub>a</sub>的最大值和最小值;δ<sub>r</sub>:方向舵偏角,<maths num="0004"><![CDATA[<math><mrow><msub><mi>&delta;</mi><mi>r</mi></msub><mo>&Element;</mo><mrow><mo>[</mo><msubsup><mi>&delta;</mi><mi>min</mi><mi>r</mi></msubsup><mo>,</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>r</mi></msubsup><mo>]</mo></mrow><mo>,</mo></mrow></math>]]></maths><img file="A200910077120C000211.GIF" wi="63" he="47" />和<img file="A200910077120C000212.GIF" wi="68" he="47" />分别是δ<sub>r</sub>的最大值和最小值;landGearSwitch:起落架收放开关,landGearSwitch∈{0,1},0表示放下起落架,1表示收起起落架;speedBreakSwitch:减速板收放开关,speedBreakSwitch∈{0,1},0表示放下减速板,1表示收起减速板;步骤一中总结的各飞行模式的控制方法如下:<maths num="0005"><![CDATA[<math><mrow><msub><mi>C</mi><msub><mi>s</mi><mn>0</mn></msub></msub><mo>:</mo><msub><mi>&delta;</mi><mi>T</mi></msub><mo>=</mo><msubsup><mi>&delta;</mi><mi>min</mi><mi>T</mi></msubsup><mo>,</mo></mrow></math>]]></maths>即僚机发动机保持关闭状态;<maths num="0006"><![CDATA[<math><mrow><msub><mi>C</mi><msub><mi>s</mi><mn>1</mn></msub></msub><mo>:</mo><msub><mi>&delta;</mi><mi>T</mi></msub><mo>=</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>T</mi></msubsup><mo>,</mo></mrow></math>]]></maths>brake=B<sub>max</sub>,即僚机油门加到最大,同时踩住刹车;<maths num="0007"><![CDATA[<math><mrow><msub><mi>C</mi><msub><mi>s</mi><mn>2</mn></msub></msub><mo>:</mo><mi>brake</mi><mo>=</mo><msub><mi>B</mi><mi>min</mi></msub><mo>,</mo></mrow></math>]]></maths>即僚机放开刹车,僚机开始滑跑;<img file="A200910077120C000216.GIF" wi="66" he="52" /><img file="A200910077120C000217.GIF" wi="67" he="51" /><img file="A200910077120C000218.GIF" wi="67" he="52" /><img file="A200910077120C000219.GIF" wi="48" he="51" />均为编队队形保持控制方法,具有相同的控制率,只是控制率参数的取值不同;<img file="A200910077120C000220.GIF" wi="77" he="51" />如果长机对僚机下达收起起落架指令,则landGearSwitch=1;如果长机对僚机下达放下起落架指令,则landGearSwitch=0;<img file="A200910077120C000221.GIF" wi="73" he="51" />如果长机对僚机下达收起减速板指令,则speedBreakSwitch=1;如果长机对僚机下达放下减速板指令,则speedBreakSwitch=0;<maths num="0008"><![CDATA[<math><mrow><msub><mi>C</mi><msub><mi>s</mi><mn>9</mn></msub></msub><mo>:</mo><msub><mi>&delta;</mi><mi>T</mi></msub><mo>=</mo><msubsup><mi>&delta;</mi><mi>min</mi><mi>T</mi></msubsup><mo>,</mo></mrow></math>]]></maths>brake=B<sub>max</sub>;<maths num="0009"><![CDATA[<math><mrow><msub><mi>C</mi><msub><mi>s</mi><mn>10</mn></msub></msub><mo>:</mo><mi>breake</mi><mo>=</mo><mn>0</mn><mo>,</mo></mrow></math>]]></maths>即僚机停止在跑道上后,放开刹车,仿真结束;其中brake为刹车输入量,break∈[B<sub>min</sub>,B<sub>max</sub>],B<sub>min</sub>和B<sub>max</sub>分别是brake的最大值和最小值;<img file="A200910077120C000224.GIF" wi="43" he="52" />表示飞行模式s<sub>i</sub>的控制方法,i∈{0,1,2,3,4,5,6,7,8,9,10},每一个飞行模式的控制方法中没有描述的其它控制量的取值与转换到该飞行模式之前的飞行模式的控制量的取值相等;所述的编队队形保持飞行模式下的控制方法<img file="A200910077120C00031.GIF" wi="70" he="54" /><img file="A200910077120C00032.GIF" wi="71" he="54" /><img file="A200910077120C00033.GIF" wi="72" he="54" /><img file="A200910077120C00034.GIF" wi="52" he="54" />的飞行控制率的设计具体如下:(a)油门控制率:<maths num="0010"><![CDATA[<math><mrow><msub><mi>&Delta;&delta;</mi><mi>T</mi></msub><mo>=</mo><msub><mi>K</mi><mi>xp</mi></msub><msub><mi>x</mi><mi>e</mi></msub><mo>+</mo><msub><mi>K</mi><mi>xd</mi></msub><msub><mover><mi>x</mi><mo>.</mo></mover><mi>e</mi></msub><mo>+</mo><msub><mi>K</mi><mi>xi</mi></msub><mo>&Integral;</mo><msub><mi>x</mi><mi>e</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0011"><![CDATA[<math><mrow><msub><mi>&delta;</mi><mi>T</mi></msub><mo>=</mo><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>&delta;</mi><mi>max</mi><mi>T</mi></msubsup><mo>,</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>T</mi></msubsup><mo>&le;</mo><msub><mi>&delta;</mi><mrow><mi>T</mi><mo>_</mo></mrow></msub><mo>+</mo><msub><mi>&Delta;&delta;</mi><mi>T</mi></msub></mtd></mtr><mtr><mtd><msub><mi>&delta;</mi><mrow><mi>T</mi><mo>_</mo></mrow></msub><mo>+</mo><msub><mi>&Delta;&delta;</mi><mi>T</mi></msub><mo>,</mo><msubsup><mi>&delta;</mi><mi>min</mi><mi>T</mi></msubsup><mo>&lt;</mo><msub><mi>&delta;</mi><mrow><mi>T</mi><mo>_</mo></mrow></msub><mo>+</mo><msub><mi>&Delta;&delta;</mi><mi>T</mi></msub><mo>&lt;</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>T</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>&delta;</mi><mi>min</mi><mi>T</mi></msubsup><mo>,</mo><msub><mi>&delta;</mi><mrow><mi>T</mi><mo>_</mo></mrow></msub><mo>+</mo><msub><mi>&Delta;&delta;</mi><mi>T</mi></msub><mo>&le;</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>T</mi></msubsup></mtd></mtr></mtable></mfenced></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中Δδ<sub>T</sub>是油门的增量,δ<sub>T_</sub>是当前时刻的油门值,x<sub>e</sub>表示当前时刻,僚机与长机质心连线矢量在僚机机体x轴的实际投影值与x<sub>command</sub>的差值,K<sub>xp</sub>是比例系数,K<sub>xd</sub>是微分系数,K<sub>xi</sub>是积分系数;(b)滚转舵控制率:<img file="A200910077120C00037.GIF" wi="1551" he="75" /><maths num="0012"><![CDATA[<math><mrow><msub><mi>&delta;</mi><mi>a</mi></msub><mo>=</mo><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>&delta;</mi><mi>max</mi><mi>a</mi></msubsup><mo>,</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>a</mi></msubsup><mo>&le;</mo><msub><mi>&delta;</mi><mrow><mi>a</mi><mo>_</mo></mrow></msub><mo>+</mo><msub><mi>&Delta;&delta;</mi><mi>a</mi></msub></mtd></mtr><mtr><mtd><msub><mi>&delta;</mi><mrow><mi>a</mi><mo>_</mo></mrow></msub><mo>+</mo><msub><mi>&Delta;&delta;</mi><mi>a</mi></msub><mo>,</mo><msubsup><mi>&delta;</mi><mi>min</mi><mi>a</mi></msubsup><mo>&lt;</mo><msub><mi>&delta;</mi><mrow><mi>a</mi><mo>_</mo></mrow></msub><mo>+</mo><msub><mi>&Delta;&delta;</mi><mi>a</mi></msub><mo>&lt;</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>a</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>&delta;</mi><mi>min</mi><mi>a</mi></msubsup><mo>,</mo><msub><mi>&delta;</mi><mrow><mi>a</mi><mo>_</mo></mrow></msub><mo>+</mo><msub><mi>&Delta;&delta;</mi><mi>a</mi></msub><mo>&le;</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>a</mi></msubsup></mtd></mtr></mtable></mfenced></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中;偏航角<img file="A200910077120C00039.GIF" wi="53" he="34" />机体轴x在水平面上的投影与地轴x<sub>g</sub>间夹角,机头右偏航为正;滚转角φ:机体轴z与通过机体轴x的铅垂面间夹角,飞机向右滚转时为正;Δδ<sub>a</sub>是滚转舵偏角的增量,δ<sub>a_</sub>是当前时刻的滚转舵偏角,y<sub>e</sub>表示当前时刻,僚机与长机质心连线矢量在僚机机体y轴的实际投影值与y<sub>command</sub>的差值,φ<sub>e</sub>表示长机滚转角与僚机滚转角的差值,<img file="A200910077120C000310.GIF" wi="38" he="35" />表示长机方向角与僚机方向角的差值,K<sub>yp</sub>,K<sub>yd</sub>,K<sub>yi</sub>分别是y<sub>e</sub>的比例系数,微分系数和积分系数,K<sub>φp</sub>,K<sub>φd</sub>分别是φ<sub>e</sub>的比例系数和微分系数,<img file="A200910077120C000311.GIF" wi="100" he="53" /><img file="A200910077120C000312.GIF" wi="71" he="52" />分别是<img file="A200910077120C000313.GIF" wi="38" he="35" />的比例系数和微分系数;(c)俯仰舵控制率<maths num="0013"><![CDATA[<math><mrow><msub><mi>&Delta;&delta;</mi><mi>e</mi></msub><mo>=</mo><msub><mi>K</mi><mi>zp</mi></msub><msub><mi>z</mi><mi>e</mi></msub><mo>+</mo><msub><mi>K</mi><mi>zd</mi></msub><msub><mover><mi>z</mi><mo>.</mo></mover><mi>e</mi></msub><mo>+</mo><msub><mi>K</mi><mi>zi</mi></msub><mo>&Integral;</mo><msub><mi>z</mi><mi>e</mi></msub><mo>+</mo><msub><mi>K</mi><mi>&theta;p</mi></msub><msub><mi>&theta;</mi><mi>e</mi></msub><mo>+</mo><msub><mi>K</mi><mi>&theta;d</mi></msub><msub><mover><mi>&theta;</mi><mo>.</mo></mover><mi>e</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0014"><![CDATA[<math><mrow><msub><mi>&delta;</mi><mi>e</mi></msub><mo>=</mo><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>&delta;</mi><mi>max</mi><mi>e</mi></msubsup><mo>,</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>e</mi></msubsup><mo>&le;</mo><msub><mi>&delta;</mi><mrow><mi>e</mi><mo>_</mo></mrow></msub><mo>+</mo><msub><mi>&Delta;&delta;</mi><mi>e</mi></msub></mtd></mtr><mtr><mtd><msub><mi>&delta;</mi><mrow><mi>e</mi><mo>_</mo></mrow></msub><mo>+</mo><msub><mi>&Delta;&delta;</mi><mi>e</mi></msub><mo>,</mo><msubsup><mi>&delta;</mi><mi>min</mi><mi>e</mi></msubsup><mo>&lt;</mo><msub><mi>&delta;</mi><mrow><mi>e</mi><mo>_</mo></mrow></msub><mo>+</mo><msub><mi>&Delta;&delta;</mi><mi>e</mi></msub><mo>&lt;</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>e</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>&delta;</mi><mi>min</mi><mi>e</mi></msubsup><mo>,</mo><msub><mi>&delta;</mi><mrow><mi>e</mi><mo>_</mo></mrow></msub><mo>+</mo><msub><mi>&Delta;&delta;</mi><mi>e</mi></msub><mo>&le;</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>e</mi></msubsup></mtd></mtr></mtable></mfenced></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中俯仰角θ:机体轴x与水平面间夹角,抬头为正;Δδ<sub>e</sub>是俯仰舵偏角的增量,δ<sub>e_</sub>是当前时刻的俯仰舵偏角,z<sub>e</sub>表示当前时刻长机高度与僚机高度的差值,θ<sub>e</sub>是长机俯仰角与僚机俯仰角的差值,K<sub>zp</sub>,K<sub>zd</sub>,K<sub>zi</sub>分别是z<sub>e</sub>的比例系数,微分系数和积分系数,K<sub>θp</sub>,K<sub>θd</sub>分别是<img file="A200910077120C000316.GIF" wi="33" he="47" />的比例系数和微分系数;(d)方向舵控制率:Δδ<sub>r</sub>=K<sub>rp</sub>r<sub>w</sub>          (7)<maths num="0015"><![CDATA[<math><mrow><msub><mi>&delta;</mi><mi>r</mi></msub><mo>=</mo><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>&delta;</mi><mi>max</mi><mi>r</mi></msubsup><mo>,</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>r</mi></msubsup><mo>&le;</mo><msub><mi>&delta;</mi><mrow><mi>r</mi><mo>_</mo></mrow></msub><mo>+</mo><msub><mi>&Delta;&delta;</mi><mi>r</mi></msub></mtd></mtr><mtr><mtd><msub><mi>&delta;</mi><mrow><mi>r</mi><mo>_</mo></mrow></msub><mo>+</mo><msub><mi>&Delta;&delta;</mi><mi>r</mi></msub><mo>,</mo><msubsup><mi>&delta;</mi><mi>min</mi><mi>r</mi></msubsup><mo>&lt;</mo><msub><mi>&delta;</mi><mrow><mi>r</mi><mo>_</mo></mrow></msub><mo>+</mo><msub><mi>&Delta;&delta;</mi><mi>r</mi></msub><mo>&lt;</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>r</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>&delta;</mi><mi>min</mi><mi>r</mi></msubsup><mo>,</mo><msub><mi>&delta;</mi><mrow><mi>r</mi><mo>_</mo></mrow></msub><mo>+</mo><msub><mi>&Delta;&delta;</mi><mi>r</mi></msub><mo>&le;</mo><msubsup><mi>&delta;</mi><mi>max</mi><mi>r</mi></msubsup></mtd></mtr></mtable></mfenced></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中Δδ<sub>r</sub>是方向舵偏角的增量,δ<sub>r_</sub>是当前时刻的方向舵偏角,r<sub>w</sub>是僚机的偏航角速度,K<sub>rp</sub>是r<sub>w</sub>的比例系数;根据控制率公式,使僚机保持编队队形。
地址 100083北京市海淀区学院路37号