CN103592946A - Active section self-adaptive longitudinal guidance method based on view acceleration measurement - Google Patents

Active section self-adaptive longitudinal guidance method based on view acceleration measurement Download PDF

Info

Publication number
CN103592946A
CN103592946A CN201310503680.4A CN201310503680A CN103592946A CN 103592946 A CN103592946 A CN 103592946A CN 201310503680 A CN201310503680 A CN 201310503680A CN 103592946 A CN103592946 A CN 103592946A
Authority
CN
China
Prior art keywords
apparent
program
axial
pitch
program angle
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
CN201310503680.4A
Other languages
Chinese (zh)
Other versions
CN103592946B (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.)
China Academy of Launch Vehicle Technology CALT
Beijing Aerospace Automatic Control Research Institute
Original Assignee
China Academy of Launch Vehicle Technology CALT
Beijing Aerospace Automatic Control Research Institute
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 China Academy of Launch Vehicle Technology CALT, Beijing Aerospace Automatic Control Research Institute filed Critical China Academy of Launch Vehicle Technology CALT
Priority to CN201310503680.4A priority Critical patent/CN103592946B/en
Publication of CN103592946A publication Critical patent/CN103592946A/en
Application granted granted Critical
Publication of CN103592946B publication Critical patent/CN103592946B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Navigation (AREA)

Abstract

一种基于视加速度测量的主动段程序自适应纵向制导方法,根据轴向视速度增量累加和自动适应标准程序角,确保在较大的动力偏差下仍能以较高的精度按程序轨迹飞行。轴向视速度增量累加和-俯仰程序角数据表,以诸元形式装订,运载器仅依靠视加速度测量量即可生成程序指令,无需其他输入条件及复杂运算;制导指令根据实际动力水平自适应变化,无需辅助导引,从而避免了导引量分配的矛盾。在实际使用中对输入输出量均采取可靠性措施确保此发明工程应用可行。

Figure 201310503680

A program-adaptive longitudinal guidance method in the active segment based on apparent acceleration measurement, accumulating and automatically adapting to the standard program angle according to the axial apparent velocity increments, to ensure that the program trajectory can still be flown with high accuracy under large dynamic deviations . Axial apparent speed incremental accumulation and -pitch program angle data table are bound in the form of various elements. The vehicle can generate program instructions only by relying on the apparent acceleration measurement, without other input conditions and complicated calculations; the guidance instructions are automatically generated according to the actual power level. Adapt to changes without auxiliary guidance, thus avoiding the contradiction of guidance amount distribution. In actual use, reliability measures are taken for the input and output quantities to ensure that the engineering application of the invention is feasible.

Figure 201310503680

Description

一种基于视加速度测量的主动段程序自适应纵向制导方法A Method of Active Segment Program Adaptive Longitudinal Guidance Based on Apparent Acceleration Measurement

技术领域technical field

本发明涉及一种制导方法,特别涉及一种主动段程序飞行自适应纵向制导方法。The invention relates to a guidance method, in particular to an adaptive longitudinal guidance method for active segment procedure flight.

背景技术Background technique

运载器飞行的特点是垂直起飞,然后按照预先计算好的飞行程序转弯飞行,制导指令通常以随时间变化的俯仰程序角序列形式给出,为便于使用,制导参数即为以时间为自变量的俯仰程序角插值表或者以时间为自变量的俯仰程序角多项式拟合系数表。箭上飞行软件基于飞行计时进行线性插值计算或多项式计算,即可实时发出制导程序指令。该方法不依靠测量设备,计算简便、可靠。如果飞行条件诸如大气状态、发动机特性、箭体结构等等都符合理想情况,则运载器在程序控制信号的作用下,将完全按照理论计算的轨迹飞行。但实际上,有许多干扰因素使飞行条件显著偏离理想情况,如发动机的秒耗量偏差、比冲偏差、起飞质量偏差、推力偏斜和横移、风等。即使运载器严格按照预定程序飞行,以上这些干扰因素都将干扰作用于运载器的力和力矩,使其飞行偏离预定轨迹。为使得运载器在纵平面内相对预定轨迹偏离在容许的范围内,还需进行法向导引。为保证飞行中的姿态稳定和姿态跟踪精度,须对导引作出限制。The characteristics of the flight of the carrier are to take off vertically, and then turn and fly according to the pre-calculated flight program. The guidance command is usually given in the form of a sequence of pitch program angles that vary with time. For ease of use, the guidance parameters are time-dependent variables The pitch program angle interpolation table or the pitch program angle polynomial fitting coefficient table with time as the independent variable. Arrow flight software performs linear interpolation calculation or polynomial calculation based on flight timing, and can issue guidance program instructions in real time. This method does not rely on measuring equipment, and the calculation is simple and reliable. If the flight conditions such as atmospheric state, engine characteristics, rocket body structure, etc. are all in line with the ideal situation, the vehicle will fly completely according to the theoretically calculated trajectory under the action of the program control signal. But in fact, there are many disturbing factors that make the flight conditions significantly deviate from the ideal situation, such as the second consumption deviation of the engine, the specific impulse deviation, the take-off mass deviation, the thrust deviation and lateral drift, and the wind. Even if the vehicle flies strictly according to the predetermined procedure, the above disturbance factors will interfere with the force and moment acting on the vehicle, causing its flight to deviate from the predetermined trajectory. In order to make the vehicle deviate from the predetermined trajectory in the vertical plane within the allowable range, normal guidance is also required. In order to ensure attitude stability and attitude tracking accuracy in flight, guidance must be restricted.

