CN114111805A - Multi-source multi-category measurement data position reference high-precision alignment method for carrier rocket - Google Patents

Multi-source multi-category measurement data position reference high-precision alignment method for carrier rocket Download PDF

Info

Publication number
CN114111805A
CN114111805A CN202111132200.9A CN202111132200A CN114111805A CN 114111805 A CN114111805 A CN 114111805A CN 202111132200 A CN202111132200 A CN 202111132200A CN 114111805 A CN114111805 A CN 114111805A
Authority
CN
China
Prior art keywords
measurement
difference
formula
radar
optical
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.)
Pending
Application number
CN202111132200.9A
Other languages
Chinese (zh)
Inventor
陈德明
王婷婷
陈永健
唐志强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
UNIT 63620 OF PLA
Original Assignee
UNIT 63620 OF 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 UNIT 63620 OF PLA filed Critical UNIT 63620 OF PLA
Priority to CN202111132200.9A priority Critical patent/CN114111805A/en
Publication of CN114111805A publication Critical patent/CN114111805A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/24Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • G01C21/165Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments

Abstract

The invention belongs to the field of space launching and the field of space measurement and control, and discloses a position reference correction method for multi-source and multi-category measurement data of a carrier rocket. According to the technical approach provided by the invention, the actual measurement precision of the measurement position deviation of the radar measurement system, the optical measurement system, the satellite navigation measurement system and the inertial navigation measurement system in the projectile coordinate system determines the correction precision of the position measurement elements of the three measurement systems, so that the problem of high-precision alignment of the multi-source multi-type measurement data position reference of the carrier rocket can be effectively solved.

Description

