Interference measurement method, interference measurement device and interference measurement system

Abstract

The present application provides an interferometric measurement method, an interferometric measurement device, and an interferometric measurement system, the method comprising: receiving a plurality of first estimated data and a plurality of second estimated data, the first estimated data being obtained by the first station performing frequency estimation and phase estimation on the first test data, the second estimated data being obtained by the second station performing frequency estimation and phase estimation on the second test data, the first test data being the observation data of the first station on the target detector, and the second test data being the observation data of the second station on the target detector; obtaining an initial theoretical delay, the initial theoretical delay being the theoretical delay at the first acquisition moment, the theoretical delay being the theoretical transmission delay of the first station relative to the second station; determining the interferometric measurement delay observation value of the first station and the second station according to each first estimated data, each second estimated data, and the initial theoretical delay. The present application solves the problem of high communication link bandwidth requirements for data transmission in interferometric measurement.

Description

Technical Field

[0001] The present application relates to the field of interference measurement technology, and in particular, to an interference measurement method, an interference measurement device, a computer-readable storage medium, a processor, and an interference measurement system. Background Art

[0002] Interferometry technology provides good lateral constraints for deep space probes with high-precision angle measurement information, which can greatly improve the accuracy of deep space probe orbit determination. Since the 1970s, the Jet Propulsion Laboratory (JPL) of the United States has applied interferometry technology to the Voyager deep space exploration mission and subsequent Mars exploration missions, significantly improving deep space exploration capabilities. The European Space Agency, the Japan Aerospace Exploration Agency (JAXA), and others are promoting the research and application of interferometry technology in deep space exploration missions. In China, with the lunar exploration project as the driving force, the China Interferometry Network (CVN) and the China Deep Space Network (CDSN) interferometry system have been built. Interferometry technology has been successfully applied to deep space exploration missions such as Chang'e 3, reentry and return tests, and Chang'e 4. The delay accuracy of differential interferometry is about 1ns, and the delay accuracy of same-beam interferometry is in the order of ps, reaching the world's advanced level. At the same time, the China Deep Space Exploration Network has successively participated in international deep space exploration missions such as Venus Express and IVS Joint Measurement, and the interferometry accuracy is comparable to that of foreign research institutions.

[0003] The CDSN interferometric measurement system has only had one baseline (JM01-KS01) for a long time, which limits its orbit determination support capability in deep space exploration missions. In 2017, the Namibia Station (NB01) and the South American Deep Space Station (NM01) were built, and CDSN has the capability of multi-station, multi-baseline, and global continuous interferometric measurement. The KS01-NM01 baseline is about 12,000 km long. Under the same interferometric measurement accuracy, the angle constraint capability of the spacecraft is stronger. At the same time, the interferometric measurement arc has been significantly improved. According to statistics, for the Chang'e-4 relay satellite, the interferometric measurement arc under the four-station condition is about 2 times higher than the JM01-KS01 baseline; for the Mars probe, the interferometric measurement arc under the four-station condition is about 3 times higher than the JM01-KS01 baseline, which lays the foundation for the subsequent improvement of the CDSN interferometric measurement system's orbit determination support capability.

[0004] Traditional correlation processing is the main mode for current deep space exploration missions. The raw data of interferometric measurements are first collected and recorded at the measuring station, and then transmitted to the relevant processing center in real time or afterwards. Finally, the relevant center outputs the observation data for orbit determination through data processing. The transmission of massive data in interferometric measurements has very high requirements on the bandwidth of the communication link. The observation data of the JM01-KS01 baseline can be transmitted through a dedicated line, which basically meets the real-time and post-transmission requirements. However, due to the influence of the communication link, the observation data of the two foreign measuring stations cannot be transmitted to the relevant processing center in real time (1 second of observation data takes at least 8 seconds to transmit), and post-transmission is also relatively difficult. This greatly limits the timeliness of interferometric measurement data processing, and directly limits the Deep Space Network interferometric measurement system's ability to support detector orbit determination during the mission.

[0005] Therefore, it is urgent to solve the problem in the prior art that the data transmission in the interference measurement process has high requirements on the communication link bandwidth.

[0006] The above information disclosed in the background technology section is only used to enhance the understanding of the background technology of the technology described in this article. Therefore, the background technology may contain certain information that does not form the prior art known in this country for those skilled in the art. Summary of the invention

[0007] The main purpose of the present application is to provide an interference measurement method, an interference measurement device, a computer-readable storage medium, a processor and an interference measurement system to solve the problem in the prior art that data transmission in the interference measurement process has high requirements on the communication link bandwidth.

[0008] According to one aspect of an embodiment of the present invention, there is provided an interference measurement method, comprising: receiving a plurality of first estimation data and a plurality of second estimation data, wherein the first estimation data is obtained by a first measuring station performing frequency estimation and phase estimation on first test data, and the second estimation data is obtained by a second measuring station performing frequency estimation and phase estimation on second test data, the first test data is observation data of a target detector by the first measuring station, and the second test data is observation data of the target detector by the second measuring station; obtaining an initial theoretical delay, wherein the initial theoretical delay is a theoretical delay at a first acquisition moment, and the theoretical delay is a theoretical transmission delay of the first measuring station relative to the second measuring station; and determining an interference measurement delay observation value of the first measuring station and the second measuring station according to each of the first estimation data, each of the second estimation data and the initial theoretical delay.

[0009] Optionally, the first estimated data includes a plurality of first Doppler frequencies and first residual phase estimation values ​​corresponding to the first acquisition moments, the second estimated data includes a plurality of second Doppler frequencies and second residual phase estimation values ​​corresponding to the second acquisition moments, and the interferometric measurement delay observation values ​​of the first measuring station and the second measuring station are determined according to each of the first estimated data, each of the second estimated data and the initial theoretical time delay, including: a first obtaining step, obtaining a first frequency coefficient, a first phase coefficient and a first radio frequency phase function according to the first Doppler frequency and the first residual phase estimation value corresponding to each of the first acquisition moments, wherein the first frequency coefficient, the first phase coefficient and the first radio frequency phase function are coefficients of the frequency function and coefficients of the phase function of the first measuring station, respectively. and a radio frequency phase function; a second obtaining step, according to the second Doppler frequency corresponding to each second acquisition moment and the second residual phase estimate, obtaining a second frequency coefficient, a second phase coefficient and a second radio frequency phase function, wherein the second frequency coefficient, the second phase coefficient and the second radio frequency phase function are respectively the coefficient of the frequency function, the coefficient of the phase function and the radio frequency phase function of the second measuring station; a first determining step, according to the first radio frequency phase function, the second radio frequency phase function, each of the first Doppler frequencies, each of the second Doppler frequencies, the first frequency coefficient, the second frequency coefficient, the initial theoretical time delay and a predetermined frequency, determining the interference measurement delay observation quantity, wherein the predetermined frequency is the downlink frequency of the first measuring station or the second measuring station.

[0010] Optionally, the frequency function of the first measuring station is a first frequency function, the phase function of the first measuring station is a first phase function, and according to the first Doppler frequency corresponding to each of the first acquisition moments and the first residual phase estimation value, a first frequency coefficient, a first phase coefficient and a first RF phase function are obtained, including: using a least squares method to establish a first initial frequency function f(t)=f dA t+f 0A And the first initial phase function The first initial frequency function represents the first Doppler frequency as a function of the first acquisition time t A The first initial phase function is obtained by integrating the first initial frequency function, and the coefficient of the first initial frequency function [f dA , f 0A ] and the coefficients of the first initial phase function are unknown; according to each of the first acquisition moments, the corresponding first Doppler frequencies and the first initial frequency function, determining the unknown coefficient in the first initial frequency as the first frequency coefficient Obtain the first frequency function; and determine, according to the first frequency coefficient and the first initial phase function, the unknown coefficient in the first initial phase function as the first phase coefficient Obtaining the first phase function; determining the first RF phase function according to the first phase coefficient and the first residual phase estimate corresponding to each of the first acquisition moments in, is the first phase coefficient, is the first residual phase estimate.