当存在严重的干扰,特别是较大的动力偏差时,制导系统设计就面临以下矛盾:一方面,为保证纵平面内飞行轨迹的跟踪精度,要增强导引;另一方面,为实现按照预定飞行轨迹的稳定飞行,要限制导引量的大小。When there are serious disturbances, especially large dynamic deviations, the design of the guidance system faces the following contradictions: on the one hand, in order to ensure the tracking accuracy of the flight trajectory in the vertical plane, it is necessary to enhance the guidance; For stable flight of the flight path, it is necessary to limit the amount of guidance.

发明内容Contents of the invention

本发明解决的技术问题是:克服现有技术的不足,提供了一种基于视加速度测量的主动段程序自适应纵向制导方法,解决了在主动段动力偏差较大的情况下,按照预定轨迹飞行的问题,干扰的适应性更强,无需导引辅助,即可实现纵向制导。The technical problem solved by the present invention is: to overcome the deficiencies of the prior art, to provide a program-adaptive longitudinal guidance method for the active section based on apparent acceleration measurement, and to solve the problem of flying according to the predetermined trajectory when the dynamic deviation of the active section is large The problem of interference is stronger, and longitudinal guidance can be achieved without guidance assistance.

本发明的技术方案是:一种基于视加速度测量的主动段程序自适应纵向制导方法,其特征在于步骤如下:The technical solution of the present invention is: an active section program adaptive longitudinal guidance method based on apparent acceleration measurement, characterized in that the steps are as follows:

1)基于视加速度理论参数生成标准程序序列;1) Generate a standard program sequence based on the theoretical parameters of apparent acceleration;

11)选取运载器轴向作为视加速度测量轴;11) Select the axial direction of the vehicle as the axis of apparent acceleration measurement;

12)将预设的时间-俯仰程序角序列代入运载器动力学模型,获得运载器的标称轨迹;根据标称轨迹,获取运载器主动段的轴向视速度增量累加和、俯仰程序角;12) Substitute the preset time-pitch program angle sequence into the dynamic model of the vehicle to obtain the nominal trajectory of the vehicle; according to the nominal trajectory, obtain the cumulative sum of the axial apparent velocity increments and the pitch program angle of the active section of the vehicle ;

13)将轴向视速度增量累加和-俯仰程序角作图,根据控制精度阈值选取数表间隔,并建立轴向视速度增量累加和-俯仰程序角数表,作为标准程序序列;13) Make a graph of the accumulated axial apparent speed increment and - pitch program angle, select the table interval according to the control accuracy threshold, and establish the accumulated axial apparent speed increment and - pitch program angle table as a standard program sequence;

14)在运载器发射准备阶段,将轴向视速度增量累加和-俯仰程序角数表以诸元形式装订;14) In the launch preparation stage of the vehicle, the axial apparent velocity incremental accumulation and -pitch program angle table are bound in the form of various elements;

2)根据测量得到的视加速度测量量,基于标准程序序列进行自适应制导;2) According to the measured apparent acceleration measurement, adaptive guidance is performed based on the standard program sequence;

21)对视加速度测量量进行合理性判别并剔除野值后,获得轴向视速度增量累加和;21) After judging the rationality of the apparent acceleration measurement and eliminating the outliers, the cumulative sum of the axial apparent velocity increments is obtained;

22)运载器飞行过程中,利用步骤21)获得的轴向视速度增量累加和,对轴向视速度增量累加和-俯仰程序角数据表进行线性插值获得俯仰程序角;22) During the flight of the vehicle, use the cumulative sum of axial apparent velocity increments obtained in step 21) to perform linear interpolation on the cumulative sum of axial apparent velocity increments - pitch program angle data table to obtain the pitch program angle;

23)根据预设的程序角变化率幅值对步骤22)获得的俯仰程序角进行合理性判别,并以判别后的俯仰程序角指令控制运载器按预定轨迹飞行。23) Judge the rationality of the pitch program angle obtained in step 22) according to the preset program angle change rate amplitude, and use the judged pitch program angle command to control the vehicle to fly according to the predetermined trajectory.

步骤13)中根据控制精度阈值还可选取分段拟合参数,并建立拟合参数表。In step 13), segment fitting parameters can also be selected according to the control accuracy threshold, and a fitting parameter table can be established.

步骤21)对视加速度测量量进行合理性判别,剔除野值后获得轴向视速度增量累加和的具体方法为:Step 21) The rationality of the apparent acceleration measurement is judged, and the specific method for obtaining the cumulative sum of axial apparent velocity increments after eliminating outliers is as follows:

211)对轴向视速度增量ΔWx1进行合理性判别,获得判别后的轴向视速度增量

Figure BDA0000400706260000021
211) Discriminate the rationality of the axial apparent velocity increment ΔW x1 , and obtain the judged axial apparent velocity increment
Figure BDA0000400706260000021

