CN112068575B - 一种采用非线性干扰观测的无人机着舰控制方法 - Google Patents

一种采用非线性干扰观测的无人机着舰控制方法 Download PDF

Info

Publication number
CN112068575B
CN112068575B CN202010384447.9A CN202010384447A CN112068575B CN 112068575 B CN112068575 B CN 112068575B CN 202010384447 A CN202010384447 A CN 202010384447A CN 112068575 B CN112068575 B CN 112068575B
Authority
CN
China
Prior art keywords
signal
angle
attack
error
aerial vehicle
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
CN202010384447.9A
Other languages
English (en)
Other versions
CN112068575A (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.)
Yantai Nanshan University
Original Assignee
Yantai Nanshan 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 Yantai Nanshan University filed Critical Yantai Nanshan University
Priority to CN202010384447.9A priority Critical patent/CN112068575B/zh
Publication of CN112068575A publication Critical patent/CN112068575A/zh
Application granted granted Critical
Publication of CN112068575B publication Critical patent/CN112068575B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/04Control of altitude or depth
    • G05D1/06Rate of change of altitude or depth
    • G05D1/0607Rate of change of altitude or depth specially adapted for aircraft
    • G05D1/0653Rate of change of altitude or depth specially adapted for aircraft during a phase of take-off or landing
    • G05D1/0676Rate of change of altitude or depth specially adapted for aircraft during a phase of take-off or landing specially adapted for landing
    • G05D1/0684Rate of change of altitude or depth specially adapted for aircraft during a phase of take-off or landing specially adapted for landing on a moving platform, e.g. aircraft carrier

Abstract

本发明是关于一种采用非线性干扰观测的无人机着舰控制方法,其一方面采用下滑速度指令与速度测量值的误差生成期望攻角指令,然后与攻角测量值比较形成攻角误差信号,再构造误差子系统的结构信号,形成攻角子系统不确定性的非线性干扰观测器,对系统不确定性进行动态观测补偿,形成油门偏转角信号。另一方面,由给定常值负的俯仰角与高度误差信号合成的俯仰角指令信号,与实际俯仰角进行对比得到俯仰角误差信号,再与俯仰角速率信号对比生成俯仰角速率误差信号,并构造俯仰角速率子系统的结构信号,生成子系统干扰的动态补偿信号,组成升降舵控制信号,完成自动着舰控制。该方法的优点是能够对系统不确定进行动态补偿,从而提供自动着舰品质。

Description