[0011] Optionally, the frequency function of the second measuring station is a second frequency function, the phase function of the second measuring station is a second phase function, and according to the second Doppler frequency corresponding to each second acquisition moment and the second residual phase estimation value, a second frequency coefficient, a second phase coefficient and a second RF phase function are obtained, including: using the least squares method to establish a second initial frequency function f(t)=f dB t+f 0B And the second initial phase function The second initial frequency function represents the second Doppler frequency as a function of the second acquisition time t B The second initial phase function is obtained by integrating the second initial frequency function, and the coefficient of the second initial frequency function [f dB , f 0B ] and the coefficients of the second initial phase function are unknown; according to each of the second acquisition moments, the corresponding second Doppler frequencies and the second initial frequency function, determining the unknown coefficient in the second initial frequency as the second frequency coefficient Obtain the second frequency function; and determine, according to the second frequency coefficient and the second initial phase function, the unknown coefficient in the second initial phase function as the second phase coefficient Obtain the second phase function; determine the second RF phase function according to the second phase coefficient and the second residual phase estimate corresponding to each of the second acquisition moments in, is the second phase coefficient, is the second residual phase estimate.

[0012] Optionally, the interference measurement delay observation amount is determined according to the first RF phase function, the second RF phase function, each first Doppler frequency, each second Doppler frequency, the first frequency coefficient, the second frequency coefficient, the initial theoretical delay and the predetermined frequency, including: setting the first phase coefficient in the first RF phase function to be the same as the second phase coefficient in the second RF phase function, and determining the residual delay estimate of the first station relative to the second station according to the set first RF phase function, the set second RF phase function, each first Doppler frequency and each second Doppler frequency; determining the compensation delay of the first station relative to the second station according to the predetermined frequency, the first frequency coefficient and the second frequency coefficient, wherein the predetermined frequency is the downlink frequency of the first station or the second station; determining the interference measurement delay observation amount according to the residual delay estimate, the compensation delay and the initial theoretical delay. Among them, τ Com is the compensation delay, is the residual delay estimate, τ Mdl (1) is the initial theoretical delay.

[0013] Optionally, setting the first phase coefficient in the first RF phase function to be the same as the second phase coefficient in the second RF phase function, and determining the residual delay estimate of the first station relative to the second station according to the set first RF phase function, the set second RF phase function, each first Doppler frequency, and each second Doppler frequency, includes: a second determination step, setting the first phase coefficient in the first RF phase function to be the same as the second phase coefficient in the second RF phase function, and determining, according to the set first RF phase function and the set second RF phase function, that the difference between the set first RF phase function and the set second RF phase function is the correlation phase between the first station and the second station; a synthesis step, when the correlation phase is within a predetermined range, performing bandwidth synthesis on the correlation phase, each first Doppler frequency, and each second Doppler frequency to obtain the residual delay estimate; and a loop step, when the correlation phase is not within the predetermined range, sequentially executing the first obtaining step, the second obtaining step, and the second determining step at least once until the correlation phase is within the predetermined range.

[0014] Optionally, determining the compensation delay of the first measuring station relative to the second measuring station according to the predetermined frequency, the first frequency coefficient and the second frequency coefficient includes: determining the compensation delay rate according to the predetermined frequency, the first frequency coefficient and the second frequency coefficient. Among them, f coef_B is the second frequency coefficient, f coef_A is the first frequency coefficient, f sky is the predetermined frequency; and according to the compensation delay rate, determines the compensation delay

[0015] According to another aspect of an embodiment of the present invention, an interference measurement device is also provided, including a receiving unit, an acquisition unit and a first determination unit, wherein the receiving unit is used to receive multiple first estimation data and multiple second estimation data, the first estimation data is obtained by the first measuring station performing frequency estimation and phase estimation on the first test data, the second estimation data is obtained by the second measuring station performing frequency estimation and phase estimation on the second test data, the first test data is the observation data of the first measuring station on the target detector, and the second test data is the observation data of the second measuring station on the target detector; the acquisition unit is used to obtain an initial theoretical delay, the initial theoretical delay is the theoretical delay at the first acquisition moment, and the theoretical delay is the theoretical transmission delay of the first measuring station relative to the second measuring station; the first determination unit is used to determine the interference measurement delay observation value of the first measuring station and the second measuring station according to each of the first estimation data, each of the second estimation data and the initial theoretical delay.

[0016] According to yet another aspect of the embodiments of the present invention, a computer-readable storage medium is provided, wherein the computer-readable storage medium includes a stored program, wherein the program executes any one of the methods described.

[0017] According to yet another aspect of an embodiment of the present invention, a processor is provided, wherein the processor is used to run a program, wherein any one of the methods is executed when the program is run.

[0018] According to another aspect of an embodiment of the present invention, an interference measurement system is also provided, comprising a first measuring station, a second measuring station and a processing center, wherein the first measuring station is used to observe a target detector to obtain first test data, and perform frequency estimation and phase estimation on the first test data to obtain first estimated data; the second measuring station is used to observe the target detector to obtain second test data, and perform frequency estimation and phase estimation on the second test data to obtain second estimated data; the processing center comprises one or more processors, a memory, a display device and one or more programs, wherein the one or more programs are stored in the memory and are configured to be executed by the one or more processors, and the one or more programs include methods for executing any one of the methods described.

[0019] According to the technical solution of the present application, in the interference measurement method, first, first estimation data is received after the first measuring station performs frequency and phase estimation on the first test data, and second estimation data is received after the second measuring station performs frequency and phase estimation on the second test data, wherein the first test data and the second test data are data obtained by the first measuring station and the second measuring station observing the same detector; then, an initial theoretical delay is obtained, and the initial theoretical delay is the theoretical transmission delay of the first measuring station relative to the second measuring station corresponding to the first acquisition moment; finally, the interference measurement delay observation value between the two measuring stations is determined according to each of the first estimation data, each of the second estimation data and the initial theoretical delay. Compared with the solution in the prior art that directly receives the observation data of the first measuring station and the second measuring station and then processes them to obtain the interference measurement delay observation quantity, due to the problem that the data transmission volume is large and the communication link bandwidth requirement is high during the data reception process, the method of the present application receives the data after the first measuring station and the second measuring station perform frequency and phase estimation. That is to say, the present application converts the transmission of original observation data into the transmission of measurement data. Compared with the existing original observation data, the data transmission volume of the present application is significantly reduced, which ensures that the bandwidth requirement for the communication link is low and the feasibility is strong. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] The drawings constituting part of the present application are used to provide a further understanding of the present application. The exemplary embodiments and descriptions of the present application are used to explain the present application and do not constitute an improper limitation on the present application. In the drawings:

[0021] Figure 1 A schematic flow chart of an interferometric measurement method according to an embodiment of the present application is shown;

[0022] Figure 2 A schematic diagram of an interferometric measurement device according to an embodiment of the present application is shown;

[0023] FIG3( a) shows a schematic diagram of differential phase simulation results of two measuring stations according to an embodiment of the present application;

[0024] FIG3( b ) is a schematic diagram showing the simulation results of differential delay of two measuring stations according to an embodiment of the present application;

[0025] Figure 4 A schematic diagram showing a delay estimation result changing with a signal-to-noise ratio under dynamic conditions according to an embodiment of the present application is shown;

[0026] Figure 5 A schematic diagram of residual phase estimation results of each channel of Chang'e-4 relay satellite JMS-KSH measured data according to an embodiment of the present application is shown;

[0027] FIG6( a) is a schematic diagram showing a result of estimating the baseband signal phase of a relay satellite signal after delay and Doppler compensation according to an embodiment of the present application;

[0028] FIG6( b ) is a schematic diagram showing the inter-station differential phase estimation result of the relay satellite signal after delay and Doppler compensation according to an embodiment of the present application;

[0029] FIG. 7( a ) is a schematic diagram showing a comparison of delay estimation results between the interference measurement method according to an embodiment of the present application and a conventional method;

[0030] FIG7( b ) is a schematic diagram showing a comparison of the geometric delay estimation results of the interference measurement method according to an embodiment of the present application and the traditional method. DETAILED DESCRIPTION

[0031] It should be noted that, in the absence of conflict, the embodiments and features in the embodiments of the present application can be combined with each other. The present application will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

[0032] In order to enable those skilled in the art to better understand the solution of the present application, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present application.

[0033] It should be noted that the terms "first", "second", etc. in the specification and claims of the present application and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequential order. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present application described here. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device comprising a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0034] It should be understood that when an element (such as a layer, film, region, or substrate) is described as being "on" another element, the element may be directly on the other element, or there may be intermediate elements. Moreover, in the specification and claims, when an element is described as being "connected" to another element, the element may be "directly connected" to the other element, or "connected" to the other element through a third element.

[0035] As mentioned in the background technology, the data transmission in the interference measurement process in the prior art has high requirements on the communication link bandwidth. In order to solve the above problem, in a typical embodiment of the present application, an interference measurement method, an interference measurement device, a computer-readable storage medium, a processor and an interference measurement system are provided.

[0036] According to an embodiment of the present application, an interference measurement method is provided.