Figure BDA0000400706260000031
Figure BDA0000400706260000031

式中

Figure BDA0000400706260000032
为合理性判别后的轴向视速度增量,为前一周期视速度增量为主动段采样周期内最大视速度增量取值;In the formula
Figure BDA0000400706260000032
is the axial apparent velocity increment after rationality judgment, is the apparent speed increment of the previous cycle is the value of the maximum apparent velocity increment within the sampling period of the active segment;

212)根据步骤211)得到的判别后的轴向视速度增量

Figure BDA0000400706260000035
获取轴向视速度增量和
Figure BDA0000400706260000036
其中累计零点为发动机点火时刻。212) According to the judged axial apparent velocity increment obtained in step 211)
Figure BDA0000400706260000035
Get the axial apparent velocity delta and
Figure BDA0000400706260000036
The accumulative zero point is the ignition time of the engine.

步骤23)中对步骤22)获得的俯仰程序角进行合理性判别的具体方法为:The specific method for judging the rationality of the pitch program angle obtained in step 22) in step 23) is as follows:

231)获得俯仰程序角增量限幅值

Figure BDA0000400706260000037
231) Obtain the limit value of pitch program angle increment
Figure BDA0000400706260000037

Figure BDA0000400706260000038
Figure BDA0000400706260000038

式中W1、W2分别为俯仰程序角速度限幅转换点对应的视速度增量累加和;

Figure BDA0000400706260000039
分别为各段俯仰程序角增量限幅值;In the formula, W1 and W2 are respectively the cumulative sum of apparent velocity increments corresponding to the angular velocity limit conversion point of the pitch program;
Figure BDA0000400706260000039
Respectively, the pitch program angle increment limit value of each segment;

232)利用俯仰程序角增量限幅值

Figure BDA00004007062600000310
对俯仰程序角进行合理性判别:232) Use pitch program angle increment limit value
Figure BDA00004007062600000310
Pitch program angle Make a reasonable judgment:

Figure BDA00004007062600000312
Figure BDA00004007062600000312

式中

Figure BDA00004007062600000313
为前一周期俯仰程序角指令,
Figure BDA00004007062600000314
为基于视加速度测量的主动段程序自适应纵向制导技术得到的程序指令,即判别后的俯仰程序角指令。In the formula
Figure BDA00004007062600000313
is the pitch program angle command in the previous cycle,
Figure BDA00004007062600000314
It is the program command obtained by the active segment program adaptive longitudinal guidance technology based on apparent acceleration measurement, that is, the pitch program angle command after discrimination.

本发明与现有技术相比的优点在于:The advantage of the present invention compared with prior art is:

(1)制导指令根据实际动力水平自适应变化,无需辅助导引,从而避免了为增强导引与导引约束的矛盾。(1) The guidance command changes adaptively according to the actual power level, without auxiliary guidance, thus avoiding the contradiction between guidance enhancement and guidance constraints.

(2)增强了对干扰的适应性,扩大了程序角的适用范围,从而保证在较大的动力偏差下仍能以较高的精度按程序轨迹飞行。(2) The adaptability to interference is enhanced, and the application range of the program angle is expanded, so as to ensure that it can still fly according to the program trajectory with high accuracy under large dynamic deviation.

(3)实现简单,轴向视速度增量累加和-俯仰程序角数据表以诸元形式装订,运载器仅依靠视加速度测量即可生成程序指令,无需其他输入条件及复杂运算。(3) The implementation is simple, the axial apparent velocity incremental accumulation and the -pitch program angle data table are bound in the form of elements, and the vehicle can generate program instructions only by relying on the apparent acceleration measurement, without other input conditions and complex calculations.

(4)在实际使用中仅依赖惯性测量装置中的加速度计,惯性测量技术的成熟可靠;对输入输出量均采取可靠性措施确保此发明工程应用可行。(4) In actual use, only the accelerometer in the inertial measurement device is relied on, and the inertial measurement technology is mature and reliable; reliability measures are taken for the input and output to ensure that the engineering application of this invention is feasible.

附图说明Description of drawings

图1为基于视加速度理论参数生成标准程序序列步骤。Figure 1 shows the sequence steps of standard program generation based on theoretical parameters of apparent acceleration.

图2为根据视加速度测量量基于标准程序序列进行自适应制导实现步骤。Figure 2 shows the implementation steps of adaptive guidance based on the standard program sequence based on the apparent acceleration measurement.

图3为某单级发动机运载器时间-俯仰程序角序列所对应的曲线图。Fig. 3 is a graph corresponding to the time-pitch program angle sequence of a single-stage engine carrier.

图4为与图3相对应的时间-轴向视速度增量累加和曲线。Fig. 4 is a time-axial apparent velocity incremental sum curve corresponding to Fig. 3 .

图5为与图3相对应的轴向视速度增量累加和-俯仰程序角曲线。Fig. 5 is the axial apparent speed incremental sum-pitch program angle curve corresponding to Fig. 3 .

具体实施方式Detailed ways

下面结合附图对本发明作进一步详细说明。The present invention will be described in further detail below in conjunction with the accompanying drawings.