Multi-source multi-category measurement data position reference high-precision alignment method for carrier rocket
Technical Field
The invention belongs to the field of space launching and the field of space measurement and control, and relates to a position reference high-precision alignment method for multi-source and multi-type measurement data of a carrier rocket.
Background
In space launch measurement data, there are three main categories of measurement data related to launch orbits: the method comprises the following steps that firstly, various external ballistic measurement data, namely measurement data of a radar measurement system (including a pulse radar and a high-precision speed measuring radar) and an optical measurement system; the second is inertial navigation ballistic data, namely output pulse data of an inertial measurement device (an accelerometer and a rate gyro) of the carrier rocket, which is downloaded by a telemetry system on the carrier rocket, and ballistic data (namely inertial navigation ballistic data) which is resolved by a navigation computer according to the pulse data, wherein the output pulse data of the inertial measurement device and the inertial navigation ballistic data have a strict corresponding relation, and can be viewed equally in application, and for convenience of application, the description in the following generally only refers to the inertial navigation ballistic data except for special description; and thirdly, satellite navigation measurement data downloaded by the carrier rocket through the rocket-borne telemetry system.
In the three types of measurement data, the measurement position of the radar measurement system is the antenna position of a transponder installed on the carrier rocket (the measurement position of the pulse radar is the antenna of the pulse radar transponder, the measurement position of the high-precision speed measuring radar is the speed measuring radar transponder), and the measurement position of the optical measurement system is the flame center position of the carrier rocket engine; the measurement position of the satellite navigation measurement system is the position of a satellite navigation information receiving antenna arranged on the carrier rocket; the measurement position of the inertial navigation measurement system is the measurement position of the inertial measurement device of the carrier rocket, namely the origin of the measurement coordinate system of the navigation platform. The inconsistency of the position reference can cause position difference and speed difference among the three types of measurement data, the existence of the difference not only directly influences the high-precision fusion of the three types of measurement data, but also inevitably influences the comprehensive analysis of the key process, the basic performance of key actions and the comprehensive performance of the flight process of the carrier rocket by utilizing inertial navigation ballistic data, outer ballistic measurement data, satellite navigation measurement data and various types of telemetering data of various systems of the carrier rocket.
In the aerospace launching task, the position alignment refers to aligning the measurement position of the radar measurement system, the measurement position of the optical measurement system and the measurement position of the satellite navigation measurement system to the measurement position of the inertial navigation measurement system, namely, correcting the position measurement data of other 3 types of measurement systems by taking the measurement position of the inertial navigation measurement system as a reference so as to keep the same position reference as the inertial navigation ballistic data. The current commonly used position alignment measure is to correct the difference between the engine flame center position measured by the optical measurement system and the measurement position of the inertial navigation measurement system, i.e. the photometric position correction. The difference between the measurement positions of the pulse radar and the speed measuring radar and the measurement position of the inertial navigation measurement system is not corrected generally, and the difference of speed measuring radar speed measurement metadata caused by the inconsistency of the measurement positions is never corrected.
The primary condition of high-precision fusion processing of multiple types of measurement data from multiple sources of a launch orbit of a carrier rocket is that the consistency of measurement positions is required. At present, no alignment principle and high-precision correction method for ensuring the consistency of position references of three types of measurement data are available.
Disclosure of Invention
The invention aims to provide a high-precision alignment method for position reference of multi-source multi-type measurement data of a carrier rocket, which utilizes the difference between the measurement position of an inertial navigation measurement system and the measurement position of other measurement systems, corrects position measurement elements corresponding to the measurement data of other measurement systems according to the measurement geometric relation, and corrects all the measurement data related to the position in the measurement data of other measurement systems to the measurement position of the inertial navigation measurement system, thereby realizing the consistency of the multi-source multi-type measurement data on the spatial position and providing technical support for the high-precision fusion processing and the fine application of the multi-source multi-type measurement data.
In order to achieve the purpose, the technical scheme of the invention is as follows:
a multi-source multi-category measurement data position reference high-precision alignment method for a carrier rocket comprises the following steps:
step one, converting the position difference of a projectile system into a navigation system position difference;
step two, converting the navigation system position difference into a geocentric system position difference;
step three, correcting position measurement data by different measurement systems;
correcting the position measurement data of the satellite navigation measurement system, and executing the step four;
for the correction of the position measurement data of the radar measurement system, executing the fifth step to the seventh step;
for correcting the position measurement data of the optical measurement system, executing the step eight to the step ten;
correcting the satellite navigation position measurement data based on the geocentric position difference;
step five, converting the position difference of the geocentric system into the position difference of a measuring system of the radar measuring system;
step six, converting the position difference of a measurement system of the radar measurement system into a position measurement element difference of the radar measurement system;
step seven, correcting the position measurement metadata of the radar measurement system based on the position measurement metadata difference of the radar measurement system;
step eight, converting the position difference of the geocentric system into the position difference of the measuring system of the optical measuring system;
ninth, the position difference of the measurement system of the optical measurement system is converted into the position measurement element difference of the optical measurement system;
and step ten, correcting the position measurement metadata of the optical measurement system based on the position measurement metadata difference of the optical measurement system.
Further, the first step further comprises:
converting the projectile system position difference d into a navigation system position difference according to the formula (1)
Figure BDA0003275548770000031
Figure BDA0003275548770000032
In the formula, TG1The direction conversion matrix for carrying a rocket projectile coordinate system to a navigation platform measurement coordinate system comprises the following specific calculation formula:
Figure BDA0003275548770000033
in the formula (I), the compound is shown in the specification,
Figure BDA0003275548770000041
for the measurement of the attitude angle of the projectile coordinate system relative to the navigation platform, this data can be obtained directly from the computer word telemetry parameters of the launch vehicle.
Further, the second step further comprises:
differentiating the position of the navigation system according to the formula (2)
Figure BDA0003275548770000042
Converted into the position difference of the geocentric system
Figure BDA0003275548770000043
Figure BDA0003275548770000044
In the formula, TOGThe direction conversion matrix from the carrier rocket navigation measurement coordinate system to the geocentric coordinate system is calculated by the following specific formula:
Figure BDA0003275548770000045
where t is the time of flight of the launch vehicle relative to the time of initiation of the navigation computation, ωeIs the angular rate of rotation of the earth, B0、L0、A0The geographical latitude, the geographical longitude and the launching azimuth of the launching point of the carrier rocket are respectively.
Further, the fourth step further comprises:
the position difference of the geocentric system obtained in the second step
Figure BDA0003275548770000046
Calculating satellite navigation position measurement data [ x ] according to equation (3)sys zs]TCorrection value of
Figure BDA0003275548770000047
Figure BDA0003275548770000048
Further, the fifth step further comprises:
the difference in the position of the earth center system according to the formula (4)
Figure BDA0003275548770000049
Converting into position difference of measurement system of radar measurement system
Figure BDA00032755487700000410
Figure BDA00032755487700000411
In the formula, TMOThe direction conversion matrix from the geocentric coordinate system to the radar measurement system measurement coordinate system is specifically calculated as follows:
Figure BDA0003275548770000051
wherein, B and L are the geographical latitude and the geographical longitude of the radar measuring system respectively.
Further, the sixth step further comprises:
measuring the position difference of the system of radar measurement according to the formula (5)
Figure BDA0003275548770000052
Converting into radar measurement system position measured element difference [ dR dA dE]T
Figure BDA0003275548770000053
In the formula, TAMThe position difference of a measuring system of the radar measuring system to a position measuring element of the radar measuring systemThe specific calculation formula of the conversion matrix of the difference is as follows:
Figure BDA0003275548770000054
wherein [ RA E ]]TThe position measuring element of the radar measuring system is shown, wherein R represents an oblique distance measuring element, A represents an azimuth angle measuring element, and E represents a high and low angle measuring element.
Further, the seventh step further comprises:
according to the position measured element difference [ dR dA dE ] of the radar measuring system obtained in the fifth step]TCalculating radar measurement system position measurement metadata [ RA E ] according to the formula (7)]TCorrection value of
Figure BDA0003275548770000055
Figure BDA0003275548770000056
Further, the step eight further comprises:
differentiating the position of the earth's center system according to the formula (8)
Figure BDA0003275548770000061
Converting into position differences of measurement systems of optical measurement systems
Figure BDA0003275548770000062
Figure BDA0003275548770000063
In the formula (I), the compound is shown in the specification,
Figure BDA0003275548770000064
the direction conversion matrix from the geocentric coordinate system to the measurement coordinate system of the optical measurement system comprises the following specific calculation formula:
Figure BDA0003275548770000065
in the formula, BGAnd LGRespectively the geographical latitude and the geographical longitude of the optical measurement system.
Further, the ninth step further comprises:
9.1 estimating the slant-distance measuring element of the optical measuring system
Figure BDA0003275548770000066
Firstly, according to the position data sequence of the carrier rocket measured by the satellite navigation measuring system within 10s before the t moment
Figure BDA0003275548770000067
The sequence number i corresponds to the time of flight t of the launch vehicleiI is 1,2, …, n, calculating the position data sequence of the carrier rocket in the optical measurement system measurement coordinate system according to the following formula:
Figure BDA0003275548770000068
in the formula (I), the compound is shown in the specification,
Figure BDA0003275548770000069
is a direction conversion matrix from a geocentric coordinate system to an optical measurement system measurement coordinate system;
secondly, according to the position data sequence of the carrier rocket in the optical measurement system measuring coordinate system
Figure BDA00032755487700000610
Calculating the slope distance sequence of the carrier rocket in the measuring coordinate system of the optical measuring system according to the following formula
Figure BDA00032755487700000611
Estimating the slant distance measuring element of the optical measuring system according to the formula (9) by using a second-order polynomial end point smoothing method
Figure BDA00032755487700000612
Figure BDA00032755487700000613
In the formula (I), the compound is shown in the specification,
Figure BDA0003275548770000071
9.2 measuring the position difference of the system by the optical measuring system according to the formula (10)
Figure BDA0003275548770000072
Conversion into position measurement element difference of optical measurement system
Figure BDA0003275548770000073
Figure BDA0003275548770000074
In the formula (I), the compound is shown in the specification,
Figure BDA0003275548770000075
the method is a conversion matrix from the position difference of a measurement system of an optical measurement system to the position measurement element difference of the optical measurement system, and the specific calculation formula is as follows:
Figure BDA0003275548770000076
wherein (A E) is an angle measurement unit of the optical measurement system.
Further, the step of decimal further comprises:
the position measured element difference of the optical measuring system obtained according to the step eight
Figure BDA0003275548770000077
Calculating an optical measuring system according to equation (11)Location measurement metadata
Figure BDA0003275548770000078
Correction value of
Figure BDA0003275548770000079
Figure BDA00032755487700000710
The invention has the advantages that:
according to the basic method of coordinate conversion and the basic characteristics of space launching measurement data, the basic fact that the positions of a satellite navigation measurement antenna and a radar responder antenna are fixed on an arrow coordinate system is fully utilized, the basic transfer relationship of the influence of the difference between the outer ballistic trajectory measurement position and the inertial navigation measurement position on various measurement system position measurement elements is established, the correction method of the satellite navigation measurement system position measurement data, the correction method of the radar measurement system position measurement element measurement data and the correction method of the optical measurement system angle measurement element measurement data are provided, and the practical problem of position consistency among various measurement data in a space launching task is well solved.
Drawings
Fig. 1 is a schematic view of a position alignment approach.
Fig. 2 is a schematic diagram of a main process of position alignment.
In fig. 2, 4 parts of the square frame mark mainly implement coordinate transformation of the measurement position difference between various types of measurement data, and mainly solve the principle and method problem of position measurement element correction; the 3 parts of the oval frame mark mainly realize the specific correction of the difference between the measurement position of the 3-type measurement system and the measurement position of the inertial navigation measurement system, and mainly solve the specific method and process problems of the correction of the position measurement element.
Detailed Description
As described above, the position measurement data of the other 3 types of measurement systems is corrected with the measurement position of the inertial navigation measurement system as a reference.
The main notations used in the present invention are described below:
Figure BDA0003275548770000081
Figure BDA0003275548770000091
the method for aligning the position reference of the multi-source multi-category measurement data of the carrier rocket at high precision comprises the following steps:
step one, converting the position difference of a projectile system into a navigation system position difference;
the position difference of the elastic system is determined according to the formula (1)
Figure BDA0003275548770000092
Converting into a navigation system position difference
Figure BDA0003275548770000093
Figure BDA0003275548770000101
In the formula, TG1The direction conversion matrix for carrying a rocket projectile coordinate system to a navigation platform measurement coordinate system comprises the following specific calculation formula:
Figure BDA0003275548770000102
in the formula (I), the compound is shown in the specification,
Figure BDA0003275548770000103
for the measurement of the attitude angle of the projectile coordinate system relative to the navigation platform, this data can be obtained directly from the computer word telemetry parameters of the launch vehicle.
Step two, converting the navigation system position difference into a geocentric system position difference;
differentiating the position of the navigation system according to the formula (2)
Figure BDA0003275548770000104
Converted into the position difference of the geocentric system
Figure BDA0003275548770000105
Figure BDA0003275548770000106
In the formula, TOGThe direction conversion matrix from the carrier rocket navigation measurement coordinate system to the geocentric coordinate system is calculated by the following specific formula:
Figure BDA0003275548770000107
where t is the time of flight of the launch vehicle relative to the time of initiation of the navigation computation, ωeIs the angular rate of rotation of the earth, B0、L0、A0The geographical latitude, the geographical longitude and the launching azimuth of the launching point of the carrier rocket are respectively.
Step one and step two are true for the satellite navigation measurement system, the radar measurement system and the optical measurement system. The correction of the position measurement data is carried out for different measurement systems, starting from the following steps:
step three, correcting position measurement data by different measurement systems;
correcting the position measurement data of the satellite navigation measurement system, and executing the step four;
for the correction of the position measurement data of the radar measurement system, executing the fifth step to the seventh step;
for correcting the position measurement data of the optical measurement system, executing the step eight to the step ten;
correcting the satellite navigation position measurement data based on the geocentric position difference;
the position difference of the geocentric system obtained in the second step
Figure BDA0003275548770000111
Calculating satellite navigation position measurement data [ x ] according to equation (3)sys zs]TCorrection value of
Figure BDA0003275548770000112
Figure BDA0003275548770000113
Step five, converting the position difference of the geocentric system into the position difference of a measuring system of the radar measuring system;
the difference in the position of the earth center system according to the formula (4)
Figure BDA0003275548770000114
Converting into position difference of measurement system of radar measurement system
Figure BDA0003275548770000115
Figure BDA0003275548770000116
In the formula, TMOThe direction conversion matrix from the geocentric coordinate system to the radar measurement system measurement coordinate system is specifically calculated as follows:
Figure BDA0003275548770000117
wherein, B and L are the geographical latitude and the geographical longitude of the radar measuring system respectively.
Step six, converting the position difference of a measurement system of the radar measurement system into a position measurement element difference of the radar measurement system;
measuring the position difference of the system of radar measurement according to the formula (5)
Figure BDA0003275548770000119
Converting into radar measurement system position measured element difference [ dR dA dE]T
Figure BDA0003275548770000118
In the formula, TAMThe method is a conversion matrix from the position difference of a radar measurement system to the position measurement element difference of the radar measurement system, and the specific calculation formula is as follows:
Figure BDA0003275548770000121
wherein [ RA E ]]TThe position measuring element of the radar measuring system is shown, wherein R represents an oblique distance measuring element, A represents an azimuth angle measuring element, and E represents a high and low angle measuring element.
Specifically, the derivation process of the formula (5) is as follows:
position coordinate vector [ x ] of carrier rocket in radar measurement coordinate systemM yM zM]TPosition measurement element [ RA E ] of radar measurement system]TThe following basic relationships exist:
Figure BDA0003275548770000122
fully differentiating the above formula to obtain
Figure BDA0003275548770000123
By performing an equivalent transformation on the above formula, then
Figure BDA0003275548770000124
To obtain the formula (5).
Step seven, correcting the position measurement metadata of the radar measurement system based on the position measurement metadata difference of the radar measurement system;
radar measurements from step fiveSystem position delta [ dR dA dE ]]TCalculating radar measurement system position measurement metadata [ RA E ] according to the formula (7)]TCorrection value of
Figure BDA0003275548770000131
Figure BDA0003275548770000132
Step eight, converting the position difference of the geocentric system into the position difference of the measuring system of the optical measuring system;
differentiating the position of the earth's center system according to the formula (8)
Figure BDA0003275548770000133
Converting into position differences of measurement systems of optical measurement systems
Figure BDA0003275548770000134
Figure BDA0003275548770000135
In the formula (I), the compound is shown in the specification,
Figure BDA0003275548770000136
the direction conversion matrix from the geocentric coordinate system to the measurement coordinate system of the optical measurement system comprises the following specific calculation formula:
Figure BDA0003275548770000137
in the formula, BGAnd LGRespectively the geographical latitude and the geographical longitude of the optical measurement system.
Ninth, the position difference of the measurement system of the optical measurement system is converted into the position measurement element difference of the optical measurement system;
the method comprises the following steps of:
9.1 estimating the slant-distance measuring element of the optical measuring system
Figure BDA0003275548770000138
For any measurement time, not recording the time t, and measuring a position data sequence of the carrier rocket by a satellite navigation measurement system within 10s before the time t
Figure BDA0003275548770000139
The sequence number i corresponds to the time of flight t of the launch vehicleiI is 1,2, …, n, calculating the position data sequence of the carrier rocket in the optical measurement system measurement coordinate system according to the following formula:
Figure BDA00032755487700001310
in the formula (I), the compound is shown in the specification,
Figure BDA00032755487700001311
is a direction conversion matrix from the geocentric coordinate system to the measurement coordinate system of the optical measurement system.
Secondly, according to the position data sequence of the carrier rocket in the optical measurement system measuring coordinate system
Figure BDA0003275548770000141
Calculating the slope distance sequence of the carrier rocket in the measuring coordinate system of the optical measuring system according to the following formula
Figure BDA0003275548770000142
Estimating the slant distance measuring element of the optical measuring system according to the formula (9) by using a second-order polynomial end point smoothing method
Figure BDA0003275548770000143
Figure BDA0003275548770000144
In the formula (I), the compound is shown in the specification,
Figure BDA0003275548770000145
9.2 measuring the position difference of the system by the optical measuring system according to the formula (10)
Figure BDA0003275548770000146
Conversion into position measurement element difference of optical measurement system
Figure BDA0003275548770000147
Figure BDA0003275548770000148
In the formula (I), the compound is shown in the specification,
Figure BDA0003275548770000149
the method is a conversion matrix from the position difference of a measurement system of an optical measurement system to the position measurement element difference of the optical measurement system, and the specific calculation formula is as follows:
Figure BDA00032755487700001410
wherein (A E) is an angle measurement unit of the optical measurement system.
Specifically, the derivation process of the formula (10) is as follows:
position coordinate vector [ x ] of carrier rocket in optical measurement system measurement coordinate systemM yM zM]TPosition measuring unit for optical measuring system
Figure BDA0003275548770000151
The following basic relationships exist:
Figure BDA0003275548770000152
fully differentiating the above formula to obtain
Figure BDA0003275548770000153
By performing an equivalent transformation on the above formula, then
Figure BDA0003275548770000154
To obtain the formula (10).
Step ten, correcting the position measurement metadata of the optical measurement system based on the position measurement metadata difference of the optical measurement system;
the position measured element difference of the optical measuring system obtained according to the step eight
Figure BDA0003275548770000155
Calculating optical measurement system position measurement metadata according to equation (11)
Figure BDA0003275548770000156
Correction value of
Figure BDA0003275548770000157
Figure BDA0003275548770000158
The position reference correction method for the multi-source multi-category measurement data of the carrier rocket provided by the invention utilizes the difference between the measurement position of the inertial navigation measurement system and the measurement position of other measurement systems, and corrects the measurement element deviation caused by correcting the measurement position deviation, thereby realizing the consistency of all position measurement data on the spatial position. According to the technical approach provided by the invention, the actual measurement precision of the measurement position deviation of the radar measurement system, the optical measurement system, the satellite navigation measurement system and the inertial navigation measurement system in the projectile coordinate system determines the correction precision of the position measurement elements of the three measurement systems, so that the problem of high-precision alignment of the multi-source multi-type measurement data position reference of the carrier rocket can be effectively solved.