[0037] Figure 1 is a flow chart of the interference measurement method according to an embodiment of the present application. Figure 1 As shown, the method comprises the following steps:

[0038] Step S101, receiving a plurality of first estimation data and a plurality of second estimation data, wherein the first estimation data is obtained by a first measuring station performing frequency estimation and phase estimation on first test data, and the second estimation data is obtained by a second measuring station performing frequency estimation and phase estimation on second test data, the first test data is observation data of a target detector by the first measuring station, and the second test data is observation data of the target detector by the second measuring station;

[0039] Step S102, obtaining an initial theoretical delay, where the initial theoretical delay is a theoretical delay at the first acquisition moment, and the theoretical delay is a theoretical transmission delay between the first measuring station and the second measuring station;

[0040] Step S103, determining the interferometric delay observation value of the first measuring station and the second measuring station according to the first estimation data, the second estimation data and the initial theoretical delay.

[0041] In the above-mentioned interference measurement method, first, first estimation data is received after the first measuring station performs frequency and phase estimation on the first test data, and second estimation data is received after the second measuring station performs frequency and phase estimation on the second test data, wherein the first test data and the second test data are data obtained by the first measuring station and the second measuring station observing the same detector; then, an initial theoretical delay is obtained, and the above-mentioned initial theoretical delay is the theoretical transmission delay of the above-mentioned first measuring station relative to the above-mentioned second measuring station corresponding to the first acquisition moment; finally, the interference measurement delay observation value between the two measuring stations is determined according to each of the above-mentioned first estimation data, each of the above-mentioned second estimation data and the above-mentioned initial theoretical delay. Compared with the solution in the prior art that directly receives the observation data of the first measuring station and the second measuring station and then processes them to obtain the interference measurement delay observation quantity, due to the problem that the data transmission volume is large and the communication link bandwidth requirement is high during the data reception process, the above method of the present application receives the data after the first measuring station and the second measuring station perform frequency and phase estimation. That is to say, the present application converts the transmission of original observation data into the transmission of measurement data. Compared with the existing original observation data, the data transmission volume of the present application is significantly reduced, ensuring that the bandwidth requirement for the communication link is low and the feasibility is strong.

[0042] It should be noted that the first measuring station can use any suitable known technology in the prior art to perform high-precision frequency and phase estimation on the first test data to obtain the first test data, and the second measuring station can do the same, which will not be elaborated in detail here.

[0043] According to a specific embodiment of the present application, the first estimated data includes a plurality of first Doppler frequencies and a first residual phase estimation value corresponding to a plurality of first acquisition moments, the second estimated data includes a plurality of second Doppler frequencies and a second residual phase estimation value corresponding to a plurality of second acquisition moments, and the interference measurement delay observation value of the first measuring station and the second measuring station is determined according to each of the first estimated data, each of the second estimated data and the initial theoretical time delay, including: a first obtaining step, according to the first Doppler frequency and the first residual phase estimation value corresponding to each of the first acquisition moments, obtaining a first frequency coefficient, a first phase coefficient and a first RF phase function, wherein the first frequency coefficient, the first phase coefficient and the first RF phase function are the coefficient of the frequency function of the first measuring station, the phase coefficient and the first phase coefficient, respectively. function coefficients and RF phase function; a second obtaining step, according to the second Doppler frequencies corresponding to each of the above-mentioned second acquisition moments and the above-mentioned second residual phase estimation value, obtaining the second frequency coefficient, the second phase coefficient and the second RF phase function, wherein the second frequency coefficient, the second phase coefficient and the second RF phase function are respectively the coefficients of the frequency function, the coefficients of the phase function and the RF phase function of the above-mentioned second measuring station; a first determining step, according to the first RF phase function, the second RF phase function, each of the above-mentioned first Doppler frequencies, each of the above-mentioned second Doppler frequencies, the above-mentioned first frequency coefficient, the above-mentioned second frequency coefficient, the above-mentioned initial theoretical time delay and a predetermined frequency, determining the above-mentioned interference measurement delay observation quantity, wherein the above-mentioned predetermined frequency is the downlink frequency of the above-mentioned first measuring station or the second measuring station. Compared with the prior art which requires multiple processing steps such as integer bit compensation, Doppler frequency compensation, decimal bit compensation, fringe rotation, etc. to obtain the above-mentioned interferometric measurement delay observation quantity, the present application can obtain the above-mentioned interferometric measurement delay observation quantity based on the Doppler frequency, frequency coefficient, phase coefficient and RF phase function of the two measuring stations, thereby ensuring that the process of obtaining the interferometric measurement delay observation quantity is relatively simple, and making the calculation efficiency of the interferometric measurement delay observation quantity higher.

[0044] Specifically, the above-mentioned frequency function is used to characterize the measured change function of the output frequency of the above-mentioned target detector with the acquisition time, the above-mentioned phase function is used to characterize the measured change function of the output phase of the above-mentioned detector with the acquisition time, and the above-mentioned RF phase function is used to characterize the measured change function of the RF phase value of the above-mentioned target detector with the acquisition time.

[0045] According to another specific embodiment of the present application, the frequency function of the first measuring station is a first frequency function, the phase function of the first measuring station is a first phase function, and the first frequency coefficient, the first phase coefficient and the first RF phase function are obtained according to the first Doppler frequency corresponding to each of the first acquisition moments and the first residual phase estimation value, including: using the least squares method to establish a first initial frequency function f(t)=f dA t+f 0A And the first initial phase function The first initial frequency function represents the first Doppler frequency as a function of the first acquisition time t A The first initial phase function is obtained by integrating the first initial frequency function. The coefficient of the first initial frequency function [f dA , f 0A ] and the coefficients of the first initial phase function are unknown; according to each of the above-mentioned first acquisition moments, the corresponding each of the above-mentioned first Doppler frequencies and the above-mentioned first initial frequency function, determine the unknown coefficient in the above-mentioned first initial frequency as the above-mentioned first frequency coefficient Obtain the first frequency function; determine the unknown coefficient in the first initial phase function as the first phase coefficient according to the first frequency coefficient and the first initial phase function Obtain the first phase function; determine the first RF phase function according to the first phase coefficient and the first residual phase estimate corresponding to each of the first acquisition moments in, is the first phase coefficient mentioned above, is the first residual phase estimation value. In this embodiment, a polynomial fitting is first performed on the Doppler frequency and phase of the first measuring station to obtain a frequency function and a phase function, and then the phase function is corrected according to the residual phase estimation value to obtain a radio frequency phase function, thereby realizing Doppler compensation. While further ensuring a low data transmission volume and a low bandwidth requirement of the communication link, the accuracy of the interference measurement is further ensured to be high.

[0046] In order to simplify the calculation, the first phase coefficient of the first measuring station is middle, Represents the initial phase, which is 0 here, that is, Then, the first RF phase function can be simplified to

[0047] In another specific embodiment of the present application, the frequency function of the second measuring station is a second frequency function, the phase function of the second measuring station is a second phase function, and the second frequency coefficient, the second phase coefficient and the second RF phase function are obtained according to the second Doppler frequency corresponding to each of the second acquisition moments and the second residual phase estimation value, including: using the least squares method to establish a second initial frequency function f(t)=f dB t+f 0B And the second initial phase function The second initial frequency function represents the second Doppler frequency as a function of the second acquisition time t B The second initial phase function is obtained by integrating the second initial frequency function. The coefficient of the second initial frequency function [f dB , f 0B ] and the coefficients of the second initial phase function are unknown; according to each of the above-mentioned second acquisition moments, the corresponding each of the above-mentioned second Doppler frequencies and the above-mentioned second initial frequency function, determine the unknown coefficient in the above-mentioned second initial frequency as the above-mentioned second frequency coefficient Obtain the second frequency function; and determine the unknown coefficient in the second initial phase function as the second phase coefficient according to the second frequency coefficient and the second initial phase function. Obtain the second phase function; determine the second RF phase function according to the second phase coefficient and the second residual phase estimate corresponding to each of the second acquisition moments in, is the second phase coefficient mentioned above, is the second residual phase estimation value. In this embodiment, a polynomial fitting is first performed on the Doppler frequency and phase of the second measuring station to obtain a frequency function and a phase function, and then the phase function is corrected according to the residual phase estimation value to obtain a radio frequency phase function, thereby realizing Doppler compensation. While further ensuring a low data transmission volume and a low bandwidth requirement of the communication link, the accuracy of the interference measurement is further ensured to be high.

[0048] In order to simplify the calculation, the second phase coefficient of the second measuring station is middle, Represents the initial phase, which is 0 here, that is, Then, the second phase function can be simplified to

