CN112782710B - Low-orbit space target tracking method and system based on interferometry system - Google Patents
Low-orbit space target tracking method and system based on interferometry system Download PDFInfo
- Publication number
- CN112782710B CN112782710B CN202011589416.3A CN202011589416A CN112782710B CN 112782710 B CN112782710 B CN 112782710B CN 202011589416 A CN202011589416 A CN 202011589416A CN 112782710 B CN112782710 B CN 112782710B
- Authority
- CN
- China
- Prior art keywords
- measurement
- target
- orbit
- station
- ambiguity
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
- G01S17/08—Systems determining position data of a target for measuring distance only
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/66—Tracking systems using electromagnetic waves other than radio waves
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
A low-orbit space target tracking method and system based on an interferometry system relate to the technical field of aerospace, and aim to establish a new active moving mode target tracking technical means through application improvement of interferometry technology and multi-base radar technology, and improve the accuracy and reliability of low-orbit non-cooperative space target tracking. The method comprises the steps of constructing a combined tracking strategy in a one-transmitting and multi-receiving mobile station distribution mode, wherein a main station transmits electromagnetic waves and measures target distance and angle, a base line formed by a secondary station and the main station carries out high-precision time difference measurement, and then the characteristics of system measurement information composition are utilized to establish a multi-method combined interference phase ambiguity resolution strategy, so that target tracking under different system constraints and tracking timeliness requirements is realized.
Description
Technical Field
The invention relates to the technical field of aerospace, in particular to a low-orbit space target tracking method and system based on an interferometry system.
Background
With the increasingly deep technical development and space application, the threat of space debris to the in-orbit spacecraft is continuously aggravated, the potential confrontation under a specific environment is gradually strengthened, the innovative application of a reliable high-precision measurement technology is urgently needed, and the tracking capability of non-cooperative space targets such as the space debris is improved. Specifically, on one hand, the efficiency and the precision of target tracking are improved, and on the other hand, the maintenance capability of the measuring system is also improved.
At present, two main systems of optics and radar are mainly used for tracking low-orbit space targets. The former mainly obtains high-precision angle measurement data, and obtains target position or track information through a plurality of devices or long-time observation, but the measurement system has larger dependence on weather conditions. The device can provide ranging, angle measuring and speed measuring information, can realize single device resolving, but the conventional radar angle measuring precision is limited, and the tracking precision is limited. In addition, existing spatial target measurement devices are typically bulky and may be less resistant to interference and sustained operation in future high-intensity applications or harsh environmental conditions.
The interferometry technique utilizes multiple stations to obtain the time difference (i.e., distance difference) between a space target and different stations through radio interferometry, which corresponds to the measurement of the projection of the spacecraft orbit along the baseline direction. The interferometry is actually equivalent to high-precision angle measurement, and is characterized in that target reflected (scattered) electromagnetic waves are obtained by using spatially separated equipment (antennas), and then a high-precision direction of arrival is calculated through relative position and phase interference processing between the equipment, so that the space target tracking measurement precision is obviously improved. However, most of the current interferometric systems are applied to the orbit measurement of high orbit space targets and deep space spacecraft, and rely on long-term accumulation. In addition, most of the existing interferometry is fixed layout, and needs external guide information with certain precision, which is not beneficial to practical application in a high-intensity mode.
Therefore, the invention gives full play to the precision advantage of interferometry, adopts a single-transmitting multi-receiving maneuvering multi-station mode, provides a tracking measurement method and a system based on an interferometry system, and provides a new technical means for non-cooperative low-orbit space target tracking.
Disclosure of Invention
The invention aims to solve the technical problem that a low-orbit non-cooperative space target tracking method and system of a one-shot and multi-shot interferometry system are provided based on an interferometry technology and a multi-base radar technology and face the tracking requirement of a near-earth space non-cooperative spacecraft, and the information acquisition requirement of situation perception tasks such as space debris monitoring is met.
In order to achieve the purpose, the invention fully utilizes the characteristics of small equipment scale and high angle measurement precision of an interference measurement system, combines the advantages of active detection and multi-station interference to design a measurement method, and achieves the purposes of reducing the system processing difficulty, improving the target tracking precision and enhancing the flexibility and reliability of the system by transmitting, receiving, separately arranging and flexibly arranging stations to realize system configuration.
The invention provides a low-orbit space target tracking method based on an interferometry system, which comprises the following steps:
step1: adopting a one-transmission and multi-reception measurement mode to construct a measurement model:
transmitting electromagnetic waves through a master station antenna and receiving target echoes to obtain target distance and angle measurement information; the plurality of secondary station antennas only receive the target echoes, and are combined with the main station to form a plurality of measurement baselines to obtain the time difference of the echoes reaching each survey station; combining a plurality of time difference measurement equations with a master station ranging equation to form a measurement model for target fine tracking;
step2: interference phase ambiguity resolution strategy:
ambiguity initial value calculation based on primary positioning information of a main station is carried out, and ambiguity estimation accuracy and reliability are improved by changing two modes of electromagnetic wave frequency and mobile station base line configuration and utilizing the characteristic that interference phase ambiguities measured by different frequencies and different base lines are relatively prime;
and 3, step 3: target tracking and positioning:
a mode of positioning first and then orbit determination is adopted, a short baseline connection line interference measurement mode is introduced, and the target real-time position and speed estimation is carried out by combining the high-precision positioning of the master station and the high-precision time difference measurement information of the baseline; combining the space target tracking data for multiple circles to realize high-precision orbit estimation; and the selectable mode of the differential observation system is constructed by adjusting the layout of the observation stations, so that the positioning and orbit determination precision is further improved.
Further, the step1 comprises the following steps:
and (3) radar measurement modeling of the master station:
adopting a conventional three-coordinate radar measurement mode, transmitting signals and receiving signalsEcho, extracting distance and angle information from a target to a master station radar antenna phase center from the echo; the modeling is as follows, the three-dimensional coordinates of the target space position at the moment of t are recorded as
Master station radar antenna phase center position
The master station measurement model is a typical mobile three-coordinate radar measurement model
Where ρ is 0 (t)、A 0 (t)、E 0 (t) are distance, azimuth and pitch measurements, respectively, all varying with time.
Further, the step1 comprises the following steps:
multi-station interferometry modeling:
all stations receive target echo at the same time, and the time difference of the signals arriving at the secondary station and the primary station is tau 10 、τ 20 、τ 30 Respectively calculating the time difference of the same signal arriving at the secondary station and the primary station as follows:
where c is the speed of light, ρ 0 (t)、ρ 1 (t)、ρ 2 (t)、ρ 3 (t) are respectively target to Master station X 0 (t) and secondary station X 1 (t)、X 1 (t)、X 3 Distance of (t) (. Rho) 1 (t)、ρ 2 (t)、ρ 3 The measurement equation of (t) is similar to the first equation of (0); recording the signal wavelength as lambda, and the phase observed quantity and the integer ambiguity obtained by the interference measurement of the main station 0 and the auxiliary station i are respectively phi i0 (t) and N (t) are
φ 10 (t)+λ·N i0 (t)＝ρ i (t)-ρ 0 (t) (0)
In summary, the baseline interferometry model is
Further, the step2 comprises the following steps:
and (3) resolving ambiguity by using a main station positioning result:
firstly: calculating and acquiring a target space position measured value X (t) through coordinate transformation;
secondly, the method comprises the following steps: according to the known site coordinates X of two stations forming a base line A (t) and X B (t) calculating a distance difference of the target at the baseline measurement;
and thirdly: according to the spectral characteristics of the measured signal, the integer constant speed is recorded, and the ambiguity in the distance difference is
And finally: analyzing the calculation precision of the integer ambiguity according to the precision of X (t), and if the error of the denominator at the right end exceeds 0.5 lambda, judging that N is equal to AB And (5) as an initial value of the ambiguity, switching to the next link, and otherwise, ending the ambiguity calculation process.
Further, the step2 comprises the following steps:
carrying out variable frequency deblurring:
firstly, selecting a transmitting frequency point f 1 And f 2 Calculating the integer ambiguity through a formula (0); because the X (t) obtained by the main station has measurement errors, the calculated ambiguity has a plurality of values
And
second, since the target location and baseline are determined, so
And
a corresponding relation exists, a combination which is consistent with the two is selected from the set of the previous step, and the corresponding is possible ambiguity estimation;
and finally, if a plurality of groups of corresponding elements exist, increasing the transmitting frequency points and guiding the public elements of the plurality of groups of sets to be unique, thereby obtaining the final ambiguity estimation value.
Further, the step2 comprises the following steps:
multi-baseline deblurring:
firstly, arranging each measuring station according to a target prior orbit, so that the lengths of all base lines are mutually prime, and increasing the non-fuzzy range;
and secondly, calculating each base line ambiguity by using X (t) obtained by the main station, and obtaining an ambiguity calculation value through a set corresponding relation.
Further, the step 3 comprises the following steps:
and (3) interference measurement:
firstly, for the measuring stations forming the connecting line measuring base line, the measuring precision of the corresponding time difference is improved by using the same time-frequency standard and eliminating modes such as common errors;
secondly, because the relative position relation of each base line is known, the base line resolving result is used as a reference, and auxiliary information is provided for ambiguity calculation of other base lines of the system;
and (3) real-time positioning:
firstly, combining a space target orbit dynamics model and the measurement model provided in the first step, and establishing a motion equation and a measurement equation of a target;
secondly, establishing a composite filter for real-time filtering of the target position and speed according to the system characteristics;
high-precision rail fixing:
after each tracking circle of the target is finished, combining the measurement data of more than three circles by adopting the measurement model provided in the step1 and combining a spacecraft dynamics constraint model to perform precise orbit determination on the target;
after each orbit determination calculation, evaluating the contribution of the current base line to the orbit calculation in time, forecasting the orbit of the next circle, adjusting the base line configuration according to the evaluation result of the base line and the forecasting result of the orbit, and serving the subsequent measurement and improvement of the orbit;
precision optimization:
under the condition that observation conditions allow, a differential interference measurement mode is introduced, so that the measurement precision is further improved, and the precision of target positioning and orbit determination is improved.
Furthermore, the composite filter comprises a main filter and an auxiliary filter, and the main filter forms a dimension-expanding combined filtering model by utilizing all information of the active measurement information of the main station and the passive measurement information of the multiple stations; and the secondary filter only utilizes the three-coordinate information of the master radar to calculate the target position information in real time.
The invention discloses a low-orbit space target tracking system based on an interferometry system, which adopts a low-orbit non-cooperative space target tracking method and comprises a system composition, a system deployment mode and a system working flow, wherein the system composition comprises a transmitting-receiving antenna system, a signal acquisition and data preprocessing system, a data transmission system and a data comprehensive processing system; the system deployment mode comprises a main station and a plurality of secondary stations, and adopts transceiving split placement and mobile vehicle-mounted; the system has the following working procedures: and adopting the strategies of multi-station cooperation, decentralized acquisition and centralized processing.
The invention can realize the following beneficial technical effects:
for tracking of low-orbit non-cooperative space targets, the existing method and system are difficult to give consideration to both precision and flexibility, and the application under the high-strength confrontation condition is limited. The invention provides a target tracking method and a target tracking system in an active moving mode based on an interferometric technical system, and the target tracking method and the target tracking system have the advantages of high tracking precision, low processing complexity, flexible deployment, low construction cost and wide application, and particularly comprise the following steps:
(1) According to the method and the system, the target high-precision distance measurement and the time difference (angle measurement) information can be simultaneously obtained, the target observation geometry can be flexibly improved by adjusting the layout of the measuring station to optimize the baseline configuration, and a high-precision high-reliability target tracking result is obtained at a lower processing cost;
(2) According to the method and the system, the flexibility of the system can be fully utilized, the algorithm complexity of phase ambiguity resolution in interferometry is reduced by optimally designing the type of the transmitted signal and the arrangement mode of the mobile baseline, and the time difference extraction precision is improved at lower processing cost;
(3) According to the method and the system, a receiving and transmitting separately-arranged mode is adopted, each survey station adopts a mobile platform, the system can freely select a deployment site and a baseline configuration, dynamic adjustment can be performed in a measurement gap or even in a measurement process, measurement accuracy and reliability requirements of different targets can be considered, and the countermeasure and the survival capability of the system can be improved.
(4) According to the method and the system, the designed interference measurement system has the advantages of small system scale, single-station measurement function and strong universality of software and hardware, and can greatly reduce the system construction cost while ensuring the target tracking precision;
(5) The method and the system can be expanded to a non-cooperative space target measuring mode, and the measurement of the non-cooperative space target is realized through receiving space signals radiated by the target, or self telemetering and load working signals and the like and through interference baselines obtained by a plurality of receiving stations.
(6) The method and the system can be directly applied to low-orbit space target tracking and can be expanded to a medium-high orbit target, and meanwhile, the methods such as interference phase solution blurring and baseline configuration design have certain reference significance for the development of the theory and technology of interferometry.
Drawings
FIG. 1 is a schematic diagram of one-shot multiple-shot spatial target interferometry based on the present invention;
FIG. 2 is a schematic diagram of differential interferometry;
FIG. 3 is a schematic diagram of the principal components and deployment of an interferometric system according to an embodiment of the invention;
fig. 4 is a schematic system flow diagram of the present invention.
Detailed Description
The invention provides a low-orbit non-cooperative space target tracking method and system based on a one-shot multi-shot interferometry system. In the following, embodiments and features of the embodiments in the present application may be combined with each other without conflict, and specific embodiments of the present invention will be described in detail.
1. Low-orbit space target tracking method based on one-shot multi-shot interferometry
The invention provides a one-transmission multi-reception space target interferometry method, taking one-transmission and three-reception as an example, and a measurement schematic diagram is shown in figure 1. Through 1 transmission master station and 3 secondary stations, combine single-station radar location and many baselines time difference measurement, realize non-cooperation space target position and the orbit calculation of high accuracy. The specific process comprises three main steps of measurement model construction, interference phase ambiguity resolution and target tracking and positioning. The receiving secondary stations can be increased or decreased according to the constraint of tracking requirement, construction cost and the like.
The method comprises the following steps: measurement model construction
Transmitting electromagnetic waves through a master station antenna and receiving target echoes to obtain target distance and angle measurement information; the plurality of secondary station antennas only receive target echoes, and are combined with the main station to form a plurality of measurement baselines to obtain the time difference of the echoes reaching each survey station; combining a plurality of time difference measurement equations with a master station ranging equation to form a measurement model for target fine tracking; the method comprises the steps of establishing a mathematical model of the relation between the space position of a space target and measurement information according to the observation geometry and the type of the measurement information in a mode that a main station transmits and a main station and an auxiliary station receive simultaneously. According to the measurement schematic diagram of fig. 1, the modeling includes two links of a master station radar measurement and a multi-station interferometry.
Step1: master station radar measurement modeling
The master station adopts a conventional three-coordinate radar measurement mode, transmits signals and receives echoes, and extracts the distance from a target to the phase center of a radar antenna of the master station from the echoesAnd divergence angle information. The steps are modeled as follows, and the three-dimensional coordinate of the target space position at the moment t is recorded as X (t) = [ X (t), y (t), z (t)] T The central position X of the radar antenna phase of the master station 0 (t)＝[x 0 (t),y 0 (t),z 0 (t)] T The master station measurement model is a typical mobile three-coordinate radar measurement model
Where ρ is 0 (t)、A 0 (t)、E 0 (t) are distance, azimuth and pitch measurements, respectively, all varying with time.
Step2: multi-station interferometric modeling
The multi-station interferometry is exemplified by a measurement system of one primary station and 3 secondary stations, and the measurement model of the procedure is given, and other cases can be similarly derived. Specifically, each station receives target echoes simultaneously, and the time difference between the arrival of the signal at the secondary station and the arrival of the signal at the primary station are respectively tau 10 、τ 20 、τ 30 The time difference of arrival of the same signal at the secondary station and the primary station is calculated as follows
Where c is the speed of light, ρ 0 (t)、ρ 1 (t)、ρ 2 (t)、ρ 3 (t) are respectively target to Master station X 0 (t) and secondary station X 1 (t)、X 1 (t)、X 3 Distance of (t) (. Rho) 1 (t)、ρ 2 (t)、ρ 3 The measurement equation of (t) is similar to the first equation of (0). Here, the primary station is used as the starting point of time difference calculation, the starting reference is changed, and similarly, tau can be calculated 01 、τ 21 And τ 31 And the number of the 4 stations is only 3 independent baselines, and the number of the corresponding independent time difference measurement information is only 3.
Recording the signal wavelength as lambda, and the phase observed quantity and the integer ambiguity obtained by the interference measurement of the main station 0 and the auxiliary station i are respectively phi i0 (t) and N (t), thenIs provided with
φ 10 (t)+λ·N i0 (t)＝ρ i (t)-ρ 0 (t) (0)
In summary, the baseline interferometry model is
Station position X in equations (0) and (0) 0 (t)、X 1 (t)、X 2 (t)、X 3 (t) can be obtained by its own positioning system, considered as a known quantity, and in order to simplify the processing difficulty and improve the phase measurement accuracy, the position of the measuring station can be fixed within the measurement time of one tracking arc of the target.
It should be noted that in the modeling equation of the step, no matter the radar measurement of the main station or the time difference measurement of the multiple stations, measurement system errors and random errors are inevitable, the specific composition analysis and processing strategy needs to be specifically analyzed in the application process, and the relevant details are not detailed in the invention.
Step two: interferometric phase ambiguity resolution strategy
The method comprises the steps of sequentially realizing ambiguity initial value resolving based on primary positioning information of a main station by adopting a multi-method combined resolving strategy, and improving ambiguity estimation precision and reliability by changing two modes of electromagnetic wave frequency and mobile station base line configuration and utilizing the characteristic of relatively prime interference phase ambiguity measured by different frequencies and different base lines; the step is to solve the integer ambiguity of the phase difference by processing the interference measurement signals of the baselines formed by the multiple stations, and provide a basis for calculating the target position or the orbit according to the time difference. As shown in fig. 2, the invention improves the accuracy and reliability of the phase ambiguity step by step through the calculation strategy of three links.
Step1: disambiguation of position results using a Master station
The step is to provide target three-coordinate measurement information by using a master station radar and acquire an initial value of the ambiguity of the interferometric measurement.
First, a target spatial position measurement value X (t) is obtained by calculation through coordinate transformation.
Secondly, according to the known site coordinates X of the two stations forming the base line A (t) and X B (t) calculating the distance difference of the target under the baseline measurement.
Thirdly, according to the frequency spectrum characteristics of the measured signal, the integer constant speed is recorded, and the ambiguity in the distance difference is
Finally, analyzing the calculation precision of the integer ambiguity according to the precision of X (t), and if the error of the denominator at the right end exceeds 0.5 lambda, judging that N is equal to N AB And (5) as an initial value of the ambiguity, switching to the next link, and otherwise, ending the ambiguity calculation process.
In order to improve the calculation accuracy and reliability of the ambiguity and reduce the uncertainty influence of the single-point coordinates of the main station, N of continuous multiple moments can be taken in one measurement AB The average or median of (t) is taken as the ambiguity calculation result.
Step2: variable frequency deblurring
The steps give full play to the advantages of active measurement, and the ambiguity resolution precision is improved by designing the frequency of a transmitting signal and combining different frequency points.
Firstly, a transmitting frequency point f is selected 1 And f 2 The integer ambiguity is calculated by equation (0). Because the X (t) obtained by the main station has measurement error, the calculated ambiguity has a plurality of values
And
second, since the target location and baseline are determined, so
And
there is a corresponding relation, and the combination of the two is selected from the set of the above steps, and the corresponding is the possible ambiguity estimation.
And finally, if a plurality of groups of corresponding elements exist, increasing the transmitting frequency points and guiding the public elements of the plurality of groups of sets to be unique, thereby obtaining the final ambiguity estimation value.
Step 3: multi-baseline deblurring
The step is to design the direction and the length of a space baseline by using the layout flexibility of the mobile measuring station, change the observation geometry and improve the ambiguity resolution precision.
Firstly, arranging each measuring station according to a target prior orbit, so that the lengths of all base lines are mutually prime numbers, and increasing the non-fuzzy range.
And secondly, calculating each baseline ambiguity by using X (t) obtained by the main station, and acquiring ambiguity calculation values by collecting corresponding relations in a Step2 similar manner. If a baseline combination is difficult to achieve accurate ambiguity fixing, the baseline configuration can be changed appropriately during the measurement.
Step three: target tracking and positioning
The method of positioning first and then fixing rails is adopted, so that the requirements of target tracking tasks with different timeliness are met; on the basis of improving the time difference measurement accuracy by introducing a short baseline connection line interference measurement mode, firstly, combining the high-accuracy positioning of a main station and the high-accuracy time difference measurement information of a baseline, and carrying out real-time position and speed estimation on a target through a real-time filtering algorithm, and secondly, combining the multi-turn tracking data of a space target, and realizing high-accuracy track estimation; in addition, the selectable mode of the differential observation system is constructed by adjusting the layout of the measuring station, so that the positioning and orbit determination precision is further improved.
The step uses the measurement information to complete the real-time position and speed and high-precision orbit calculation of the space target. As shown in fig. 3, the present invention adopts a method of positioning first and then tracking to meet different timeliness task requirements, and mainly comprises 4 main links, as follows.
Step1: interferometric measurement
The invention obtains the multi-station time difference based on the interference measurement as the main measurement information of target tracking. The step is that in the measuring process, individual measuring stations are selected to form a base line below 1km, a short base line connection interference measuring mode is established, and the measuring precision is improved.
Firstly, for the stations forming the connecting line measurement base line, the measurement precision of the corresponding time difference is improved by using the same time-frequency standard and eliminating modes such as common errors.
Secondly, because the relative position relation of each base line is known, the base line resolving result is used as a reference, and auxiliary information is provided for ambiguity calculation of other base lines of the system.
Step2: real-time positioning
The method mainly aims at the target tracking requirement of a single circle, and is characterized in that the active measurement information of a main station and the passive measurement information of multiple stations are comprehensively utilized to calculate the position and the speed of a target in real time and provide real-time positioning information of the target.
Firstly, combining a space target orbit dynamics model and the measurement model provided in the first step, establishing a motion equation and a measurement equation of the target. The orbit dynamics model can refer to relevant books, and the description is omitted here.
And secondly, establishing a composite filter for real-time filtering of the target position and speed according to the system characteristics. The composite filter comprises a main filter and an auxiliary filter, wherein the main filter forms an extended-dimension combined filter model by utilizing all information of active measurement information of a main station and passive measurement information of a plurality of stations, namely ranging, angle measurement and time difference information, and the calculation of the position and the speed of a target is carried out, and the result is output as positioning information. The secondary filter only utilizes the three-coordinate information of the radar of the main station to calculate the target position information in real time, on one hand, an initial value is provided in the initialization stage of the main filter, on the other hand, the secondary filter is used as the check information of the main filter in the normal operation period, and abnormal phenomena such as filter divergence and the like caused by the abnormality of an interferometric system are avoided.
Step 3: high-precision orbit determination
The method mainly aims at the requirement of target high-precision tracking, and combines the measurement information of a plurality of circles to realize target track estimation. This step involves two main links.
Firstly, after each tracking circle of the target is finished, the measurement model provided in the first step is adopted, measurement data above three circles are combined, and the spacecraft dynamics constraint model is combined to perform precise orbit determination of the target. The calculation adopts a statistical orbit determination method, and the initial value of the orbit is obtained from Step 2.
And secondly, after each orbit determination calculation, the contribution of the current base line to the orbit calculation is evaluated in time, the orbit of the next circle is forecasted, and the base line configuration is adjusted according to the base line evaluation and orbit forecasting results, so that the subsequent orbit measurement and improvement are provided.
Step 4: precision optimization
The step is to introduce a differential interference measurement mode under the condition that observation conditions allow, so that the measurement precision is further improved, and the precision of target positioning and orbit determination is improved.
According to the observation capability of the observation station, if reference sources (such as a radio frequency source, a navigation satellite and the like) which are viewed by the observation station in common exist in the observation sight direction of the target in a certain period of time, as shown in fig. 4, the observation of the reference sources and the target by the observation station is combined to form a differential observation system, so that the common system error influence of clock error between stations and the like is eliminated.
2. Low-orbit space target tracking system based on one-shot multi-shot interferometry
The embodiment of the invention also provides a design scheme of the low-orbit non-cooperative space target tracking system based on the one-shot multiple-shot interferometry system, and the detailed design steps of the system composition and the working flow are as follows.
The method comprises the following steps: design of system components
As an implementation mode, the system is composed of a transmitting-receiving antenna system, a signal acquisition and data preprocessing system, a data transmission system, a data comprehensive processing system and the like.
The receiving and transmitting antenna system obtains original signals for distance measurement, angle measurement and time difference measurement in a multi-station combined mode;
the signal acquisition and processing system completes the acquisition and pretreatment of the original data required by single-station positioning and interferometry;
the data transmission system selects a flexibly adopted information transmission strategy according to different distance base lengths and system constraints, and is responsible for collecting the acquired signals to the data comprehensive processing system;
the data comprehensive processing system comprises three parts of signal interference processing, real-time positioning and high-precision orbit determination resolving, and is used for respectively finishing the acquisition of time difference measurement data, the real-time filtering resolving of target position and speed and the high-precision orbit determination of a target by combining multiple circles of tracking data.
Step two: system deployment style design
The system comprises a main station and a plurality of secondary stations, and is deployed in the modes of transceiving split placement, mobile vehicle-mounted, flexible communication, decentralized processing and the like, and the deployment principles are as follows:
a receiving-transmitting integrated antenna is additionally installed at the primary station, a receiving antenna is additionally installed at the secondary station, and a plurality of receiving antennas can be additionally installed at one survey station to form a short baseline for connection measurement (as shown in fig. 1 for the primary station and the secondary station 1);
a signal acquisition and preprocessing system is distributed in each station, a main station acquisition system finishes distance measurement, angle measurement signal acquisition and preprocessing and acquisition and preprocessing of signals before interference, and a secondary station acquisition system finishes acquisition and preprocessing of signals before interference;
the data transmission system is responsible for collecting the collected signals to the data comprehensive processing system, and communication means such as direct-connected optical fibers, special networks or encrypted 5G are flexibly adopted according to different distance baselines as required, as shown in FIG. 1;
the data integrated processing system completes signal interference processing and target positioning solution, and can be deployed in a separate processing center or arranged in a master station.
Step three: system workflow design
The system adopts the strategy of multi-station cooperation, decentralized acquisition and centralized processing, and the working process can be briefly expressed as follows:
flexibly setting the layout of each station according to the target prior information, and guiding a master station radar to autonomously track a target;
after the master station captures the target and obtains the positioning information, guiding the secondary station to track;
collecting the measurement information of the primary station and the secondary station, and extracting the high-precision multi-station time difference;
the target high-precision real-time positioning is realized by using the distance measurement information of the master station and the time difference information of multiple stations;
evaluation of the contribution of the base line to the track, improvement of the station layout;
and combining the observation data of multiple circles to perform target precise orbit determination.
In summary, the invention provides a low-orbit space target tracking method and system based on an interferometric system, and the method comprises the following steps: constructing a measurement model, an interference phase ambiguity resolution strategy and target tracking and positioning by adopting a one-transmission and multi-reception measurement mode; the system comprises a system composition, a system deployment mode and a system working process, wherein the system composition comprises a receiving and transmitting antenna system, a signal acquisition and data preprocessing system, a data transmission system and a data comprehensive processing system; the system deployment mode comprises a main station and a plurality of secondary stations, and adopts transceiving split placement and mobile vehicle-mounted; the system has the following working procedures: and adopting the strategies of multi-station cooperation, decentralized acquisition and centralized processing. According to the invention, through application improvement of an interferometric measurement technology and a multi-base radar technology, a new active moving mode target tracking technical means is established, the precision and reliability of low-orbit non-cooperative space target tracking are improved, a multi-method combined interferometric phase ambiguity resolution strategy is established by using the characteristics of system measurement information composition, and target tracking under different system constraints and tracking timeliness requirements is realized.
In the previous description, numerous specific details were set forth in order to provide a thorough understanding of the present invention. The foregoing description is only a preferred embodiment of the invention, which can be embodied in many different forms than described herein, and therefore the invention is not limited to the specific embodiments disclosed above. And that those skilled in the art may, using the methods and techniques disclosed above, make numerous possible variations and modifications to the disclosed embodiments, or modify equivalents thereof, without departing from the scope of the claimed embodiments. Any simple modification, equivalent change and modification of the above embodiments according to the technical essence of the present invention are within the scope of the technical solution of the present invention.
Claims (7)
1. A low-orbit space target tracking method based on an interferometry system is characterized in that: the method comprises the following steps:
step1: adopting a one-transmission and multi-reception measurement mode to construct a measurement model:
transmitting electromagnetic waves through a master station antenna and receiving target echoes to obtain target distance and angle measurement information; a plurality of secondary station antennas only receive target echoes, and are combined with the main station to form a plurality of measuring baselines, and the time difference of the echoes reaching each measuring station is obtained; combining a plurality of time difference measurement equations with a master station ranging equation to form a measurement model for target precise tracking;
step2: interference phase ambiguity resolution strategy:
ambiguity initial value calculation based on primary positioning information of a main station is carried out, and ambiguity estimation accuracy and reliability are improved by changing two modes of electromagnetic wave frequency and baseline configuration of mobile station arrangement and utilizing the characteristic of interference phase ambiguity reciprocity of different frequency and different baseline measurement;
and step 3: target tracking and positioning:
a mode of positioning first and then orbit determination is adopted, a short baseline connection line interference measurement mode is introduced, and target real-time position and speed estimation is carried out by combining the high-precision positioning of the main station and the high-precision time difference measurement information of the baseline; tracking the measurement data for multiple circles in combination with the space target to realize high-precision track estimation; and the selectable mode of the differential observation system is constructed by adjusting the layout of the observation stations, so that the positioning and orbit determination precision is further improved.
2. The low-orbit spatial target tracking method according to claim 1, characterized in that: the step2 further comprises the following steps:
and (3) resolving ambiguity by using a main station positioning result:
firstly: calculating and acquiring a target space position measurement value X (t) through coordinate transformation;
secondly, the method comprises the following steps: according toKnown site coordinates X of two stations forming a base line A (t) and X B (t) calculating a distance difference of the target at the baseline measurement;
and thirdly: recording based on spectral characteristics of the measured signal
To get the integer uniform velocity, the ambiguity in the distance difference is
And finally: analyzing the calculation precision of the integer ambiguity according to the precision of X (t), and if the error of the denominator at the right end exceeds 0.5 lambda, judging that N is equal to AB And (5) as an initial value of the ambiguity, switching to the next link, and otherwise, ending the ambiguity calculation process.
3. The low-orbit spatial target tracking method according to claim 2, characterized in that: the step2 further comprises the following steps:
and (3) variable frequency deblurring:
firstly, selecting a transmitting frequency point f 1 And f 2 Calculating the integer ambiguity through a formula (15); because the X (t) obtained by the main station has measurement error, the calculated ambiguity has a plurality of values
And
；
second, since the target location and baseline are determined, so
And
there is a correspondence, from the set of above stepsSelecting a combination with the consistency of the two, wherein the corresponding combination is an alternative ambiguity estimation value;
and finally, if a plurality of groups of corresponding elements exist, increasing the transmitting frequency points until the common elements of the plurality of groups of sets are unique, thereby obtaining the final ambiguity estimation value.
4. The low-orbit spatial target tracking method according to claim 3, characterized in that: the step2 further comprises the following steps:
multi-baseline deblurring:
firstly, arranging each measuring station according to a target prior orbit, so that the lengths of all base lines are mutually prime, and increasing the non-fuzzy range;
and secondly, calculating each base line ambiguity by using X (t) obtained by the main station, and obtaining an ambiguity calculation value through a set corresponding relation.
5. The low-orbit spatial target tracking method according to claim 1, characterized in that: the step 3 further comprises the following steps:
and (3) interference measurement:
firstly, for the stations forming the connecting line measurement base line, the measurement precision of the corresponding time difference is improved by using the same time-frequency standard and eliminating a common error mode;
secondly, because the relative position relation of each base line is known, the base line resolving result is used as a reference, and auxiliary information is provided for ambiguity calculation of other base lines of the system;
and (3) real-time positioning:
firstly, combining a space target orbit dynamics model and the measurement model provided in the first step, and establishing a motion equation and a measurement equation of a target;
secondly, establishing a composite filter for real-time filtering of the target position and speed according to the system characteristics;
high-precision rail fixing:
after the measurement data of each circle of the target is tracked, combining the measurement data of more than three circles by adopting the measurement model provided in the step1 and combining a spacecraft dynamics constraint model to perform precise orbit determination of the target;
after each orbit determination calculation, evaluating the contribution of the current base line to the orbit calculation in time, forecasting the orbit of the next circle, adjusting the configuration of the base line according to the evaluation result of the base line and the forecast result of the orbit, and serving the subsequent orbit measurement;
precision optimization:
under the condition that observation conditions allow, a differential interference measurement mode is introduced, so that the measurement accuracy is further improved, and the accuracy of target positioning and orbit determination is improved.
6. The low-orbit spatial target tracking method according to claim 5, characterized in that: the composite filter comprises a main filter and an auxiliary filter, and the main filter forms a combined filtering model for dimension expansion by utilizing all information of the active measurement information of the main station and the passive measurement information of the multiple stations; and the secondary filter only utilizes the three-coordinate information of the master station radar to calculate the target position information in real time.
7. A low-orbit space target tracking system based on an interferometry system, which adopts the low-orbit space target tracking method of any one of claims 1 to 6, and is characterized in that: the system comprises a system composition, a system deployment mode and a system working process, wherein the system composition comprises a receiving and transmitting antenna system consisting of a main station antenna and a secondary station antenna, a signal acquisition and data preprocessing system, a data transmission system and a data comprehensive processing system; the system deployment mode comprises a main station and a plurality of secondary stations, and adopts transceiving split placement and mobile vehicle-mounted; the system has the following working procedures: and adopting the strategies of multi-station cooperation, decentralized acquisition and centralized processing.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011589416.3A CN112782710B (en) | 2020-12-29 | 2020-12-29 | Low-orbit space target tracking method and system based on interferometry system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011589416.3A CN112782710B (en) | 2020-12-29 | 2020-12-29 | Low-orbit space target tracking method and system based on interferometry system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112782710A CN112782710A (en) | 2021-05-11 |
| CN112782710B true CN112782710B (en) | 2022-10-14 |
Family
ID=75753142
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011589416.3A Active CN112782710B (en) | 2020-12-29 | 2020-12-29 | Low-orbit space target tracking method and system based on interferometry system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112782710B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106507955B (en) * | 2007-05-10 | 2011-06-29 | 中国科学院国家天文台 | Two-stage distributed collaboration source localization method based on Pseudo Code Spread Spectrum |
| CN103968844A (en) * | 2013-02-04 | 2014-08-06 | 上海新跃仪表厂 | Large ellipse maneuverable spacecraft autonomous navigation method based on low-orbit platform tracking measurement |
| CN104309817A (en) * | 2014-10-11 | 2015-01-28 | 中国科学院国家授时中心 | Multiple parallel address receiver-based Beidou navigation satellite region orbit determination method |
| CN108061908A (en) * | 2017-12-11 | 2018-05-22 | 山东科技大学 | SLR station three-dimensional coordinate geometry calculation methods based on the spaceborne GNSS technologies of low orbit satellite |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11073608B2 (en) * | 2018-10-11 | 2021-07-27 | Raytheon Company | Resolving radar angle ambiguities using a multiple hypothesis tracker |
-
2020
- 2020-12-29 CN CN202011589416.3A patent/CN112782710B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106507955B (en) * | 2007-05-10 | 2011-06-29 | 中国科学院国家天文台 | Two-stage distributed collaboration source localization method based on Pseudo Code Spread Spectrum |
| CN103968844A (en) * | 2013-02-04 | 2014-08-06 | 上海新跃仪表厂 | Large ellipse maneuverable spacecraft autonomous navigation method based on low-orbit platform tracking measurement |
| CN104309817A (en) * | 2014-10-11 | 2015-01-28 | 中国科学院国家授时中心 | Multiple parallel address receiver-based Beidou navigation satellite region orbit determination method |
| CN108061908A (en) * | 2017-12-11 | 2018-05-22 | 山东科技大学 | SLR station three-dimensional coordinate geometry calculation methods based on the spaceborne GNSS technologies of low orbit satellite |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112782710A (en) | 2021-05-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104316903B (en) | A kind of three station positioning using TDOA performance test appraisal procedures | |
| JP2020509386A (en) | Method and apparatus for capturing surroundings | |
| CN103746757B (en) | A kind of single star interference source localization method based on satellite multi-beam antenna | |
| Aernouts et al. | TDAoA: A combination of TDoA and AoA localization with LoRaWAN | |
| CN109782289B (en) | Underwater vehicle positioning method based on baseline geometric structure constraint | |
| CN105549005A (en) | Dynamic target direction of arrive tracking method based on mesh dividing | |
| CN103698743A (en) | Ionospheric-reflection-based time difference of arrival positioning method for shortwave radiation source | |
| KR101953185B1 (en) | Multistatic passive coherent location method for estimating target location using fm network | |
| CN102879784A (en) | Unified imaging method for synthetic aperture radar (SAR) in four modes | |
| CN110954865A (en) | Short wave time difference positioning method based on ionosphere information | |
| US20190204433A1 (en) | Method of tracking target by using 2d radar with sensor | |
| CN101561499A (en) | Single-station Doppler distance-measuring and positioning method | |
| Aernouts et al. | Combining TDoA and AoA with a particle filter in an outdoor LoRaWAN network | |
| CN110554373A (en) | Interferometric time measuring and ranging method | |
| CN104267420A (en) | Satellite-borne three-dimensional moving object positioning method, device and system | |
| CN115113155A (en) | Airborne distributed aperture coherent synthetic radar testing and evaluating method | |
| CN112782710B (en) | Low-orbit space target tracking method and system based on interferometry system | |
| CN103616669B (en) | In two-dimensional frequency, channel-equalization method is carried out to non-homogeneous scene | |
| CN113933876B (en) | Multi-star communication time difference positioning data fusion processing method | |
| KR101887877B1 (en) | Multistatic passive coherent location system and method for estimating target location using fm network | |
| Burns et al. | IFSAR for the rapid terrain visualization demonstration | |
| CN113203985B (en) | Direct positioning method for shortwave same-frequency signals | |
| Mikhalev et al. | Passive emitter geolocation using agent-based data fusion of AOA, TDOA and FDOA measurements | |
| Li et al. | Performance Analysis of Collaborative Positioning Method in Automated Car Driving. | |
| Kelner et al. | The empirical verification of the location method based on the Doppler effect |
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 |