以某单级发动机运载器为例:Take a single-stage engine carrier as an example:

1)基于视加速度理论参数生成标准程序序列:1) Generate a standard program sequence based on the theoretical parameters of apparent acceleration:

11)根据运载器动力配置与惯性测量敏感元件的安装方式,选取运载器轴向作为视加速度测量轴。在主动段,运载器轴向上推力起主要作用,所以运载器主动段的轴向视速度增量累加和也呈单调递增趋势,完全可作为程序角插值自变量,即由视速度增量累加和确定唯一的程序角。11) According to the power configuration of the vehicle and the installation method of the inertial measurement sensitive components, the axial direction of the vehicle is selected as the apparent acceleration measurement axis. In the active section, the axial upward thrust of the vehicle plays a major role, so the cumulative sum of the axial apparent velocity increments of the active section of the vehicle also shows a monotonous increasing trend, which can be used as the independent variable of the program angle interpolation, that is, the cumulative sum of the apparent velocity increments and identify unique program angles.

12)如图3所示,以某单级发动机运载器为例,其时间-俯仰程序角序列所对应的曲线如图所示。横坐标T表示时间,纵坐标FICX表示俯仰程序角;由图3可以看出运载器起飞后进行迅速转弯,保持一段时间的定轴飞行后,程序角缓慢下降;12) As shown in Figure 3, taking a single-stage engine carrier as an example, the curve corresponding to its time-pitch program angle sequence is shown in the figure. The abscissa T represents time, and the ordinate FICX represents the pitch program angle; it can be seen from Figure 3 that the vehicle makes a rapid turn after take-off, and after maintaining a fixed-axis flight for a period of time, the program angle decreases slowly;

根据预先设计好的时间-俯仰程序角,代入运载器动力学模型进行仿真,获得运载器的标称轨迹,并进一步获得运载器主动段的轴向视速度增量累加和Wk、俯仰程序角

Figure BDA0000400706260000041
画出轴向视速度增量累加和-俯仰程序角曲线如图5,记录Wk序列。According to the pre-designed time-pitch program angle, it is substituted into the dynamic model of the vehicle for simulation, and the nominal trajectory of the vehicle is obtained, and the axial apparent velocity incremental sum W k and the pitch program angle of the active section of the vehicle are further obtained
Figure BDA0000400706260000041
Draw the axial apparent speed incremental accumulation and -pitch program angle curve as shown in Figure 5, record W k , sequence.

13)对于步骤12)获得的Wk

Figure BDA0000400706260000043
序列,可采用分段细化刻度生成
Figure BDA0000400706260000051
数据表,作为标准程序序列,也可根据Wk进行分段多项式拟合,生成拟合参数表。这里以前者为例进行详细说明。观察图5曲线变化规律,在W1、W2处程序角变化率变化显著,可取分段数为3。根据控制精度规划各段插值刻度。图5中,在视速度增量累加和小于W1时,曲线非线性明显,应细化插值刻度,最终建立数据表
Figure BDA0000400706260000053
13) For W k obtained in step 12),
Figure BDA0000400706260000043
sequence, which can be generated using a piecewise refinement scale
Figure BDA0000400706260000051
data sheet, as a standard program sequence, is also available according to W k for Perform piecewise polynomial fitting and generate a fitting parameter table. Here we take the former as an example to describe in detail. Observing the change law of the curve in Figure 5, the program angle change rate changes significantly at W1 and W2, and the number of segments is 3. Plan the interpolation scale of each segment according to the control accuracy. In Fig. 5, when the cumulative sum of apparent velocity increments is less than W1, the curve is obviously non-linear, and the interpolation scale should be refined to finally establish a data table
Figure BDA0000400706260000053

14)在运载器发射准备阶段,将设计好的轴向视速度增量累加和-俯仰程序角数表以诸元形式装订。14) In the launch preparation stage of the vehicle, the designed axial apparent velocity incremental accumulation and -pitch program angle number tables are bound in the form of various elements.

2)根据测量得到的视加速度测量量,基于标准程序序列进行自适应制导;2) According to the measured apparent acceleration measurement, adaptive guidance is performed based on the standard program sequence;

21)飞行器起控之后开始对视速度增量ΔWx1进行合理性判别,剔除野值,计算经合理性判别后的轴向视速度增量

Figure BDA0000400706260000054
其合理性判别方法如下:21) After the aircraft starts to control, it starts to judge the rationality of the apparent speed increment ΔW x1 , eliminates outliers, and calculates the axial apparent speed increment after the rationality judgment
Figure BDA0000400706260000054
The method of judging its rationality is as follows:

&Delta;&Delta; WW &OverBar;&OverBar; xx 11 == &Delta;&Delta; WW xx 11 ,, 00 << &Delta;&Delta; WW xx 11 << &Delta;&Delta; WW &OverBar;&OverBar; MAXMAX &Delta;&Delta; WW &OverBar;&OverBar; xx 11 __ ,, elseelse -- -- -- (( 11 ))

即起控后满足合理性条件

