CN102591212B - 一种时变测量延迟输出信号飞行器纵向运动状态观测方法 - Google Patents

一种时变测量延迟输出信号飞行器纵向运动状态观测方法 Download PDF

Info

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
Application number
CN 201210051956
Other languages
English (en)
Other versions
CN102591212A (zh
Inventor
刘金琨
贺庆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN 201210051956 priority Critical patent/CN102591212B/zh
Publication of CN102591212A publication Critical patent/CN102591212A/zh
Application granted granted Critical
Publication of CN102591212B publication Critical patent/CN102591212B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明一种时变测量延迟输出信号飞行器纵向运动状态观测方法,它有四大步骤:步骤一、获取飞行器在俯仰平面内的动力学模型;步骤二、设计观测器;步骤三、调节观测参数,查看观测效果;步骤四、设计结束。该观测方法根据包含有一定时变延迟后的输出和系统的模型,来重构出系统的当前时刻的状态,从而实现对系统状态的观测。它不仅可以实现对飞行器纵向运动动力学模型的观测,还可以实现对其他控制对象的观测,因此,该发明方法在飞行控制技术领域里具有较好的实用价值和良好的应用前景。

Description

一种时变测量延迟输出信号飞行器纵向运动状态观测方法
(一)技术领域
本发明涉及一种时变测量延迟输出信号飞行器纵向运动状态观测方法,它是针对包含时变测量延迟输出的飞行器纵向动力学模型,设计了一种状态观测的方法,属于飞行控制技术领域。
(二)背景技术
对飞行器的研究一直以来都是研究者所共同关注的一个热点问题,而飞行器控制系统是保证飞行器在空间中稳定飞行的一个必不可少的环节。在飞行器控制系统设计的过程中需要考虑两个方面的因素:稳定性和快速性。稳定性是保证飞行器能够在空间中稳定地飞行,而快速性是要求飞行器能够在高速飞行的过程中对于一些突变情况迅速做出正确的反应。
飞行器在空间运动的过程中,需要对其姿态角(航迹倾角、攻角、俯仰角等)进行控制,而控制飞行器姿态角是通过计算控制律来实现的。在计算控制律之前需要提前获取系统的状态。因此,获取飞行器在当前时刻的飞行状态是实现对飞行器控制的一个重要环节。
但是对于通常的飞行器而言,由于其定位和测量状态的仪器都是通过GPS或者一些惯性的测量单元来完成的。而在从卫星到飞行器的传输测量信号的过程中,不可避免的会存在一定的延迟,这会给之后的控制律的设计带来非常的不便。而由于测量延迟在实际传输过程中,并不一定是常值,其随时间变化的特性,给观测器的设计带来了进一步的困难。
(三)发明内容
1.发明目的
针对以上所提到的问题,本发明提供了一种时变测量延迟输出信号飞行器纵向运动状态观测方法。该观测方法能够根据包含有一定时变延迟后的输出和系统的模型,来重构出系统的当前时刻的状态,从而实现对系统状态的观测。本发明不仅可以实现对飞行器纵向运动动力学模型的观测,还可以实现对其他控制对象的观测,因此,该发明方法具有一定的适用性。
2.技术方案
为达到上述目的,本发明针对包含时变测量延迟输出的飞行器纵向动力学模型,利用测得的带有延迟和噪声的航迹倾角信号,给出了一种状态观测的方法。
本发明一种时变测量延迟输出信号飞行器纵向运动状态观测方法,该方法包括以下步骤:
步骤1:获取飞行器在俯仰平面内的动力学模型
本发明所针对的观测对象是俯仰平面内运动的飞行器。首先获取纵向运动的动力学模型。飞行器在俯仰平面内的结构示意图如图1所示。在该平面中,飞行器受到来自发动机的推力、空气施加的升力和阻力以及俯仰力矩。通过对所受到的力进行分析,可以得到飞行器在纵向运动中的动力学方程,对其进行线性化处理后,可以得到如下所示的动力学模型表达式:
γ · = L ‾ α α - g V T cos γ + L ‾ o α · = q + g V T cos γ - L ‾ o - L ‾ α α ψ · = q q · = M o + M δ δ - - - ( 1 )
并且有
L ‾ o = L o m V T , L ‾ α = L α m V T
式中,γ、α、ψ分别表示飞行器航迹倾角、攻角和俯仰角且有γ=ψ-α;q为俯仰角变化率;VT为航速;m和g分别为飞行器质量和重力加速度;Lα表示升力曲线斜率;Lo表示其他对升力的影响因素;Mδ表示控制俯仰力矩;Mo表示其他来源力矩,通常由公式Mo=Mαα+Mqq近似,δ为舵面偏角。
定义x1=γ,x2=ψ,x3=q,可将飞行器纵向模型写成状态空间的表达形式:
x · 1 = a 1 x 2 - a 1 x 1 + a 2 - a 3 cos x 1 x · 2 = x 3 x · 3 = b 1 u + b 2 ( x 2 - x 1 ) + b 3 x 3 - - - ( 2 )
其中,系统参数
Figure BDA0000139852150000025
Figure BDA0000139852150000026
a3=g/VT,b1=Mδ>0,b2=Mα,b3=Mq
此外,由于该系统模型中的输出y=x1(t),包含有一定的延迟环节,可将时变延时输出表达为:
y ‾ = x 1 ( t - δ ( t ) )
其中δ(t)为随时间变化的延迟时间,满足|δ(t)|≤1.5,其变化过程仿真图如图6所示。
步骤2:设计观测器
根据上个步骤中获得的动力学模型,对其进行分析,便可设计出满足要求的观测器结构。在保证观测误差收敛的同时,还需要保证观测速度能快速达到设计要求。该步骤是本发明方法中最重要的一个步骤。
观测器结构中包含了模型估计环节和增益调节环节,其数学表达式如下所示:
x ^ · 1 = a 1 x ^ 2 - a 1 x ^ 1 + a 2 - a 3 cos x ^ 1 + G ( 1 ) x ^ · 2 = x ^ 3 + G ( 2 ) x ^ · 3 = b 1 u + b 2 ( x ^ 2 - x ^ 1 ) + b 3 x ^ 3 + G ( 3 ) - - - ( 3 )
其中G(·)为设计的观测器调节增益环节,其表达式为:
G ( · ) = Q - 1 K ( y ‾ - x ^ 1 ( t - δ ( t ) ) )
其中
Q = 1 0 0 M 1 a 1 0 M 2 M 3 a 1
M 1 = - a 1 + a 3 sin x ^ 1 M 3 = - a 1 2 + a 1 a 3 sin x ^ 1
M 2 = a 1 2 - 2 a 1 a 3 sin x ^ 1 - a 1 a 3 x ^ 1 cos x ^ 1 + a 1 a 3 x ^ 2 cos x ^ 2 + a 2 a 3 cos x ^ 1 + a 3 2 sin 2 x ^ 1 - a 3 2 co s 2 x ^ 1
而K=[k1k2k3]T为欲配置极点的配置参数矩阵。若期望配置的极点位置为λ1=-e1,λ2=-e2,λ3=-e3,ei>0,则矩阵K的计算公式为:
k 1 = e 1 + e 2 + e 3 k 2 = e 1 e 2 + e 2 e 3 + e 3 e 1 k 3 = e 1 e 2 e 3 - - - ( 4 )
对上述提出的观测器进行Lyapunov-Razumikhin收敛性分析。首先根据式(2)和式(3),整理出观测误差ξ状态方程的表达式,然后定义Lyapunov函数V(ξ),利用Razumikhin理论,对其进行收敛性分析。采用观测器式(3),可以证明出
Figure BDA0000139852150000038
其中ω>0,即可保证所设计的观测器的收敛性,同时也说明了所提出的观测方法的有效性。但是由于所设计的观测器中包含有很多的参数,为了使观测效果达到最好,即观测器的响应速度最快和观测误差超调量最小,需要对观测器中的参数进行调整。
步骤3:调节观测参数,查看观测效果
按照上一步中所设计的观测器结构,对观测器中的某些参数进行调节,并查看观测效果是否满足要求,如果不满足则继续调节参数,一直到符合要求为止。
在仿真中,选择飞行器系统的参数为:
Figure BDA0000139852150000041
Figure BDA0000139852150000042
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
以上参数均是飞行器本身所固有的参数,不需要做调整。而在仿真过程中,主要调节的参数是欲配置的极点坐标位置。经过多调节后,最终所选择欲配置的极点坐标的绝对值大小为:
e 1 = 0.261 e 2 = 0.262 e 3 = 0.263
此外,在仿真中还加入噪声环节,所加入的噪声的幅值为0.001。
步骤4:设计结束
整个设计过程重点考虑了两个方面的要求,一个是所设计的观测器的收敛性问题;另外一个是要保证观测效果能够满足发明要求。第一步中,需要确定飞行器在俯仰平面内的动力学模型;第二步中,需要设计出满足要求的观测器;第三步,是对之前所设计出的观测器进行参数调整。
3.优点及功效
本发明方法针对输出中存在时变延迟的飞行器纵向运动模型,利用测得的带有延迟和噪声的航迹倾角信号,实现对状态变量的估计和观测,解决了由于测量、信号传输过程或设备等因素造成的信号延时问题,从而有利于控制系统的设计。此外,如果在输出量中存在一定的测量噪声,采用本方法,能从含有噪声的延迟输出信号中提取出待观测信号的实际值,具有一定的抗噪能力。
(四)附图说明
图1本发明中飞行器纵向模型示意图
图2本发明观测系统仿真示意图
图3.1本发明航迹倾角跟踪效果仿真示意图
图3.2本发明俯仰角跟踪效果仿真示意图
图3.3本发明俯仰角速率跟踪效果仿真示意图
图4本发明加入噪声的延迟输出仿真示意图
图5本发明延迟输出信号与观测信号仿真示意图
图6本发明时变延迟时间的仿真图
图7本发明实施流程框图
图中符号说明如下:
γ、α、ψ分别表示飞行器航迹倾角、攻角和俯仰角且有γ=ψ-α;VT为航速;δ为舵面偏角。
(五)具体实施方式
下面将结合附图和实施实例对本发明做进一步的详细说明。
见图7,本发明一种时变测量延迟输出信号飞行器纵向运动状态观测方法,该方法包括以下步骤:
步骤1:获取飞行器在俯仰平面内的动力学模型
本发明所针对的观测对象是,在俯仰平面内运动的飞行器。飞行器在俯仰平面内的结构示意图如图1所示,在该平面中,飞行器受到来自发动机的推力、空气施加的升力和阻力以及俯仰力矩。通过对所受到的力进行分析,可以得到飞行器在纵向运动中的动力学方程,对其进行线性化处理后,可以得到如下所示的动力学模型表达式:
γ · = L ‾ α α - g V T cos γ + L ‾ o α · = q + g V T cos γ - L ‾ o - L ‾ α α ψ · = q q · = M o + M δ δ - - - ( 1 )
并且有
L ‾ o = L o m V T , L ‾ α = L α m V T
式中,γ、α、ψ分别表示飞行器航迹倾角、攻角和俯仰角且有γ=ψ-α;q为俯仰角变化率;VT为航速;m和g分别为飞行器质量和重力加速度;Lα表示升力曲线斜率;Lo表示其他对升力的影响因素;Mδ表示控制俯仰力矩;Mo表示其他来源力矩,通常由公式Mo=Mαα+Mqq近似,δ为舵面偏角。
定义x1=γ,x2=ψ,x3=q,可将飞行器纵向模型写成状态空间的表达形式:
x · 1 = a 1 x 2 - a 1 x 1 + a 2 - a 3 cos x 1 x · 2 = x 3 x · 3 = b 1 u + b 2 ( x 2 - x 1 ) + b 3 x 3 - - - ( 2 )
其中,系统参数
Figure BDA0000139852150000062
Figure BDA0000139852150000063
a3=g/VT,b1=Mδ>0,b2=Mα,b3=Mq
此外,由于该系统模型中的输出y=x1(t),包含有一定的延迟环节,可将时变延时输出表达为:
y ‾ = x 1 ( t - δ ( t ) )
其中δ(t)为随时间变化的延迟时间,满足|δ(t)|≤1.5,其变化过程仿真图如图6所示。而图2为本发明观测系统仿真示意图。
步骤2:设计观测器
根据上个步骤中获得的动力学模型,对其进行分析,可以设计出满足要求的观测器结构。观测器结构中包含了模型估计环节和增益调节环节,其数学表达式如下所示:
x ^ · 1 = a 1 x ^ 2 - a 1 x ^ 1 + a 2 - a 3 cos x ^ 1 + G ( 1 ) x ^ · 2 = x ^ 3 + G ( 2 ) x ^ · 3 = b 1 u + b 2 ( x ^ 2 - x ^ 1 ) + b 3 x ^ 3 + G ( 3 ) - - - ( 3 )
其中G(·)为设计的观测器调节增益环节,其表达式为:
G ( · ) = Q - 1 K ( y ‾ - x ^ 1 ( t - δ ( t ) ) )
其中
Q = 1 0 0 M 1 a 1 0 M 2 M 3 a 1
M 1 = - a 1 + a 3 sin x ^ 1 M 3 = - a 1 2 + a 1 a 3 sin x ^ 1
M 2 = a 1 2 - 2 a 1 a 3 sin x ^ 1 - a 1 a 3 x ^ 1 cos x ^ 1 + a 1 a 3 x ^ 2 cos x ^ 2 + a 2 a 3 cos x ^ 1 + a 3 2 sin 2 x ^ 1 - a 3 2 co s 2 x ^ 1
而K=[k1k2k3]T为欲配置极点的配置参数矩阵。若期望配置的极点位置为λ1=-e1,λ2=-e2,λ3=-e3,ei>0,则矩阵K的计算公式为:
k 1 = e 1 + e 2 + e 3 k 2 = e 1 e 2 + e 2 e 3 + e 3 e 1 k 3 = e 1 e 2 e 3 - - - ( 4 )
对上述提出的观测器进行Lyapunov-Razumikhin收敛性分析。首先根据式(2)和式(3),整理出观测误差ξ状态方程的表达式,然后定义Lyapunov函数V(ξ),利用Razumikhin理论,对其进行收敛性分析。采用观测器式(3),可以证明出
Figure BDA0000139852150000072
其中ω>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
以上参数均是飞行器本身所固有的参数,不需要做调整。而在仿真过程中,主要调节的参数是欲配置的极点坐标位置。经过多调节后,最终所选择欲配置的极点坐标的绝对值大小为:
e 1 = 0.261 e 2 = 0.262 e 3 = 0.263
此外,在仿真中还加入噪声环节,所加入的噪声的幅值为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:获取飞行器在俯仰平面内的动力学模型
对飞行器在俯仰平面内的受力情况分析,飞行器受到来自发动机的推力、空气施加的升力和阻力以及俯仰力矩;从而得到飞行器在纵向运动中的动力学方程,对其进行线性化处理后,得到如下所示的动力学模型表达式:
γ · = L ‾ α α - g V T cos γ + L ‾ o α · = q + g V T cos γ - L ‾ o - L ‾ α α ψ · = q q · = M o + M δ δ - - - ( 1 )
并且有
L ‾ o = L o mV T , L ‾ α = L α mV T
式中,γ、α、ψ分别表示飞行器航迹倾角、攻角和俯仰角且有γ=ψ-α;q为俯仰角变化率;VT为航速;m和g分别为飞行器质量和重力加速度;Lα表示升力曲线斜率;Lo表示其他对升力的影响因素;Mδ表示控制俯仰力矩;Mo表示其他来源力矩,通常由公式Mo=Mαα+Mqq近似,δ为舵面偏角;
定义x1=γ,x2=ψ,x3=q,将飞行器纵向模型写成状态空间的表达形式:
x · 1 = a 1 x 2 - a 1 x 1 + a 2 - a 3 cos x 1 x · 2 = x 3 x · 3 = b 1 u + b 2 ( x 2 - x 1 ) + b 3 x 3 - - - ( 2 ) 其中,系统参数
Figure FDA00003361299900014
a3=g/VT,b1=Mδ>0,b2=Mα,b3=Mq
此外,由于该系统模型中的输出y=x1(t),包含有一定的延迟环节,将时变延时输出表达为:
y ‾ = x 1 ( t - δ ( t ) )
其中δ(t)为随时间变化的延迟时间,满足|δ(t)|≤1.5;
步骤2:设计观测器
根据上个步骤中获得的动力学模型,对其进行分析,便设计出满足要求的观测器结构;在保证观测误差收敛的同时,还需要保证观测速度能快速达到设计要求;观测器结构中包含了模型估计环节和增益调节环节,其数学表达式如下所示:
x ^ . 1 = a 1 x ^ 2 - a 1 x ^ 1 + a 2 - a 3 cos x ^ 1 + G ( 1 ) x ^ . 2 = x ^ 3 + G ( 2 ) x ^ . 3 = b 1 u + b 2 ( x ^ 2 - x ^ 1 ) + b 3 x ^ 3 + G ( 3 ) - - - ( 3 )
其中G(.)为设计的观测器调节增益环节,其表达式为:
G ( . ) = Q - 1 K ( y ‾ - x ^ 1 ( t - δ ( t ) ) )
其中,G(1)、G(2)、G(3)分别为 G ( . ) = Q - 1 K ( y ‾ - x ^ 1 ( t - δ ( t ) ) ) 的第1、2、3项元素;
其中,
Q = 1 0 0 M 1 a 1 0 M 2 M 3 a 1
M 1 = - a 1 + a 3 sin x ^ 1 M 3 = - a 1 2 + a 1 a 3 sin x ^ 1
M 2 = a 1 2 - 2 a 1 a 3 sin x ^ 1 - a 1 a 3 x ^ 1 cos x ^ 1 + a 1 a 3 x ^ 2 cos x ^ 2 + a 2 a 3 cos x ^ 1 + a 3 2 sin 2 x ^ 1 - a 3 2 cos 2 x ^ 1
而K=[k1 k2 k3]T为欲配置极点的配置参数矩阵;若期望配置的极点位置为λ1=-e12=-e23=-e3,ei>0,则矩阵K的计算公式为:
k 1 = e 1 + e 2 + e 3 k 2 = e 1 e 2 + e 2 e 3 + e 3 e 1 k 3 = e 1 e 2 e 3 - - - ( 4 )
对上述提出的观测器进行Lyapunov-Razumikhin收敛性分析;首先根据式(2)和式(3),整理出观测误差ξ状态方程的表达式,然后定义Lyapunov函数V(ξ),利用Razumikhin理论,对其进行收敛性分析;采用观测器式(3),证明出
Figure FDA00003361299900028
其中ω>0,即可保证所设计的观测器的收敛性,同时也说明了所提出的观测方法的有效性;但是由于所设计的观测器中包含有很多的参数,为了使观测效果达到最好,即观测器的响应速度最快和观测误差超调量最小,需要对观测器中的参数进行调整;
步骤3:调节观测参数,查看观测效果
按照上一步中所设计的观测器结构,对观测器中的参数进行调节,并查看观测效果是否满足要求,如果不满足则继续调节参数,一直到符合要求为止;
在仿真中,选择飞行器系统的参数为:
Figure FDA00003361299900031
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
以上参数均是飞行器本身所固有的参数,不需要做调整;而在仿真过程中,调节的参数是欲配置的极点坐标位置,经过多调节后,最终所选择欲配置的极点坐标的绝对值大小为:
e 1 = 0.261 e 2 = 0.262 e 3 = 0.263
此外,在仿真中还加入噪声环节,所加入的噪声的幅值为0.001;
步骤4:设计结束
整个设计过程重点考虑了两个方面的要求,一个是所设计的观测器的收敛性问题;另外一个是要保证观测效果能够满足发明要求;第一步中,需要确定飞行器在俯仰平面内的动力学模型;第二步中,需要设计出满足要求的观测器;第三步,是对之前所设计出的观测器进行参数调整。
CN 201210051956 2012-03-01 2012-03-01 一种时变测量延迟输出信号飞行器纵向运动状态观测方法 Expired - Fee Related CN102591212B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (1)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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