CN113448346B - Trajectory optimization method based on telemetering information - Google Patents

Trajectory optimization method based on telemetering information Download PDF

Info

Publication number
CN113448346B
CN113448346B CN202010230978.2A CN202010230978A CN113448346B CN 113448346 B CN113448346 B CN 113448346B CN 202010230978 A CN202010230978 A CN 202010230978A CN 113448346 B CN113448346 B CN 113448346B
Authority
CN
China
Prior art keywords
trajectory
angle
telemetering
acceleration
calculating
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.)
Active
Application number
CN202010230978.2A
Other languages
Chinese (zh)
Other versions
CN113448346A (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.)
No63729 Troops Pla
Original Assignee
No63729 Troops Pla
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 No63729 Troops Pla filed Critical No63729 Troops Pla
Priority to CN202010230978.2A priority Critical patent/CN113448346B/en
Publication of CN113448346A publication Critical patent/CN113448346A/en
Application granted granted Critical
Publication of CN113448346B publication Critical patent/CN113448346B/en
Active 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/10Simultaneous control of position or course in three dimensions
    • G05D1/107Simultaneous control of position or course in three dimensions specially adapted for missiles

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)
  • Navigation (AREA)
  • Complex Calculations (AREA)

Abstract

The utility model discloses a trajectory preferred method based on telemetering information, which is characterized by comprising the following steps of: (1) a correlation model of the telemetering angle and the telemetering speed; (2) calculating a correlation model of the remote vision acceleration and the external acceleration; (3) defining a window function; (4) outer trajectory and telemetering information windowing processing; (5) calculating the error of the same telemetering information in each outer ballistic window; (6) calculating the difference value of any trajectory j and the telemetering angle; (7) calculating the correlation difference between the remote vision acceleration and the external acceleration; (8) calculating the error of the same telemetering information in each outer ballistic window; (9) and (4) selecting an optimal trajectory. The trajectory optimization method based on the telemetering information can optimize the external trajectory through credible telemetering parameters, and improves trajectory optimization precision and stability.

Description