Claims (10)

1. A multi-source multi-category measurement data position reference high-precision alignment method for a carrier rocket is characterized by comprising the following steps:
step one, converting the position difference of a projectile system into a navigation system position difference;
step two, converting the navigation system position difference into a geocentric system position difference;
step three, correcting position measurement data by different measurement systems;
correcting the position measurement data of the satellite navigation measurement system, and executing the step four;
for the correction of the position measurement data of the radar measurement system, executing the fifth step to the seventh step;
for correcting the position measurement data of the optical measurement system, executing the step eight to the step ten;
correcting the satellite navigation position measurement data based on the geocentric position difference;
step five, converting the position difference of the geocentric system into the position difference of a measuring system of the radar measuring system;
step six, converting the position difference of a measurement system of the radar measurement system into a position measurement element difference of the radar measurement system;
step seven, correcting the position measurement metadata of the radar measurement system based on the position measurement metadata difference of the radar measurement system;
step eight, converting the position difference of the geocentric system into the position difference of the measuring system of the optical measuring system;
ninth, the position difference of the measurement system of the optical measurement system is converted into the position measurement element difference of the optical measurement system;
and step ten, correcting the position measurement metadata of the optical measurement system based on the position measurement metadata difference of the optical measurement system.
2. The launch vehicle multi-source multi-category measurement data position reference high precision alignment method of claim 1, wherein said step one further comprises:
the position difference of the elastic system is determined according to the formula (1)
Figure FDA0003275548760000011
Converting into a navigation system position difference
Figure FDA0003275548760000012
Figure FDA0003275548760000021
In the formula, TG1The direction conversion matrix for carrying a rocket projectile coordinate system to a navigation platform measurement coordinate system comprises the following specific calculation formula:
Figure FDA0003275548760000022
in the formula (I), the compound is shown in the specification,
Figure FDA0003275548760000023
for the measurement of the attitude angle of the projectile coordinate system relative to the navigation platform, this data can be obtained directly from the computer word telemetry parameters of the launch vehicle.
3. The launch vehicle multi-source multi-category measurement data position reference high precision alignment method of claim 2, wherein said step two further comprises:
differentiating the position of the navigation system according to the formula (2)
Figure FDA0003275548760000024
Converted into the position difference of the geocentric system
Figure FDA0003275548760000025
Figure FDA0003275548760000026
In the formula, TOGThe direction conversion matrix from the carrier rocket navigation measurement coordinate system to the geocentric coordinate system is calculated by the following specific formula:
Figure FDA0003275548760000027
where t is the time of flight of the launch vehicle relative to the time of initiation of the navigation computation, ωeIs the angular rate of rotation of the earth, B0、L0、A0The geographical latitude, the geographical longitude and the launching azimuth of the launching point of the carrier rocket are respectively.
4. The launch vehicle multi-source multi-category measurement data position reference high precision alignment method of claim 3, wherein said step four further comprises:
the position difference of the geocentric system obtained in the second step
Figure FDA0003275548760000028
Calculating satellite navigation position measurement data [ x ] according to equation (3)s ys zs]TCorrection value of
Figure FDA0003275548760000029
Figure FDA0003275548760000031
5. The launch vehicle multi-source multi-category measurement data position reference high precision alignment method of claim 4, wherein said step five further comprises:
the difference in the position of the earth center system according to the formula (4)
Figure FDA0003275548760000032
Conversion into a measurement system position of a radar measurement systemSet difference
Figure FDA0003275548760000033
Figure FDA0003275548760000034
In the formula, TMOThe direction conversion matrix from the geocentric coordinate system to the radar measurement system measurement coordinate system is specifically calculated as follows:
Figure FDA0003275548760000035
wherein, B and L are the geographical latitude and the geographical longitude of the radar measuring system respectively.
6. The launch vehicle multi-source multi-category measurement data position reference high precision alignment method of claim 5, wherein said sixth step further comprises:
measuring the position difference of the system of radar measurement according to the formula (5)
Figure FDA0003275548760000036
Converting into radar measurement system position measured element difference [ dR dA dE]T
Figure FDA0003275548760000037
In the formula, TAMThe method is a conversion matrix from the position difference of a radar measurement system to the position measurement element difference of the radar measurement system, and the specific calculation formula is as follows:
Figure FDA0003275548760000041
wherein [ RA E ]]TPosition measuring cell for radar measuring system, whichAnd R in the middle represents an oblique distance measuring element, A represents an azimuth angle measuring element, and E represents a high and low angle measuring element.
7. The launch vehicle multi-source multi-category measurement data position reference high precision alignment method of claim 6, wherein said seventh step further comprises:
according to the position measured element difference [ dR dA dE ] of the radar measuring system obtained in the fifth step]TCalculating radar measurement system position measurement metadata [ RA E ] according to the formula (7)]TCorrection value of
Figure FDA0003275548760000042
Figure FDA0003275548760000043
8. The launch vehicle multi-source multi-category measurement data position reference high accuracy alignment method of claim 7, wherein said step eight further comprises:
differentiating the position of the earth's center system according to the formula (8)
Figure FDA0003275548760000044
Converting into position differences of measurement systems of optical measurement systems
Figure FDA0003275548760000045
Figure FDA0003275548760000046
In the formula (I), the compound is shown in the specification,
Figure FDA0003275548760000047
the direction conversion matrix from the geocentric coordinate system to the measurement coordinate system of the optical measurement system comprises the following specific calculation formula:
Figure FDA0003275548760000048
in the formula, BGAnd LGRespectively the geographical latitude and the geographical longitude of the optical measurement system.
9. The launch vehicle multi-source multi-category measurement data position reference high precision alignment method of claim 8, wherein said ninth step further comprises:
9.1 estimating the slant-distance measuring element of the optical measuring system
Figure FDA0003275548760000051
Firstly, according to the position data sequence of the carrier rocket measured by the satellite navigation measuring system within 10s before the t moment
Figure FDA0003275548760000052
The sequence number i corresponds to the time of flight t of the launch vehicleiI is 1,2, …, n, calculating the position data sequence of the carrier rocket in the optical measurement system measurement coordinate system according to the following formula:
Figure FDA0003275548760000053
in the formula (I), the compound is shown in the specification,
Figure FDA0003275548760000054
is a direction conversion matrix from a geocentric coordinate system to an optical measurement system measurement coordinate system;
secondly, according to the position data sequence of the carrier rocket in the optical measurement system measuring coordinate system
Figure FDA0003275548760000055
Calculating the slope distance sequence of the carrier rocket in the measuring coordinate system of the optical measuring system according to the following formula
Figure FDA0003275548760000056
Estimating the slant distance measuring element of the optical measuring system according to the formula (9) by using a second-order polynomial end point smoothing method
Figure FDA0003275548760000057
Figure FDA0003275548760000058
In the formula (I), the compound is shown in the specification,
Figure FDA0003275548760000059
9.2 measuring the position difference of the system by the optical measuring system according to the formula (10)
Figure FDA00032755487600000510
Conversion into position measurement element difference of optical measurement system
Figure FDA00032755487600000511
Figure FDA00032755487600000512
In the formula (I), the compound is shown in the specification,
Figure FDA0003275548760000061
the method is a conversion matrix from the position difference of a measurement system of an optical measurement system to the position measurement element difference of the optical measurement system, and the specific calculation formula is as follows:
Figure FDA0003275548760000062
wherein (A E) is an angle measurement unit of the optical measurement system.
10. The high-precision alignment method for position references of multi-source multi-class measurement data of a launch vehicle of claim 9, wherein said step of decimal further comprises:
the position measured element difference of the optical measuring system obtained according to the step eight
Figure FDA0003275548760000063
Calculating optical measurement system position measurement metadata according to equation (11)
Figure FDA0003275548760000064
Correction value of
Figure FDA0003275548760000065
Figure FDA0003275548760000066
CN202111132200.9A 2021-09-23 2021-09-23 Multi-source multi-category measurement data position reference high-precision alignment method for carrier rocket Pending CN114111805A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111132200.9A CN114111805A (en) 2021-09-23 2021-09-23 Multi-source multi-category measurement data position reference high-precision alignment method for carrier rocket

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111132200.9A CN114111805A (en) 2021-09-23 2021-09-23 Multi-source multi-category measurement data position reference high-precision alignment method for carrier rocket