[0049] In order to further achieve the acquisition of observation results with the same physical meaning, consistent results and comparable accuracy as the traditional correlation processing method, according to an embodiment of the present application, the interference measurement delay observation is determined according to the first RF phase function, the second RF phase function, each of the first Doppler frequencies, each of the second Doppler frequencies, the first frequency coefficient, the second frequency coefficient, the initial theoretical delay and the predetermined frequency, including: setting the first phase coefficient in the first RF phase function to be the same as the second phase coefficient in the second RF phase function, and determining the residual delay estimate of the first station relative to the second station according to the set first RF phase function, the set second RF phase function, each of the first Doppler frequencies and each of the second Doppler frequencies; determining the compensation delay of the first station relative to the second station according to the predetermined frequency, the first frequency coefficient and the second frequency coefficient, wherein the predetermined frequency is the downlink frequency of the first station or the second station; determining the interference measurement delay observation according to the residual delay estimate, the compensation delay and the initial theoretical delay. Among them, τ Com To compensate for the above delay, is the residual delay estimate, τ Mdl (1) is the initial theoretical delay. First, assume that the second phase coefficient of the second measuring station is equal to the first phase coefficient of the first measuring station, obtain the residual delay estimate of the first measuring station relative to the second measuring station, and then calculate the compensation delay to compensate for the part where the second phase coefficient is not equal to the first phase coefficient. Finally, according to the residual delay estimate, the compensation delay and the initial theoretical delay, determine the interferometric delay observation, which further ensures that a more accurate interferometric delay observation can be obtained.

[0050] Specifically, the first phase coefficient in the first RF phase function is set to be the same as the second phase coefficient in the second RF phase function, and the residual delay estimate of the first station relative to the second station is determined based on the set first RF phase function, the set second RF phase function, the first Doppler frequencies and the second Doppler frequencies, including: a second determination step, setting the first phase coefficient in the first RF phase function to be the same as the second phase coefficient in the second RF phase function, and determining, based on the set first RF phase function and the set second RF phase function, that the difference between the set first RF phase function and the set second RF phase function is the correlation phase between the first station and the second station, that is, the correlation phase The integration step, when the above-mentioned correlation phase is within the predetermined range, the bandwidth integration of the above-mentioned correlation phase, each of the above-mentioned first Doppler frequencies and each of the above-mentioned second Doppler frequencies is performed to obtain the above-mentioned residual delay estimation; the loop step, when the above-mentioned correlation phase is not within the predetermined range, the above-mentioned first obtaining step, the above-mentioned second obtaining step and the above-mentioned second determining step are sequentially performed at least once until the above-mentioned correlation phase is within the predetermined range. Theoretically, after the delay model compensation, the correlation phase of the two stations It should be relatively stable. In the case of instability, in order to further ensure the high accuracy of the interferometric delay observation, the above-mentioned cycle steps should be repeated. The changing trend of the data can be modeled to achieve iterative processing and model updating.

[0051] In addition, the bandwidth integration method can be any suitable method in the prior art for performing bandwidth integration processing on the frequencies and related phases of the two measuring stations. Those skilled in the art can flexibly set it according to actual conditions. The predetermined range is [-π, π].

[0052] In order to further realize the compensation of the residual delay estimation, thereby further ensuring that the observation results with the same physical meaning, consistent results and comparable accuracy as the traditional correlation processing method are obtained, in another specific embodiment of the present application, according to the predetermined frequency, the first frequency coefficient and the second frequency coefficient, the compensation delay of the first measuring station relative to the second measuring station is determined, including: according to the predetermined frequency, the first frequency coefficient and the second frequency coefficient, determining the compensation delay rate Among them, f coef_B is the second frequency coefficient mentioned above, f coef_A is the first frequency coefficient mentioned above, f sky is the predetermined frequency; according to the compensation delay rate, the compensation delay is determined

[0053] In a specific embodiment, the down-converted signal is assumed to be consistent with the downlink carrier frequency of the signal. Since the target detector moves relative to the first measuring station and the second measuring station, the frequency of the down-converted signal changes with the actual situation. The signal can be modeled as:

[0054]

[0055] Wherein, A and B represent the first and second stations of the interferometric measurement, respectively; i represents the number of DOR beacons in the downlink signal. The S band generally has 3 DOR signals (one of which is a carrier), and the X band generally has 5 DOR signals (one of which is a carrier); a i , b i are the DOR signal amplitudes of the first measuring station A and the second measuring station B respectively; are the initial phases of the i-th DOR signal of the first station A and the second station, respectively; τ g is the delay of the second station B relative to the first station A, which varies with time. The goal of interferometry is to estimate τ g .

[0056] Assume that the total data processing time is T, the preset processing cycle is Tp, and the number of processing cycles is N Tp , then the data collection time of the first station A is t A It can be expressed as t A =[0:(N Tp -1)]×Tp+0.5×Tp, for the second measuring station B, take the initial theoretical delay τ Mdl (1) By correcting the acquisition time scale of the second measuring station B, the delay model compensation is realized, and the data acquisition time t of the second measuring station B is obtained. B t B =[0:(N Tp -1)]×Tp+0.5×Tp+τ Mdl (1) = t + τ Mdl (1) By correcting the fitting time scale with the delay model, the delay model compensation is achieved.

[0057] The method further includes: determining the theoretical time delay corresponding to the plurality of acquisition moments according to the first estimation data and the plurality of estimation data, and establishing the theoretical time delay τ by using the least square method. Mdl The theoretical time delay function of (i) τ(i) = a1t 2 +a2t+a3, where i satisfies 1≤i≤N Tp The coefficients in the theoretical delay function are unknown. Substitute the theoretical delays corresponding to multiple acquisition moments into the theoretical delay function to calculate the coefficients of the theoretical delay function.

[0058] It should be noted that the steps shown in the flowcharts of the accompanying drawings can be executed in a computer system such as a set of computer executable instructions, and that, although a logical order is shown in the flowcharts, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0059] The embodiment of the present application also provides an interference measurement device. It should be noted that the interference measurement device of the embodiment of the present application can be used to execute the interference measurement method provided in the embodiment of the present application. The interference measurement device provided in the embodiment of the present application is introduced below.

[0060] Figure 2 Schematic diagram of an interferometric measurement device according to an embodiment of the present application. Figure 2As shown, the device includes a receiving unit 10, an acquiring unit 20 and a first determining unit 30, wherein the receiving unit 10 is used to receive a plurality of first estimation data and a plurality of second estimation data, wherein the first estimation data is obtained by the first station performing frequency estimation and phase estimation on the first test data, and the second estimation data is obtained by the second station performing frequency estimation and phase estimation on the second test data, the first test data is the observation data of the target detector by the first station, and the second test data is the observation data of the target detector by the second station; the acquiring unit 20 is used to acquire an initial theoretical delay, wherein the initial theoretical delay is the theoretical delay at the first acquisition moment, and the theoretical delay is the transmission delay of the first station relative to the second station in theory; the first determining unit 30 is used to determine the interferometric measurement delay observation value of the first station and the second station according to each of the first estimation data, each of the second estimation data and the initial theoretical delay.

[0061] In the above-mentioned interference measurement device, first estimation data after the first measuring station performs frequency and phase estimation on the first test data, and second estimation data after the second measuring station performs frequency and phase estimation on the second test data are received by the receiving unit, wherein the first test data and the second test data are data obtained by the first measuring station and the second measuring station observing the same detector; the initial theoretical delay is obtained by the acquisition unit, and the above-mentioned initial theoretical delay is the theoretical transmission delay of the above-mentioned first measuring station relative to the above-mentioned second measuring station corresponding to the first acquisition moment; and the interference measurement delay observation value between the two measuring stations is determined by the first determination unit according to each of the above-mentioned first estimation data, each of the above-mentioned second estimation data and the above-mentioned initial theoretical delay. Compared with the solution in the prior art that directly receives the observation data of the first measuring station and the second measuring station and then processes them to obtain the interference measurement delay observation quantity, due to the problem that the data transmission volume is large and the communication link bandwidth requirement is high during the data reception process, the above-mentioned device of the present application receives the data after the first measuring station and the second measuring station perform frequency and phase estimation. That is to say, the present application converts the transmission of original observation data into the transmission of measurement data. Compared with the existing original observation data, the data transmission volume of the present application is significantly reduced, ensuring that the bandwidth requirement for the communication link is low and the feasibility is strong.

[0062] It should be noted that the first measuring station can use any suitable known technology in the prior art to perform high-precision frequency and phase estimation on the first test data to obtain the first test data, and the second measuring station can do the same, which will not be elaborated in detail here.

