发明名称 高速飞行器快速在线综合识别方法
摘要 本发明公开了一种高速飞行器快速在线综合识别方法,用于解决现有基于特征模型的高超声速飞行器鲁棒控制方法控制精度低的技术问题。技术方案是首先建立高速飞行器小扰动线性化模型,再建立高速飞行器纵向运动的定义特征状态量和特征模型,然后对高速飞行器的特征状态量进行在线快速综合识别。由于该方法以传感器可测姿态所描述的系统瞬时动力学方程为系统特征模型,对系统特征模型进行小扰动线性化,获得可用于控制的线性扰动模型,将扰动量表示为真实量与平衡点的差,并通过定义新的结构参数,将扰动模型变形为真实物理量所描述的瞬时方程,将结构参数定义为特征状态量,将瞬时方程定义为特征模型,提高了高速飞行器的控制精度。
申请公布号 CN103197545A 申请公布日期 2013.07.10
申请号 CN201310059417.0 申请日期 2013.02.25
申请人 西北工业大学 发明人 周军;林鹏;王霄婷;王立祺;葛振振
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 1.一种高速飞行器快速在线综合识别方法,其特征在于包括以下步骤:步骤一、将高速飞行器视为刚体,其纵向运动模型是<img file="FDA00002853851700011.GIF" wi="1334" he="416" />分离长短周期运动,并采用小扰动线性化的方法来获取短时姿态低阶扰动模型;若取弹道某点V<sub>0</sub>,θ<sub>0</sub>,h<sub>0</sub>处的线性化参考点为α<sub>0</sub>及其配平舵偏δ<sub>z0</sub>,俯仰角速度<img file="FDA000028538517000111.GIF" wi="251" he="77" />则高速飞行器小扰动线性化模型是<maths num="0001"><![CDATA[<math><mrow><mfenced open='' close='}'><mtable><mtr><mtd><mover><mi>&alpha;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><msub><mi>&omega;</mi><mi>z</mi></msub></mtd></mtr><mtr><mtd><msub><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mi>z</mi></msub><mo>=</mo><msub><mi>a</mi><mrow><mi>pd</mi><mn>1</mn></mrow></msub><msub><mi>&omega;</mi><mi>z</mi></msub><mo>+</mo><msub><mi>a</mi><mrow><mi>pd</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>&alpha;</mi><mo>-</mo><msub><mi>&alpha;</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>b</mi><mi>pd</mi></msub><mrow><mo>(</mo><msub><mi>&delta;</mi><mi>z</mi></msub><mo>-</mo><msub><mi>&delta;</mi><mrow><mi>z</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,<img file="FDA00002853851700014.GIF" wi="354" he="70" />为阻尼力矩所产生的角加速度对角速度的导数;<img file="FDA00002853851700015.GIF" wi="347" he="69" />为稳定力矩所产生的角加速度对扰动攻角的导数,代表单位扰动攻角产生角加速度的能力;<img file="FDA00002853851700016.GIF" wi="337" he="70" />为操纵力矩所对应的角加速度对扰动舵偏的导数,代表单位扰动舵偏产生角加速度的能力;步骤二、建立高速飞行器纵向运动的特征状态量为<maths num="0002"><![CDATA[<math><mrow><mfenced open='' close='}'><mtable><mtr><mtd><msub><mi>a</mi><mrow><mi>p</mi><mn>2</mn></mrow></msub><mo>=</mo><msub><mi>a</mi><mrow><mi>pd</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&alpha;</mi><mn>0</mn></msub><mo>/</mo><mi>&alpha;</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>b</mi><mi>p</mi></msub><mo>=</mo><msub><mi>b</mi><mi>pd</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&delta;</mi><mrow><mi>z</mi><mn>0</mn></mrow></msub><mo>/</mo><msub><mi>&delta;</mi><mi>z</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,a<sub>p2</sub>是稳定力矩与实际攻角的比值,代表单位实际攻角产生角加速度的能力;b<sub>p</sub>是操纵力矩与实际舵偏的比值,代表单位实际舵偏产生角加速度的能力;基于式(2)和式(3),得到高速飞行器纵向运动的特征模型为<maths num="0003"><![CDATA[<math><mrow><mfenced open='' close='}'><mtable><mtr><mtd><mover><mi>&alpha;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><msub><mi>&omega;</mi><mi>z</mi></msub></mtd></mtr><mtr><mtd><msub><mover><mi>&omega;</mi><mo>&CenterDot;</mo></mover><mi>z</mi></msub><mo>=</mo><msub><mi>a</mi><mrow><mi>pd</mi><mn>1</mn></mrow></msub><msub><mi>&omega;</mi><mi>z</mi></msub><mo>+</mo><msub><mi>a</mi><mrow><mi>p</mi><mn>2</mn></mrow></msub><mi>&alpha;</mi><mo>+</mo><msub><mi>b</mi><mi>p</mi></msub><msub><mi>&delta;</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤三、特征状态量在线快速综合识别步骤为:(a)采集传感器数据,即飞行姿态矢量<img file="FDA00002853851700019.GIF" wi="408" he="70" />舵偏矢量<img file="FDA000028538517000110.GIF" wi="286" he="86" />飞行弹道状态矢量Γ=[Q,V,H]<sup>T</sup>;(b)用采集的传感器数据,按照下式计算a<sub>pd1</sub>,<maths num="0004"><![CDATA[<math><mrow><msub><mover><mi>a</mi><mo>^</mo></mover><mrow><mi>pd</mi><mn>1</mn></mrow></msub><mo>=</mo><msub><mover><mi>k</mi><mo>^</mo></mover><mi>z</mi></msub><mi>Q</mi><msup><mi>SL</mi><mn>2</mn></msup><mo>/</mo><msub><mi>J</mi><mi>z</mi></msub><mi>V</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,<img file="FDA00002853851700022.GIF" wi="41" he="68" />是阻尼力矩系数<img file="FDA00002853851700023.GIF" wi="69" he="60" />的近似值,根据气动阻尼特性将其拟合为可测物理量的函数用于在线计算,或者根据<img file="FDA00002853851700024.GIF" wi="69" he="60" />的变化范围选择为随某可测状态分段变化的常值;对于欠阻尼受控对象,<img file="FDA00002853851700025.GIF" wi="41" he="68" />作定常处理;(c)计算配平系数μ=-a<sub>p2</sub>/b<sub>p</sub>;在采样周期足够小的情况下,假设参数a<sub>p2</sub>,b<sub>p</sub>变化可忽略不计,配平系数估算式为<maths num="0005"><![CDATA[<math><mrow><mover><mi>&mu;</mi><mo>^</mo></mover><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><msub><mi>x</mi><mn>4</mn></msub><mo>-</mo><msub><mover><mi>a</mi><mo>^</mo></mover><mn>1</mn></msub><msub><mi>x</mi><mn>3</mn></msub><mo>)</mo></mrow><msub><mi>u</mi><mn>1</mn></msub><mo>-</mo><mrow><mo>(</mo><msub><mi>x</mi><mn>3</mn></msub><mo>-</mo><msub><mover><mi>a</mi><mo>^</mo></mover><mn>1</mn></msub><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow><msub><mi>u</mi><mn>2</mn></msub></mrow><mrow><mrow><mo>(</mo><msub><mi>x</mi><mn>4</mn></msub><mo>-</mo><msub><mover><mi>a</mi><mo>^</mo></mover><mn>1</mn></msub><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><mrow><mo>(</mo><msub><mi>x</mi><mn>3</mn></msub><mo>-</mo><msub><mover><mi>a</mi><mo>^</mo></mover><mn>1</mn></msub><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow><msub><mi>x</mi><mn>2</mn></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>(d)用配平系数估值来计算特征状态量<maths num="0006"><![CDATA[<math><mrow><mfenced open='' close='}'><mtable><mtr><mtd><msub><mover><mi>a</mi><mo>^</mo></mover><mrow><mi>p</mi><mn>2</mn></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>x</mi><mn>3</mn></msub><mo>-</mo><msub><mover><mi>a</mi><mo>^</mo></mover><mn>1</mn></msub><msub><mi>x</mi><mn>2</mn></msub></mrow><mrow><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mi>u</mi><mn>1</mn></msub><mo>/</mo><mover><mi>&mu;</mi><mo>^</mo></mover></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mover><mi>b</mi><mo>^</mo></mover><mi>p</mi></msub><mo>=</mo><mo>-</mo><mfrac><msub><mover><mi>a</mi><mo>^</mo></mover><mrow><mi>p</mi><mn>2</mn></mrow></msub><mover><mi>&mu;</mi><mo>^</mo></mover></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow><mo>.</mo></mrow></math>]]></maths>
地址 710072 陕西省西安市友谊西路127号