Figure BDA0000400706260000056
取当前采样周期的轴向视速度增量;否则取前一周期视速度增量
Figure BDA0000400706260000057
其中
Figure BDA0000400706260000058
为主动段采样周期内最大视速度增量取值,避免加表输出的野值问题。轴向视速度增量和WNx1累计零点为发动机点火时刻。计算方法为:That is, the condition of rationality is satisfied after control
Figure BDA0000400706260000056
Take the axial apparent velocity increment of the current sampling period; otherwise, take the apparent velocity increment of the previous period
Figure BDA0000400706260000057
in
Figure BDA0000400706260000058
The value of the maximum apparent velocity increment in the sampling period of the active segment is used to avoid the problem of outliers output by adding tables. The axial apparent velocity increment and the cumulative zero point of W Nx1 are the ignition moment of the engine. The calculation method is:

WW Nxnx 11 == &Sigma;&Delta;&Sigma;&Delta; WW &OverBar;&OverBar; xx 11 ;;

22)运载器飞行过程中,利用数据表

Figure BDA00004007062600000510
俯仰程序角按轴向视速度累加和WNx1分段线性插值。根据WNx1插值俯仰程序角,其公式如下22) During the flight of the carrier, use the data sheet
Figure BDA00004007062600000510
The pitch program angle is accumulated according to the axial apparent velocity and W Nx1 piecewise linear interpolation. Interpolate the pitch program angle according to W Nx1 , the formula is as follows

23)根据运载器性能、设计要求及轨迹形式对步骤22)获得的俯仰程序角

Figure BDA00004007062600000512
进行合理性判别,即对程序角变化率进行限幅。如图4所示,为与图3相对应的时间-轴向视速度增量累加和曲线。图4为本发明由时间-俯仰程序角数据表到轴向视速度增量累加和-俯仰程序角数据表的转换方法,根据设好的时间-俯仰程序角,获得运载器的标称轨迹,并记录相应时刻的运载器轴向视速度增量累加和W。T表示时间,W表示视速度增量累加和。t1、t2为程序角变化率变化显著时刻,也为程序转弯段、定轴段、调姿段切换时刻,相应视速度增量累加和分别为W1、W2。三段俯仰程序角增量限幅值分别为其数值由运载器性能及设计要求决定。23) According to the vehicle performance, design requirements and trajectory form, adjust the pitch program angle obtained in step 22)
Figure BDA00004007062600000512
To judge the rationality, that is, to limit the rate of change of the program angle. As shown in FIG. 4 , it is the time-axis apparent velocity incremental sum curve corresponding to FIG. 3 . Fig. 4 is the conversion method of the present invention from the time-pitch program angle data table to the axial apparent velocity incremental accumulation and-pitch program angle data table, according to the set time-pitch program angle, the nominal track of the carrier is obtained, And record the accumulative sum W of the axial apparent velocity increment of the vehicle at the corresponding moment. T represents time, and W represents the cumulative sum of apparent velocity increments. t 1 and t 2 are the moments when the rate of change of program angle changes significantly, and also the switching moments of program turning section, fixed axis section, and attitude adjustment section, and the cumulative sums of corresponding apparent speed increments are W1 and W2 respectively. The three-stage pitch program angle increment limit values are respectively Its value is determined by the vehicle performance and design requirements.

俯仰程序角增量限幅值

Figure BDA0000400706260000062
计算公式为Pitch program angle increment limit value
Figure BDA0000400706260000062
The calculation formula is

Figure BDA0000400706260000063
Figure BDA0000400706260000063

利用

Figure BDA0000400706260000064
对俯仰程序角
Figure BDA0000400706260000065
进行合理性判别:use
Figure BDA0000400706260000064
Pitch program angle
Figure BDA0000400706260000065
Make a reasonable judgment:

Figure BDA0000400706260000066
Figure BDA0000400706260000066

为前一周期俯仰程序角指令,

Figure BDA0000400706260000068
即为基于视加速度测量的主动段程序自适应纵向制导技术的出的程序指令,运载器在此指令导引下按照预定轨迹飞行。 is the pitch program angle command in the previous cycle,
Figure BDA0000400706260000068
It is the program command issued by the active segment program adaptive longitudinal guidance technology based on apparent acceleration measurement, and the vehicle will fly according to the predetermined trajectory under the guidance of this command.

本发明也可应用于具有多级发动机的运载器,针对不同发动机特性,基于视加速度理论参数设计生成相应的标准程序序列,从而实现自适应制导。其设计过程同上所述。The invention can also be applied to a vehicle with multi-stage engines, and according to different engine characteristics, corresponding standard program sequences are designed and generated based on theoretical parameters of apparent acceleration, so as to realize self-adaptive guidance. Its design process is the same as above.

本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。The content that is not described in detail in the description of the present invention belongs to the well-known technology of those skilled in the art.

Claims (4)