[0063] According to a specific embodiment of the present application, the first estimation data includes a first Doppler frequency and a first residual phase estimation value corresponding to a plurality of first acquisition moments, the second estimation data includes a second Doppler frequency and a second residual phase estimation value corresponding to a plurality of second acquisition moments, and the first determination unit includes a first obtaining module, a second obtaining module and a first determination module, wherein the first obtaining module is used for the first obtaining step, and obtains a first frequency coefficient, a first phase coefficient and a first RF phase function according to the first Doppler frequency and the first residual phase estimation value corresponding to each of the first acquisition moments, wherein the first frequency coefficient, the first phase coefficient and the first RF phase function are respectively the coefficient of the frequency function, the coefficient of the phase function and the RF phase function of the first measuring station; the second The obtaining module is used in the second obtaining step, and obtains the second frequency coefficient, the second phase coefficient and the second RF phase function according to the second Doppler frequency corresponding to each of the second acquisition moments and the second residual phase estimate, wherein the second frequency coefficient, the second phase coefficient and the second RF phase function are respectively the coefficient of the frequency function, the coefficient of the phase function and the RF phase function of the second measuring station; the first determining module is used in the first determining step, and determines the interference measurement delay observation according to the first RF phase function, the second RF phase function, each of the first Doppler frequencies, each of the second Doppler frequencies, the first frequency coefficient, the second frequency coefficient, the initial theoretical time delay and the predetermined frequency, wherein the predetermined frequency is the downlink frequency of the first measuring station or the second measuring station. Compared with the prior art which requires multiple processing steps such as integer bit compensation, Doppler frequency compensation, decimal bit compensation, fringe rotation, etc. to obtain the above-mentioned interferometric measurement delay observation quantity, the present application can obtain the above-mentioned interferometric measurement delay observation quantity based on the Doppler frequency, frequency coefficient, phase coefficient and RF phase function of the two measuring stations, thereby ensuring that the process of obtaining the interferometric measurement delay observation quantity is relatively simple, and making the calculation efficiency of the interferometric measurement delay observation quantity higher.

[0064] Specifically, the above-mentioned frequency function is used to characterize the measured change function of the output frequency of the above-mentioned target detector with the acquisition time, the above-mentioned phase function is used to characterize the measured change function of the output phase of the above-mentioned detector with the acquisition time, and the above-mentioned RF phase function is used to characterize the measured change function of the RF phase value of the above-mentioned target detector with the acquisition time.

[0065] According to another specific embodiment of the present application, the frequency function of the first measuring station is a first frequency function, the phase function of the first measuring station is a first phase function, the first obtaining module includes a first establishing submodule, a first determining submodule, a second determining submodule and a third determining submodule, and the first establishing submodule is used to establish a first initial frequency function f(t)=f by using the least squares method. dA t+f 0A And the first initial phase function The first initial frequency function represents the first Doppler frequency as a function of the first acquisition time t A The first initial phase function is obtained by integrating the first initial frequency function. The coefficient of the first initial frequency function [f dA , f 0A ] and the coefficients of the first initial phase function The first determination submodule is used to determine the unknown coefficient in the first initial frequency as the first frequency coefficient according to each of the first acquisition moments, the corresponding first Doppler frequencies and the first initial frequency function. The first frequency function is obtained; the second determination submodule is used to determine the unknown coefficient in the first initial phase function as the first phase coefficient according to the first frequency coefficient and the first initial phase function. The first phase function is obtained; the third determination submodule is used to determine the first RF phase function according to the first phase coefficient and the first residual phase estimation value corresponding to each of the first acquisition moments. in, is the first phase coefficient mentioned above, is the first residual phase estimation value. In this embodiment, a polynomial fitting is first performed on the Doppler frequency and phase of the first measuring station to obtain a frequency function and a phase function, and then the phase function is corrected according to the residual phase estimation value to obtain a radio frequency phase function, thereby realizing Doppler compensation. While further ensuring a low data transmission volume and a low bandwidth requirement of the communication link, the accuracy of the interference measurement is further ensured to be high.

[0066] In order to simplify the calculation, the first phase coefficient of the first measuring station is middle, Represents the initial phase, which is 0 here, that is, Then, the first RF phase function can be simplified to

[0067] In another specific embodiment of the present application, the frequency function of the second measuring station is a second frequency function, the phase function of the second measuring station is a second phase function, and the second obtaining module includes a second establishing submodule, a fourth determining submodule, a fifth determining submodule and a sixth determining submodule, wherein the second establishing submodule is used to establish a second initial frequency function f(t)=f by using the least squares method. dB t+f 0B And the second initial phase function The second initial frequency function represents the second Doppler frequency as a function of the second acquisition time t B The second initial phase function is obtained by integrating the second initial frequency function. The coefficient of the second initial frequency function [f dB , f 0B ] and the coefficients of the second initial phase function The fourth determination submodule is used to determine the unknown coefficient in the second initial frequency as the second frequency coefficient according to each of the second acquisition moments, the corresponding second Doppler frequencies and the second initial frequency function. The above-mentioned second frequency function is obtained; the above-mentioned fifth determination submodule determines the unknown coefficient in the above-mentioned second initial phase function as the above-mentioned second phase coefficient according to the above-mentioned second frequency coefficient and the above-mentioned second initial phase function The sixth determination submodule is used to determine the second RF phase function according to the second phase coefficient and the second residual phase estimation value corresponding to each of the second acquisition moments. in, is the second phase coefficient mentioned above, is the second residual phase estimation value. In this embodiment, a polynomial fitting is first performed on the Doppler frequency and phase of the second measuring station to obtain a frequency function and a phase function, and then the phase function is corrected according to the residual phase estimation value to obtain a radio frequency phase function, thereby realizing Doppler compensation. While further ensuring a low data transmission volume and a low bandwidth requirement of the communication link, the accuracy of the interference measurement is further ensured to be high.

[0068] In order to simplify the calculation, the second phase coefficient of the second measuring station is middle, Represents the initial phase, which is 0 here, that is, Then, the second phase function can be simplified to

[0069] In order to further achieve the acquisition of observation results with the same physical meaning, consistent results and comparable accuracy as the traditional correlation processing method, according to an embodiment of the present application, the above-mentioned first determination module includes a setting submodule, a seventh determination submodule and an eighth determination submodule, wherein the above-mentioned setting submodule is used to set the above-mentioned first phase coefficient in the above-mentioned first RF phase function to be the same as the above-mentioned second phase coefficient in the above-mentioned second RF phase function, and determine the residual delay estimate of the above-mentioned first station relative to the above-mentioned second station according to the set first RF phase function, the set second RF phase function, each of the above-mentioned first Doppler frequencies and each of the above-mentioned second Doppler frequencies; the above-mentioned seventh determination submodule is used to determine the compensation delay of the first station relative to the above-mentioned second station according to a predetermined frequency, the above-mentioned first frequency coefficient and the above-mentioned second frequency coefficient, wherein the above-mentioned predetermined frequency is the downlink frequency of the above-mentioned first station or the second station; the above-mentioned eighth determination submodule is used to determine the above-mentioned interference measurement delay observation according to the above-mentioned residual delay estimate, the above-mentioned compensation delay and the above-mentioned initial theoretical delay. Among them, τ Com To compensate for the above delay, is the residual delay estimate, τ Mdl (1) is the initial theoretical delay. First, assume that the second phase coefficient of the second measuring station is equal to the first phase coefficient of the first measuring station, obtain the residual delay estimate of the first measuring station relative to the second measuring station, and then calculate the compensation delay to compensate for the part where the second phase coefficient is not equal to the first phase coefficient. Finally, according to the residual delay estimate, the compensation delay and the initial theoretical delay, determine the interferometric delay observation, which further ensures that a more accurate interferometric delay observation can be obtained.

[0070] Specifically, the setting submodule is also used in the second determination step to set the first phase coefficient in the first RF phase function to be the same as the second phase coefficient in the second RF phase function, and according to the set first RF phase function and the set second RF phase function, determine that the difference between the set first RF phase function and the set second RF phase function is the correlation phase between the first measuring station and the second measuring station, that is, the correlation phase The setting submodule is also used in the synthesis step. When the correlation phase is within the predetermined range, the bandwidth synthesis is performed on the correlation phase, each of the first Doppler frequencies and each of the second Doppler frequencies to obtain the residual delay estimation. The setting submodule is also used in the loop step. When the correlation phase is not within the predetermined range, the first obtaining step, the second obtaining step and the second determining step are sequentially performed at least once until the correlation phase is within the predetermined range. Theoretically, after the delay model compensation, the correlation phase of the two stations It should be relatively stable. In the case of instability, in order to further ensure the high accuracy of the interferometric delay observation, the above-mentioned cycle steps should be repeated. The changing trend of the data can be modeled to achieve iterative processing and model updating.