Trajectory optimization method based on telemetering information
Technical Field
The utility model relates to the technical field of aerospace outer trajectory processing, in particular to a trajectory optimization method based on telemetering information.
Background
In recent years, various new technologies, new states, and new requirements have been increasingly demanded for the accuracy of the center trajectory of the finger control. The command center has the advantages that the number of real-time trajectories formed by measurement information based on optical equipment, pulse radar equipment, speed measuring radar equipment, remote measuring equipment and the like is large, the total types of errors contained in all trajectories are large, and meanwhile, due to the fact that no unified quality evaluation standard exists, the difficulty in selecting the good trajectory and the bad trajectory is large, and an abnormal trajectory data correction method is lacked.
At present, the method for fusing the trajectory of the command center mainly comprises an optimization fusion algorithm and a weighting fusion algorithm, wherein the method comprises the steps of firstly determining trajectory sequencing by a manual optimization method to form a trajectory optimization sequence frame; secondly, selecting at most n ballistic trajectory parameters at the current moment from ballistic trajectory optimization sequence frames, carrying out difference with the parameters of the previous period to obtain absolute values, and sequencing the absolute values to remove the maximum value and the minimum value; and finally, determining the calculation results of the optimal fusion and the weighted fusion according to the number of the residual trajectories. This preferred method of ballistic trajectory has several major limitations: the selection is mainly carried out by depending on experience and manual operation, and the correctness and the reliability of the selection cannot be completely ensured; the preferred result mainly considers the external information and ignores the remote information.
Disclosure of Invention
The utility model aims to solve the problems, and provides a method based on telemetering information and capable of solving the problem of difficult ballistic preference outside a command center.
In order to achieve the purpose, the utility model adopts the following technical scheme, which comprises the following steps:
(1) correlation model of telemetry angle and telemetry speed: outer ballistic velocity of (v)x,vy,vz) Phi, psi, gamma, alpha and beta are respectively telemetered pitch angle, yaw angle, roll angle, attack angle and sideslip angle, and the calculation formula is as follows:
Figure GDA0002663939120000021
(2) calculating a correlation model of the remote vision acceleration and the external acceleration: (a)x,ay,az) For external acceleration measurement, (g)x,gy,gz) Is the component of gravitational acceleration in the launch coordinate system,
Figure GDA0002663939120000022
to look at the component of the acceleration in the transmit coordinate system, the calculation formula is as follows:
Figure GDA0002663939120000023
(3) defining a window function: t is the width of the window function, k is a positive integer, k Δ T represents that the function can be shifted along the time axis by an integer multiple of the time interval, and the calculation formula is as follows:
Figure GDA0002663939120000024
(4) outer trajectory, telemetry information windowing: multiplying the window function defined in step (3) with all outer ballistic, telemetry information;
(5) and calculating the error of the same telemetering information in each outer ballistic window: based on the results of the external ballistic and telemetry information windowing processing in the step (4), when k is equal to 0, the difference between the n external ballistic and telemetry angle functions is as follows:
Figure GDA0002663939120000025
in the above formula, f1(v),f2(v)…fn(v) And f1(a),f2(a)…fn(a) Respectively an external speed measuring function and an external acceleration measuring function which participate in the preference, g (phi, gamma, alpha, beta) and g (a)w,g0) Respectively a function of telemetric angle and a function of telemetric apparent acceleration, T0Is the zero point of takeoff, and delta T is the sampling interval;
(6) calculating the difference between any trajectory j and the telemetering angle as
Figure GDA0002663939120000026
Figure GDA0002663939120000027
In the above formula, M and P are respectively
Figure GDA0002663939120000028
The mean and variance of (A) are calculated as follows
Figure GDA0002663939120000031
(7) Calculating the correlation difference between the remote vision acceleration and the external acceleration: when k is 0, the difference between the n ballistic trajectories and the telemetry angle function is calculated as
Figure GDA0002663939120000032
The difference between any one trajectory j and the telemetry angle is as follows
Figure GDA0002663939120000033
Wherein M and P are each independently
Figure GDA0002663939120000034
The mean and variance of (A) are calculated as follows
Figure GDA0002663939120000035
(8) And (3) calculating the error of the same telemetric information in each outer ballistic window:
Figure GDA0002663939120000036
(9) selecting an optimal trajectory:
selecting an optimal trajectory based on the error calculation results of the steps (5) to (8) as follows:
p=min[ΔΨ1,ΔΨ2,…,ΔΨn]
where min [ ] indicates the minimum error.
The utility model has the beneficial effects that: the utility model provides a trajectory optimization method based on telemetering information, the telemetering information is used as a main parameter for representing the working state of an internal system on an arrow and has a certain coupling relation with external measuring information, the external measuring trajectory can be optimized through credible telemetering parameters, trajectory optimization precision and stability are improved, an association model of key telemetering parameters and external measuring speed parameters is established, a sliding window stepping trajectory optimization method based on the telemetering information is provided, and the external trajectory is optimized.
Drawings
Figure 1 is a schematic view of external ballistic windowing.
Fig. 2 is a diagram of telemetry information windowing.
Fig. 3 is a telemetry calculation result and each external velocity function curve, wherein a diagram (a) is a telemetry angle function curve, and diagrams (b) to (f) are external velocity function curves corresponding to trajectories 1 to 5, respectively.
Detailed Description
The following examples are presented to enable those skilled in the art to more fully understand the present invention and are not intended to limit the utility model in any way.
The utility model is further described with reference to fig. 1 and 2 and is implemented by the following technical solutions.
Taking a rocket as an example, the telemetry T0 time, the external T0 time and the time interval delta T0And measurement data such as a pitch angle phi, a yaw angle psi, a roll angle gamma, a flight time t, an apparent acceleration component, an aiming azimuth angle, a launching point latitude, an outer trajectory and the like are input.
Step 1), calculating according to the formula in the step 1) to obtain telemetering angle function curves at the left end and the right end of an equation and each external speed function curve;
fig. 3 shows the telemetry calculation result and each external velocity function curve, where fig. 3(a) is a telemetry angle function curve, and fig. 3(b) -3 (f) are external velocity function curves corresponding to trajectory 1-trajectory 5, respectively.
Step 2), calculating according to the formula in the step 2) to obtain an external acceleration curve;
step 3), giving specific definition of the window function;
step 4), windowing external trajectory and telemetering information;
step 5), calculating the error between each trajectory and the telemetering standard curve;
and 6) selecting the trajectory with the minimum error as the optimal trajectory.
Those skilled in the art will appreciate that the above embodiments are merely exemplary embodiments and that various changes, substitutions, and alterations can be made without departing from the spirit and scope of the utility model.