一种采用非线性干扰观测的无人机着舰控制方法
技术领域
本发明涉及无人机自控着舰控制领域,具体而言,涉及一种采用非线性干扰观测的无人机着舰控制方法。
背景技术
飞行器在着舰过程中,由于甲板长度有限,导致对着舰控制的精度要求非常高。由于飞行器在运动过程中的攻角、俯仰角、速度之间的耦合铰链比较严重,引起对着舰控制的精度要求不仅仅体现在外在的高度控制精度上,而且对飞行器内在的攻角、姿态角速率、姿态角的控制都必须按照预定的期望规律进行精确控制。同时,着舰过程中存在中比较严重的风干扰、海浪起伏干扰等,具体作用在飞行器上则体现为力与力矩的干扰两大方面。上述力与力矩的干扰比较复杂,而且无法进行预先估计与准确测量,因此对干扰进行动态补偿,不失为提高着舰过程品质的一个可选可行的手段。基于上述背景原因,本发明提出了一种对攻角子系统与俯仰角速率子系统分别进行非线性干扰观测器设计的方法,能够对系统响应的力与力矩的不确定性进行动态估计与补偿,从而提高了无人机着舰过程的动态品质,也使得本发明具有很高的工程应用价值。
需要说明的是,在上述背景技术部分发明的信息仅用于加强对本发明的背景的理解,因此可以包括不构成对本领域普通技术人员已知的现有技术的信息。
发明内容
本发明的目的在于提供一种采用非线性干扰观测的无人机着舰控制方法,进而至少在一定程度上克服由于相关技术的限制和缺陷而导致的自动着舰过程中未知不确定干扰的所引起的着舰效果不佳的问题。
步骤S10,根据飞行任务需要,设置指数下滑期望速度信号,在无人机上安装HPS-1H型速度传感装置,测量无人机速度信号,两者进行比较生成速度误差信号,并由速度误差信号生成期望攻角信号;
步骤S20,在无人机上安装SMV-1型攻角传感器,测量无人机攻角信号,并与期望攻角信号进行比较,得到攻角误差信号,并对攻角误差信号进行非线性变换与积分,得到得到攻角误差非线性信号与其积分信号;
步骤S30,在无人机上安装BWD-VG300型惯性测量设备,测量无人机俯仰角速率信号,并根据无人机升降舵偏转角信号,以及所述的攻角误差非线性信号与非线性积分信号,构造形成攻角子系统非线性干扰观测器,得到攻角子系统干扰观测估计补偿信号;
步骤S40,根据所述的攻角子系统干扰观测估计补偿信号、攻角误差信号、攻角误差非线性信号、攻角误差积分信号、攻角误差非线性积分信号与攻角子系统结构信号,进行线性叠加,得到无人机的油门偏转角信号;
步骤S50,在无人机上安装JC-KYW28A型无线电高度表,对无人机高度进行测量,并与甲板高度以及甲板运动补偿信号进行对比,得到高度误差信号,并进行非线性积分得到高度误差非线性积分信号,然后通过高度误差与积分信号生成俯仰角指令信号;
步骤S60,根据所述的俯仰角指令信号并与俯仰角测量信号进行比较得到俯仰角误差信号,并进行非线性变换得到俯仰角非线性误差信号,构建无人机俯仰角速率的期望信号,并与BWD-VG300型惯性测量设备测量的无人机俯仰角速率信号进行对比与非线性积分,得到俯仰角速率误差信号与俯仰角速率误差非线性积分信号;
步骤S70,根据所述的无人机速度测量信号,攻角测量信号,无人机俯仰角速率信号与无人机油门偏角信号,构建俯仰角速度子系统结构信号,然后根据升降舵信号以及俯仰角速率期望信号的导数与俯仰角速率误差信息,构造俯仰角速率子系统的非线性干扰观测器;
步骤S80,根据所述的俯仰角速率误差信号与俯仰角速率误差非线性积分信号,俯仰角速率子系统结构信号以及俯仰角速率子系统的干扰观测估计补偿信号,构造无人机升降舵信号,实现无人机的着舰控制。
在本发明的一种示例实施例中,根据飞行任务需要,设置指数下滑期望速度信号,在无人机上安装HPS-1H型速度传感装置,测量无人机速度信号,两者进行比较生成速度误差信号,并由速度误差信号生成期望攻角信号包括:
Figure GDA0003685762720000031
ev=v-vd
Figure GDA0003685762720000032
其中vd为无人机着舰控制的期望速度信号,v0为无人机的初始速度,vf为无人机的末段着舰速度。τ1为速度下滑快慢的时间调节参数。
Figure GDA0003685762720000033
为指数函数,其中t为时间信号。v为无人机上安装HPS-1H型速度传感装置测量得到的无人机飞速速度,ev为速度误差信号,αd为期望攻角信号,α0为常值指令信号,kv、εi为常值参数,其详细设计见后文案例实施。
在本发明的一种示例实施例中,测量无人机攻角信号并与期望攻角信号进行比较,得到攻角误差信号,并对攻角误差信号进行非线性变换与积分,得到非线性误差信号与其积分信号包括:
eα=α-αd
Figure GDA0003685762720000034
Figure GDA0003685762720000035
其中α为在无人机上安装SMV-1型攻角传感器测量所得的攻角信号,αd为期望攻角信号,eα为攻角误差信号,eα1为攻角误差非线性变换信号,s1为攻角误差非线性积分信号,dt表示对时间信号的积分,ε1、ε2为常值参数信号,其详细设计见后文案例实施。
在本发明的一种示例实施例中,测量无人机的俯仰角速率,并根据无人机升降舵偏转角信号,构造形成攻角子系统非线性干扰观测器,得到攻角子系统干扰观测估计补偿信号包括:
Figure GDA0003685762720000041
Figure GDA0003685762720000042
δp1=a21v+a22α+a23ωz+a24δz+a25eα1+a26s1
其中ωz为在无人机上安装BWD-VG300型惯性测量设备测量得到的无人机俯仰角速率信号,δz为无人机升降舵偏转信号,δp为无人机油门偏角信号,eα1为攻角误差非线性信号,s1为攻角误差非线性积分信号,a2i(i=1,2,3,4,5,6)为常值比例参数信号,其详细设计见后文案例实施。δp1为攻角子系统结构信号,
Figure GDA0003685762720000043
为攻角子系统非线性干扰观测器的状态变量,其初始值选取为0,
Figure GDA0003685762720000047
为期望攻角信号的导数信号。
Figure GDA0003685762720000044
为攻角子系统非线性估计补偿信号。
在本发明的一种示例实施例中,根据所述的攻角子系统干扰观测估计补偿信号、攻角误差信号、攻角误差非线性信号、攻角误差积分信号、攻角误差非线性积分信号与攻角子系统结构信号,进行线性叠加,得到无人机的油门偏转角信号包括:
s2=∫eαdt;
Figure GDA0003685762720000045
其中
Figure GDA0003685762720000046
为攻角子系统干扰观测估计补偿信号、eα为攻角误差信号、eα1为攻角非线性误差信号、s2为攻角误差积分信号、s1为攻角误差非线性积分信号,δp1为攻角子系统结构信号,δp为无人机的油门偏转角信号,k1i(i=1,2,3,4)为常值比例参数信号,其详细设计见后文案例实施。
在本发明的一种示例实施例中,由高度误差信号进行非线性积分得到高度误差非线性积分信号,然后通过高度误差与积分信号生成俯仰角指令信号包括:
eh=h-ha-hb
Figure GDA0003685762720000051
Figure GDA0003685762720000052
其中h为在无人机上安装JC-KYW28A型无线电高度表测量得到的高度信号,ha为甲板高度信号,hb为甲板运动补偿信号,eh为高度误差信号,其中甲板运动补偿信号hb的设计已有成熟的公开方法,非本发明重点内容,不在此不再详细累述。sh为高度误差非线性积分信号,dt表示对时间信号的积分,θ1为俯仰角指令信号,εh、k21、k22、ha、θ1a为常值参数,其详细设计见后文案例实施。
在本发明的一种示例实施例中,根据俯仰角误差信号并进行非线性变换得到俯仰角非线性误差信号,构建无人机俯仰角速率的期望信号,并与无人机俯仰角速率信号进行对比与非线性积分,得到俯仰角速率误差信号与俯仰角速率误差非线性积分信号包括:
eθ=θ-θd
Figure GDA0003685762720000053
Figure GDA0003685762720000054
Figure GDA0003685762720000055
Figure GDA0003685762720000056
其中θd为所述的俯仰角指令信号,θ与ωz为采用BWD-VG300型惯性测量设备得到的俯仰角信号与俯仰角速率信号,eθ为俯仰角误差信号,eθ1为俯仰角误差非线性信号。
Figure GDA0003685762720000061
为俯仰角速率指令信号,τa、τb、k31、k32为常值参数信号,其详细设计见后文案例实施。eω为俯仰角速率误差信号,sω为俯仰角速率误差非线性积分信号。
在本发明的一种示例实施例中,根据所述的无人机速度测量信号,攻角测量信号,俯仰角速度测量信号与无人机油门偏角信号,构建俯仰角速度子系统结构信号,然后根据升降舵信号以及俯仰角速率期望信号的导数与俯仰角速率误差信息,构造俯仰角速率子系统的非线性干扰观测器包括:
δz1=a41v+a42α+a43ωz+a44δp+a45eω+a46sω
Figure GDA0003685762720000062
Figure GDA0003685762720000063
其中v为无人机速度的测量信号、α为攻角测量信号、ωz为俯仰角速率信号、δp为无人机油门偏角信号、eω为俯仰角速率误差信号、sω为俯仰角速率误差非线性积分信号,δz1为俯仰角速率子系统结构信号,a41、a42、a43、a44、a45与a46为常值参数信号,其详细设计见后文案例实施。δz为升降舵信号,
Figure GDA0003685762720000064
为俯仰角速率期望值导数信号,
Figure GDA0003685762720000065
为俯仰角速率子系统非线性干扰观测器的状态变量,其初始值选取为0,
Figure GDA0003685762720000066
为俯仰角速率子系统的干扰观测估计补偿信号。
在本发明的一种示例实施例中,根据所述的俯仰角速率误差信号与俯仰角速度误差非线性积分信号,俯仰角速率子系统结构信号以及俯仰角速率子系统的干扰观测估计补偿信号,构造无人机升降舵信号包括:
Figure GDA0003685762720000067
其中eω为俯仰角速率误差信号、sω为俯仰角速率误差非线性积分信号、δz1为俯仰角速率子系统结构信号、δz为升降舵信号,k41、k42、k43为常值参数信号,其详细设计见后文案例实施。
最后,将所述的升降舵信号输送给无人机升降舵系统,即可驱动无人机在甲板上自动着舰。
有益效果
本发明提供的一种采用非线性干扰观测的无人机着舰控制方法,其优点在于,一方面采用了攻角误差以及误差非线性积分信号构成了攻角子系统结构信号,并组成非线性干扰观测器,从而实现了对攻角子系统的干扰与不确定性的动态补偿。另一方面,根据着舰过程中的高度误差信号生成俯仰角速率误差信号,并进行非线性积分驱动生成俯仰角速率子系统的结构信号,并组成子系统的干扰观测器,对俯仰角速率子系统的干扰与不确定性进行动态补偿。从而实现了无人机油门控制与升降舵之间的解耦,尤其是实现了两者的抗干扰动态补偿,从而提高了着舰过程的抗干扰能力,提高了无人机的着舰品质,也使得本发明方法具有很高的工程应用价值。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本发明。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本发明的实施例,并与说明书一起用于解释本发明的原理。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明提供的一种采用非线性干扰观测的无人机着舰控制方法的流程图;
图2是本发明实施例所提供方法的无人机期望攻角信号曲线(单位:度);
图3是本发明实施例所提供方法的无人机攻角信号曲线(单位:度);
图4是本发明实施例所提供方法的无人机攻角误差信号曲线(单位:度);
图5是本发明实施例所提供方法的无人机攻角误差非线性积分信号曲线(无单位);
图6是本发明实施例所提供方法的无人机攻角子系统非线性估计补偿信号曲线(无单位);
图7是本发明实施例所提供方法的无人机油门偏转角信号曲线(无单位);
图8是本发明实施例所提供方法的无人机俯仰角指令信号曲线(单位:度);
图9是本发明实施例所提供方法的无人机俯仰角误差信号曲线(单位:度);
图10是本发明实施例所提供方法的无人机俯仰角速率信号曲线(单位:弧度/秒);
图11是本发明实施例所提供方法的无人机俯仰角速率误差信号曲线(单位:弧度/秒);
图12是本发明实施例所提供方法的无人机升降舵信号(无单位);
图13是本发明实施例所提供方法的无人机俯仰角速率子系统的干扰观测估计补偿信号(单位:度);
图14是本发明实施例所提供方法的无人机高度信号与高度指令信号对比曲线(单位:米)。
具体实施方式
现在将参考附图基础上更全面地描述示例实施方式。然而,示例实施方式能够以多种形式实施,且不应被理解为限于在此阐述的范例;相反,提供这些实施方式使得本发明将更加全面和完整,并将示例实施方式的构思全面地传达给本领域的技术人员。所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施方式中。在下面的描述中,提供许多具体细节从而给出对本发明的实施方式的充分理解。然而,本领域技术人员将意识到,可以实践本发明的技术方案而省略所述特定细节中的一个或更多,或者可以采用其它的方法、组元、装置、步骤等。在其它情况下,不详细示出或描述公知技术方案以避免喧宾夺主而使得本发明的各方面变得模糊。
本发明提供了一种采用非线性干扰观测的无人机着舰控制方法,其首先由指数下滑的期望速度指令与实际速度测量值对比得到的速度误差信号进行非线性变换得到期望期望攻角信号。再由攻角测量信号与期望期望攻角信号进行对比与积分生成攻角误差与非线性积分信号。再测量无人机的俯仰角速率,组成攻角子系统的结构信号,并由油门偏转角与升降舵信号,组成攻角子系统的干扰观测器,得到攻角子系统不确定性的动态观测补偿信号。其次,由高度误差信号与常值负的俯仰角指令信号切换生成期望俯仰角指令信号,再与俯仰角的测量值进行比较得到俯仰角误差信号,并进行非线性变换得到俯仰角非线性误差信号。再综合俯仰角误差与非线性误差信号形成俯仰角速率期望信号,并与俯仰角速率的测量值比较与积分形成俯仰角速率误差信号与积分信号,然后构造俯仰角速率子系统的结构信号,组成俯仰角速率子系统的干扰观测器,并得到俯仰角速率子系统的不确定性与干扰的动态补偿信号,最终与俯仰角速率误差与误差积分信号,组成无人机升降舵的控制规律,实现无人机的自动着舰控制。
下面,将结合附图对本发明的一种采用非线性干扰观测的无人机着舰控制方法进行进一步的解释以及说明。参考图1所示,该一种采用非线性干扰观测的无人机着舰控制方法包括以下步骤:
步骤S10,根据飞行任务需要,设置指数下滑期望速度信号,在无人机上安装HPS-1H型速度传感装置,测量无人机速度信号,两者进行比较生成速度误差信号,并由速度误差信号生成期望攻角信号。
具体的,首先,根据无人机着舰控制的飞行任务需要,设置指数下滑的期望速度信号,记作vd,其计算方式如下:
Figure GDA0003685762720000101
其中v0为无人机的初始速度,vf为无人机的末段着舰速度。τ1为速度下滑快慢的时间调节参数。
Figure GDA0003685762720000102
为指数函数,其中t为时间信号。
其次,在无人机上安装HPS-1H型速度传感装置,测量无人机的飞速速度,记作v,并与期望速度信号进行比较,得到速度误差信号,记作ev,其计算方式如下:
ev=v-vd
最后,根据速度误差信号,进行非线性变换,得到期望攻角信号,记作αd,其计算方式如下:
Figure GDA0003685762720000103
其中α0为常值指令信号,kv、εi为常值参数,其详细设计见后文案例实施。
步骤S20,在无人机上安装SMV-1型攻角传感器,测量无人机攻角信号,并与期望攻角信号进行比较,得到攻角误差信号,并对攻角误差信号进行非线性变换与积分,得到攻角误差非线性信号与其积分信号。
具体的,首先在无人机上安装SMV-1型攻角传感器,测量无人机的攻角,记作α。然后,对所述的期望攻角信号αd进行比较,得到误差信号,记作eα,其计算方式如下:
eα=α-αd
其次,对攻角误差信号进行非线性变换,得到误差非线性变换信号,记作eα1,其计算方式如下:
Figure GDA0003685762720000104
其中ε1为常值参数信号,其详细设计见后文案例实施。
最后,针对所述的攻角误差信号进行非线性积分,得到攻角误差非线性积分信号,记作s1,其计算方式如下:
Figure GDA0003685762720000111
其中dt表示对时间信号的积分,ε2为常值参数信号,其详细设计见后文案例实施。
步骤S30,在无人机上安装BWD-VG300型惯性测量设备,测量无人机俯仰角速率信号,并根据无人机升降舵偏转角信号,以及所述的攻角误差非线性信号与非线性积分信号,构造形成攻角子系统非线性干扰观测器,得到攻角子系统干扰观测估计补偿信号。
具体的,首先在无人机上安装BWD-VG300型惯性测量设备,测量无人机的俯仰角速率信号,记作ωz。然后,引入无人机升降舵偏转信号,记作δz;引入无人机油门偏角信号δp,记作。
最后,根据所述的攻角误差非线性信号eα1与误差非线性积分信号s1,构造如下的攻角子系统非线性干扰观测器:
Figure GDA0003685762720000112
Figure GDA0003685762720000113
δp1=a21v+a22α+a23ωz+a24δz+a25eα1+a26s1
其中a2i(i=1,2,3,4,5,6)为常值比例参数信号,其详细设计见后文案例实施。δp1为攻角子系统结构信号,
Figure GDA0003685762720000114
为攻角子系统非线性干扰观测器的状态变量,其初始值选取为0,
Figure GDA0003685762720000115
为期望攻角信号的导数信号。
Figure GDA0003685762720000116
为攻角子系统非线性估计补偿信号。
步骤S40,根据所述的攻角子系统干扰观测估计补偿信号、攻角误差信号、攻角误差非线性信号、攻角误差积分信号、攻角误差非线性积分信号与攻角子系统结构信号,进行线性叠加,得到无人机的油门偏转角信号。
具体的,根据所述的攻角子系统干扰观测估计补偿信号
Figure GDA0003685762720000117
攻角误差信号eα、攻角非线性误差信号eα1、攻角误差积分信号s2、攻角误差非线性积分信号s1与攻角子系统结构信号δp1,进行如下的线性组合,得到无人机的油门偏转角信号,记作δp,其计算方式如下:
s2=∫eαdt;
Figure GDA0003685762720000121
其中k1i(i=1,2,3,4)为常值比例参数信号,其详细设计见后文案例实施。
步骤S50,在无人机上安装JC-KYW28A型无线电高度表,对无人机高度进行测量,并与甲板高度以及甲板运动补偿信号进行对比,得到高度误差信号,并进行非线性积分得到高度误差非线性积分信号,然后通过高度误差与积分信号生成俯仰角指令信号。
具体的,首先,在无人机上安装JC-KYW28A型无线电高度表,对无人机高度进行测量,记作h;其次,设置甲板高度为ha,甲板运动补偿信号为hb,然后与无人机高度信号进行对比得到误差信号,记作eh,其计算方式如下:
eh=h-ha-hb
其中甲板运动补偿信号hb的设计已有成熟的公开方法,非本发明重点内容,不在此不再详细累述。
然后,对高度误差信号进行非线性积分,得到高度误差非线性积分信号,记作sh,其计算方式如下:
Figure GDA0003685762720000122
其中dt表示对时间信号的积分,εh为常值参数信号,其详细设置见后文案例实施。
最后,根据所述的高度误差信号与高度误差非线性积分信号,设计俯仰角指令信号,记作θ1,其计算方式如下:
Figure GDA0003685762720000131
其中k21、k22、ha、θ1a为常值参数,其详细设计见后文案例实施。
步骤S60,根据所述的俯仰角指令信号并与俯仰角测量信号进行比较得到俯仰角误差信号,并进行非线性变换得到俯仰角非线性误差信号,构建无人机俯仰角速率的期望信号,并与BWD-VG300型惯性测量设备测量的无人机俯仰角速率信号进行对比与非线性积分,得到俯仰角速率误差信号与俯仰角速率误差非线性积分信号。
具体的,首先,根据所述的俯仰角指令信号θd,与所述的BWD-VG300型惯性测量设备得到的俯仰角信号θ进行比较,得到俯仰角误差信号,记作eθ,其计算方式如下:
eθ=θ-θd
其次,根据所述的俯仰角误差信信号,进行非线性变换,得到俯仰角误差非线性信号,记作eθ1,其计算方式如下:
Figure GDA0003685762720000132
其中τa为常值参数,其详细设计见后文案例实施。
再次,根据俯仰角误差信号eθ与俯仰角误差非线性信号eθ1进行线性组合,得到俯仰角速率指令信号,记作
Figure GDA0003685762720000133
其计算方式如下:
Figure GDA0003685762720000134
其中k31、k32为常值参数信号,其详细设计见后文案例实施。
然后,根据所述的俯仰角速率指令信号
Figure GDA0003685762720000135
与采用BWD-VG300型惯性测量设备测量的俯仰角速率信号ωz进行比较,得到俯仰角速率误差信号记作eω,其计算方式如下:
Figure GDA0003685762720000136
最后,对俯仰角速率误差信号进行积分,得到俯仰角速率误差非线性积分信号sω,其计算方式如下:
Figure GDA0003685762720000141
其中τb为常值参数,其详细设计见后文案例实施。
步骤S70,根据所述的无人机速度测量信号,攻角测量信号,无人机俯仰角速率信号与无人机油门偏角信号,构建俯仰角速度子系统结构信号,然后根据升降舵信号以及俯仰角速率期望信号的导数与俯仰角速率误差信息,构造俯仰角速率子系统的非线性干扰观测器。
具体的,根据所述的无人机速度测量信号v、攻角测量信号α、俯仰角速率信号ωz、无人机油门偏角信号δp、俯仰角速率误差信号eω与俯仰角速率误差非线性积分信号sω,进行比例叠加,构造俯仰角速率子系统结构信号,记作δz1,其计算方式如下:
δz1=a41v+a42α+a43ωz+a44δp+a45eω+a46sω
其中a41、a42、a43、a44、a45与a46为常值参数信号,其详细设计见后文案例实施。
其次,根据所述的俯仰角速率子系统结构信号δz1、升降舵信号δz与俯仰角速率期望值导数信号
Figure GDA0003685762720000142
构建俯仰角速率子系统的非线性干扰观测器如下:
Figure GDA0003685762720000143
Figure GDA0003685762720000144
其中
Figure GDA0003685762720000145
为俯仰角速率子系统非线性干扰观测器的状态变量,其初始值选取为0,
Figure GDA0003685762720000146
为俯仰角速率子系统的干扰观测估计补偿信号。
步骤S80,根据所述的俯仰角速率误差信号与俯仰角速率误差非线性积分信号,俯仰角速率子系统结构信号以及俯仰角速率子系统的干扰观测估计补偿信号,构造无人机升降舵信号,实现无人机的着舰控制。
具体的,针对所述的俯仰角速率误差信号eω、俯仰角速率误差非线性积分信号sω、俯仰角速率子系统结构信号δz1、进行线性叠加,得到升降舵信号,记作δz,其计算方式如下:
Figure GDA0003685762720000151
其中k41、k42、k43为常值参数信号,其详细设计见后文案例实施。
最后,将所述的升降舵信号输送给无人机升降舵系统,即可驱动无人机在甲板上自动着舰。
案例实施与计算机仿真模拟结果分析
为验证本发明所提供方法的正确性与有效性,特提供如下案例仿真进行模拟。
在步骤S10中,选取τ1=0.5、α0=-3.5、kv=1、εi=15,得到期望攻角信号如图2所示。
在步骤S20中,安装SMV-1型攻角传感器得到攻角信号曲线如图3所示,选取ε2=0.005,得到攻角误差信号如图4所示,攻角误差非线性积分信号如图5所示。
在步骤S30中,选取a21=-0.13,a22=-10,a23=29,a24=-2.1,a25=-130,a26=-16,得到攻角子系统非线性估计补偿信号如图6所示。
在步骤S40中,选取k11=-40,k12=-60,k13=-40,k14=-15,得到油门偏转角信号如图7所示。
在步骤S50中,选取k21=-0.12、k22=-0.002、ha=8、θ1a=-0.1,得到俯仰角指令信号如图8所示。
在步骤S60中,选取τa=0.005,τb=0.2,得到俯仰角误差信号如图9所示,俯仰角速率信号如图10所示,俯仰角速率误差信号如图11所示。
在步骤S70中,选取a41=0.0016、a42=-0.6425、a43=-0.1402、a44=0.0012、a45=2、a46=0.2,得到俯仰角速率子系统的干扰观测估计补偿信号如图12所示。
在步骤S80中,选取k41=6、k42=0.15、k43=0.5,得到升降舵信号如图13所示。最终得到无人机高度信号与高度指令信号对比曲线如图14所示。
图2、图3与图4分别给出了期望攻角曲线、攻角信号与攻角误差信号,可以看出攻角能够快速跟踪期望信号,攻角误差能够快速收敛到0。图6与图12为本发明所提供的观测器分别对攻角子系统与俯仰角速率子系统的系统不确定性进行观测得到的补偿信号。可以看出观测器所得到的补偿信号能够随着甲板的起伏波动而动态补偿,从而实现了整个无人机着舰控制系统对系统的外干扰的实时动态补偿。图7与图13给出了无人机的油门偏转角与升降舵信号,可以看出其均满足工程实际要求。图14给出了自动着舰高度信号与高度指令信号对比曲线,可以看出无人机末端的高度能够随着舰艇高度起伏而自动适应跟踪。综上所述,本发明所提供的非线性干扰观测补偿的无人机着舰控制方法,能够实现无人的自动着舰高度控制,尤其是其能够实时动态地对外部环境干扰进行补偿,从而提供无人机自动着舰的控制品质,也使得本发明具有很高的工程实用价值。
本领域技术人员在考虑说明书及实践这类的发明后,将容易想到本发明的其他实施例。本申请旨在涵盖本发明的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未指明的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本发明的真正范围和精神由权利要求指出。