[0071] In addition, the bandwidth integration method can be any suitable method in the prior art for performing bandwidth integration processing on the frequencies and related phases of the two measuring stations. Those skilled in the art can flexibly set it according to actual conditions. The predetermined range is [-π, π].

[0072] In order to further realize the compensation of the residual delay estimation, thereby further ensuring that the observed quantity results with the same physical meaning, consistent results and comparable accuracy as the traditional correlation processing method are obtained, in another specific embodiment of the present application, the seventh determination submodule is also used to determine the compensation delay rate according to the predetermined frequency, the first frequency coefficient and the second frequency coefficient. Among them, f coef_B is the second frequency coefficient mentioned above, f coef_A is the first frequency coefficient mentioned above, f sky is the predetermined frequency; according to the compensation delay rate, the compensation delay is determined

[0073] In a specific embodiment, the down-converted signal is assumed to be consistent with the downlink carrier frequency of the signal. Since the target detector moves relative to the first measuring station and the second measuring station, the frequency of the down-converted signal changes with the actual situation. The signal can be modeled as:

[0074]

[0075] Wherein, A and B represent the first and second stations of the interferometric measurement, respectively; i represents the number of DOR beacons in the downlink signal. The S band generally has 3 DOR signals (one of which is a carrier), and the X band generally has 5 DOR signals (one of which is a carrier); a i , b i are the DOR signal amplitudes of the first measuring station A and the second measuring station B respectively; are the initial phases of the i-th DOR signal of the first station A and the second station, respectively; τ g is the delay of the second station B relative to the first station A, which varies with time. The goal of interferometry is to estimate τ g .

[0076] Assume that the total data processing time is T, the preset processing cycle is Tp, and the number of processing cycles is N Tp , then the data collection time of the first station A is t A It can be expressed as t A =[0:(N Tp -1)]×Tp+0.5×Tp, for the second measuring station B, take the initial theoretical delay τ Mdl (1) By correcting the acquisition time scale of the second measuring station B, the delay model compensation is realized, and the data acquisition time t of the second measuring station B is obtained. B t B =[0:(N Tp -1)]×Tp+0.5×Tp+τ Mdl (1) = t + τ Mdl (1) By correcting the fitting time scale with the delay model, the delay model compensation is achieved.

[0077] The device further includes a second determining unit, which is used to determine the theoretical time delay corresponding to the plurality of acquisition moments according to the first estimated data and the plurality of estimated data, and to establish the theoretical time delay τ by using the least square method. Mdl (i) The theoretical time delay function τ(t) = a1t 2 +a2t+a3, where i satisfies 1≤i≤N Tp The coefficients in the theoretical delay function are unknown. Substitute the theoretical delays corresponding to multiple acquisition moments into the theoretical delay function to calculate the coefficients of the theoretical delay function.

[0078] The interference measurement device includes a processor and a memory. The receiving unit, the acquisition unit, the first determination unit, etc. are all stored in the memory as program units. The processor executes the program units stored in the memory to implement corresponding functions.

[0079] The processor includes a kernel, and the kernel retrieves the corresponding program unit from the memory. One or more kernels can be set, and the kernel parameters can be adjusted to obtain a first test data, a second test data and a processing center. The first test station is used to observe the target detector, obtain first test data, and perform frequency estimation and phase estimation on the first test data to obtain first estimation data; the second test station is used to observe the target detector, obtain second test data, and perform frequency estimation and phase estimation on the second test data to obtain second estimation data; the processing center includes one or more processors, a memory, a display device and one or more programs, wherein the one or more programs are stored in the memory and are configured to be executed by the one or more processors, and the one or more programs include a method for executing any of the above methods. .

[0080] The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM, and the memory includes at least one memory chip.

[0081] An embodiment of the present invention provides a computer-readable storage medium on which a program is stored. When the program is executed by a processor, the above-mentioned interference measurement method is implemented.

[0082] An embodiment of the present invention provides a processor, and the processor is used to run a program, wherein the interference measurement method is executed when the program is run.

[0083] An embodiment of the present invention provides a device, the device comprising a processor, a memory, and a program stored in the memory and executable on the processor, and when the processor executes the program, at least the following steps are implemented:

[0084] Step S101, receiving a plurality of first estimation data and a plurality of second estimation data, wherein the first estimation data is obtained by a first measuring station performing frequency estimation and phase estimation on first test data, and the second estimation data is obtained by a second measuring station performing frequency estimation and phase estimation on second test data, the first test data is observation data of a target detector by the first measuring station, and the second test data is observation data of the target detector by the second measuring station;

[0085] Step S102, obtaining an initial theoretical delay, where the initial theoretical delay is a theoretical delay at the first acquisition moment, and the theoretical delay is a theoretical transmission delay between the first measuring station and the second measuring station;

[0086] Step S103, determining the interferometric delay observation value of the first measuring station and the second measuring station according to the first estimation data, the second estimation data and the initial theoretical delay.

[0087] The devices in this article can be servers, PCs, PADs, mobile phones, etc.

[0088] The present application also provides a computer program product, which, when executed on a data processing device, is suitable for executing a program for initializing at least the following method steps:

[0089] Step S101, receiving a plurality of first estimation data and a plurality of second estimation data, wherein the first estimation data is obtained by a first measuring station performing frequency estimation and phase estimation on first test data, and the second estimation data is obtained by a second measuring station performing frequency estimation and phase estimation on second test data, the first test data is observation data of a target detector by the first measuring station, and the second test data is observation data of the target detector by the second measuring station;

[0090] Step S102, obtaining an initial theoretical delay, where the initial theoretical delay is a theoretical delay at the first acquisition moment, and the theoretical delay is a theoretical transmission delay between the first measuring station and the second measuring station;

[0091] Step S103, determining the interferometric delay observation value of the first measuring station and the second measuring station according to the first estimation data, the second estimation data and the initial theoretical delay.

[0092] According to another typical embodiment of the present application, an interference measurement system is also provided, including a first measuring station, a second measuring station and a processing center, wherein the first measuring station is used to observe the target detector to obtain first test data, and perform frequency estimation and phase estimation on the first test data to obtain first estimation data; the second measuring station is used to observe the target detector to obtain second test data, and perform frequency estimation and phase estimation on the second test data to obtain second estimation data; the processing center includes one or more processors, a memory, a display device and one or more programs, wherein the one or more programs are stored in the memory and are configured to be executed by the one or more processors, and the one or more programs include programs for executing any one of the above methods.

[0093] In the above-mentioned interference measurement system, the first measuring station and the second measuring station respectively perform frequency and phase estimation on the test data obtained by the observation target detector, and send the estimated data to the processing center after obtaining the estimated data, and the above-mentioned processing center is used to execute any of the above-mentioned methods. Compared with the scheme in the prior art that directly receives the observation data of the first measuring station and the second measuring station, and then processes it to obtain the interference measurement delay observation, due to the problem that the data transmission volume is large and the communication link bandwidth requirement is high during the data reception process, in the above-mentioned system of the present application, the processor center receives the data after the first measuring station and the second measuring station perform frequency and phase estimation, that is, in the system of the present application, the first measuring station and the second measuring station convert the transmission of the original observation data into the transmission of the measurement data. Compared with the existing original observation data, the data transmission volume of the present application is significantly reduced, which ensures that the bandwidth requirement for the communication link is low and the feasibility is strong.

[0094] The data transmission requirements of the interferometric measurement system proposed in the present application are significantly reduced. In the interferometric measurement system, the first measuring station and the second measuring station can directly perform high-precision frequency and phase estimation after receiving and recording the original measurement data at the measuring station, and transmit the results to the processing center without transmitting the original observation data. That is, the interferometric measurement system proposed in the present application converts the transmission of original observation data into measurement information transmission, and the amount of data is reduced by several orders of magnitude, and can even be converted into telemetry data transmission, which significantly reduces the data transmission link requirements. On the other hand, in the interferometric measurement system of the present application, the data processing efficiency of the processing center is relatively high.

[0095] In order to enable those skilled in the art to more clearly understand the technical solution of the present application, it will be explained below in conjunction with specific embodiments.

[0096] Example 1

