CN111650947A - A non-linear control method for stratospheric airship height - Google Patents
A non-linear control method for stratospheric airship height Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000010354 integration Effects 0.000 claims abstract description 8
- 238000013016 damping Methods 0.000 abstract description 5
- 230000008569 process Effects 0.000 abstract description 2
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 abstract 1
- 230000006641 stabilisation Effects 0.000 abstract 1
- 238000011105 stabilization Methods 0.000 abstract 1
- 230000009466 transformation Effects 0.000 abstract 1
- 230000007704 transition Effects 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 5
- 230000003044 adaptive effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000009194 climbing Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/04—Control of altitude or depth
- G05D1/042—Control of altitude or depth specially adapted for aircraft
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64B—LIGHTER-THAN AIR AIRCRAFT
- B64B1/00—Lighter-than-air aircraft
- B64B1/06—Rigid airships; Semi-rigid airships
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/08—Control of attitude, i.e. control of roll, pitch, or yaw
- G05D1/0808—Control 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
本发明公开了一种平流层飞艇高度非线性控制方法,属于飞艇飞行控制领域。其仅需测量高度信号,并与期望高度指令进行比较形成高度误差信号。同时采用非线性变换与积分形成终端滑模信号。同时构建差分错位器,得到终端滑模信号的差分错位信号,为系统提供阻尼。然后构建干扰观测器,对高度通道的非线性项与不确定性进行观测估计,最终由终端滑模信号、干扰观测信号、差分错位信号形成最终的飞艇俯仰角期望信号,由飞艇姿态稳定系统进行跟踪,最终实现平流层飞艇的高度精准平稳控制。该方法的优点在于对大高度信号与小高度信号均具有很好的适应性,同时高度控制平稳,过度过程动态特性好,非常适合平流层飞艇的高度控制。
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.
Description
技术领域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-hd;e h =hh d ;
sef=∫ef1dt;s ef =∫e f1 dt;
s1=c5eh+c6sef;s 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)/T1;s 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)/T2;s 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;
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:
其中z为干扰观测器的状态变量,其初始值为零,即z(1)=0。为干扰观测信号为,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. 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:
其中s1为终端滑模信号,sd为差分错位信号,为干扰观测信号,θd为最终的平流层飞艇俯仰角期望信号,其初始值选为0,k3、k4、k5、k7为常值参数,其详细设置见后文案例实施。where s 1 is the terminal sliding mode signal, s d is the differential dislocation signal, 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-hd。Specifically, 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:
其中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+c6sef;s 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)/T1;s 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)/T2;s 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:
步骤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。其次,设置干扰观测信号为其计算方式如下: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 It is calculated as follows:
最后,按照如下干扰观测器来递推求取干扰观测器的状态变量z:Finally, the state variable z of the disturbance observer is calculated recursively according to the following disturbance observer:
其中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以及干扰观测信号设计最终的平流层飞艇俯仰角期望信号,记作θd,其初始值选为0,后续计算如下:Specifically, according to the terminal sliding mode signal s 1 , the differential dislocation signal s d and the interference observation signal 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:
其中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中,设置期望高度指令,记作采用高度表测量平流层飞艇的飞行高度如图2所示。In step S10, a desired height command is set, denoted as 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)
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)
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)
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 |
-
2020
- 2020-07-06 CN CN202010642425.8A patent/CN111650947B/en not_active Expired - Fee Related
Patent Citations (5)
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)
Title |
---|
肖鸣超 等: "一种平流层飞艇高度控制方法", 《电光与控制》 * |
Cited By (2)
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 |