CN102591212B - 一种时变测量延迟输出信号飞行器纵向运动状态观测方法 - Google Patents
一种时变测量延迟输出信号飞行器纵向运动状态观测方法 Download PDFInfo
- Publication number
- CN102591212B CN102591212B CN 201210051956 CN201210051956A CN102591212B CN 102591212 B CN102591212 B CN 102591212B CN 201210051956 CN201210051956 CN 201210051956 CN 201210051956 A CN201210051956 A CN 201210051956A CN 102591212 B CN102591212 B CN 102591212B
- Authority
- CN
- China
- Prior art keywords
- observer
- aircraft
- observation
- parameter
- overbar
- 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.)
- Expired - Fee Related
Links
Images
Abstract
本发明一种时变测量延迟输出信号飞行器纵向运动状态观测方法,它有四大步骤:步骤一、获取飞行器在俯仰平面内的动力学模型;步骤二、设计观测器;步骤三、调节观测参数,查看观测效果;步骤四、设计结束。该观测方法根据包含有一定时变延迟后的输出和系统的模型,来重构出系统的当前时刻的状态,从而实现对系统状态的观测。它不仅可以实现对飞行器纵向运动动力学模型的观测,还可以实现对其他控制对象的观测,因此,该发明方法在飞行控制技术领域里具有较好的实用价值和良好的应用前景。
Description
(一)技术领域
本发明涉及一种时变测量延迟输出信号飞行器纵向运动状态观测方法,它是针对包含时变测量延迟输出的飞行器纵向动力学模型,设计了一种状态观测的方法,属于飞行控制技术领域。
(二)背景技术
对飞行器的研究一直以来都是研究者所共同关注的一个热点问题,而飞行器控制系统是保证飞行器在空间中稳定飞行的一个必不可少的环节。在飞行器控制系统设计的过程中需要考虑两个方面的因素:稳定性和快速性。稳定性是保证飞行器能够在空间中稳定地飞行,而快速性是要求飞行器能够在高速飞行的过程中对于一些突变情况迅速做出正确的反应。
飞行器在空间运动的过程中,需要对其姿态角(航迹倾角、攻角、俯仰角等)进行控制,而控制飞行器姿态角是通过计算控制律来实现的。在计算控制律之前需要提前获取系统的状态。因此,获取飞行器在当前时刻的飞行状态是实现对飞行器控制的一个重要环节。
但是对于通常的飞行器而言,由于其定位和测量状态的仪器都是通过GPS或者一些惯性的测量单元来完成的。而在从卫星到飞行器的传输测量信号的过程中,不可避免的会存在一定的延迟,这会给之后的控制律的设计带来非常的不便。而由于测量延迟在实际传输过程中,并不一定是常值,其随时间变化的特性,给观测器的设计带来了进一步的困难。
(三)发明内容
1.发明目的
针对以上所提到的问题,本发明提供了一种时变测量延迟输出信号飞行器纵向运动状态观测方法。该观测方法能够根据包含有一定时变延迟后的输出和系统的模型,来重构出系统的当前时刻的状态,从而实现对系统状态的观测。本发明不仅可以实现对飞行器纵向运动动力学模型的观测,还可以实现对其他控制对象的观测,因此,该发明方法具有一定的适用性。
2.技术方案
为达到上述目的,本发明针对包含时变测量延迟输出的飞行器纵向动力学模型,利用测得的带有延迟和噪声的航迹倾角信号,给出了一种状态观测的方法。
本发明一种时变测量延迟输出信号飞行器纵向运动状态观测方法,该方法包括以下步骤:
步骤1:获取飞行器在俯仰平面内的动力学模型
本发明所针对的观测对象是俯仰平面内运动的飞行器。首先获取纵向运动的动力学模型。飞行器在俯仰平面内的结构示意图如图1所示。在该平面中,飞行器受到来自发动机的推力、空气施加的升力和阻力以及俯仰力矩。通过对所受到的力进行分析,可以得到飞行器在纵向运动中的动力学方程,对其进行线性化处理后,可以得到如下所示的动力学模型表达式:
并且有
式中,γ、α、ψ分别表示飞行器航迹倾角、攻角和俯仰角且有γ=ψ-α;q为俯仰角变化率;VT为航速;m和g分别为飞行器质量和重力加速度;Lα表示升力曲线斜率;Lo表示其他对升力的影响因素;Mδ表示控制俯仰力矩;Mo表示其他来源力矩,通常由公式Mo=Mαα+Mqq近似,δ为舵面偏角。
定义x1=γ,x2=ψ,x3=q,可将飞行器纵向模型写成状态空间的表达形式:
此外,由于该系统模型中的输出y=x1(t),包含有一定的延迟环节,可将时变延时输出表达为:
其中δ(t)为随时间变化的延迟时间,满足|δ(t)|≤1.5,其变化过程仿真图如图6所示。
步骤2:设计观测器
根据上个步骤中获得的动力学模型,对其进行分析,便可设计出满足要求的观测器结构。在保证观测误差收敛的同时,还需要保证观测速度能快速达到设计要求。该步骤是本发明方法中最重要的一个步骤。
观测器结构中包含了模型估计环节和增益调节环节,其数学表达式如下所示:
其中G(·)为设计的观测器调节增益环节,其表达式为:
其中
而K=[k1k2k3]T为欲配置极点的配置参数矩阵。若期望配置的极点位置为λ1=-e1,λ2=-e2,λ3=-e3,ei>0,则矩阵K的计算公式为:
对上述提出的观测器进行Lyapunov-Razumikhin收敛性分析。首先根据式(2)和式(3),整理出观测误差ξ状态方程的表达式,然后定义Lyapunov函数V(ξ),利用Razumikhin理论,对其进行收敛性分析。采用观测器式(3),可以证明出其中ω>0,即可保证所设计的观测器的收敛性,同时也说明了所提出的观测方法的有效性。但是由于所设计的观测器中包含有很多的参数,为了使观测效果达到最好,即观测器的响应速度最快和观测误差超调量最小,需要对观测器中的参数进行调整。
步骤3:调节观测参数,查看观测效果
按照上一步中所设计的观测器结构,对观测器中的某些参数进行调节,并查看观测效果是否满足要求,如果不满足则继续调节参数,一直到符合要求为止。
a1=1,a2=-0.1,a3=g/Vt
b1=1,b2=0.1,b3=-0.02
以上参数均是飞行器本身所固有的参数,不需要做调整。而在仿真过程中,主要调节的参数是欲配置的极点坐标位置。经过多调节后,最终所选择欲配置的极点坐标的绝对值大小为:
此外,在仿真中还加入噪声环节,所加入的噪声的幅值为0.001。
步骤4:设计结束
整个设计过程重点考虑了两个方面的要求,一个是所设计的观测器的收敛性问题;另外一个是要保证观测效果能够满足发明要求。第一步中,需要确定飞行器在俯仰平面内的动力学模型;第二步中,需要设计出满足要求的观测器;第三步,是对之前所设计出的观测器进行参数调整。
3.优点及功效
本发明方法针对输出中存在时变延迟的飞行器纵向运动模型,利用测得的带有延迟和噪声的航迹倾角信号,实现对状态变量的估计和观测,解决了由于测量、信号传输过程或设备等因素造成的信号延时问题,从而有利于控制系统的设计。此外,如果在输出量中存在一定的测量噪声,采用本方法,能从含有噪声的延迟输出信号中提取出待观测信号的实际值,具有一定的抗噪能力。
(四)附图说明
图1本发明中飞行器纵向模型示意图
图2本发明观测系统仿真示意图
图3.1本发明航迹倾角跟踪效果仿真示意图
图3.2本发明俯仰角跟踪效果仿真示意图
图3.3本发明俯仰角速率跟踪效果仿真示意图
图4本发明加入噪声的延迟输出仿真示意图
图5本发明延迟输出信号与观测信号仿真示意图
图6本发明时变延迟时间的仿真图
图7本发明实施流程框图
图中符号说明如下:
γ、α、ψ分别表示飞行器航迹倾角、攻角和俯仰角且有γ=ψ-α;VT为航速;δ为舵面偏角。
(五)具体实施方式
下面将结合附图和实施实例对本发明做进一步的详细说明。
见图7,本发明一种时变测量延迟输出信号飞行器纵向运动状态观测方法,该方法包括以下步骤:
步骤1:获取飞行器在俯仰平面内的动力学模型
本发明所针对的观测对象是,在俯仰平面内运动的飞行器。飞行器在俯仰平面内的结构示意图如图1所示,在该平面中,飞行器受到来自发动机的推力、空气施加的升力和阻力以及俯仰力矩。通过对所受到的力进行分析,可以得到飞行器在纵向运动中的动力学方程,对其进行线性化处理后,可以得到如下所示的动力学模型表达式:
并且有
式中,γ、α、ψ分别表示飞行器航迹倾角、攻角和俯仰角且有γ=ψ-α;q为俯仰角变化率;VT为航速;m和g分别为飞行器质量和重力加速度;Lα表示升力曲线斜率;Lo表示其他对升力的影响因素;Mδ表示控制俯仰力矩;Mo表示其他来源力矩,通常由公式Mo=Mαα+Mqq近似,δ为舵面偏角。
定义x1=γ,x2=ψ,x3=q,可将飞行器纵向模型写成状态空间的表达形式:
此外,由于该系统模型中的输出y=x1(t),包含有一定的延迟环节,可将时变延时输出表达为:
其中δ(t)为随时间变化的延迟时间,满足|δ(t)|≤1.5,其变化过程仿真图如图6所示。而图2为本发明观测系统仿真示意图。
步骤2:设计观测器
根据上个步骤中获得的动力学模型,对其进行分析,可以设计出满足要求的观测器结构。观测器结构中包含了模型估计环节和增益调节环节,其数学表达式如下所示:
其中G(·)为设计的观测器调节增益环节,其表达式为:
其中
而K=[k1k2k3]T为欲配置极点的配置参数矩阵。若期望配置的极点位置为λ1=-e1,λ2=-e2,λ3=-e3,ei>0,则矩阵K的计算公式为:
对上述提出的观测器进行Lyapunov-Razumikhin收敛性分析。首先根据式(2)和式(3),整理出观测误差ξ状态方程的表达式,然后定义Lyapunov函数V(ξ),利用Razumikhin理论,对其进行收敛性分析。采用观测器式(3),可以证明出其中ω>0,即可保证所设计的观测器的收敛性,同时也说明了所提出的观测方法的有效性。但是由于所设计的观测器中包含有很多的参数,为了使观测效果达到最好,即观测器的响应速度最快和观测误差超调量最小,需要对观测器中的参数进行调整。
步骤3:调节观测参数,查看观测效果
按照上一步中所设计的观测器,对观测器中的某些参数进行调节,并查看观测效果是否满足要求,主要查看观测器的响应速度和观测误差的超调量。如果不满足则继续调节参数,直到控制效果达到要求;如果控制效果能够满足要求,则保留该组参数,一直到设计结束。在仿真中,选择飞行器系统的参数为: Mα=0.1,Mq=-0.02,Mδ=1.0.设定稳定航速为VT=200m/s,g取9.8m/s2,即:
a1=1,a2=-0.1,a3=g/Vt
b1=1,b2=0.1,b3=-0.02
以上参数均是飞行器本身所固有的参数,不需要做调整。而在仿真过程中,主要调节的参数是欲配置的极点坐标位置。经过多调节后,最终所选择欲配置的极点坐标的绝对值大小为:
此外,在仿真中还加入噪声环节,所加入的噪声的幅值为0.001,加入噪声之后的测量输出信号如图4所示。仿真最终得到的结果图,如图3.1、图3.2、图3.3、图4、图5、图6所示。
步骤4:设计结束
通过以上几个步骤设计之后,可以得到满意的观测效果,也意味着观测方法的有效性。通过仿真,可以进一步验证所设计的观测方法的正确性。由仿真结果图3.1、图3.2、图3.3可以看出,所设计的观测器能够很好的实现对飞行器三个状态量(飞行器航迹倾角和俯仰角以及俯仰角变化率)的观测,可以根据延迟后的输出来估计出系统在当前时刻的状态;图4展示了噪声和时变延迟对系统输出的影响;图5对比了系统实际输出和通过观测器观测之后的输出之间的关系;而图6给出了系统时变延迟环节的变化过程。
综上所述,对于包含有时变延时输出的飞行器纵向运动模型而言,使用本发明提出的观测状态的方法,能够快速实现对其状态的高精度观测。此外,由于在测量输出的过程中不可避免的会引入测量噪声,而本方法也考虑了含有噪声时的结果,并且得到了很好的仿真结果,进一步说明了本发明方法的有效性和抗噪性。
Claims (1)
1.一种时变测量延迟输出信号飞行器纵向运动状态观测方法,其特征在于:该方法包括以下步骤:
步骤1:获取飞行器在俯仰平面内的动力学模型
对飞行器在俯仰平面内的受力情况分析,飞行器受到来自发动机的推力、空气施加的升力和阻力以及俯仰力矩;从而得到飞行器在纵向运动中的动力学方程,对其进行线性化处理后,得到如下所示的动力学模型表达式:
并且有
式中,γ、α、ψ分别表示飞行器航迹倾角、攻角和俯仰角且有γ=ψ-α;q为俯仰角变化率;VT为航速;m和g分别为飞行器质量和重力加速度;Lα表示升力曲线斜率;Lo表示其他对升力的影响因素;Mδ表示控制俯仰力矩;Mo表示其他来源力矩,通常由公式Mo=Mαα+Mqq近似,δ为舵面偏角;
定义x1=γ,x2=ψ,x3=q,将飞行器纵向模型写成状态空间的表达形式:
此外,由于该系统模型中的输出y=x1(t),包含有一定的延迟环节,将时变延时输出表达为:
其中δ(t)为随时间变化的延迟时间,满足|δ(t)|≤1.5;
步骤2:设计观测器
根据上个步骤中获得的动力学模型,对其进行分析,便设计出满足要求的观测器结构;在保证观测误差收敛的同时,还需要保证观测速度能快速达到设计要求;观测器结构中包含了模型估计环节和增益调节环节,其数学表达式如下所示:
其中G(.)为设计的观测器调节增益环节,其表达式为:
其中,G(1)、G(2)、G(3)分别为 的第1、2、3项元素;
其中,
而K=[k1 k2 k3]T为欲配置极点的配置参数矩阵;若期望配置的极点位置为λ1=-e1,λ2=-e2,λ3=-e3,ei>0,则矩阵K的计算公式为:
对上述提出的观测器进行Lyapunov-Razumikhin收敛性分析;首先根据式(2)和式(3),整理出观测误差ξ状态方程的表达式,然后定义Lyapunov函数V(ξ),利用Razumikhin理论,对其进行收敛性分析;采用观测器式(3),证明出其中ω>0,即可保证所设计的观测器的收敛性,同时也说明了所提出的观测方法的有效性;但是由于所设计的观测器中包含有很多的参数,为了使观测效果达到最好,即观测器的响应速度最快和观测误差超调量最小,需要对观测器中的参数进行调整;
步骤3:调节观测参数,查看观测效果
按照上一步中所设计的观测器结构,对观测器中的参数进行调节,并查看观测效果是否满足要求,如果不满足则继续调节参数,一直到符合要求为止;
a1=1,a2=-0.1,a3=g/Vt
b1=1,b2=0.1,b3=-0.02
以上参数均是飞行器本身所固有的参数,不需要做调整;而在仿真过程中,调节的参数是欲配置的极点坐标位置,经过多调节后,最终所选择欲配置的极点坐标的绝对值大小为:
此外,在仿真中还加入噪声环节,所加入的噪声的幅值为0.001;
步骤4:设计结束
整个设计过程重点考虑了两个方面的要求,一个是所设计的观测器的收敛性问题;另外一个是要保证观测效果能够满足发明要求;第一步中,需要确定飞行器在俯仰平面内的动力学模型;第二步中,需要设计出满足要求的观测器;第三步,是对之前所设计出的观测器进行参数调整。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201210051956 CN102591212B (zh) | 2012-03-01 | 2012-03-01 | 一种时变测量延迟输出信号飞行器纵向运动状态观测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201210051956 CN102591212B (zh) | 2012-03-01 | 2012-03-01 | 一种时变测量延迟输出信号飞行器纵向运动状态观测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102591212A CN102591212A (zh) | 2012-07-18 |
CN102591212B true CN102591212B (zh) | 2013-09-18 |
Family
ID=46480047
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201210051956 Expired - Fee Related CN102591212B (zh) | 2012-03-01 | 2012-03-01 | 一种时变测量延迟输出信号飞行器纵向运动状态观测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102591212B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102915036B (zh) * | 2012-07-26 | 2014-12-03 | 北京航空航天大学 | 一种参数不确定性飞行器倾斜角控制系统极限环抑制方法 |
CN102880182B (zh) * | 2012-09-12 | 2015-01-14 | 北京航空航天大学 | 一种存在网络随机延迟的微小型无人飞行器控制方法 |
CN102929143B (zh) * | 2012-10-10 | 2015-05-06 | 西北工业大学 | 飞行器时间滞后模型的控制器设计方法 |
CN102929135B (zh) * | 2012-10-10 | 2015-05-27 | 西北工业大学 | 飞行器时间滞后时变模型的控制器设计方法 |
CN102929144B (zh) * | 2012-10-10 | 2015-07-29 | 西北工业大学 | 飞行器区间模型的控制器设计方法 |
CN102929142B (zh) * | 2012-10-10 | 2015-04-01 | 西北工业大学 | 飞行器不确定时变模型的控制器设计方法 |
CN104038132B (zh) * | 2014-06-12 | 2016-08-24 | 国电科学技术研究院 | 具有时变测量延迟输出和噪声的伺服电机的状态观测方法 |
CN104808493B (zh) * | 2015-04-21 | 2017-08-25 | 国电科学技术研究院 | 一种基于延时观测器的汽轮发电机主汽门开度预测控制方法 |
CN105093932B (zh) * | 2015-06-23 | 2016-06-01 | 北京航天自动控制研究所 | 一种确定lpv变增益控制器的鲁棒性的方法 |
CN104991566B (zh) * | 2015-07-07 | 2016-06-08 | 北京航天自动控制研究所 | 一种用于高超声速飞行器的参数不确定性lpv系统建模方法 |
CN107367941B (zh) * | 2017-08-25 | 2020-05-19 | 中国人民解放军海军航空工程学院 | 高超声速飞行器攻角观测方法 |
CN109062237B (zh) * | 2018-09-17 | 2021-07-20 | 南京航空航天大学 | 一种无人倾转旋翼机自抗扰姿态控制方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6298318B1 (en) * | 1998-07-01 | 2001-10-02 | Ching-Fang Lin | Real-time IMU signal emulation method for test of Guidance Navigation and Control systems |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2010325107B2 (en) * | 2009-09-09 | 2016-08-25 | Aerovironment, Inc. | Systems and devices for remotely operated unmanned aerial vehicle report-suppressing launcher with portable RF transparent launch tube |
-
2012
- 2012-03-01 CN CN 201210051956 patent/CN102591212B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6298318B1 (en) * | 1998-07-01 | 2001-10-02 | Ching-Fang Lin | Real-time IMU signal emulation method for test of Guidance Navigation and Control systems |
Non-Patent Citations (8)
Title |
---|
刘强等.高速声速飞行器的滑模观测器设计.《宇航学报》.2004,第25卷(第6期),第588-592页. |
基于观测器的飞行器时变控制效益估计算法;池沛等;《控制与决策》;20081031;第23卷(第10期);第1186-1191页 * |
孟中杰等.弹性高超声速飞行器建模及精细姿态控制.《宇航学报》.2011,第32卷(第8期),第1683-1687页. |
弹性高超声速飞行器建模及精细姿态控制;孟中杰等;《宇航学报》;20110831;第32卷(第8期);第1683-1687页 * |
池沛等.基于观测器的飞行器时变控制效益估计算法.《控制与决策》.2008,第23卷(第10期),第1186-1191页. |
邓以高等.飞行器姿态控制方法综述.《战术导弹控制技术》.2006,(第2期),第7-12页. |
飞行器姿态控制方法综述;邓以高等;《战术导弹控制技术》;20061231(第2期);第7-12页 * |
高速声速飞行器的滑模观测器设计;刘强等;《宇航学报》;20041130;第25卷(第6期);第588-592页 * |
Also Published As
Publication number | Publication date |
---|---|
CN102591212A (zh) | 2012-07-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102591212B (zh) | 一种时变测量延迟输出信号飞行器纵向运动状态观测方法 | |
CN102749851B (zh) | 一种挠性高超声速飞行器的精细抗干扰跟踪控制器 | |
CN104443427B (zh) | 飞行器颤振预测系统及方法 | |
Ryan et al. | Current and future research in active control of lightweight, flexible structures using the X-56 aircraft | |
Wang et al. | Robust nonlinear flight control of a high-performance aircraft | |
CN102654772B (zh) | 一种基于控制力受限情况下飞行器航迹倾角反演控制方法 | |
CN108536020A (zh) | 一种针对垂直起降重复使用运载器的模型参考自适应滑模控制方法 | |
CN103324083B (zh) | 基于鲁棒观测器的非线性船舶运动控制方法 | |
CN104102128A (zh) | 一种适用于小型无人飞行器的抗干扰姿态控制方法 | |
CN103838145A (zh) | 基于级联观测器的垂直起降飞机鲁棒容错控制系统及方法 | |
EP2551737B1 (en) | Method and apparatus for minimizing dynamic structural loads of an aircraft | |
Hjartarson et al. | LPV aeroservoelastic control using the LPVTools toolbox | |
CN103587723A (zh) | 一种再入初始段解析式纵向在线轨迹设计及跟踪方法 | |
CN108153323B (zh) | 一种高空无人飞行器高精度再入制导方法 | |
CN103869701A (zh) | 基于姿态序列解算的飞行器新型实时制导方法 | |
CN104950901A (zh) | 无人直升机姿态误差有限时间收敛非线性鲁棒控制方法 | |
Williams | Three-dimensional aircraft terrain-following via real-time optimal control | |
CN102540882A (zh) | 一种基于最小参数学习法的飞行器航迹倾角控制方法 | |
CN106019937A (zh) | 一种混杂系统的抗干扰控制方法 | |
Suh et al. | Virtual Deformation Control of the X 56A Model with Simulated Fiber Optic Sensors | |
CN103921957B (zh) | 一种探月飞船跳跃式再入的跃起点能量管理方法 | |
Matt et al. | Frequency domain system identification of a small flying-wing uas | |
Lee et al. | Simulation of pilot workload for a helicopter operating in a turbulent ship airwake | |
Takase et al. | Aircraft gust alleviation preview control with a discrete-time LPV model | |
Grauer et al. | Aeroelastic modeling of X-56A stiff-wing configuration flight test data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130918 Termination date: 20150301 |
|
EXPY | Termination of patent right or utility model |