[0097] The interferometric measurement method of the present invention is used to calculate the interferometric measurement delay observation quantity. The sampling frequency is 50kHz (the sampling frequency commonly used in the interferometric measurement complex sampling mode), the carrier signal frequency is 12.3421kHz, the DOR signal and the carrier are spaced 1kHz apart; the initial phase is 0.5537rad; the delay rate is 10 -4(corresponding to a Doppler frequency of about 1Hz), data duration of 100s, integration time of 1s; the second station B is delayed by 2.4ms relative to the first station A; the relative residual delay of the two stations is estimated to be 2.4ms-100ns. The simulation results when the signal-to-noise ratio is 20dB are shown in Figure 3(a) and Figure 3(b). It can be seen from Figure 3(a) that the differential phase trends of the three channels are consistent and relatively stable. The differential phase trend reflects the influence of residual Doppler. Figure 3(b) shows the delay estimation results under the current conditions. It can be seen that the estimated value is consistent with the true value trend, and the deviation between the two is within 10ns; the statistical average deviation of the delay estimation is about 0.1671ns, and the random difference is about 4.6760ns.

[0098] Figure 4 A schematic diagram of the delay estimation results under dynamic conditions as the signal-to-noise ratio changes is given. It can be seen that as the signal-to-noise ratio changes, the deviation between the estimated value and the true value is on the order of 10ns; when the signal-to-noise ratio is greater than 12dB, the deviation between the two is less than 2ns, and the random error is less than 40ns.

[0099] Example 2

[0100] The processing and verification were carried out using the CE04 relay satellite Jiamusi (JMS)-Kashgar (KSH) deep space station data on a certain day in September 2018. Figure 5 The residual phase estimation of each channel at the two stations is given. It can be seen that the residual phase trend of each channel at the JMS station is relatively stable, and the random accuracy of the residual phase estimation is about 20mrad (20.4mrad, 11.8mrad, 17.3mrad); the residual phase of each channel at the KSH station has a small step at 160s, and the random accuracy of the residual phase estimation is about 80mrad (84.4mrad, 66.8mrad, 71.2mrad).

[0101] Figure 6(a) and Figure 6(b) show the differential phase estimation between stations after delay and Doppler compensation, where Figure 6(a) is the phase estimation of the baseband signal obtained by fitting based on phase modeling; Figure 6(b) is the phase estimation normalized to the interval [-π, π] after phase differentiation. As can be seen from Figure 6(b), the differential phases of the three channels of the two stations are relatively stable, but there is still a slight trend of change. After fitting and statistics, the random accuracy of the differential phases of the three channels are approximately 0.0518 rad, 0.0641 rad, and 0.0799 rad, respectively.

[0102] Figures 7(a) and 7(b) show the delay estimation comparison between the interference measurement method (StaCorr) of the present application and the traditional method (FXCorr). As can be seen from Figure 7(a), the random jitter of the residual delay of the traditional correlation processing is small, and the random accuracy is about 0.2859ns; the residual delay of the interference measurement method of the present application has an obvious step phenomenon, which is mainly affected by the phase estimation performance (as shown in Figure 6(b)), and the random accuracy is about 0.9109ns; the delay estimation performance can be improved by improving the phase estimation algorithm in the future. Figure 7(b) is a comparison of the geometric delay estimation and the residual delay estimation and the traditional correlation processing results. It can be seen that the geometric delay estimation value obtained by this method differs from the traditional correlation processing result by about 5ns, which includes the difference in link delay and device delay (the data source of the traditional correlation processing and the correlation processing method proposed in this application is inconsistent).

[0103] For 300s observation, the raw data of traditional correlation processing is about 2.3GB, and the estimated data after single-station processing by the interferometric measurement method proposed in this application is about 40KB, and the data transmission volume can be reduced by about 5 orders of magnitude, which is significantly reduced; the processing center implements correlation processing based on the estimated data of the measuring station, and obtains the same observation amount as the traditional correlation processing method, which can be directly used for the detector to determine the orbit. Therefore, the above-mentioned interferometric measurement method proposed in this application can provide the same observation amount as the traditional correlation processing while significantly reducing the data transmission requirements of the interferometric measurement system, and is expected to improve the feasibility of the interferometric measurement of the overseas baseline of my country's deep space network.

[0104] In the above embodiments of the present invention, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0105] In the several embodiments provided in this application, it should be understood that the disclosed technical content can be implemented in other ways. Among them, the device embodiments described above are only schematic. For example, the division of the above-mentioned units can be a logical function division. There may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of units or modules, which can be electrical or other forms.

[0106] The units described above as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple units. Some or all of the units may be selected according to actual needs to achieve the purpose of the present embodiment.

[0107] In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0108] If the above-mentioned integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions for a computer device (which can be a personal computer, a server or a network device, etc.) to perform all or part of the steps of the above-mentioned methods of each embodiment of the present invention. The aforementioned storage medium includes: U disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), mobile hard disk, magnetic disk or optical disk and other media that can store program codes.

[0109] From the above description, it can be seen that the above embodiments of the present application achieve the following technical effects:

[0110] 1) In the above-mentioned interference measurement method of the present application, firstly, first estimation data after the first measuring station performs frequency and phase estimation on the first test data, and second estimation data after the second measuring station performs frequency and phase estimation on the second test data, are received, wherein the first test data and the second test data are data obtained by the first measuring station and the second measuring station observing the same detector; then, an initial theoretical delay is obtained, and the above-mentioned initial theoretical delay is the theoretical transmission delay of the above-mentioned first measuring station relative to the above-mentioned second measuring station corresponding to the first acquisition moment; finally, the interference measurement delay observation value between the two measuring stations is determined according to each of the above-mentioned first estimation data, each of the above-mentioned second estimation data and the above-mentioned initial theoretical delay. Compared with the solution in the prior art that directly receives the observation data of the first measuring station and the second measuring station and then processes them to obtain the interference measurement delay observation quantity, due to the problem that the data transmission volume is large and the communication link bandwidth requirement is high during the data reception process, the above method of the present application receives the data after the first measuring station and the second measuring station perform frequency and phase estimation. That is to say, the present application converts the transmission of original observation data into the transmission of measurement data. Compared with the existing original observation data, the data transmission volume of the present application is significantly reduced, ensuring that the bandwidth requirement for the communication link is low and the feasibility is strong.

[0111] 2) In the above-mentioned interference measurement device of the present application, a receiving unit receives first estimation data after the first measuring station performs frequency and phase estimation on the first test data, and receives second estimation data after the second measuring station performs frequency and phase estimation on the second test data, wherein the first test data and the second test data are data obtained by the first measuring station and the second measuring station observing the same detector; an initial theoretical delay is obtained by an acquisition unit, and the above-mentioned initial theoretical delay is the theoretical transmission delay of the above-mentioned first measuring station relative to the above-mentioned second measuring station corresponding to the first acquisition moment; and a first determination unit determines the interference measurement delay observation value between the two measuring stations according to each of the above-mentioned first estimation data, each of the above-mentioned second estimation data and the above-mentioned initial theoretical delay. Compared with the solution in the prior art that directly receives the observation data of the first measuring station and the second measuring station and then processes them to obtain the interference measurement delay observation quantity, due to the problem that the data transmission volume is large and the communication link bandwidth requirement is high during the data reception process, the above-mentioned device of the present application receives the data after the first measuring station and the second measuring station perform frequency and phase estimation. That is to say, the present application converts the transmission of original observation data into the transmission of measurement data. Compared with the existing original observation data, the data transmission volume of the present application is significantly reduced, ensuring that the bandwidth requirement for the communication link is low and the feasibility is strong.

[0112] 3) In the above-mentioned interference measurement system of the present application, the first measuring station and the second measuring station respectively perform frequency and phase estimation on the test data obtained by the observation target detector, and send the estimated data to the processing center after obtaining the estimated data, and the above-mentioned processing center is used to execute any of the above-mentioned methods. Compared with the scheme in the prior art that directly receives the observation data of the first measuring station and the second measuring station, and then processes it to obtain the interference measurement delay observation, due to the problem that the data transmission volume is large and the communication link bandwidth requirement is high during the data reception process, in the above-mentioned system of the present application, the processor center receives the data after the first measuring station and the second measuring station perform frequency and phase estimation, that is, in the system of the present application, the first measuring station and the second measuring station convert the transmission of the original observation data into the transmission of the measurement data. Compared with the existing original observation data, the data transmission volume of the present application is significantly reduced, which ensures that the bandwidth requirement for the communication link is low and the feasibility is strong.