Claims (9)

1.一种采用非线性干扰观测的无人机着舰控制方法,其特征在于,包括以下步骤:
步骤S10,根据飞行任务需要,设置指数下滑期望速度信号,在无人机上安装HPS-1H型速度传感装置,测量无人机速度信号,两者进行比较生成速度误差信号,并由速度误差信号生成期望攻角信号;
步骤S20,在无人机上安装SMV-1型攻角传感器,测量无人机攻角信号,并与期望攻角信号进行比较,得到攻角误差信号,并对攻角误差信号进行非线性变换与积分,得到攻角误差非线性信号与其积分信号;
步骤S30,在无人机上安装BWD-VG300型惯性测量设备,测量无人机俯仰角速率信号,并根据无人机升降舵偏转角信号,以及所述的攻角误差非线性信号与非线性积分信号,构造形成攻角子系统非线性干扰观测器,得到攻角子系统干扰观测估计补偿信号;
步骤S40,根据所述的攻角子系统干扰观测估计补偿信号、攻角误差信号、攻角误差非线性信号、攻角误差积分信号、攻角误差非线性积分信号与攻角子系统结构信号,进行线性叠加,得到无人机的油门偏转角信号;
步骤S50,在无人机上安装JC-KYW28A型无线电高度表,对无人机高度进行测量,并与甲板高度以及甲板运动补偿信号进行对比,得到高度误差信号,并进行非线性积分得到高度误差非线性积分信号,然后通过高度误差与积分信号生成俯仰角指令信号;
步骤S60,根据所述的俯仰角指令信号并与俯仰角测量信号进行比较得到俯仰角误差信号,并进行非线性变换得到俯仰角非线性误差信号,构建无人机俯仰角速率的期望信号,并与BWD-VG300型惯性测量设备测量的无人机俯仰角速率信号进行对比与非线性积分,得到俯仰角速率误差信号与俯仰角速率误差非线性积分信号;
步骤S70,根据所述的无人机速度测量信号、攻角测量信号、无人机俯仰角速率信号与无人机油门偏角信号,构建俯仰角速度子系统结构信号,然后根据升降舵信号以及俯仰角速率期望信号的导数与俯仰角速率误差信息,构造俯仰角速率子系统的非线性干扰观测器;
步骤S80,根据所述的俯仰角速率误差信号与俯仰角速率误差非线性积分信号,俯仰角速率子系统结构信号以及俯仰角速率子系统的干扰观测估计补偿信号,构造无人机升降舵信号,实现无人机的着舰控制。
2.根据权利要求1所述的一种采用非线性干扰观测的无人机着舰控制方法,其特征在于,根据飞行任务需要,设置指数下滑期望速度信号,测量无人机速度信号,两者进行比较生成速度误差信号,并由速度误差信号生成期望攻角信号包括:
Figure FDA0003685762710000021
ev=v-vd
Figure FDA0003685762710000022
其中vd为无人机着舰控制的期望速度信号,v0为无人机的初始速度,vf为无人机的末段着舰速度;τ1为速度下滑快慢的时间调节参数;
Figure FDA0003685762710000023
为指数函数,其中t为时间信号;v为无人机上安装HPS-1H型速度传感装置测量得到的无人机飞速速度,ev为速度误差信号,αd为期望攻角信号,α0为常值指令信号,kv、εi为常值参数。
3.根据权利要求1所述的一种采用非线性干扰观测的无人机着舰控制方法,其特征在于,测量无人机攻角信号并与期望攻角信号进行比较,得到攻角误差信号,并对攻角误差信号进行非线性变换与积分,得到非线性误差信号与其积分信号包括:
eα=α-αd
Figure FDA0003685762710000024
Figure FDA0003685762710000031
其中α为在无人机上安装SMV-1型攻角传感器测量所得的攻角信号,αd为期望攻角信号,eα为攻角误差信号,eα1为攻角误差非线性变换信号,s1为攻角误差非线性积分信号,dt表示对时间信号的积分,ε1、ε2为常值参数信号。
4.根据权利要求1所述的一种采用非线性干扰观测的无人机着舰控制方法,其特征在于,测量无人机的俯仰角速率,并根据无人机升降舵偏转角信号,构造形成攻角子系统非线性干扰观测器,得到攻角子系统干扰观测估计补偿信号包括:
Figure FDA0003685762710000032
Figure FDA0003685762710000033
δp1=a21v+a22α+a23ωz+a24δz+a25eα1+a26s1
其中ωz为在无人机上安装BWD-VG300型惯性测量设备测量得到的无人机俯仰角速率信号,δz为无人机升降舵偏转信号,δp为无人机油门偏角信号,eα1为攻角误差非线性信号,s1为攻角误差非线性积分信号,a2i(i=1,2,3,4,5,6)为常值比例参数信号;δp1为攻角子系统结构信号,
Figure FDA0003685762710000034
为攻角子系统非线性干扰观测器的状态变量,其初始值选取为0,
Figure FDA0003685762710000035
为期望攻角信号的导数信号;
Figure FDA0003685762710000036
为攻角子系统非线性估计补偿信号。
5.根据权利要求1所述的一种采用非线性干扰观测的无人机着舰控制方法,其特征在于,根据所述的攻角子系统干扰观测估计补偿信号、攻角误差信号、攻角误差非线性信号、攻角误差积分信号、攻角误差非线性积分信号与攻角子系统结构信号,进行线性叠加,得到无人机的油门偏转角信号包括:
s2=∫eαdt;
Figure FDA0003685762710000037
其中
Figure FDA0003685762710000047
为攻角子系统干扰观测估计补偿信号、eα为攻角误差信号、eα1为攻角非线性误差信号、s2为攻角误差积分信号、s1为攻角误差非线性积分信号,δp1为攻角子系统结构信号,δp为无人机的油门偏转角信号,k1i(i=1,2,3,4)为常值比例参数信号。
6.根据权利要求1所述的一种采用非线性干扰观测的无人机着舰控制方法,其特征在于,由高度误差信号进行非线性积分得到高度误差非线性积分信号,然后通过高度误差与积分信号生成俯仰角指令信号包括:
eh=h-ha-hb
Figure FDA0003685762710000041
Figure FDA0003685762710000042
其中h为在无人机上安装JC-KYW28A型无线电高度表测量得到的高度信号,ha为甲板高度信号,hb为甲板运动补偿信号,eh为高度误差信号,sh为高度误差非线性积分信号,dt表示对时间信号的积分,θ1为俯仰角指令信号,εh、k21、k22、ha、θ1a为常值参数。
7.根据权利要求1所述的一种采用非线性干扰观测的无人机着舰控制方法,其特征在于,根据俯仰角误差信号并进行非线性变换得到俯仰角非线性误差信号,构建无人机俯仰角速率的期望信号以及俯仰角速率误差信号与俯仰角速率误差非线性积分信号包括:
eθ=θ-θd
Figure FDA0003685762710000043
Figure FDA0003685762710000044
Figure FDA0003685762710000045
Figure FDA0003685762710000046
其中θd为所述的俯仰角指令信号,θ与ωz为采用BWD-VG300型惯性测量设备得到的俯仰角信号与俯仰角速率信号,eθ为俯仰角误差信号,eθ1为俯仰角误差非线性信号;
Figure FDA0003685762710000051
为俯仰角速率指令信号,τa、τb、k31、k32为常值参数信号;eω为俯仰角速率误差信号,sω为俯仰角速率误差非线性积分信号。
8.根据权利要求1所述的一种采用非线性干扰观测的无人机着舰控制方法,其特征在于,根据所述的无人机速度测量信号,攻角测量信号,俯仰角速度测量信号与无人机油门偏角信号,构建俯仰角速度子系统结构信号与俯仰角速率子系统的非线性干扰观测器包括:
δz1=a41v+a42α+a43ωz+a44δp+a45eω+a46sω
Figure FDA0003685762710000052
Figure FDA0003685762710000053
其中v为无人机速度的测量信号、α为攻角测量信号、ωz为俯仰角速率信号、δp为无人机油门偏角信号、eω为俯仰角速率误差信号、sω为俯仰角速率误差非线性积分信号,δz1为俯仰角速率子系统结构信号,a41、a42、a43、a44、a45与a46为常值参数信号;δz为升降舵信号,
Figure FDA0003685762710000054
为俯仰角速率期望值导数信号,
Figure FDA0003685762710000055
为俯仰角速率子系统非线性干扰观测器的状态变量,其初始值选取为0,
Figure FDA0003685762710000056
为俯仰角速率子系统的干扰观测估计补偿信号。
9.根据权利要求1所述的一种采用非线性干扰观测的无人机着舰控制方法,其特征在于,根据所述的俯仰角速率误差信号与俯仰角速率误差非线性积分信号,俯仰角速率子系统结构信号以及俯仰角速率子系统的干扰观测估计补偿信号,构造无人机升降舵信号包括:
Figure FDA0003685762710000057
其中eω为俯仰角速率误差信号、sω为俯仰角速率误差非线性积分信号、δz1为俯仰角速率子系统结构信号、δz为升降舵信号,k41、k42、k43为常值参数信号。
CN202010384447.9A 2020-05-07 2020-05-07 一种采用非线性干扰观测的无人机着舰控制方法 Expired - Fee Related CN112068575B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010384447.9A CN112068575B (zh) 2020-05-07 2020-05-07 一种采用非线性干扰观测的无人机着舰控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010384447.9A CN112068575B (zh) 2020-05-07 2020-05-07 一种采用非线性干扰观测的无人机着舰控制方法