Publications (1)

Publication Number Publication Date
CN114111805A true CN114111805A (en) 2022-03-01

Family

ID=80441267

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111132200.9A Pending CN114111805A (en) 2021-09-23 2021-09-23 Multi-source multi-category measurement data position reference high-precision alignment method for carrier rocket

Country Status (1)

Country Link
CN (1) CN114111805A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116068549A (en) * 2023-03-06 2023-05-05 中国人民解放军32035部队 On-orbit spacecraft instantaneous attitude inversion method based on ground-based light thunder image fusion

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11271088A (en) * 1998-03-24 1999-10-05 Chubu Electric Power Co Inc Automatic surveying device
JP2010261842A (en) * 2009-05-08 2010-11-18 Ihi Corp Inertial navigation device, flying object, and navigation data calculation method
CN103913163A (en) * 2014-03-24 2014-07-09 中国人民解放军63680部队 Method for calculating A-E-C triaxial antenna coordinate of shipborne satellite communication earth station
CN104021311A (en) * 2014-06-23 2014-09-03 宫志华 Hermite function constraint-based data fusion calculation method
CN106052717A (en) * 2016-07-15 2016-10-26 中国西安卫星测控中心 Method for calibrating radio wave refraction correction effects by virtue of precision trajectory
CN109781116A (en) * 2018-11-16 2019-05-21 中国西安卫星测控中心 Error auto-correction fusion and positioning method based on active sensor mean iterative
CN110954092A (en) * 2019-11-28 2020-04-03 上海航天控制技术研究所 Collaborative navigation method based on relative measurement information assistance
CN111027204A (en) * 2019-12-05 2020-04-17 中国人民解放军63620部队 Method for fusion processing of space emission light, thunder and remote and navigation satellite measurement data
CN112255647A (en) * 2020-10-29 2021-01-22 中国人民解放军63921部队 Optimization method for trajectory measurement information post fusion processing
US20210080287A1 (en) * 2019-09-18 2021-03-18 Harbin Engineering University Method for initial alignment of radar assisted airborne strapdown inertial navigation system

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11271088A (en) * 1998-03-24 1999-10-05 Chubu Electric Power Co Inc Automatic surveying device
JP2010261842A (en) * 2009-05-08 2010-11-18 Ihi Corp Inertial navigation device, flying object, and navigation data calculation method
CN103913163A (en) * 2014-03-24 2014-07-09 中国人民解放军63680部队 Method for calculating A-E-C triaxial antenna coordinate of shipborne satellite communication earth station
CN104021311A (en) * 2014-06-23 2014-09-03 宫志华 Hermite function constraint-based data fusion calculation method
CN106052717A (en) * 2016-07-15 2016-10-26 中国西安卫星测控中心 Method for calibrating radio wave refraction correction effects by virtue of precision trajectory
CN109781116A (en) * 2018-11-16 2019-05-21 中国西安卫星测控中心 Error auto-correction fusion and positioning method based on active sensor mean iterative
US20210080287A1 (en) * 2019-09-18 2021-03-18 Harbin Engineering University Method for initial alignment of radar assisted airborne strapdown inertial navigation system
CN110954092A (en) * 2019-11-28 2020-04-03 上海航天控制技术研究所 Collaborative navigation method based on relative measurement information assistance
CN111027204A (en) * 2019-12-05 2020-04-17 中国人民解放军63620部队 Method for fusion processing of space emission light, thunder and remote and navigation satellite measurement data
CN112255647A (en) * 2020-10-29 2021-01-22 中国人民解放军63921部队 Optimization method for trajectory measurement information post fusion processing

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李冬;孙杰;: "弹道跟踪数据融合处理的快速算法", 飞行器测控学报, no. 04, pages 69 - 74 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116068549A (en) * 2023-03-06 2023-05-05 中国人民解放军32035部队 On-orbit spacecraft instantaneous attitude inversion method based on ground-based light thunder image fusion
CN116068549B (en) * 2023-03-06 2023-06-06 中国人民解放军32035部队 On-orbit spacecraft instantaneous attitude inversion method based on ground-based light thunder image fusion