[0113] The above description is only the preferred embodiment of the present application and is not intended to limit the present application. For those skilled in the art, the present application may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. An interferometric measurement method, characterized in that: include: Receive a plurality of first estimation data and a plurality of second estimation data, wherein the first estimation data is obtained by the first measuring station performing frequency estimation and phase estimation on the first test data, the second estimation data is obtained by the second measuring station performing frequency estimation and phase estimation on the second test data, the first test data is observation data of the target detector by the first measuring station, and the second test data is observation data of the target detector by the second measuring station; Acquire an initial theoretical delay, where the initial theoretical delay is a theoretical delay at a first acquisition moment, and the theoretical delay is a theoretical transmission delay of the first measuring station relative to the second measuring station; An observed interferometric delay of the first measuring station and the second measuring station is determined according to each of the first estimation data, each of the second estimation data and the initial theoretical delay.

2. The method according to claim 1, characterized in that The first estimation data includes a plurality of first Doppler frequencies and first residual phase estimation values ​​corresponding to the first acquisition moments, the second estimation data includes a plurality of second Doppler frequencies and second residual phase estimation values ​​corresponding to the second acquisition moments, and determining the interferometric delay observation value of the first measuring station and the second measuring station according to each of the first estimation data, each of the second estimation data and the initial theoretical delay, including: A first obtaining step, obtaining a first frequency coefficient, a first phase coefficient and a first radio frequency phase function according to the first Doppler frequency and the first residual phase estimation value corresponding to each of the first acquisition moments, wherein the first frequency coefficient, the first phase coefficient and the first radio frequency phase function are respectively a coefficient of a frequency function, a coefficient of a phase function and a radio frequency phase function of the first measuring station; A second obtaining step, obtaining a second frequency coefficient, a second phase coefficient and a second radio frequency phase function according to the second Doppler frequency and the second residual phase estimation value corresponding to each of the second acquisition moments, wherein the second frequency coefficient, the second phase coefficient and the second radio frequency phase function are respectively a coefficient of a frequency function, a coefficient of a phase function and a radio frequency phase function of the second measuring station; The first determination step is to determine the interference measurement delay observation value according to the first RF phase function, the second RF phase function, each of the first Doppler frequencies, each of the second Doppler frequencies, the first frequency coefficient, the second frequency coefficient, the initial theoretical delay and a predetermined frequency, wherein the predetermined frequency is a downlink frequency of the first measuring station or the second measuring station.

3. The method according to claim 2, characterized in that The frequency function of the first measuring station is a first frequency function, the phase function of the first measuring station is a first phase function, and a first frequency coefficient, a first phase coefficient, and a first radio frequency phase function are obtained according to the first Doppler frequency and the first residual phase estimation value corresponding to each of the first acquisition moments, including: The first initial frequency function f(t)=f is established by the least square method. dA t+f 0A And the first initial phase function The first initial frequency function represents the first Doppler frequency as a function of the first acquisition time t A The first initial phase function is obtained by integrating the first initial frequency function, and the coefficient of the first initial frequency function [f dA , f 0A ] and the coefficients of the first initial phase function All unknown; According to each of the first acquisition moments, the corresponding first Doppler frequencies and the first initial frequency function, determining the unknown coefficient in the first initial frequency as the first frequency coefficient obtaining the first frequency function; According to the first frequency coefficient and the first initial phase function, determine the unknown coefficient in the first initial phase function as the first phase coefficient obtaining the first phase function; Determine the first RF phase function according to the first phase coefficient and the first residual phase estimate corresponding to each of the first acquisition moments in, is the first phase coefficient, is the first residual phase estimate.

4. The method according to claim 2, characterized in that: The frequency function of the second measuring station is a second frequency function, the phase function of the second measuring station is a second phase function, and a second frequency coefficient, a second phase coefficient and a second radio frequency phase function are obtained according to the second Doppler frequency and the second residual phase estimation value corresponding to each second acquisition moment, including: The second initial frequency function f(t)=f is established by the least square method. dB t+f 0B And the second initial phase function The second initial frequency function represents the second Doppler frequency as a function of the second acquisition time t B The second initial phase function is obtained by integrating the second initial frequency function, and the coefficient of the second initial frequency function [f dB , f 0B ] and the coefficients of the second initial phase function All unknown; According to each of the second acquisition moments, the corresponding second Doppler frequencies and the second initial frequency function, determine the unknown coefficient in the second initial frequency as the second frequency coefficient obtaining the second frequency function; According to the second frequency coefficient and the second initial phase function, determine the unknown coefficient in the second initial phase function as the second phase coefficient obtaining the second phase function; Determine the second RF phase function according to the second phase coefficient and the second residual phase estimate corresponding to each second acquisition moment in, is the second phase coefficient, is the second residual phase estimate.

5. The method according to any one of claims 2 to 4, characterized in that Determining the interferometric delay observation amount according to the first radio frequency phase function, the second radio frequency phase function, each of the first Doppler frequencies, each of the second Doppler frequencies, the first frequency coefficient, the second frequency coefficient, the initial theoretical delay, and a predetermined frequency, comprises: Setting the first phase coefficient in the first radio frequency phase function to be the same as the second phase coefficient in the second radio frequency phase function, and determining a residual delay estimate of the first station relative to the second station according to the set first radio frequency phase function, the set second radio frequency phase function, each of the first Doppler frequencies, and each of the second Doppler frequencies; Determine a compensation delay of the first station relative to the second station according to a predetermined frequency, the first frequency coefficient, and the second frequency coefficient, wherein the predetermined frequency is a downlink frequency of the first station or the second station; Determine the interferometric delay observation according to the residual delay estimate, the compensation delay and the initial theoretical delay Among them, τ Con is the compensation delay, is the residual delay estimate, τ Mdl (1) is the initial theoretical delay.

6. The method according to any one of claims 2 to 4, characterized in that The method further comprises: setting the first phase coefficient in the first radio frequency phase function to be the same as the second phase coefficient in the second radio frequency phase function, and determining a residual delay estimate of the first station relative to the second station according to the set first radio frequency phase function, the set second radio frequency phase function, each of the first Doppler frequencies, and each of the second Doppler frequencies, including: a second determination step, setting the first phase coefficient in the first radio frequency phase function to be the same as the second phase coefficient in the second radio frequency phase function, and determining, based on the set first radio frequency phase function and the set second radio frequency phase function, a difference between the set first radio frequency phase function and the set second radio frequency phase function as the correlation phase between the first measuring station and the second measuring station; a synthesis step, when the correlation phase is within a predetermined range, performing bandwidth synthesis on the correlation phase, each of the first Doppler frequencies and each of the second Doppler frequencies to obtain the residual delay estimate; A loop step, when the relevant phase is not within the predetermined range, sequentially performing the first obtaining step, the second obtaining step and the second determining step at least once until the relevant phase is within the predetermined range.

7. The method according to any one of claims 2 to 4, characterized in that Determining a compensation delay of a first measuring station relative to the second measuring station according to a predetermined frequency, the first frequency coefficient, and the second frequency coefficient includes: Determine the compensation delay rate according to the predetermined frequency, the first frequency coefficient and the second frequency coefficient Among them, f coef_B is the second frequency coefficient, f coef_A is the first frequency coefficient, f sky is the predetermined frequency; Determine the compensation delay according to the compensation delay rate 8. An interferometric measurement device, characterized in that: include: a receiving unit, configured to receive a plurality of first estimation data and a plurality of second estimation data, wherein the first estimation data is obtained by the first measuring station performing frequency estimation and phase estimation on the first test data, and the second estimation data is obtained by the second measuring station performing frequency estimation and phase estimation on the second test data, the first test data is observation data of the target detector by the first measuring station, and the second test data is observation data of the target detector by the second measuring station; An acquisition unit, configured to acquire an initial theoretical delay, wherein the initial theoretical delay is a theoretical delay at a first acquisition moment, and the theoretical delay is a theoretical transmission delay of the first measuring station relative to the second measuring station; The first determining unit is used to determine the interferometric delay observation value of the first measuring station and the second measuring station according to each of the first estimation data, each of the second estimation data and the initial theoretical delay.

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium includes a stored program, wherein the program executes the method of any one of claims 1 to 7.

10. A processor, characterized in that: The processor is used to run a program, wherein the program executes the method according to any one of claims 1 to 7 when running.

11. An interferometric measurement system, characterized in that: include: A first measuring station is used to observe the target detector to obtain first test data, and perform frequency estimation and phase estimation on the first test data to obtain first estimated data; A second measuring station is used to observe the target detector to obtain second test data, and perform frequency estimation and phase estimation on the second test data to obtain second estimated data; A processing center, comprising one or more processors, a memory, a display device and one or more programs, wherein the one or more programs are stored in the memory and are configured to be executed by the one or more processors, and the one or more programs include methods for executing any one of claims 1 to 7.

Images