Publications (2)

Publication Number Publication Date
CN112068575A CN112068575A (zh) 2020-12-11
CN112068575B true CN112068575B (zh) 2022-08-23

Family

ID=73656152

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010384447.9A Expired - Fee Related CN112068575B (zh) 2020-05-07 2020-05-07 一种采用非线性干扰观测的无人机着舰控制方法

Country Status (1)

Country Link
CN (1) CN112068575B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114047776B (zh) * 2021-11-27 2023-10-13 山东创惠电子科技有限责任公司 一种基于双时变增益自适应的飞行器攻角跟踪控制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1901153A1 (en) * 2006-09-12 2008-03-19 OFFIS e.V. Control system for unmanned 4-rotor-helicopter
CN106406102A (zh) * 2016-12-20 2017-02-15 哈尔滨工业大学 一种含干扰观测器的高超声速飞行器跟踪控制方法
CN107942651A (zh) * 2017-10-20 2018-04-20 南京航空航天大学 一种近空间飞行器控制系统
CN109144084A (zh) * 2018-07-11 2019-01-04 哈尔滨工业大学 一种基于固定时间收敛观测器的垂直起降重复使用运载器姿态跟踪控制方法
CN109917800A (zh) * 2019-03-19 2019-06-21 中国科学院深圳先进技术研究院 无人机控制方法、装置、计算机设备及存储介质

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1901153A1 (en) * 2006-09-12 2008-03-19 OFFIS e.V. Control system for unmanned 4-rotor-helicopter
CN106406102A (zh) * 2016-12-20 2017-02-15 哈尔滨工业大学 一种含干扰观测器的高超声速飞行器跟踪控制方法
CN107942651A (zh) * 2017-10-20 2018-04-20 南京航空航天大学 一种近空间飞行器控制系统
CN109144084A (zh) * 2018-07-11 2019-01-04 哈尔滨工业大学 一种基于固定时间收敛观测器的垂直起降重复使用运载器姿态跟踪控制方法
CN109917800A (zh) * 2019-03-19 2019-06-21 中国科学院深圳先进技术研究院 无人机控制方法、装置、计算机设备及存储介质

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Robust Backstepping Control for Flight Motion Simulator Based on Nonlinear Disturbance Oberver;ZH Liu,et al.;《Journal of System Simulation》;20081231;第5354-5357页 *
考虑输入饱和的制导控制一体化设计;王松艳 等;《航空学报》;20170112(第10期);第320897-1-320897-11页 *