Similar Documents

Publication Publication Date Title
CN110487301B (en) Initial alignment method of radar-assisted airborne strapdown inertial navigation system
CN103217159B (en) A kind of SINS/GPS/ polarized light integrated navigation system modeling and initial alignment on moving base method
US6876926B2 (en) Method and system for processing pulse signals within an inertial navigation system
CN104165640B (en) Near-space missile-borne strap-down inertial navigation system transfer alignment method based on star sensor
CN110926468B (en) Communication-in-motion antenna multi-platform navigation attitude determination method based on transfer alignment
CN101788296B (en) SINS/CNS deep integrated navigation system and realization method thereof
CN100356139C (en) Miniature assembled gesture measuring system for mini-satellite
CN105910602B (en) A kind of Combinated navigation method
CN111121766B (en) Astronomical and inertial integrated navigation method based on starlight vector
CN109633724B (en) Passive target positioning method based on single-satellite and multi-ground-station combined measurement
CN110672128B (en) Starlight/inertia combined navigation and error online calibration method
CN105184002A (en) Simulation analysis method for pointing angle of data transmission antenna
CN105115508A (en) Post data-based rotary guided projectile quick air alignment method
US3214575A (en) Celestial-inertial navigation system
CN105180728A (en) Front data based rapid air alignment method of rotary guided projectiles
CN102679999A (en) Star sensor installation error four-position calibrating and compensating method
CN114111805A (en) Multi-source multi-category measurement data position reference high-precision alignment method for carrier rocket
CN112284412B (en) Ground static alignment method for avoiding precision reduction caused by singular Euler transformation
CN116481564A (en) Polar region double-inertial navigation collaborative calibration method based on Psi angle error correction model
CN114167459B (en) Implementation method of high-precision intelligent remote sensing small satellite
CN113884704B (en) High-precision alignment method for speed references of multi-source multi-class measurement data of carrier rocket
CN115950419A (en) Combined navigation method, device and system for subminiature unmanned aerial vehicle
CN114111804A (en) Multi-source multi-category measurement data time zero high-precision alignment method for carrier rocket
CN111090830B (en) On-orbit light pressure identification method for high-orbit non-cooperative target
CN109506645B (en) Star sensor mounting matrix ground accurate measurement method

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