1.一种基于视加速度测量的主动段程序自适应纵向制导方法,其特征在于步骤如下:1. an active segment program adaptive longitudinal guidance method based on apparent acceleration measurement, characterized in that the steps are as follows: 1)基于视加速度理论参数生成标准程序序列;1) Generate a standard program sequence based on the theoretical parameters of apparent acceleration; 11)选取运载器轴向作为视加速度测量轴;11) Select the axial direction of the vehicle as the axis of apparent acceleration measurement; 12)将预设的时间-俯仰程序角序列代入运载器动力学模型,获得运载器的标称轨迹;根据标称轨迹,获取运载器主动段的轴向视速度增量累加和、俯仰程序角;12) Substitute the preset time-pitch program angle sequence into the dynamic model of the vehicle to obtain the nominal trajectory of the vehicle; according to the nominal trajectory, obtain the cumulative sum of the axial apparent velocity increments and the pitch program angle of the active section of the vehicle ; 13)将轴向视速度增量累加和-俯仰程序角作图,根据控制精度阈值选取数表间隔,并建立轴向视速度增量累加和-俯仰程序角数表,作为标准程序序列;13) Make a graph of the accumulated axial apparent speed increment and - pitch program angle, select the table interval according to the control accuracy threshold, and establish the accumulated axial apparent speed increment and - pitch program angle table as a standard program sequence; 14)在运载器发射准备阶段,将轴向视速度增量累加和-俯仰程序角数表以诸元形式装订;14) In the launch preparation stage of the vehicle, the axial apparent velocity incremental accumulation and -pitch program angle table are bound in the form of various elements; 2)根据测量得到的视加速度测量量,基于标准程序序列进行自适应制导;2) According to the measured apparent acceleration measurement, adaptive guidance is performed based on the standard program sequence; 21)对视加速度测量量进行合理性判别并剔除野值后,获得轴向视速度增量累加和;21) After judging the rationality of the apparent acceleration measurement and eliminating the outliers, the cumulative sum of the axial apparent velocity increments is obtained; 22)运载器飞行过程中,利用步骤21)获得的轴向视速度增量累加和,对轴向视速度增量累加和-俯仰程序角数据表进行线性插值获得俯仰程序角;22) During the flight of the vehicle, use the cumulative sum of axial apparent velocity increments obtained in step 21) to perform linear interpolation on the cumulative sum of axial apparent velocity increments - pitch program angle data table to obtain the pitch program angle; 23)根据预设的程序角变化率幅值对步骤22)获得的俯仰程序角进行合理性判别,并以判别后的俯仰程序角指令控制运载器按预定轨迹飞行。23) Judge the rationality of the pitch program angle obtained in step 22) according to the preset program angle change rate amplitude, and use the judged pitch program angle command to control the vehicle to fly according to the predetermined trajectory. 2.根据权利要求1所述的一种基于视加速度测量的主动段程序自适应纵向制导方法,其特征在于:步骤13)中根据控制精度阈值还可选取分段拟合参数,并建立拟合参数表。2. An active segment program adaptive longitudinal guidance method based on apparent acceleration measurement according to claim 1, characterized in that: in step 13), segmental fitting parameters can also be selected according to the control accuracy threshold, and the fitting is established Parameters Table. 3.根据权利要求1所述的一种基于视加速度测量的主动段程序自适应纵向制导方法,其特征在于:步骤21)对视加速度测量量进行合理性判别,剔除野值后获得轴向视速度增量累加和的具体方法为:3. A method for adaptive longitudinal guidance based on apparent acceleration measurement in the active segment according to claim 1, characterized in that: step 21) judge the rationality of the measured amount of apparent acceleration, and obtain the axial visual field after eliminating outliers The specific method of accumulating the speed increment is as follows: 211)对轴向视速度增量ΔWx1进行合理性判别,获得判别后的轴向视速度增量
Figure FDA0000400706250000021
211) Discriminate the rationality of the axial apparent velocity increment ΔW x1 , and obtain the judged axial apparent velocity increment
Figure FDA0000400706250000021
Figure FDA0000400706250000022
Figure FDA0000400706250000022
式中
Figure FDA0000400706250000023
为合理性判别后的轴向视速度增量,
Figure FDA0000400706250000024
为前一周期视速度增量
Figure FDA0000400706250000025
为主动段采样周期内最大视速度增量取值;
In the formula
Figure FDA0000400706250000023
is the axial apparent velocity increment after rationality judgment,
Figure FDA0000400706250000024
is the apparent speed increment of the previous cycle
Figure FDA0000400706250000025
is the value of the maximum apparent velocity increment within the sampling period of the active segment;
212)根据步骤211)得到的判别后的轴向视速度增量
Figure FDA0000400706250000026
获取轴向视速度增量和
Figure FDA0000400706250000027
其中累计零点为发动机点火时刻。
212) According to the judged axial apparent velocity increment obtained in step 211)
Figure FDA0000400706250000026
Get the axial apparent velocity delta and
Figure FDA0000400706250000027
The accumulative zero point is the ignition time of the engine.
4.根据权利要求1所述的一种基于视加速度测量的主动段程序自适应纵向制导方法,其特征在于:步骤23)中对步骤22)获得的俯仰程序角进行合理性判别的具体方法为:4. A method for adaptive longitudinal guidance based on apparent acceleration measurement in the active segment according to claim 1, characterized in that: in step 23), the specific method for judging the rationality of the pitch program angle obtained in step 22) is as follows : 231)获得俯仰程序角增量限幅值
Figure FDA0000400706250000028
231) Obtain the limit value of pitch program angle increment
Figure FDA0000400706250000028
式中W1、W2分别为俯仰程序角速度限幅转换点对应的视速度增量累加和;
Figure FDA00004007062500000210
分别为各段俯仰程序角增量限幅值;
In the formula, W1 and W2 are respectively the cumulative sum of apparent velocity increments corresponding to the angular velocity limit conversion point of the pitch program;
Figure FDA00004007062500000210
Respectively, the pitch program angle increment limit value of each segment;
232)利用俯仰程序角增量限幅值对俯仰程序角
Figure FDA00004007062500000212
进行合理性判别:
232) Use pitch program angle increment limit value Pitch program angle
Figure FDA00004007062500000212
Make a reasonable judgment:
Figure FDA00004007062500000213
Figure FDA00004007062500000213
式中为前一周期俯仰程序角指令,为基于视加速度测量的主动段程序自适应纵向制导技术得到的程序指令,即判别后的俯仰程序角指令。In the formula is the pitch program angle command in the previous cycle, It is the program command obtained by the active segment program adaptive longitudinal guidance technology based on apparent acceleration measurement, that is, the pitch program angle command after discrimination.
CN201310503680.4A 2013-10-23 2013-10-23 A kind of powered phase program self-adaptation longitudinal guidance method measured based on apparent acceleration Active CN103592946B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310503680.4A CN103592946B (en) 2013-10-23 2013-10-23 A kind of powered phase program self-adaptation longitudinal guidance method measured based on apparent acceleration

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310503680.4A CN103592946B (en) 2013-10-23 2013-10-23 A kind of powered phase program self-adaptation longitudinal guidance method measured based on apparent acceleration