Also Published As

Publication number Publication date
CN112068575A (zh) 2020-12-11

Similar Documents

Publication Publication Date Title
CN111538346B (zh) 一种四旋翼直升飞机的干扰观测补偿飞行控制方法
CN102393641B (zh) 基于甲板运动补偿的舰载机自动着舰引导控制方法
US8209069B1 (en) Aircraft sensor anticipator system
CN113219838B (zh) 一种基于自适应时变带宽观测器的三轴惯性稳定平台高精度控制方法
CN110794864B (zh) 基于姿态角速率与攻角测量的飞行器稳定控制方法
WO2022105408A1 (zh) 基于径向基神经网络的自动驾驶仪及其解耦控制方法
CN111061286B (zh) 一种滤波微分提供阻尼实现飞行器侧向过载控制的方法
CN111459184A (zh) 一种采用分段攻角指令的无人飞行器自动着舰控制方法
CN112068575B (zh) 一种采用非线性干扰观测的无人机着舰控制方法
Niermeyer et al. Open-loop quadrotor flight dynamics identification in frequency domain via closed-loop flight testing
CN105865441A (zh) 一种针对多源干扰系统的复合分层自适应滤波器
CN111813133B (zh) 一种基于相对精密单点定位的无人机舰船自主着陆方法
CN111399530A (zh) 一种基于反解传函的小型飞行器攻角滑模跟踪方法
CN109976380A (zh) 基于卡尔曼滤波估计的隔离度辨识校正方法及系统
Su et al. Anti-disturbance rapid vibration suppression of the flexible aerial refueling hose
Matt et al. Frequency domain system identification of a small flying-wing uas
Adams et al. Modeling and input shaping control of a micro coaxial radio-controlled helicopter carrying a suspended load
Seher-Weiss et al. EC135 system identification for model following control and turbulence modeling
CN113848730A (zh) 一种基于内模原理和自适应时变带宽观测器的三轴惯性稳定平台高精度控制方法
CN111650947B (zh) 一种平流层飞艇高度非线性控制方法
CN109101035A (zh) 一种用于高空滑翔uuv纵平面弹道控制的方法
CN112197768B (zh) 一种测量侧向过载的飞行器反演干扰观测转弯控制方法
Matt et al. Longitudinal System Identification for a Small Flying-wing UAS
Tunik et al. Onboard stabilization of the information and measurement system in the double gimbal suspension of aircraft
Mulder et al. Aerodynamic hinge moment coefficient estimation using automatic fly-by-wire control inputs

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: 20220823