Claims (1)

1. A trajectory optimization method based on telemetering information is characterized by comprising the following steps:
(1) correlation model of telemetry angle and telemetry speed: an outer ballistic velocity of
Figure DEST_PATH_IMAGE002
),
Figure DEST_PATH_IMAGE004
Figure DEST_PATH_IMAGE006
Figure DEST_PATH_IMAGE008
Figure DEST_PATH_IMAGE010
Figure DEST_PATH_IMAGE012
Respectively telemetering a pitch angle, a yaw angle, a rolling angle, an attack angle and a sideslip angle, and calculating formulas are as follows:
Figure DEST_PATH_IMAGE014
(2) calculating a correlation model of the remote vision acceleration and the external acceleration: (
Figure DEST_PATH_IMAGE016
) For external measurement of acceleration: (
Figure DEST_PATH_IMAGE018
) As a component of gravitational acceleration in the launch coordinate system, ((ii))
Figure DEST_PATH_IMAGE020
) To look at the component of the acceleration in the transmit coordinate system, the calculation formula is as follows:
Figure DEST_PATH_IMAGE022
(3) defining a window function:
Figure DEST_PATH_IMAGE024
is the width of the window function, k is a positive integer,
Figure DEST_PATH_IMAGE026
the expression function may be shifted along the time axis by an integer multiple of the time interval, and is calculated as follows:
Figure DEST_PATH_IMAGE028
(4) outer trajectory, telemetry information windowing: multiplying the window function defined in step (3) with all outer ballistic, telemetry information;
(5) and calculating the error of the same telemetering information in each outer ballistic window: windowing the processing result based on the external trajectory and the telemetric information in the step (4), when
Figure DEST_PATH_IMAGE030
Then, the difference between the n outer trajectories and the telemetry angle function is:
Figure DEST_PATH_IMAGE032
in the above-mentioned formula,
Figure DEST_PATH_IMAGE034
and
Figure DEST_PATH_IMAGE036
respectively an external speed measurement function and an external acceleration measurement function which participate in the preference,
Figure DEST_PATH_IMAGE038
and
Figure DEST_PATH_IMAGE040
respectively a function of telemetric angle and a function of telemetric apparent acceleration, T0In order to take-off the zero point,
Figure DEST_PATH_IMAGE042
is the sampling interval;
(6) calculating the difference between any trajectory j and the telemetering angle as
Figure DEST_PATH_IMAGE044
In the above formula
Figure DEST_PATH_IMAGE046
And
Figure DEST_PATH_IMAGE048
are respectively as
Figure DEST_PATH_IMAGE050
The mean and variance of (A) are calculated as follows
Figure DEST_PATH_IMAGE052
(7) Calculating the correlation difference between the remote vision acceleration and the external acceleration: when in use
Figure 638928DEST_PATH_IMAGE030
Then, the difference value of n ballistic trajectories and the telemetering angle function is calculated as
Figure DEST_PATH_IMAGE054
The difference between any one trajectory j and the telemetry angle is as follows
Figure DEST_PATH_IMAGE056
Wherein
Figure 87227DEST_PATH_IMAGE046
And
Figure 132543DEST_PATH_IMAGE048
are respectively as
Figure DEST_PATH_IMAGE058
The mean and variance of (A) are calculated as follows
Figure DEST_PATH_IMAGE060
(8) And (3) calculating the error of the same telemetric information in each outer ballistic window:
Figure DEST_PATH_IMAGE062
(9) selecting an optimal trajectory:
selecting an optimal trajectory based on the error calculation results of the steps (5) to (8) as follows:
Figure DEST_PATH_IMAGE064
where min [ ] indicates the minimum error.
CN202010230978.2A 2020-03-27 2020-03-27 Trajectory optimization method based on telemetering information Active CN113448346B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010230978.2A CN113448346B (en) 2020-03-27 2020-03-27 Trajectory optimization method based on telemetering information

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010230978.2A CN113448346B (en) 2020-03-27 2020-03-27 Trajectory optimization method based on telemetering information

Publications (2)

Publication Number Publication Date
CN113448346A CN113448346A (en) 2021-09-28
CN113448346B true CN113448346B (en) 2022-04-26

Family

ID=77807897

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010230978.2A Active CN113448346B (en) 2020-03-27 2020-03-27 Trajectory optimization method based on telemetering information

Country Status (1)

Country Link
CN (1) CN113448346B (en)

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0692425A1 (en) * 1994-07-14 1996-01-17 Orbital Sciences Corporation Method and system for formationkeeping between orbiting spacecraft by varying their ballistic coefficients
US5631654A (en) * 1996-02-05 1997-05-20 The Regents Of The University Of California Ballistic projectile trajectory determining system
CN104571125B (en) * 2014-12-18 2016-01-27 北京控制工程研究所 A kind ofly utilize the multiple control method returning condition of normal trajectory reply
CN104950301A (en) * 2015-04-09 2015-09-30 北京寰亚翔宇科技有限责任公司 Exterior trajectory measurement method
CN107101649B (en) * 2017-05-25 2019-08-23 北京航天自动控制研究所 A kind of in-orbit error separating method of spacecraft Guidance instrumentation
CN110674567A (en) * 2019-08-23 2020-01-10 中国人民解放军63729部队 Rocket power situation judgment method based on external acceleration
CN110672060B (en) * 2019-08-23 2022-09-16 中国人民解放军63729部队 Rocket attitude angle condition judgment method based on external measurement speed

Also Published As

Publication number Publication date
CN113448346A (en) 2021-09-28

Similar Documents

Publication Publication Date Title
CN109343550B (en) Spacecraft angular velocity estimation method based on rolling time domain estimation
CN110320925B (en) Four-rotor aircraft safety control method based on high-order disturbance observer
CN103712598B (en) Attitude determination method of small unmanned aerial vehicle
CN110926460A (en) Uwb positioning abnormal value processing method based on IMU
CN112148036B (en) Bilateral tracking control method of fixed time estimator of networked robot system
CN111123700B (en) Constraint full-course satisfied optimal control system for obstacle-detouring flight of hypersonic aircraft
CN110595434B (en) Quaternion fusion attitude estimation method based on MEMS sensor
CN114527795B (en) Aircraft intelligent control method based on increment online learning
CN113448346B (en) Trajectory optimization method based on telemetering information
CN112000127A (en) Reverse-step-method-based aircraft lateral combined control method
CN111045447A (en) High-precision hypersonic aircraft trajectory optimization multi-scale optimal control system
CN114964226A (en) Noise adaptive strong tracking extended Kalman filter four-rotor attitude resolving method
CN115129072A (en) Terminal sliding mode control method under position tracking deviation constraint of fixed wing unmanned aerial vehicle
CN116923730B (en) Spacecraft attitude active fault-tolerant control method with self-adjusting preset performance constraint
CN111413996B (en) Four-rotor performance-guaranteeing trajectory tracking control method based on event-triggered ESO
CN111061294B (en) Trajectory optimization self-adaptive optimal control system for nonstationary hypersonic aircraft
CN112506209B (en) Reentry vehicle prediction control method based on self-adaptive prediction period
CN110285711B (en) Guidance ammunition flight attitude correction simulation system and method based on information physical system
CN113110428A (en) Carrier-based aircraft landing fixed time trajectory tracking method based on limited backstepping control
CN114111804A (en) Multi-source multi-category measurement data time zero high-precision alignment method for carrier rocket
CN116700306B (en) Integrated guidance control method for strapdown guided aircraft
CN115615437B (en) Factor graph integrated navigation method
CN113552895B (en) Outer trajectory break point correction method based on remote vision acceleration
CN111338364B (en) High-precision controller for optimizing trajectory of hypersonic aerocraft with quick response
CN114578838B (en) Reinforced learning active disturbance rejection attitude control method suitable for aircrafts of various configurations

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