Publications (2)

Publication Number Publication Date
CN103592946A true CN103592946A (en) 2014-02-19
CN103592946B CN103592946B (en) 2016-02-10

Family

ID=50083137

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310503680.4A Active CN103592946B (en) 2013-10-23 2013-10-23 A kind of powered phase program self-adaptation longitudinal guidance method measured based on apparent acceleration

Country Status (1)

Country Link
CN (1) CN103592946B (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103863558A (en) * 2014-03-10 2014-06-18 北京航天自动控制研究所 Powered phase shutdown method for vehicle
CN105973516A (en) * 2015-12-11 2016-09-28 北京强度环境研究所 Pulsation thrust method for identification of solid rocket engine
CN107063244A (en) * 2017-04-14 2017-08-18 北京航天自动控制研究所 A kind of aircraft flight process analogy method
CN108536020A (en) * 2018-07-17 2018-09-14 哈尔滨工业大学 A kind of model reference adaptive sliding model control method for VTOL Reusable Launch Vehicles
CN109541941A (en) * 2018-11-19 2019-03-29 哈尔滨工业大学 A kind of adaptive augmentation anti-interference fault-tolerance approach for VTOL vehicle boosting flight
CN109798902A (en) * 2019-03-11 2019-05-24 北京星际荣耀空间科技有限公司 One kind being suitable for carrier rocket and enters the orbit modified interative guidance method
CN109857140A (en) * 2019-01-30 2019-06-07 北京星际荣耀空间科技有限公司 Carrier rocket pitch program angle calculation method, system, equipment and storage medium
CN110989669A (en) * 2019-12-11 2020-04-10 西安智翔防务技术有限公司 Online self-adaptive guidance algorithm for active section of multistage boosting gliding aircraft
WO2022095643A1 (en) * 2020-11-03 2022-05-12 蓝箭航天空间科技股份有限公司 Self-adaptive iterative guidance method and device for aerospace vehicle

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4275605A (en) * 1978-12-13 1981-06-30 Rockwell International Corporation Acceleration sensitive gyroscope stabilized platform
RU2102704C1 (en) * 1994-08-03 1998-01-20 Фирма "Саратовгазприборавтоматика" Device for determining and recording of geometrical parameters of pipelines
CN101723096A (en) * 2009-12-24 2010-06-09 北京航天自动控制研究所 Control method for reducing grand gesture perturbation under interative guidance
US20110181466A1 (en) * 2010-01-22 2011-07-28 Astrium Limited Receiver and method for authenticating satellite signals
US20120316824A1 (en) * 2010-01-15 2012-12-13 United States Government, As Represented By The Secretary Of The Navy Aerial bogey discrimination technique

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4275605A (en) * 1978-12-13 1981-06-30 Rockwell International Corporation Acceleration sensitive gyroscope stabilized platform
RU2102704C1 (en) * 1994-08-03 1998-01-20 Фирма "Саратовгазприборавтоматика" Device for determining and recording of geometrical parameters of pipelines
CN101723096A (en) * 2009-12-24 2010-06-09 北京航天自动控制研究所 Control method for reducing grand gesture perturbation under interative guidance
US20120316824A1 (en) * 2010-01-15 2012-12-13 United States Government, As Represented By The Secretary Of The Navy Aerial bogey discrimination technique
US20110181466A1 (en) * 2010-01-22 2011-07-28 Astrium Limited Receiver and method for authenticating satellite signals

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103863558B (en) * 2014-03-10 2015-12-30 北京航天自动控制研究所 A kind of powered-flight phase closedown method of vehicle
CN103863558A (en) * 2014-03-10 2014-06-18 北京航天自动控制研究所 Powered phase shutdown method for vehicle
CN105973516A (en) * 2015-12-11 2016-09-28 北京强度环境研究所 Pulsation thrust method for identification of solid rocket engine
CN105973516B (en) * 2015-12-11 2019-04-05 北京强度环境研究所 A method of the pulsation thrust of solid propellant rocket for identification
CN107063244B (en) * 2017-04-14 2019-07-12 北京航天自动控制研究所 A kind of aircraft flight process analogy method
CN107063244A (en) * 2017-04-14 2017-08-18 北京航天自动控制研究所 A kind of aircraft flight process analogy method
CN108536020A (en) * 2018-07-17 2018-09-14 哈尔滨工业大学 A kind of model reference adaptive sliding model control method for VTOL Reusable Launch Vehicles
CN109541941A (en) * 2018-11-19 2019-03-29 哈尔滨工业大学 A kind of adaptive augmentation anti-interference fault-tolerance approach for VTOL vehicle boosting flight
CN109541941B (en) * 2018-11-19 2020-07-21 哈尔滨工业大学 Self-adaptive amplification anti-interference fault-tolerant method for active section flight of vertical take-off and landing carrier
CN109857140A (en) * 2019-01-30 2019-06-07 北京星际荣耀空间科技有限公司 Carrier rocket pitch program angle calculation method, system, equipment and storage medium
CN109798902A (en) * 2019-03-11 2019-05-24 北京星际荣耀空间科技有限公司 One kind being suitable for carrier rocket and enters the orbit modified interative guidance method
CN109798902B (en) * 2019-03-11 2020-09-22 北京星际荣耀空间科技有限公司 Iterative guidance method suitable for carrier rocket in-orbit correction
CN110989669A (en) * 2019-12-11 2020-04-10 西安智翔防务技术有限公司 Online self-adaptive guidance algorithm for active section of multistage boosting gliding aircraft
WO2022095643A1 (en) * 2020-11-03 2022-05-12 蓝箭航天空间科技股份有限公司 Self-adaptive iterative guidance method and device for aerospace vehicle
GB2616530A (en) * 2020-11-03 2023-09-13 Landspace Science & Tech Co Ltd Self-adaptive iterative guidance method and device for aerospace vehicle
GB2616530B (en) * 2020-11-03 2024-07-17 Landspace Science & Tech Co Ltd Self-adaptive iterative guidance method and device for aerospace vehicle

Also Published As

Publication number Publication date
CN103592946B (en) 2016-02-10

Similar Documents

Publication Publication Date Title
CN103592946B (en) A kind of powered phase program self-adaptation longitudinal guidance method measured based on apparent acceleration
CN108536020B (en) A Model Reference Adaptive Sliding Mode Control Method for VTOL Reusable Vehicles
CN105045101B (en) A kind of mechanical arm servo-drive system full-order sliding mode control method based on extended state observer
CN104950898B (en) A kind of full rank non-singular terminal Sliding Mode Attitude control method of reentry vehicle
CN105202972B (en) Multi-missile cooperative engagement guidance method based on model predictive control technique
CN107908109B (en) A Trajectory Optimization Controller for Reentry Segment of Hypersonic Vehicle Based on Orthogonal Configuration Optimization
CN104682805B (en) Full-order sliding mode variable structure position servo control method for permanent magnet synchronous motor based on extended state observer
CN107102542B (en) A Robust Adaptive Non-Singular Terminal Sliding Mode Control Method for Automatic Train Running
CN107544245B (en) Discrete repetitive control method for motor servo system using first-order inertial attraction law with disturbance expansion compensation
CN107797448B (en) Motor position discrete repetitive control method adopting disturbance expansion compensation
CN101723096B (en) Control method for reducing grand gesture perturbation under iteration guidance
CN105549381A (en) Discrete repetition control method based on attraction law
CN105005311B (en) A kind of aircraft pitch channel angle of attack synchronization and tracking control method
CN103708045B (en) The on-line parameter discrimination method that a kind of lunar exploration airship great-jump-forward reenters
CN107491080A (en) A kind of Spacecraft Attitude Control tracked based on nonlinear feedback and differential
CN111007720A (en) Nonlinear self-healing control method for multi-sensor faults in hypersonic aircraft
CN104217041A (en) Multi-constrained online Gauss pseudo spectrum type terminal guidance method
CN103400035A (en) High-reliability method for rapidly forecasting rolling dynamic derivative of aircraft
CN103224023B (en) Phase plane self-adaptation control method based on characteristic model
CN105501467A (en) Carrier rocket exhaustion shutdown control method based on zero-range orientation
CN114020019A (en) Guidance method and device for aircraft
CN115079565B (en) Variable-coefficient guidance method, device and aircraft with drop angle constraints
CN115795687A (en) An adaptive design method for angular velocity limiter of air-breathing hypersonic vehicle
CN112034879A (en) Standard trajectory tracking guidance method based on height-range ratio
CN116562052B (en) Lateral winding flight method considering available overload constraint

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant