CN107367941A - 基于非线性增益的高超声速飞行器攻角观测方法 - Google Patents
基于非线性增益的高超声速飞行器攻角观测方法 Download PDFInfo
- Publication number
- CN107367941A CN107367941A CN201710743894.7A CN201710743894A CN107367941A CN 107367941 A CN107367941 A CN 107367941A CN 201710743894 A CN201710743894 A CN 201710743894A CN 107367941 A CN107367941 A CN 107367941A
- Authority
- CN
- China
- Prior art keywords
- mover
- mrow
- msub
- angle
- hypersonic aircraft
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种基于非线性增益的高超声速飞行器攻角观测方法,通过测量高超声速飞行器弹道倾角与飞行器姿态角速率,构造攻角观测系统,通过弹道倾角与弹道倾角估计值的误差、以及姿态角速度与姿态角速率估计值的误差,构造一种非线性变增益误差规律,调节攻角观测值,通过攻角观测值来进行飞行器受力与力矩的估计,为观测器系统提供力与力矩信息,实现攻角观测值对攻角实际值的精准估计。该方法不同于一般飞行器攻角观测方法在于不需要测量飞行器的高度与速度,观测系统所需的信息较少,而且由于采用了变增益误差调节规律,攻角观测的响应速度很快,而且精度较高。
Description
技术领域
本发明属于控制技术领域,具体而言,涉及一种基于非线性增益的高超声速飞行器攻角观测方法。
背景技术
高超声速飞行器具有飞行速度快、费用低的特点,可确保廉价进入太空,也可作为全球快速运输机与空间打击武器,其鲜明的军事价值使它备受世界各国的高度重视。
高超声速飞行器与高超工程研究是目前我国集中力量并重点关注的重大项目之一,其难度与意义不亚于飞船上天与登月计划。高超声速控制问题也是近十年来国内外控制领域的热点与难点问题,而高超声速飞行器的攻角由于工程测量比较困难,而且测量精度不高,因此其攻角观测的方法非常有经济与军事价值。当然,由于高超声速运动的复杂性,其攻角准确观测与估计的复杂性也毋庸置疑。
发明内容
为了达到上述目的,本发明提供一种基于非线性增益的高超声速飞行器攻角观测方法,进而至少在一定程度上克服由于相关技术的限制和缺陷而导致的一个或者多个问题。
本发明所采用的技术方案是,一种基于非线性增益的高超声速飞行器攻角观测方法,按照以下步骤进行:
步骤一,测量弹道倾角以及姿态角速度,并根据测量到的弹道倾角以及姿态角速度构建一观测系统;
步骤二,根据所述观测系统观测高超声速飞行器的攻角。
进一步的,所述步骤一包括:
首先,采用惯导系统测量高超声速飞行器的弹道倾角γ以及姿态角速度q;
其次,根据弹道倾角γ以及姿态角速度q构建观测如下观测系统:
其中,α为高超声速飞行器的攻角;为弹道倾角观测值,为姿态角速度观测值,为攻角观测值;
为的导数,为的导数,为的导数;
为估计误差;并且有:
为高超声速飞行器质量的估计值;g为重力加速度,取值为9.8;
为高超声速飞行器的速度估计值;为高超声速飞行器的升力估计值;
为高超声速飞行器的推力估计值;为高超声速飞行器的俯仰力矩估计值;
为高超声速飞行器的转动惯量估计值;
k22、k23、ε1、k32、k33、ε2、k42、k43、ε3为控制参数;
其中,为k24的导数,k24b与k24a为正数;
其中,为k34的导数,k34b与k34a为正数;
其中,为k44的导数,k44b与k44a为正数。
进一步的,高超声速飞行器推力升力与力矩的计算方法包括:
首先,高超声速飞行器的升力估计值的计算方法包括:
其中,S为飞行器特征面积;ρ为大气密度;为高超声速飞行器的速度估计值;
其次,高超声速飞行器的推力估计值的计算方法包括:
其中:S为飞行器特征面积;ρ为大气密度;为高超声速飞行器的速度估计值;
β为发动机的油门供油因子;
最后,高超声速飞行器的俯仰力矩估计值的计算方法包括:
其中,为高超声速飞行器的特征因子;δ为高超声速飞行器的俯仰舵偏角;
S为飞行器特征面积;ρ为大气密度;
为姿态角速度观测值;为攻角观测值;为高超声速飞行器的速度估计值。
进一步的,各所述控制参数的取值为:
k22=300、k33=500、ε1=0.5、ε2=0.5、ε3=0.5、k23=0.2、k32=0.1、k42=0.066、k43=4.544。
本发明提供了一种基于非线性增益的高超声速飞行器攻角观测方法,对高超声速飞行器构造观测系统,尤其是提供了一种变增益非线性的误差调节规律,能够有效地改善攻角观测的快速性与准确性。该方法不同于一般飞行器攻角观测方法在于不需要测量飞行器的高度与速度,观测系统所需的信息较少,而且由于采用了变增益误差调节规律,攻角观测的响应速度很快,而且精度较高。因此,本发明所提供的方法具有很高的工程应用价值与经济价值。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明提供的一种基于非线性变增益的高超声速飞行器攻角观测技术原理框图。
图2是本发明实施例所提供方法的高超声速飞行器的攻角曲线。
图3是本发明实施例所提供方法的高超声速飞行器的攻角观测曲线。
图4是本发明实施例所提供方法的高超声速飞行器的攻角观测误差曲线。
图5是本发明实施例所提供方法的高超声速飞行器的弹道倾角曲线。
图6是本发明实施例所提供方法的高超声速飞行器的弹道倾角估计曲线。
图7是本发明实施例所提供方法的高超声速飞行器的弹道倾角观测误差曲线。
图8是本发明实施例所提供方法的高超声速飞行器的姿态角速率曲线。
图9是本发明实施例所提供方法的高超声速飞行器的姿态角速率观测曲线。
图10是本发明实施例所提供方法的高超声速飞行器的姿态角速率观测误差曲线。
图11是本发明实施例所提供方法的攻角与攻角观测比较曲线。
图12是本发明实施例所提供方法的弹道倾角与弹道倾角观测比较曲线。
图13是本发明实施例所提供方法的姿态角速率与姿态角速率观测比较曲线。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明公开了一种基于非线性增益的高超声速飞行器攻角观测方法,参考图1所示,可以通过测量高超声速飞行器弹道倾角与飞行器姿态角速率,构造攻角观测系统,通过弹道倾角与弹道倾角估计值的误差、以及姿态角速度与姿态角速率估计值的误差,构造一种非线性变增益误差规律,调节攻角观测值,通过攻角观测值来进行飞行器受力与力矩的估计,为观测器系统提供力与力矩信息,实现攻角观测值对攻角实际值的精准估计。
进一步的,上述一种基于非线性增益的高超声速飞行器攻角观测方法可以按照以下步骤进行:
步骤一,测量弹道倾角以及姿态角速度,并根据测量到的弹道倾角以及姿态角速度构建一观测系统。
首先,采用惯导系统测量高超声速飞行器的弹道倾角γ以及姿态角速度q;
其次,根据弹道倾角γ以及姿态角速度q构建观测如下观测系统:
其中,α为高超声速飞行器的攻角;为弹道倾角观测值,为姿态角速度观测值,为攻角观测值;
为的导数,为的导数,为的导数;
为估计误差;并且有:
为高超声速飞行器质量的估计值;g为重力加速度,取值为9.8;
为高超声速飞行器的速度估计值;为高超声速飞行器的升力估计值;
为高超声速飞行器的推力估计值;为高超声速飞行器的俯仰力矩估计值;
为高超声速飞行器的转动惯量估计值;
k22、k23、ε1、k32、k33、ε2、k42、k43、ε3为控制参数;
其中,为k24的导数,k24b与k24a为正数;
其中,为k34的导数,k34b与k34a为正数;
其中,为k44的导数,k44b与k44a为正数。
步骤二,根据所述观测系统对高超声速飞行器的攻角进行观测。
进一步的,高超声速飞行器推力升力与力矩的计算方法可以包括如下:
首先,高超声速飞行器的升力估计值的计算方法可以包括:
其中,S为飞行器特征面积;ρ为大气密度;为高超声速飞行器的速度估计值。
其次,高超声速飞行器的推力估计值的计算方法可以包括:
其中:S为飞行器特征面积;ρ为大气密度;为高超声速飞行器的速度估计值;
β为发动机的油门供油因子。
最后,高超声速飞行器的俯仰力矩估计值的计算方法可以包括:
其中,为高超声速飞行器的特征因子;δ为高超声速飞行器的俯仰舵偏角;
S为飞行器特征面积;ρ为大气密度;
为姿态角速度观测值;为攻角观测值;
为高超声速飞行器的速度估计值。
步骤三:高超声速飞行器俯仰通道计算机仿真模型的建立
为了确保上述步骤一以及步骤二中控制器的参数选取合理,可用通过计算机模拟仿真的手段进行编程,从而模拟被控对象高超声速飞行器俯仰通道的运动特性,从而方便进行观测器参数调整。在此以与前文观测器中力与力矩估算相匹配的一类高超声速飞行器模型为例示范说明,其俯仰通道可以采用如下微分方程建模表示:
其中,γ为飞行器的弹道倾角,为弹道倾角的导数,为飞行器姿态角速率的导数,为飞行器攻角的导数,Iyy为飞行器的转动惯量;
V为高超声速飞行器速度;
L、T、My(α)分别为升力、推力与俯仰力矩,其它未解释参量见前文说明。α采取PID控制规律对期望攻角进行跟踪,详见案例实施。
步骤四:观测器控制参数调试
首先将上述模型中的常参数进行数值设定,必须符合飞行器模型的物理意义,详细设定参见后文案例实施。
其次,进行观测器15个控制参数的调试。15个控制参数如下:
k22、k23、ε1、k32、k33、ε2、k42、k43、ε3、k24b、k24a、k34b、k34a、k44b、k44a。
调试分四步,第一步对k22、k33参数进行调试,其选取为较大的正数,具体选值见案例实施。
第二步进行ε1、ε2、ε3参数调试,其选取为较小的正数,一般小于1,具体选值见案例实施。
第三步进行k24b、k24a、k34b、k34a、k44b、k44a,上述参数可选取为较大的正数,具体选值见案例实施。
第四步进行k23、k32、k42、k43参数的调试,该参数正负不定,具体选值见案例实施。
最后,通过上述调试好的参数,进行仿真分析,根据攻角观测曲线与原系统的吻合程度,进行进一步的调整,值得得到较好的观测效果,从而确定观测器参数。具体仿真曲线见案例实施的结果。
通过上述四步,即实现了本发明所提供的非线性攻角观测方法。
本发明提供了一种测量弹道倾角与姿态角速率的方法,对高超声速飞行器构造观测系统,尤其是提供了一种变增益非线性的误差调节规律,能够有效地改善攻角观测的快速性与准确性。因此,本发明所提供的方法具有很高的工程应用价值与经济价值。
在步骤一中,设定超声速飞行器的速度估计值
高超声速飞行器的转动惯量估计值
高超声速飞行器质量的估计值
在步骤二中,大气密度ρ选取为ρ=0.0125。
发动机的油门供油因子β选取为β=0.8。
为高超声速飞行器的俯仰力矩估计值。
为高超声速飞行器的特征因子,选取为
S为高超声速飞行器的特征因子的特征面积,选取为S=334.7。
在步骤三中,设置飞行器的俯仰舵偏角δ的PID控制律构造如下:
从而实现对期望攻角4度的跟踪。
高超声速飞行器的速度可选与估计值相同,V=4602.5。
高超声速飞行器的转动惯量可选与估计值相同Iyy=6.77*105。
高超声速飞行器质量可选与估计值相同,m=4352.3。
步骤四中设定控制参数如下:
k22=300、k33=500、ε1=0.5、ε2=0.5、ε3=0.5、k24b=2、k24a=30、k34b=2、k34a=40、k44b=2、k44a=50,k23=0.2、k32=0.1、k42=0.066、k43=4.544。具体仿真曲线见图2至图13。
由仿真曲线可知,在4s后,攻角观测值即能较好地跟踪攻角真实值,而弹道倾角与姿态角速度也在一定的延迟之后,能够跟踪真实值。而且弹道倾角的估计值与姿态角速度的估计值,本可以直接采用测量值代替,但此时则完全不存在弹道倾角与姿态角速率的估计误差。该误差的存在与动态响应并消失主要是为了提高攻角观测的响应效果。
可见,采用本发明所提供的方法,信号恢复误差在可接受的范围内。能够通过阀值判断来实现。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
本说明书中的各个实施例均采用相关的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。
以上所述仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等,均包含在本发明的保护范围内。
Claims (4)
1.一种基于非线性增益的高超声速飞行器攻角观测方法,其特征在于,按照以下步骤进行:
步骤一,测量弹道倾角以及姿态角速度,并根据测量到的弹道倾角以及姿态角速度构建一观测系统;
步骤二,根据所述观测系统观测高超声速飞行器的攻角。
2.根据权利要求1所述的高超声速飞行器攻角观测方法,其特征在于,所述步骤一包括:
首先,采用惯导系统测量高超声速飞行器的弹道倾角γ以及姿态角速度q;
其次,根据弹道倾角γ以及姿态角速度q构建观测如下观测系统:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mover>
<mi>&gamma;</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mover>
<mi>m</mi>
<mo>^</mo>
</mover>
<mover>
<mi>V</mi>
<mo>^</mo>
</mover>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mover>
<mi>L</mi>
<mo>^</mo>
</mover>
<mo>+</mo>
<mover>
<mi>T</mi>
<mo>^</mo>
</mover>
<mi>sin</mi>
<mover>
<mi>&alpha;</mi>
<mo>^</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mover>
<mi>V</mi>
<mo>^</mo>
</mover>
</mfrac>
<mi>g</mi>
<mi> </mi>
<mi>cos</mi>
<mover>
<mi>&gamma;</mi>
<mo>^</mo>
</mover>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>22</mn>
</msub>
<mover>
<mi>&gamma;</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>23</mn>
</msub>
<mover>
<mi>q</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>24</mn>
</msub>
<mfrac>
<msup>
<mover>
<mi>&gamma;</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>3</mn>
</mrow>
</msup>
<mrow>
<msup>
<mover>
<mi>&gamma;</mi>
<mo>~</mo>
</mover>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mover>
<mi>q</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mover>
<mi>M</mi>
<mo>^</mo>
</mover>
<mi>y</mi>
</msub>
<mrow>
<mo>(</mo>
<mover>
<mi>&alpha;</mi>
<mo>^</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mover>
<mi>I</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msub>
</mfrac>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>32</mn>
</msub>
<mover>
<mi>&gamma;</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>33</mn>
</msub>
<mover>
<mi>q</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>34</mn>
</msub>
<mfrac>
<msup>
<mover>
<mi>q</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>3</mn>
</mrow>
</msup>
<mrow>
<msup>
<mover>
<mi>q</mi>
<mo>~</mo>
</mover>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mover>
<mi>&alpha;</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mover>
<mi>q</mi>
<mo>^</mo>
</mover>
<mo>-</mo>
<mover>
<mover>
<mi>&gamma;</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>42</mn>
</msub>
<mover>
<mi>&gamma;</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>43</mn>
</msub>
<mover>
<mi>q</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>44</mn>
</msub>
<mfrac>
<msup>
<mover>
<mi>&alpha;</mi>
<mo>~</mo>
</mover>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>3</mn>
</mrow>
</msup>
<mrow>
<msup>
<mover>
<mi>&alpha;</mi>
<mo>~</mo>
</mover>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<mn>3</mn>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
其中,α为高超声速飞行器的攻角;为弹道倾角观测值,为姿态角速度观测值,为攻角观测值;
为的导数,为的导数,为的导数;
为估计误差;并且有:
为高超声速飞行器质量的估计值;g为重力加速度,取值为9.8;
为高超声速飞行器的速度估计值;为高超声速飞行器的升力估计值;
为高超声速飞行器的推力估计值;为高超声速飞行器的俯仰力矩估计值;
为高超声速飞行器的转动惯量估计值;
k22、k23、ε1、k32、k33、ε2、k42、k43、ε3为控制参数;
其中,为k24的导数,k24b与k24a为正数;
其中,为k34的导数,k34b与k34a为正数;
其中,为k44的导数,k44b与k44a为正数。
3.根据权利要求2所述的高超声速飞行器攻角观测方法,其特征在于,高超声速飞行器推力升力与力矩的计算方法包括:
首先,高超声速飞行器的升力估计值的计算方法包括:
其中,S为飞行器特征面积;ρ为大气密度;为高超声速飞行器的速度估计值;
其次,高超声速飞行器的推力估计值的计算方法包括:
其中:S为飞行器特征面积;ρ为大气密度;为高超声速飞行器的速度估计值;
β为发动机的油门供油因子;
最后,高超声速飞行器的俯仰力矩估计值的计算方法包括:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>M</mi>
<mo>^</mo>
</mover>
<mi>y</mi>
</msub>
<mrow>
<mo>(</mo>
<mover>
<mi>&alpha;</mi>
<mo>^</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mover>
<mi>q</mi>
<mo>&OverBar;</mo>
</mover>
<mi>S</mi>
<mover>
<mi>c</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&lsqb;</mo>
<msub>
<mover>
<mi>C</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>&alpha;</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>C</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>&delta;</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>C</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>C</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>&alpha;</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mn>0.035</mn>
<msup>
<mover>
<mi>&alpha;</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>0.036617</mn>
<mover>
<mi>&alpha;</mi>
<mo>^</mo>
</mover>
<mo>+</mo>
<mn>5.326</mn>
<mo>*</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>6</mn>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>C</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mover>
<mi>c</mi>
<mo>&OverBar;</mo>
</mover>
<mover>
<mi>q</mi>
<mo>^</mo>
</mover>
</mrow>
<mrow>
<mn>2</mn>
<mover>
<mi>V</mi>
<mo>^</mo>
</mover>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>6.79</mn>
<msup>
<mover>
<mi>&alpha;</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>0.3015</mn>
<mover>
<mi>&alpha;</mi>
<mo>^</mo>
</mover>
<mo>-</mo>
<mn>0.2289</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>C</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>M</mi>
<mi>&delta;</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0.0292</mn>
<mrow>
<mo>(</mo>
<mi>&delta;</mi>
<mo>-</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
其中,为高超声速飞行器的特征因子;δ为高超声速飞行器的俯仰舵偏角;
S为飞行器特征面积;ρ为大气密度;
为姿态角速度观测值;为攻角观测值;为高超声速飞行器的速度估计值。
4.根据权利要求2所述的高超声速飞行器攻角观测方法,其特征在于,各所述控制参数的取值为:
k22=300、k33=500、ε1=0.5、ε2=0.5、ε3=0.5、k23=0.2、k32=0.1、k42=0.066、k43=4.544。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710743894.7A CN107367941B (zh) | 2017-08-25 | 2017-08-25 | 高超声速飞行器攻角观测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710743894.7A CN107367941B (zh) | 2017-08-25 | 2017-08-25 | 高超声速飞行器攻角观测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107367941A true CN107367941A (zh) | 2017-11-21 |
CN107367941B CN107367941B (zh) | 2020-05-19 |
Family
ID=60312136
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710743894.7A Expired - Fee Related CN107367941B (zh) | 2017-08-25 | 2017-08-25 | 高超声速飞行器攻角观测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107367941B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108490788A (zh) * | 2018-05-08 | 2018-09-04 | 中国人民解放军海军航空大学 | 一种基于双干扰观测的飞行器俯仰通道反演控制方法 |
CN110672060A (zh) * | 2019-08-23 | 2020-01-10 | 中国人民解放军63729部队 | 一种基于外测速度的箭上姿态角情况判决方法 |
CN110889170A (zh) * | 2019-11-29 | 2020-03-17 | 中国人民解放军海军工程大学 | 飞行器大角度攻击目标的落角与攻角的估计方法 |
CN111122901A (zh) * | 2019-12-31 | 2020-05-08 | 中国航空工业集团公司西安飞机设计研究所 | 一种攻角误差自动检测方法及其修正系统 |
CN111427267A (zh) * | 2020-04-01 | 2020-07-17 | 山东创惠电子科技有限责任公司 | 一种采用力与力矩自适应估计的高速飞行器攻角跟踪方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1746523A2 (en) * | 2005-07-22 | 2007-01-24 | Lockheed Martin Corporation | Method of designing a hypersonic inlet for an aircraft and system therefor |
CN102591212A (zh) * | 2012-03-01 | 2012-07-18 | 北京航空航天大学 | 一种时变测量延迟输出信号飞行器纵向运动状态观测方法 |
CN103592847A (zh) * | 2013-10-30 | 2014-02-19 | 天津大学 | 一种基于高增益观测器的高超声速飞行器非线性控制方法 |
CN104199286A (zh) * | 2014-07-15 | 2014-12-10 | 北京航空航天大学 | 一种基于滑模干扰观测器的飞行器递阶动态逆控制方法 |
CN106406102A (zh) * | 2016-12-20 | 2017-02-15 | 哈尔滨工业大学 | 一种含干扰观测器的高超声速飞行器跟踪控制方法 |
-
2017
- 2017-08-25 CN CN201710743894.7A patent/CN107367941B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1746523A2 (en) * | 2005-07-22 | 2007-01-24 | Lockheed Martin Corporation | Method of designing a hypersonic inlet for an aircraft and system therefor |
CN102591212A (zh) * | 2012-03-01 | 2012-07-18 | 北京航空航天大学 | 一种时变测量延迟输出信号飞行器纵向运动状态观测方法 |
CN103592847A (zh) * | 2013-10-30 | 2014-02-19 | 天津大学 | 一种基于高增益观测器的高超声速飞行器非线性控制方法 |
CN104199286A (zh) * | 2014-07-15 | 2014-12-10 | 北京航空航天大学 | 一种基于滑模干扰观测器的飞行器递阶动态逆控制方法 |
CN106406102A (zh) * | 2016-12-20 | 2017-02-15 | 哈尔滨工业大学 | 一种含干扰观测器的高超声速飞行器跟踪控制方法 |
Non-Patent Citations (3)
Title |
---|
CHENGBIN LIAN,等: "The Design of Nonlinear Flight Control System of Hypersonic Cruise Vehicle Based on Nonlinear Disturbance Observer", 《2013 THIRD INTERNATIONAL CONFERENCE ON INSTRUMENTATION, MEASUREMENT, COMPUTER, COMMUNICATION AND CONTROL》 * |
于进勇: "一类高超声速飞行器自适应容错控制器设计", 《中国科学:信息科学》 * |
卜祥伟: "基于非线性干扰观测器的高超声速飞行器滑模反演控制", 《控制理论与应用》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108490788A (zh) * | 2018-05-08 | 2018-09-04 | 中国人民解放军海军航空大学 | 一种基于双干扰观测的飞行器俯仰通道反演控制方法 |
CN108490788B (zh) * | 2018-05-08 | 2021-06-01 | 中国人民解放军海军航空大学 | 一种基于双干扰观测的飞行器俯仰通道反演控制方法 |
CN110672060A (zh) * | 2019-08-23 | 2020-01-10 | 中国人民解放军63729部队 | 一种基于外测速度的箭上姿态角情况判决方法 |
CN110672060B (zh) * | 2019-08-23 | 2022-09-16 | 中国人民解放军63729部队 | 一种基于外测速度的箭上姿态角情况判决方法 |
CN110889170A (zh) * | 2019-11-29 | 2020-03-17 | 中国人民解放军海军工程大学 | 飞行器大角度攻击目标的落角与攻角的估计方法 |
CN110889170B (zh) * | 2019-11-29 | 2022-09-30 | 中国人民解放军海军工程大学 | 飞行器大角度攻击目标的落角与攻角的估计方法 |
CN111122901A (zh) * | 2019-12-31 | 2020-05-08 | 中国航空工业集团公司西安飞机设计研究所 | 一种攻角误差自动检测方法及其修正系统 |
CN111427267A (zh) * | 2020-04-01 | 2020-07-17 | 山东创惠电子科技有限责任公司 | 一种采用力与力矩自适应估计的高速飞行器攻角跟踪方法 |
CN111427267B (zh) * | 2020-04-01 | 2022-08-30 | 山东创惠电子科技有限责任公司 | 一种采用力与力矩自适应估计的高速飞行器攻角跟踪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107367941B (zh) | 2020-05-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107367941A (zh) | 基于非线性增益的高超声速飞行器攻角观测方法 | |
CN104777844B (zh) | 一种高超声速临近空间飞行器航迹跟踪方法 | |
CN102540882B (zh) | 一种基于最小参数学习法的飞行器航迹倾角控制方法 | |
CN106896722B (zh) | 采用状态反馈与神经网络的高超飞行器复合控制方法 | |
CN109635494A (zh) | 一种飞行试验与地面仿真气动力数据综合建模方法 | |
CN107844128B (zh) | 一种基于复合比例导引的高超声速飞行器巡航段制导方法 | |
CN106527122A (zh) | 平流层飞艇定高飞行非线性pid控制方法 | |
CN105005099B (zh) | 一种基于捷联惯导与飞行控制系统的大气参数解算方法 | |
CN103744057A (zh) | 基于输出相关自适应卡尔曼滤波的弹道轨迹形成方法 | |
CN110187713A (zh) | 一种基于气动参数在线辨识的高超声速飞行器纵向控制方法 | |
CN105242679A (zh) | 一种四旋翼飞行器的控制系统设计方法 | |
CN103837321B (zh) | 一种水面飞行器实机稳定性试验方法 | |
Box et al. | Stochastic six-degree-of-freedom flight simulator for passively controlled high-power rockets | |
CN105180728A (zh) | 基于前数据的旋转制导炮弹快速空中对准方法 | |
CN105136423A (zh) | 考虑摩擦力的自由振动动导数试验的数据分析方法 | |
CN108318219B (zh) | 自由飞试验相似律重力补偿设计方法 | |
CN103744058A (zh) | 基于指数加权衰减记忆滤波的弹道轨迹形成方法 | |
CN109190248A (zh) | 一种用于滑翔飞行器的滑翔射程解析方法及解析系统 | |
CN103983274B (zh) | 一种适用于低精度无方位基准双轴转位设备的惯性测量单元标定方法 | |
CN106774385B (zh) | 一种采用自适应变结构的飞艇定点悬停控制方法 | |
CN106484969B (zh) | 一种大包线强机动飞行器动力学高精度仿真方法 | |
CN107894240A (zh) | 一种用于水下无人航行器在极区导航的初始粗对准方法 | |
Li et al. | Simulation method for wind tunnel based virtual flight testing | |
CN103994775B (zh) | 一种适用于低精度有方位基准双轴转位设备的惯性测量单元标定方法 | |
CN112947522A (zh) | 一种基于有限时间观测器的硬式空中加油姿态控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200519 Termination date: 20210825 |
|
CF01 | Termination of patent right due to non-payment of annual fee |