CN111650947A - A non-linear control method for stratospheric airship height - Google Patents

A non-linear control method for stratospheric airship height Download PDF

Info

Publication number
CN111650947A
CN111650947A CN202010642425.8A CN202010642425A CN111650947A CN 111650947 A CN111650947 A CN 111650947A CN 202010642425 A CN202010642425 A CN 202010642425A CN 111650947 A CN111650947 A CN 111650947A
Authority
CN
China
Prior art keywords
signal
height
sliding mode
stratospheric airship
differential
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
Application number
CN202010642425.8A
Other languages
Chinese (zh)
Other versions
CN111650947B (en
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.)
Shanghai Jiao Tong University
Original Assignee
Shanghai Jiao Tong 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 Shanghai Jiao Tong University filed Critical Shanghai Jiao Tong University
Priority to CN202010642425.8A priority Critical patent/CN111650947B/en
Publication of CN111650947A publication Critical patent/CN111650947A/en
Application granted granted Critical
Publication of CN111650947B publication Critical patent/CN111650947B/en
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, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/04Control of altitude or depth
    • G05D1/042Control of altitude or depth specially adapted for aircraft
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64BLIGHTER-THAN AIR AIRCRAFT
    • B64B1/00Lighter-than-air aircraft
    • B64B1/06Rigid airships; Semi-rigid airships
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/08Control of attitude, i.e. control of roll, pitch, or yaw
    • G05D1/0808Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Mechanical Engineering (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明公开了一种平流层飞艇高度非线性控制方法,属于飞艇飞行控制领域。其仅需测量高度信号,并与期望高度指令进行比较形成高度误差信号。同时采用非线性变换与积分形成终端滑模信号。同时构建差分错位器,得到终端滑模信号的差分错位信号,为系统提供阻尼。然后构建干扰观测器,对高度通道的非线性项与不确定性进行观测估计,最终由终端滑模信号、干扰观测信号、差分错位信号形成最终的飞艇俯仰角期望信号,由飞艇姿态稳定系统进行跟踪,最终实现平流层飞艇的高度精准平稳控制。该方法的优点在于对大高度信号与小高度信号均具有很好的适应性,同时高度控制平稳,过度过程动态特性好,非常适合平流层飞艇的高度控制。

Figure 202010642425

The invention discloses a non-linear control method for the height of a stratospheric airship, belonging to the field of airship flight control. It only needs to measure the height signal and compare it with the desired height command to form the height error signal. At the same time, nonlinear transformation and integration are used to form terminal sliding mode signal. At the same time, a differential dislocation device is constructed to obtain the differential dislocation signal of the terminal sliding mode signal, which provides damping for the system. Then a disturbance observer is constructed to observe and estimate the nonlinear term and uncertainty of the altitude channel. Finally, the final desired signal of airship pitch angle is formed by the terminal sliding mode signal, disturbance observation signal and differential dislocation signal, which is carried out by the airship attitude stabilization system. Tracking, and finally achieve highly precise and stable control of the stratospheric airship. The advantage of this method is that it has good adaptability to both high-altitude signals and small-altitude signals, and at the same time, the altitude control is stable, and the dynamic characteristics of the transition process are good, which is very suitable for the altitude control of stratospheric airships.

Figure 202010642425

Description

一种平流层飞艇高度非线性控制方法A non-linear control method for stratospheric airship height

技术领域technical field

本发明属于飞行器飞行制导领域,尤其是涉及平流层飞艇的高度非线性控制方法。The invention belongs to the field of aircraft flight guidance, in particular to a highly nonlinear control method of a stratospheric airship.

背景技术Background technique

平流层飞艇由于飞行速度较低,同时飞行气动外形庞大导致控制延迟大,姿态响应慢的特点。同时飞艇的任务平台一般要求器爬升过程非常平稳,有利于其载物的安全。一般飞行器的高度控制为例保证高度控制具有足够的稳定裕度与平滑性,一般需要测量高度信号的微分,使得控制成本增加。同时高速运动飞行器的高度控制方案无法直接照搬到平流层飞艇的控制,主要原因是其速度太快导致高度控制在无微分信号作用下,平滑性无法保证。The stratospheric airship has the characteristics of low flight speed and large flight aerodynamic shape, which leads to large control delay and slow attitude response. At the same time, the mission platform of the airship generally requires that the climbing process of the aircraft is very stable, which is conducive to the safety of its cargo. Take the altitude control of general aircraft as an example to ensure that the altitude control has sufficient stability margin and smoothness. Generally, it is necessary to measure the differential of the altitude signal, which increases the control cost. At the same time, the altitude control scheme of the high-speed moving aircraft cannot be directly copied to the control of the stratospheric airship. The main reason is that its speed is too fast, and the altitude control cannot be guaranteed under the action of the differential signal.

基于以上背景技术,本发明采用了一种基于终端滑模的非线性控制方法,仅需要测量高度信号,通过错位差分方法提供阻尼信号,并采用干扰观测的方法提供强自适应能力,使得最终高度控制具有很好的鲁棒性,高度控制效果也具有很好的动态特性。Based on the above background technology, the present invention adopts a nonlinear control method based on terminal sliding mode, which only needs to measure the height signal, provide damping signal through the dislocation difference method, and use the interference observation method to provide strong adaptive ability, so that the final height The control has good robustness, and the height control effect also has good dynamic characteristics.

需要说明的是,在上述背景技术部分发明的信息仅用于加强对本发明的背景的理解,因此可以包括不构成对本领域普通技术人员已知的现有技术的信息。It should be noted that the information disclosed in the above Background section is only for enhancing understanding of the background of the invention, and therefore may include information that does not form the prior art known to a person of ordinary skill in the art.

发明内容SUMMARY OF THE INVENTION

本发明的目的在于提供一种平流层飞艇高度非线性控制方法,进而至少在一定程度上克服传统高度控制方法带来的平流层飞艇高度控制平稳性与鲁棒性不足的问题。The purpose of the present invention is to provide a non-linear control method for the height of a stratospheric airship, thereby at least to a certain extent overcome the problems of insufficient stability and robustness of the height control of the stratospheric airship caused by the traditional height control method.

本发明提供了一种平流层飞艇高度非线性控制方法,包括以下步骤:The invention provides a non-linear control method for the height of a stratospheric airship, comprising the following steps:

步骤S10:采用高度表测量平流层飞艇的飞行高度,根据控制任务设置期望高度指令,并进行比较得到高度误差;Step S10: using an altimeter to measure the flight altitude of the stratospheric airship, setting a desired altitude command according to the control task, and comparing to obtain the altitude error;

步骤S20:根据所述的高度误差信号,设计终端型非线性项,并进行积分,得到非线性积分信号,然后构成终端滑模信号;Step S20: Design a terminal nonlinear term according to the height error signal, and perform integration to obtain a nonlinear integrated signal, and then form a terminal sliding mode signal;

步骤S30:根据终端滑模信号,设计差分错位器,得到终端滑模信号的差分错位信号;Step S30: designing a differential dislocation device according to the terminal sliding mode signal to obtain a differential dislocation signal of the terminal sliding mode signal;

步骤S40:根据所述的高度误差信号,设计高度通道干扰观测器,得到干扰观测信号;Step S40: Design an altitude channel interference observer according to the altitude error signal to obtain an interference observation signal;

步骤S50:根据终端滑模信号与差分错位信号以及干扰观测信号,设计最终的平流层飞艇俯仰角期望信号,并输送给飞艇姿态控制系统,实现高度控制。Step S50: Design the final expected signal of the pitch angle of the stratospheric airship according to the terminal sliding mode signal, the differential dislocation signal and the interference observation signal, and send it to the airship attitude control system to realize altitude control.

在本发明的一种示例实施例中,设计终端型非线性项,并进行积分,得到非线性积分信号,然后构成终端滑模信号包括:In an exemplary embodiment of the present invention, the terminal type nonlinear term is designed and integrated to obtain a nonlinear integral signal, and then forming a terminal sliding mode signal includes:

eh=h-hde h =hh d ;

Figure BDA0002571678630000021
Figure BDA0002571678630000021

sef=∫ef1dt;s ef =∫e f1 dt;

s1=c5eh+c6sefs 1 =c 5 e h +c 6 s ef ;

其中eh为高度误差信号,hd为根据平流层飞艇的飞行任务设置的期望高度指令,h为平流层飞艇的实际高度,ef1为终端型非线性项,c1、c2、c3、c4为常值参数,其详细设计见后文案例实施。sef为终端型非线性项的积分,s1为终端滑模信号,c5、c6为常值参数,其详细设计见后文案例实施。where e h is the height error signal, h d is the desired height command set according to the flight mission of the stratospheric airship, h is the actual height of the stratospheric airship, e f1 is the terminal nonlinear term, c 1 , c 2 , c 3 , c 4 is a constant value parameter, and its detailed design is shown in the following case implementation. s ef is the integral of the terminal nonlinear term, s 1 is the terminal sliding mode signal, and c 5 and c 6 are constant value parameters. The detailed design is shown in the following case implementation.

在本发明的一种示例实施例中,设计差分错位器,得到终端滑模信号的差分错位信号包括:In an exemplary embodiment of the present invention, designing a differential dislocation device to obtain the differential dislocation signal of the terminal sliding mode signal includes:

sd1=(s1-y1)/T1s d1 =(s 1 -y 1 )/T 1 ;

n≥1时,y1(n+1)=y1(n)+sd1T;When n≥1, y 1 (n+1)=y 1 (n)+s d1 T;

sd2=(s1-y2)/T2s d2 =(s 1 -y 2 )/T 2 ;

n≥1时,y2(n+1)=y2(n)+sd2T;When n≥1, y 2 (n+1)=y 2 (n)+s d2 T;

Figure BDA0002571678630000031
Figure BDA0002571678630000031

s1为所述的终端滑模信号,y1、y2分别为第一个与第二个差分器的输出信号,其初始值均为0,y1(1)=0,y2(1)=0。sd1、sd2分别为第一个与第二个差分器的增长信号,T1、T2为时间常数,其详细设计见后文案例实施。y1(n)为第一个差分器的输出信号y1的第n个数据,T为数据间的时间间隔,其详细设计见后文案例实施。y2(n)为第二个差分器的输出信号y2的第n个数据。sd即为所求的差分错位信号。s 1 is the terminal sliding mode signal, y 1 and y 2 are the output signals of the first and second differentiators respectively, and their initial values are both 0, y 1 (1)=0, y 2 (1 )=0. s d1 and s d2 are the growth signals of the first and second differentiators, respectively, and T 1 and T 2 are time constants. The detailed design is shown in the following case implementation. y 1 (n) is the nth data of the output signal y 1 of the first differentiator, and T is the time interval between the data. The detailed design is shown in the following case implementation. y 2 (n) is the n-th data of the output signal y 2 of the second differentiator. s d is the required differential dislocation signal.

在本发明的一种示例实施例中,设计高度通道干扰观测器,得到干扰观测信号包括:In an exemplary embodiment of the present invention, designing a height channel interference observer, and obtaining the interference observation signal includes:

Figure BDA0002571678630000032
Figure BDA0002571678630000032

Figure BDA0002571678630000033
Figure BDA0002571678630000033

其中z为干扰观测器的状态变量,其初始值为零,即z(1)=0。

Figure BDA0002571678630000034
为干扰观测信号为,k1、k2、k6为常值参数,其详细设计见后文案例实施,θd为俯仰角指令信号,其详细设计见下一步设计,其初始值可以选取为0,即θd(1)=0。其中T为数据间的时间间隔,与差分器的设置相同即可。Among them, z is the state variable of the disturbance observer, and its initial value is zero, that is, z(1)=0.
Figure BDA0002571678630000034
In order to interfere with the observation signal, k 1 , k 2 , and k 6 are constant value parameters, and the detailed design is shown in the following case implementation, θ d is the pitch angle command signal, and the detailed design is shown in the next design. 0, that is, θ d (1)=0. Among them, T is the time interval between data, which is the same as the setting of the differentiator.

在本发明的一种示例实施例中,根据终端滑模信号与差分错位信号以及干扰观测信号,设计最终的平流层飞艇俯仰角期望信号包括:In an exemplary embodiment of the present invention, according to the terminal sliding mode signal, the differential dislocation signal and the interference observation signal, designing the final desired signal of the pitch angle of the stratospheric airship includes:

Figure BDA0002571678630000035
Figure BDA0002571678630000035

其中s1为终端滑模信号,sd为差分错位信号,

Figure BDA0002571678630000036
为干扰观测信号,θd为最终的平流层飞艇俯仰角期望信号,其初始值选为0,k3、k4、k5、k7为常值参数,其详细设置见后文案例实施。where s 1 is the terminal sliding mode signal, s d is the differential dislocation signal,
Figure BDA0002571678630000036
In order to interfere with the observation signal, θ d is the final expected signal of the pitch angle of the stratospheric airship, and its initial value is selected as 0, and k 3 , k 4 , k 5 , and k 7 are constant parameters.

在此基础上,让飞行器姿态控制系统跟踪最终得到的平流层飞艇俯仰角期望信号,即可消除高度控制误差,实现平流层飞艇的高度控制。On this basis, let the aircraft attitude control system track the final desired signal of the pitch angle of the stratospheric airship, which can eliminate the altitude control error and realize the altitude control of the stratospheric airship.

有益效果beneficial effect

本发明提供一种平流层飞艇高度非线性控制方法,其优点之一在于采用了一类新型的终端滑模将高度误差信号进行组合,使得最终的高度控制方法对大信号与小信号都具有很好的适应能力。同时又采用了干扰观测器,对高度通道中高度信号与俯仰角信号之间存在的三角变化非线性与不确定性进行估计,使得最终的高度非线性控制方法具有很好的鲁棒性。最后,该高度非线性方法仅需要测量高度信号,采用错位差分方法提供阻尼,而不需要测量速度信号提供阻尼,因此其具有物理意义明确,实施简单,信号平滑的特点,最终使得高度变化十分平缓,非常适合平流层飞艇的实际应用需求。The present invention provides a non-linear control method for the height of a stratospheric airship. One of its advantages lies in that a new type of terminal sliding mode is used to combine the height error signals, so that the final height control method has very good effects on both large and small signals. good adaptability. At the same time, a disturbance observer is used to estimate the nonlinearity and uncertainty of the triangular variation between the height signal and the pitch angle signal in the height channel, so that the final height nonlinear control method has good robustness. Finally, this highly nonlinear method only needs to measure the height signal, and uses the dislocation differential method to provide damping, but does not need to measure the velocity signal to provide damping, so it has the characteristics of clear physical meaning, simple implementation, and smooth signal, and finally makes the height change very gentle. , which is very suitable for the practical application needs of stratospheric airships.

应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本发明。It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the invention.

附图说明Description of drawings

此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本发明的实施例,并与说明书一起用于解释本发明的原理。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the invention and together with the description serve to explain the principles of the invention. Obviously, the drawings in the following description are only some embodiments of the present invention, and for those of ordinary skill in the art, other drawings can also be obtained from these drawings without creative effort.

图1是本发明提供的一种平流层飞艇高度非线性控制方法的设计实施流程图;Fig. 1 is the design implementation flow chart of a kind of stratospheric airship height nonlinear control method provided by the present invention;

图2是本发明实施例所提供方法的平流层飞艇的飞行高度曲线(单位:米);Fig. 2 is the flight height curve (unit: meter) of the stratospheric airship of the method provided by the embodiment of the present invention;

图3是本发明实施例所提供方法的非线性积分信号曲线(无单位);3 is a nonlinear integral signal curve (unitless) of a method provided by an embodiment of the present invention;

图4是本发明实施例所提供方法的终端滑模信号(无单位);4 is a terminal sliding mode signal (unitless) of a method provided by an embodiment of the present invention;

图5是本发明实施例所提供方法的差分错位信号曲线(无单位);5 is a differential misalignment signal curve (unitless) of a method provided by an embodiment of the present invention;

图6是本发明实施例所提供方法的干扰观测信号曲线(无单位);Fig. 6 is the interference observation signal curve (unitless) of the method provided by the embodiment of the present invention;

图7是本发明实施例所提供方法的平流层飞艇俯仰角期望信号(单位:度);FIG. 7 is a desired signal of the pitch angle of the stratospheric airship (unit: degree) of the method provided by the embodiment of the present invention;

图8是本发明实施例所提供方法的平流层飞艇的俯仰舵偏角信号曲线(单位:度);8 is a pitch rudder declination signal curve (unit: degree) of a stratospheric airship according to a method provided by an embodiment of the present invention;

图9本发明实施例所提供方法的平流层飞艇的俯仰角信号曲线(单位:度);Fig. 9 pitch angle signal curve (unit: degree) of the stratospheric airship according to the method provided by the embodiment of the present invention;

具体实施方式Detailed ways

现在将参考附图更全面地描述示例实施方式。然而,示例实施方式能够以多种形式实施,且不应被理解为限于在此阐述的范例;相反,提供这些实施方式使得本发明将更加全面和完整,并将示例实施方式的构思全面地传达给本领域的技术人员。所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施方式中。在下面的描述中,提供许多具体细节从而给出对本发明的实施方式的充分理解。然而,本领域技术人员将意识到,可以实践本发明的技术方案而省略所述特定细节中的一个或更多,或者可以采用其它的方法、组元、装置、步骤等。在其它情况下,不详细示出或描述公知技术方案以避免喧宾夺主而使得本发明的各方面变得模糊。Example embodiments will now be described more fully with reference to the accompanying drawings. Example embodiments, however, can be embodied in various forms and should not be construed as limited to the examples set forth herein; rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the concept of example embodiments to those skilled in the art. The described features, structures, or characteristics may be combined in any suitable manner in one or more embodiments. In the following description, numerous specific details are provided in order to give a thorough understanding of the embodiments of the present invention. However, those skilled in the art will appreciate that the technical solutions of the present invention may be practiced without one or more of the specific details, or other methods, components, devices, steps, etc. may be employed. In other instances, well-known solutions have not been shown or described in detail to avoid obscuring aspects of the present invention.

本发明一种平流层飞艇高度非线性控制方法,其主要通过高度误差信号进行非线性积分组成终端滑模信号,同时采用差分错位方法得到终端滑模信号的错位差分,从而为系统提供阻尼信号。然后根据干扰观测器对高度通道不方便测量计算的非线性与不确定性进行估计观测,形成最终的高度非线性控制方法。由于终端滑模带来的强鲁棒性以及干扰观测器带来的自适应能力,其特别适合平流层飞艇的高度控制。The present invention is a non-linear control method for the height of a stratospheric airship, which mainly forms a terminal sliding mode signal by nonlinear integration of the height error signal, and at the same time adopts a differential dislocation method to obtain the dislocation difference of the terminal sliding mode signal, thereby providing a damping signal for the system. Then, according to the disturbance observer, the nonlinearity and uncertainty of the height channel inconvenient to measure and calculate are estimated and observed, and the final height nonlinear control method is formed. Due to the strong robustness brought by the terminal sliding mode and the adaptive capability brought by the disturbance observer, it is especially suitable for the altitude control of stratospheric airships.

以下,将结合附图对本发明实例实施例中提及的一种平流层飞艇高度非线性控制方法进行解释以及说明。参考图1所示,一种平流层飞艇高度非线性控制方法可以包括以下步骤:Hereinafter, a method for nonlinear control of the height of a stratospheric airship mentioned in the example embodiments of the present invention will be explained and described with reference to the accompanying drawings. Referring to Figure 1, a method for nonlinear control of the height of a stratospheric airship may include the following steps:

步骤S10:采用高度表测量平流层飞艇的飞行高度,根据控制任务设置期望高度指令,并进行比较得到高度误差。Step S10: Measure the flight height of the stratospheric airship by using an altimeter, set a desired height command according to the control task, and compare to obtain the height error.

具体的,首先根据平流层飞艇的飞行任务,设置期望高度指令,记作hd。其次通过高度表或其它高度测量装置,测量平流层飞艇的高度,记作h。最后,通过比较得到高度误差信号,记作eh。其定义为:eh=h-hdSpecifically, first, according to the flight mission of the stratospheric airship, a desired altitude command is set, denoted as h d . Next, measure the height of the stratospheric airship through an altimeter or other height measuring device, denoted as h. Finally, the height error signal is obtained by comparison, denoted as e h . It is defined as: e h =hh d .

步骤S20:根据所述的高度误差信号,设计终端型非线性项,并进行积分,得到非线性积分信号,然后构成终端滑模信号。Step S20: Design a terminal type nonlinear term according to the height error signal, and perform integration to obtain a nonlinear integral signal, and then form a terminal sliding mode signal.

具体的,首先,根据所述的高度误差信号,设计如下的终端型非线性项,记作ef1,其计算方式如下:Specifically, first, according to the height error signal, the following terminal nonlinear term is designed, denoted as e f1 , and its calculation method is as follows:

Figure BDA0002571678630000061
Figure BDA0002571678630000061

其中c1、c2、c3、c4为常值参数,其详细设计见后文案例实施。Among them, c 1 , c 2 , c 3 , and c 4 are constant value parameters, and the detailed design is shown in the following case implementation.

其次,根据所述的终端非线性项进行积分,得到终端积分信号,记作sef,其积分运算如下:Secondly, perform integration according to the terminal nonlinear term to obtain the terminal integration signal, denoted as s ef , and the integration operation is as follows:

sef=∫ef1dt;s ef =∫e f1 dt;

最后,根据所述的高度误差信号与终端积分信号,组合成终端滑模信号,记作s1,其计算方式如下:Finally, according to the height error signal and the terminal integral signal, a terminal sliding mode signal is combined into a terminal sliding mode signal, denoted as s 1 , and its calculation method is as follows:

s1=c5eh+c6sefs 1 =c 5 e h +c 6 s ef ;

其中c5、c6为常值参数,其详细设计见后文案例实施。Among them, c 5 and c 6 are constant value parameters, and the detailed design is shown in the following case implementation.

步骤S30:根据终端滑模信号,设计差分错位器,得到终端滑模信号的差分错位信号。Step S30: Design a differential dislocation device according to the terminal sliding mode signal to obtain a differential dislocation signal of the terminal sliding mode signal.

具体的,首先根据所述的终端滑模信号s1,构建第一个差分器,设置第一个差分器的输出信号为y1,设置其初始值为0,即y1(1)=0,求取第一个差分器的增长信号为sd1,其计算方式为:Specifically, first, build a first differentiator according to the terminal sliding mode signal s 1 , set the output signal of the first differentiator as y 1 , and set its initial value to 0, that is, y 1 (1)=0 , find the growth signal of the first differentiator as s d1 , the calculation method is:

sd1=(s1-y1)/T1s d1 =(s 1 -y 1 )/T 1 ;

其中T1为时间常数,其详细设计见后文案例实施。然后按照如下递推公式计算第一个差分器的输出信号y1的第n个数据y1(n)Among them, T 1 is the time constant, and its detailed design is shown in the following case implementation. Then calculate the nth data y 1 (n) of the output signal y 1 of the first differentiator according to the following recursive formula

n≥1时,y1(n+1)=y1(n)+sd1T;When n≥1, y 1 (n+1)=y 1 (n)+s d1 T;

其中T为数据间的时间间隔,其详细设计见后文案例实施。Among them, T is the time interval between data, and its detailed design is shown in the implementation of the following case.

其次,首先根据所述的终端滑模信号s1,构建第二个差分器,设置第二个差分器的输出信号为y2,设置其初始值为0,即y2(1)=0,求取第二个差分器的增长信号为sd2,其计算方式为:Secondly, first, build a second differentiator according to the terminal sliding mode signal s 1 , set the output signal of the second differentiator as y 2 , and set its initial value to 0, that is, y 2 (1)=0, Find the growth signal of the second differentiator as s d2 , which is calculated as:

sd2=(s1-y2)/T2s d2 =(s 1 -y 2 )/T 2 ;

其中T2为时间常数,其详细设计见后文案例实施。然后按照如下递推公式计算第二个差分器的输出信号y2的第n个数据y2(n)Among them, T 2 is the time constant, and its detailed design is shown in the following case implementation. Then calculate the nth data y 2 (n) of the output signal y 2 of the second differentiator according to the following recursive formula

n≥1时,y2(n+1)=y2(n)+sd2T;When n≥1, y 2 (n+1)=y 2 (n)+s d2 T;

其中T为数据间的时间间隔,和第一个差分器的设置相同。Where T is the time interval between data, which is the same as the setting of the first differentiator.

最后,根据所述的第一个差分器的输出信号y1与第二个差分器的输出信号y2求解差分错位信号,记作sd,其计算方式如下:Finally, the differential dislocation signal is solved according to the output signal y 1 of the first differentiator and the output signal y 2 of the second differentiator, denoted as s d , and the calculation method is as follows:

Figure BDA0002571678630000071
Figure BDA0002571678630000071

步骤S40:根据所述的高度误差信号,设计高度通道干扰观测器,得到干扰观测信号。Step S40: Design an altitude channel interference observer according to the altitude error signal to obtain an interference observation signal.

具体的,首先设置干扰观测器状态变量z的初始值为零,即z(1)=0。其次,设置干扰观测信号为

Figure BDA0002571678630000072
其计算方式如下:Specifically, first, the initial value of the disturbance observer state variable z is set to zero, that is, z(1)=0. Second, set the interference observation signal as
Figure BDA0002571678630000072
It is calculated as follows:

Figure BDA0002571678630000073
Figure BDA0002571678630000073

最后,按照如下干扰观测器来递推求取干扰观测器的状态变量z:Finally, the state variable z of the disturbance observer is calculated recursively according to the following disturbance observer:

Figure BDA0002571678630000074
Figure BDA0002571678630000074

其中k1、k2、k6为常值参数,其详细设计见后文案例实施,θd为俯仰角指令信号,其详细设计见下一步设计,其初始值可以选取为0,即θd(1)=0。其中T为数据间的时间间隔,与差分器的设置相同即可。Among them, k 1 , k 2 , and k 6 are constant value parameters, and the detailed design is shown in the implementation of the following case . (1)=0. Among them, T is the time interval between data, which is the same as the setting of the differentiator.

步骤S50:根据终端滑模信号与差分错位信号以及干扰观测信号,设计最终的平流层飞艇俯仰角期望信号,并输送给飞艇姿态控制系统,实现高度控制。Step S50: Design the final expected signal of the pitch angle of the stratospheric airship according to the terminal sliding mode signal, the differential dislocation signal and the interference observation signal, and send it to the airship attitude control system to realize altitude control.

具体的,根据终端滑模信号s1与差分错位信号sd以及干扰观测信号

Figure BDA0002571678630000081
设计最终的平流层飞艇俯仰角期望信号,记作θd,其初始值选为0,后续计算如下:Specifically, according to the terminal sliding mode signal s 1 , the differential dislocation signal s d and the interference observation signal
Figure BDA0002571678630000081
The final expected signal of the pitch angle of the stratospheric airship is designed, denoted as θ d , and its initial value is selected as 0, and the subsequent calculation is as follows:

Figure BDA0002571678630000082
Figure BDA0002571678630000082

其中k3、k4、k5、k7为常值参数,其详细设置见后文案例实施。Among them, k 3 , k 4 , k 5 , and k 7 are constant value parameters, and the detailed settings are shown in the following case implementation.

然后,将该俯仰角期望信号输送给平流层飞艇俯仰通道姿态控制系统,进行姿态控制,使得平流层飞艇的实际俯仰角跟踪所述的俯仰角期望指令θd,即可实现对平流层飞艇高度的控制,其中平流层飞艇姿态控制系统比较复杂,也有相应的公开设计方法,非本专利保护的内容,在此不再累述。Then, the desired pitch angle signal is sent to the pitch channel attitude control system of the stratospheric airship, and attitude control is performed, so that the actual pitch angle of the stratospheric airship tracks the desired pitch angle command θ d , so that the height of the stratospheric airship can be adjusted. Among them, the stratospheric airship attitude control system is relatively complex, and there are corresponding disclosed design methods, and the content not protected by this patent will not be repeated here.

案例实施与计算机仿真模拟结果分析Case implementation and analysis of computer simulation results

为验证本发明例所提供方法的有效性,进行以下的案例仿真。In order to verify the effectiveness of the method provided by the example of the present invention, the following case simulation is carried out.

在步骤S10中,设置期望高度指令,记作

Figure BDA0002571678630000083
采用高度表测量平流层飞艇的飞行高度如图2所示。In step S10, a desired height command is set, denoted as
Figure BDA0002571678630000083
Using an altimeter to measure the flight height of the stratospheric airship is shown in Figure 2.

在步骤S20中,选取c1=0.01、c2=0.01、c3=0.01、c4=3,得到非线性积分信号如图3所示。最后选取c5=0.08、c6=0.5,构成终端滑模信号如图4所示。In step S20, c 1 =0.01, c 2 =0.01, c 3 =0.01, and c 4 =3 are selected to obtain a nonlinear integral signal as shown in FIG. 3 . Finally, c 5 =0.08 and c 6 =0.5 are selected to form a terminal sliding mode signal as shown in FIG. 4 .

在步骤S30中,选取T1=10,T2=20,T=0.2,得到终端滑模信号的差分错位信号如图5所示。In step S30, T 1 =10, T 2 =20, and T = 0.2 are selected, and the differential dislocation signal of the terminal sliding mode signal is obtained as shown in FIG. 5 .

在步骤S40中,选取k1=0.05、k2=0.01、k3=0.02,设计高度通道干扰观测器,得到干扰观测信号如图6所示。In step S40, k 1 =0.05, k 2 =0.01, k 3 =0.02 are selected to design a height channel interference observer, and the interference observation signal is obtained as shown in FIG. 6 .

在步骤S50中,选取k3=15、k7=200、k4=0.02、k5=0.02,得到最终的平流层飞艇俯仰角期望信号如图7所示。最终平流层飞艇的俯仰舵偏角如图8所示,平流层飞艇的俯仰角如图9所示。In step S50, k 3 =15, k 7 =200, k 4 =0.02, and k 5 =0.02 are selected to obtain the final desired signal of the pitch angle of the stratospheric airship as shown in FIG. 7 . The final pitch angle of the stratospheric airship is shown in Figure 8, and the pitch angle of the stratospheric airship is shown in Figure 9.

由图2可以看出最终平流层飞艇的高度控制能够很好地完成大小高度指令的跟踪控制任务,飞行平稳。其中800秒前为小高度飞行任务,800秒后为大高度飞行任务。而由图7与图9可以看出,飞艇的俯仰角能够平稳跟踪俯仰角指令任务,同时由图8可以看出飞艇的最大俯仰舵偏角指令不大于30度,因此能够满足工程控制需要。由图4与图5以及图6可以看出,干扰观测信号与滑模信号以及差分错误信号的数量级比较一致,同时由于信号本身比较大,因此相应的控制参数选取也都比较小,从而是比较合理的。综上所述,可以看出本发明例所提供的方法具有信号匹配合理,物理意义明确,同时对平流层飞艇具有高度大小信号均比较平稳准确跟踪的特点,从而具有很高的实用价值。It can be seen from Figure 2 that the altitude control of the final stratospheric airship can well complete the tracking control task of the size and altitude commands, and the flight is stable. Before 800 seconds, it is a low-altitude flight mission, and after 800 seconds, it is a high-altitude flight mission. It can be seen from Figure 7 and Figure 9 that the pitch angle of the airship can smoothly track the pitch angle command task. At the same time, it can be seen from Figure 8 that the maximum pitch angle command of the airship is not greater than 30 degrees, so it can meet the needs of engineering control. It can be seen from Figure 4, Figure 5 and Figure 6 that the magnitude of the interference observation signal is relatively consistent with the sliding mode signal and the differential error signal. At the same time, because the signal itself is relatively large, the corresponding control parameters are also relatively small. reasonable. To sum up, it can be seen that the method provided by the example of the present invention has reasonable signal matching, clear physical meaning, and at the same time, the stratospheric airship has the characteristics of relatively stable and accurate tracking of both height and size signals, so it has high practical value.

本领域技术人员在考虑说明书及实践这里发明的发明后,将容易想到本发明的其他实施例。本申请旨在涵盖本发明的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未发明的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本发明的真正范围和精神由权利要求指出。Other embodiments of the invention will readily occur to those skilled in the art upon consideration of the specification and practice of the invention disclosed herein. This application is intended to cover any variations, uses or adaptations of the invention which follow the general principles of the invention and which include common knowledge or conventional techniques in the art to which the invention is not invented . The specification and examples are to be regarded as exemplary only, with the true scope and spirit of the invention being indicated by the claims.

Claims (5)

1.一种平流层飞艇高度非线性控制方法,其特征在于,包括以下步骤:1. a stratospheric airship height nonlinear control method, is characterized in that, comprises the following steps: 步骤S10:采用高度表测量平流层飞艇的飞行高度,根据控制任务设置期望高度指令,并进行比较得到高度误差;Step S10: using an altimeter to measure the flight altitude of the stratospheric airship, setting a desired altitude command according to the control task, and comparing to obtain the altitude error; 步骤S20:根据所述的高度误差信号,设计终端型非线性项,并进行积分,得到非线性积分信号,然后构成终端滑模信号;Step S20: Design a terminal nonlinear term according to the height error signal, and perform integration to obtain a nonlinear integrated signal, and then form a terminal sliding mode signal; 步骤S30:根据终端滑模信号,设计差分错位器,得到终端滑模信号的差分错位信号;Step S30: designing a differential dislocation device according to the terminal sliding mode signal to obtain a differential dislocation signal of the terminal sliding mode signal; 步骤S40:根据所述的高度误差信号,设计高度通道干扰观测器,得到干扰观测信号;Step S40: Design an altitude channel interference observer according to the altitude error signal to obtain an interference observation signal; 步骤S50:根据终端滑模信号与差分错位信号以及干扰观测信号,设计最终的平流层飞艇俯仰角期望信号,并输送给飞艇姿态控制系统,实现高度控制。Step S50: Design the final expected signal of the pitch angle of the stratospheric airship according to the terminal sliding mode signal, the differential dislocation signal and the interference observation signal, and send it to the airship attitude control system to realize altitude control. 2.根据权利要求1所述的所述的一种平流层飞艇高度非线性控制方法,其特征在于,设计终端型非线性项,并进行积分,得到非线性积分信号,然后构成终端滑模信号包括:2 . The method for nonlinear control of stratospheric airship height according to claim 1 , wherein a terminal type nonlinear term is designed and integrated to obtain a nonlinear integral signal, and then a terminal sliding mode signal is formed. 3 . include: eh=h-hde h =hh d ;
Figure FDA0002571678620000011
Figure FDA0002571678620000011
sef=∫ef1dt;s ef =∫e f1 dt; s1=c5eh+c6sefs 1 =c 5 e h +c 6 s ef ; 其中eh为高度误差信号,hd为根据平流层飞艇的飞行任务设置的期望高度指令,h为平流层飞艇的实际高度,ef1为终端型非线性项,c1、c2、c3、c4、c5、c6为常值参数。sef为终端型非线性项的积分,s1为终端滑模信号。where e h is the height error signal, h d is the desired height command set according to the flight mission of the stratospheric airship, h is the actual height of the stratospheric airship, e f1 is the terminal nonlinear term, c 1 , c 2 , c 3 , c 4 , c 5 , and c 6 are constant value parameters. s ef is the integral of the terminal nonlinear term, and s 1 is the terminal sliding mode signal.
3.根据权利要求1所述的所述的一种平流层飞艇高度非线性控制方法,其特征在于,设计差分错位器,得到终端滑模信号的差分错位信号包括:3. the described a kind of stratospheric airship height nonlinear control method according to claim 1, is characterized in that, designing differential dislocation device, the differential dislocation signal that obtains terminal sliding mode signal comprises: sd1=(s1-y1)/T1s d1 =(s 1 -y 1 )/T 1 ; n≥1时,y1(n+1)=y1(n)+sd1T;When n≥1, y 1 (n+1)=y 1 (n)+s d1 T; sd2=(s1-y2)/T2s d2 =(s 1 -y 2 )/T 2 ; n≥1时,y2(n+1)=y2(n)+sd2T;When n≥1, y 2 (n+1)=y 2 (n)+s d2 T;
Figure FDA0002571678620000021
Figure FDA0002571678620000021
s1为所述的终端滑模信号,y1、y2分别为第一个与第二个差分器的输出信号,其初始值均为0,y1(1)=0,y2(1)=0。sd1、sd2分别为第一个与第二个差分器的增长信号,T1、T2为时间常数。y1(n)为第一个差分器的输出信号y1的第n个数据,T为数据间的时间间隔。y2(n)为第二个差分器的输出信号y2的第n个数据。sd即为所求的差分错位信号。s 1 is the terminal sliding mode signal, y 1 and y 2 are the output signals of the first and second differentiators respectively, and their initial values are both 0, y 1 (1)=0, y 2 (1 )=0. s d1 and s d2 are the increasing signals of the first and second differentiators, respectively, and T 1 and T 2 are time constants. y 1 (n) is the nth data of the output signal y 1 of the first differentiator, and T is the time interval between the data. y 2 (n) is the n-th data of the output signal y 2 of the second differentiator. s d is the required differential dislocation signal.
4.根据权利要求1所述的所述的一种平流层飞艇高度非线性控制方法,其特征在于,设计高度通道干扰观测器,得到干扰观测信号包括:4. the described a kind of stratospheric airship height nonlinear control method according to claim 1, is characterized in that, designing height channel disturbance observer, obtains disturbance observation signal comprises:
Figure FDA0002571678620000022
Figure FDA0002571678620000022
Figure FDA0002571678620000023
Figure FDA0002571678620000023
其中z为干扰观测器的状态变量,其初始值为零,即z(1)=0。
Figure FDA0002571678620000024
为干扰观测信号为,k1、k2、k6为常值参数,θd为俯仰角指令信号,其详细设计见下一步设计,其初始值可以选取为0,即θd(1)=0。其中T为数据间的时间间隔,与差分器的设置相同即可。
Among them, z is the state variable of the disturbance observer, and its initial value is zero, that is, z(1)=0.
Figure FDA0002571678620000024
In order to interfere with the observation signal, k 1 , k 2 , and k 6 are constant parameters, and θ d is the pitch angle command signal. For the detailed design, see the next design. The initial value can be selected as 0, that is, θ d (1)= 0. Among them, T is the time interval between data, which is the same as the setting of the differentiator.
5.根据权利要求1所述的所述的一种平流层飞艇高度非线性控制方法,其特征在于,根据终端滑模信号与差分错位信号以及干扰观测信号,设计最终的平流层飞艇俯仰角期望信号包括:5. the described a kind of stratospheric airship height nonlinear control method according to claim 1, is characterized in that, according to terminal sliding mode signal and differential dislocation signal and interference observation signal, design final stratospheric airship pitch angle expectation Signals include:
Figure FDA0002571678620000031
Figure FDA0002571678620000031
其中s1为终端滑模信号,sd为差分错位信号,
Figure FDA0002571678620000032
为干扰观测信号,θd为最终的平流层飞艇俯仰角期望信号,其初始值选为0。k3、k4、k5、k7为常值参数。
where s 1 is the terminal sliding mode signal, s d is the differential dislocation signal,
Figure FDA0002571678620000032
In order to interfere with the observation signal, θ d is the final expected signal of the pitch angle of the stratospheric airship, and its initial value is selected as 0. k 3 , k 4 , k 5 , and k 7 are constant value parameters.
CN202010642425.8A 2020-07-06 2020-07-06 Stratospheric airship height nonlinear control method Expired - Fee Related CN111650947B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010642425.8A CN111650947B (en) 2020-07-06 2020-07-06 Stratospheric airship height nonlinear control method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010642425.8A CN111650947B (en) 2020-07-06 2020-07-06 Stratospheric airship height nonlinear control method

Publications (2)

Publication Number Publication Date
CN111650947A true CN111650947A (en) 2020-09-11
CN111650947B CN111650947B (en) 2021-09-14

Family

ID=72347731

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010642425.8A Expired - Fee Related CN111650947B (en) 2020-07-06 2020-07-06 Stratospheric airship height nonlinear control method

Country Status (1)

Country Link
CN (1) CN111650947B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112904871A (en) * 2021-01-22 2021-06-04 北京电子工程总体研究所 Capsule height control system based on propeller and pitching cable composite control

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105652880A (en) * 2016-02-24 2016-06-08 中国人民解放军海军航空工程学院 Non-linear anti-saturation height instruction generating method for aircraft large airspace flight
CN106406333A (en) * 2016-12-05 2017-02-15 烟台南山学院 Integral-type-terminal-sliding-form-based method for tracking pitching angle of stratospheric airship
CN106681337A (en) * 2017-01-05 2017-05-17 烟台南山学院 Height keeping flight control method for stratospheric airship based on odd-order sliding mode
CN107368085A (en) * 2017-08-29 2017-11-21 中国人民解放军国防科技大学 Model prediction-based method for controlling height of stratospheric airship in wind field
CN111309040A (en) * 2020-03-02 2020-06-19 中国人民解放军海军航空大学 A method for controlling longitudinal pitch angle of aircraft using simplified fractional derivative

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105652880A (en) * 2016-02-24 2016-06-08 中国人民解放军海军航空工程学院 Non-linear anti-saturation height instruction generating method for aircraft large airspace flight
CN106406333A (en) * 2016-12-05 2017-02-15 烟台南山学院 Integral-type-terminal-sliding-form-based method for tracking pitching angle of stratospheric airship
CN106681337A (en) * 2017-01-05 2017-05-17 烟台南山学院 Height keeping flight control method for stratospheric airship based on odd-order sliding mode
CN107368085A (en) * 2017-08-29 2017-11-21 中国人民解放军国防科技大学 Model prediction-based method for controlling height of stratospheric airship in wind field
CN111309040A (en) * 2020-03-02 2020-06-19 中国人民解放军海军航空大学 A method for controlling longitudinal pitch angle of aircraft using simplified fractional derivative

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
肖鸣超 等: "一种平流层飞艇高度控制方法", 《电光与控制》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112904871A (en) * 2021-01-22 2021-06-04 北京电子工程总体研究所 Capsule height control system based on propeller and pitching cable composite control
CN112904871B (en) * 2021-01-22 2023-02-28 北京电子工程总体研究所 Capsule height control system based on propeller and pitching cable composite control

Also Published As

Publication number Publication date
CN111650947B (en) 2021-09-14

Similar Documents

Publication Publication Date Title
CN109144084B (en) An Attitude Tracking Control Method for Vertical Takeoff and Landing Reusable Vehicle Based on Fixed Time Convergence Observer
CN108490788B (en) A Pitch Channel Inversion Control Method for Aircraft Based on Double Jamming Observation
CN110794864B (en) Aircraft stability control method based on attitude angle rate and attack angle measurement
CN109358634A (en) A robust adaptive control method for hypersonic aircraft
CN109240335B (en) An approach and landing guidance method for an aerospace vehicle
CN110889170B (en) Method for estimating falling angle and attack angle of large-angle attack target of aircraft
CN111208844B (en) A Nonlinear Control Method for Longitudinal Overload of Aircraft Based on Overload Angular Velocity
CN111538236B (en) A Longitudinal Overload Control Method for Aircraft Based on Fractional Approximate Differential Realization of Damping
CN106527122A (en) Stratospheric airship fixed-height flight nonlinear PID control method
CN111309040A (en) A method for controlling longitudinal pitch angle of aircraft using simplified fractional derivative
CN114004023A (en) A method for identifying aerodynamic parameters of aircraft based on cyclic neural network
CN111381598A (en) Height control method of slow unmanned aerial vehicle adopting height and acceleration double errors
CN105759609A (en) Carrier-based aircraft autonomous landing method based on explicitly nonlinear model predictive control
CN108549398B (en) Quadrotor flight control method based on fractional saturation function power switching law
CN107102547A (en) A kind of RLV landing phase Guidance Law acquisition methods based on sliding mode control theory
CN112068444A (en) Aircraft attack angle control method adopting nonlinear self-adaptive sliding mode
CN107315419B (en) A three-stage control strategy to realize long-distance fixed-point hovering of airship
CN111650947A (en) A non-linear control method for stratospheric airship height
CN106681337A (en) Height keeping flight control method for stratospheric airship based on odd-order sliding mode
CN106774385A (en) A kind of dirigible spot hover control method of use adaptive variable structure
CN119511829A (en) A disturbance rejection control method for the longitudinal model of a ground effect vehicle in the ground effect region
CN112068575B (en) Unmanned aerial vehicle carrier landing control method adopting nonlinear interference observation
CN104571100B (en) A non-minimum phase hypersonic vehicle control method
CN111650838B (en) Aircraft pitching channel stabilizing method adopting adaptive observer
CN105259762A (en) Graphical automatic optimization setting method of flight control system control parameter

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210914