Interference processing method and system suitable for detectors with limited prior orbit accuracy

Abstract

The present invention relates to the field of radio interference measurement technology, and specifically discloses an interference processing method and system suitable for use when the a priori orbit accuracy of the detector is limited, the method comprising: obtaining a polynomial model of the residual differential delay rate sequence of the target detector's acquisition signal data in multiple integration periods; performing two integration processes on the polynomial model of the residual differential delay rate sequence to obtain and calculate the residual delay and residual delay rate based on the polynomial model of the residual delay sequence; based on the residual delay and the residual delay rate, and in combination with the a priori delay model of the target detector, determining the final delay observation and the final delay rate observation. The present invention can improve the observation accuracy when the a priori orbit accuracy of the detector is limited, and overcomes the problem of low observation accuracy caused by large residual delay noise, fuzzy delay rate, and large interference fringe noise.

Description

Background Art

[0002] Very Long Baseline Interferometry (VLBI) is a technical means to obtain high-precision angular position information of the detector, supporting the orbit determination of the detector. In the traditional correlation processing process, it is necessary to apply information such as the detector's prior ephemeris. Due to the influence of various factors such as sparse tracking arc segments and detector maneuvering control, the accuracy of the detector's prior predicted orbit may be limited. At this time, the correlation processing guidance model directly calculated using the prior ephemeris has poor accuracy, resulting in the inability to effectively process data to obtain observations.

[0003] Therefore, it is urgent to provide a technical solution to solve the above problems. Summary of the invention

[0004] In order to solve the above technical problems, the present invention provides an interference processing method and system suitable for use when the accuracy of a priori orbit of the detector is limited.

[0005] In a first aspect, the present invention provides an interference processing method suitable for use with a detector with limited prior orbit accuracy. The technical solution of the method is as follows:

[0006] Obtaining a polynomial model of a residual differential delay rate sequence of the collected signal data of the target detector in multiple integration periods;

[0007] Performing two integration processes on the polynomial model of the residual differential delay rate sequence to obtain and calculate the residual delay and the residual delay rate according to the polynomial model of the residual delay sequence;

[0008] Based on the residual delay and the residual delay rate, and in combination with a priori delay model of the target detector, a final delay observation and a final delay rate observation are determined.

[0009] The beneficial effects of an interference processing method of the present invention applicable to detectors with limited prior orbit accuracy are as follows:

[0010] The method of the present invention can improve the observation accuracy when the detector's prior orbit accuracy is limited, and overcomes the problem of low observation accuracy caused by large residual delay noise, ambiguity in the delay rate, and large interference fringe noise.

[0011] On the basis of the above scheme, an interference processing method of the present invention applicable to the case where the accuracy of the detector prior orbit is limited can also be improved as follows.

[0012] In an optional manner, the method further includes:

[0013] A residual time delay rate sequence of the target detector's collected signal data in multiple integration periods is constructed, and the residual time delay rate sequence is differentially processed to obtain the residual differential time delay rate sequence, and a polynomial model of the residual differential time delay rate sequence is constructed.

[0014] In an optional manner, the step of performing two integration processes on the polynomial model of the residual differential delay rate sequence to obtain the polynomial model of the residual delay sequence includes:

[0015] The polynomial model of the residual differential delay rate sequence is integrated to obtain the polynomial model of the residual delay rate sequence, and the polynomial model of the residual delay rate sequence is integrated to obtain the polynomial model of the residual delay sequence.

[0016] In an optional manner, the step of determining a final delay observation and a final delay rate observation based on the residual delay and the residual delay rate and in combination with a priori delay model of the target detector includes:

[0017] The final delay observation is determined based on the a priori delay model and the residual delay, and the final delay rate observation is determined based on the a priori delay model and the residual delay rate.

[0018] In a second aspect, the present invention provides an interference processing system suitable for use with a detector under limited a priori orbit accuracy, and the technical solution of the system is as follows:

[0019] It includes: a data acquisition module, a first processing module and a second processing module;

[0020] The data acquisition module is used to: obtain a polynomial model of a residual differential delay rate sequence of the target detector's collected signal data in multiple integration periods;

[0021] The first processing module is used to: perform two integration processes on the polynomial model of the residual differential delay rate sequence to obtain and calculate the residual delay and the residual delay rate according to the polynomial model of the residual delay sequence;

[0022] The second processing module is used to determine a final delay observation and a final delay rate observation based on the residual delay and the residual delay rate and in combination with a priori delay model of the target detector.

[0023] The beneficial effects of an interference processing system of the present invention applicable to detectors with limited prior orbit accuracy are as follows:

[0024] The system of the present invention can improve the observation accuracy when the detector's prior orbit accuracy is limited, and overcomes the problem of low observation accuracy caused by large residual delay noise, ambiguity in the delay rate, and large interference fringe noise.

[0025] On the basis of the above scheme, the interference processing system of the present invention suitable for the detector with limited prior orbit accuracy can also be improved as follows.

[0026] In an optional manner, the method further includes: a building module; the building module is used to:

[0027] A residual time delay rate sequence of the target detector's collected signal data in multiple integration periods is constructed, and the residual time delay rate sequence is differentially processed to obtain the residual differential time delay rate sequence, and a polynomial model of the residual differential time delay rate sequence is constructed.

[0028] In an optional manner, the step of performing two integration processes on the polynomial model of the residual differential delay rate sequence in the first processing module to obtain the polynomial model of the residual delay sequence includes:

[0029] The polynomial model of the residual differential delay rate sequence is integrated to obtain the polynomial model of the residual delay rate sequence, and the polynomial model of the residual delay rate sequence is integrated to obtain the polynomial model of the residual delay sequence.

[0030] In an optional manner, the second processing module is specifically used to:

[0031] The final delay observation is determined based on the a priori delay model and the residual delay, and the final delay rate observation is determined based on the a priori delay model and the residual delay rate.

[0032] In a third aspect, a technical solution of an electronic device of the present invention is as follows:

[0033] It includes a memory, a processor, and a program stored in the memory and running on the processor. When the processor executes the program, it implements the steps of the interference processing method of the present invention that is applicable to the detector with limited prior orbit accuracy.

[0034] In a fourth aspect, the present invention provides a computer-readable storage medium having the following technical solution:

[0035] Instructions are stored in the computer-readable storage medium. When the computer-readable storage medium reads the instructions, the computer-readable storage medium executes the steps of the interference processing method of the present invention that is applicable to the detector with limited prior orbit accuracy.

[0036] The above description is only an overview of the technical solution of the present invention. In order to more clearly understand the technical means of the present invention, it can be implemented according to the contents of the specification. In order to make the above and other purposes, features and advantages of the present invention more obvious and easy to understand, the specific implementation methods of the present invention are listed below. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] The accompanying drawings are only used to illustrate the embodiments and are not to be considered as limiting the present invention. In addition, the same reference symbols are used to represent the same components throughout the accompanying drawings. In the accompanying drawings:

[0038] Figure 1 It is a flow chart of an embodiment of an interference processing method applicable to a detector with limited a priori orbit accuracy according to the present invention;

[0039] Figure 2 Schematic diagram of residual delay and delay rate obtained for interference processing;

[0040] Figure 3 Schematic diagram of the cross power spectrum and interference fringes obtained for the interference processing;

[0041] Figure 4 is a schematic diagram of residual delay and residual delay rate of the present invention;

[0042] Figure 5 is a schematic diagram of the cross-correlation spectrum and interference fringes of the present invention;

[0043] Figure 6 is a phase diagram of a high-order harmonic signal of the present invention;

[0044] Figure 7 It is a structural schematic diagram of an embodiment of an interference processing system applicable to a detector with limited a priori orbit accuracy according to the present invention;

[0045] Figure 8 The figure is a schematic structural diagram of an embodiment of an electronic device of the present invention. DETAILED DESCRIPTION

[0046] Exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. Although exemplary embodiments of the present invention are shown in the drawings, it should be understood that the present invention can be implemented in various forms and should not be limited to the embodiments set forth herein.

[0047] Figure 1 FIG. 2 is a flow chart showing an embodiment of an interference processing method provided by the present invention that is applicable to a detector with limited prior orbit accuracy. Figure 1 As shown, the following steps are included:

[0048] S1. Obtain a polynomial model of the residual differential delay rate sequence of the collected signal data of the target detector in multiple integration periods.

[0049] The target detector is: In this embodiment, there is a detector with limited a priori orbit accuracy. The collected signal data is used to reflect the energy strength of the microwave signal emitted by the target detector. The type of detector is not limited, such as satellites flying around the moon, Mars, etc., which can be selected according to actual conditions and are not limited here.

[0050] It should be noted that in this embodiment, the collected signal data is subjected to traditional correlation processing to obtain a residual delay sequence τ(n) and a residual delay rate sequence τ′(n) of multiple integration periods N. Among them, traditional correlation processing is a classic algorithm for interference processing, and the main steps are: ① integer bit compensation; ② fringe rotation; ③ FFT transformation; ④ fractional bit compensation; ⑤ coherent accumulation; ⑥ bandwidth synthesis.

[0051] S2. Performing two integration processes on the polynomial model of the residual differential delay rate sequence to obtain and calculate the residual delay and the residual delay rate based on the polynomial model of the residual delay sequence. Specifically:

[0052] 1) constructing a residual time delay rate sequence of the target detector's acquisition signal data in multiple integration periods, performing differential processing on the residual time delay rate sequence to obtain the residual differential time delay rate sequence, and constructing a polynomial model of the residual differential time delay rate sequence.

[0053] Wherein, the residual delay rate sequence τ′(n), n=1,2...N, N is the number of integration periods. The residual differential delay rate sequence τ″(n): n=1,2...,N-1, T is the integration period set during the correlation processing. The polynomial model of the residual differential delay rate sequence is: k=polyfit(T,τ″,M), T is the time vector; τ″=[τ″(1),τ″(2),...,τ″(N-1)], τ″ is the residual differential delay rate vector, τ″(N-1) is the N-1th residual differential delay rate; M is the polynomial order, k=[k M ,k M-1 ,...,k 0 ] is the polynomial coefficient vector, k M is the Mth-order polynomial coefficient of the residual differential delay rate sequence, and polyfit() represents the least squares fitting calculation.

[0054] 2) Integrating the polynomial model of the residual differential delay rate sequence to obtain the polynomial model of the residual delay rate sequence, and integrating the polynomial model of the residual delay rate sequence to obtain the polynomial model of the residual delay sequence.

[0055] Among them, the polynomial model of the residual delay rate sequence τ′(n) is:

[0056] p=[p M+1 ,...,p 2 ,p 1 ,p 0 ]

[0057]

[0058]

[0059] polyval() represents polynomial calculation, n 0 and n1 represents a data interval with a lower noise level in the residual delay rate sequence τ′(n). M+1 is the M+1th order coefficient of the polynomial of the residual delay rate sequence, p 0 is the polynomial 0th-order coefficient (constant term) of the residual delay rate sequence. The polynomial model of the residual delay sequence τ(n) is:

[0060] q=[q M+2 ,q M+1 ,...q 2 ,q 1 ,q 0 ]

[0061]

[0062]

[0063] q M+2 is the M+2th order coefficient of the polynomial of the residual delay sequence, q 0 is the 0th order coefficient (constant term) of the polynomial of the residual delay sequence.

[0064] 3) Calculating the residual delay and the residual delay rate according to the polynomial model of the residual delay sequence.

[0065] Among them, obtain the relevant processing guidance model: new (t)=Ψ(t)+polyval(q,t); Ψ(t) is the prior delay model, new (t) = Ψ(t) + polyval(q, t) to perform traditional correlation processing to obtain the residual delay and residual delay rate.

[0066] S3. Based on the residual delay and the residual delay rate, and in combination with a priori delay model of the target detector, determine a final delay observation and a final delay rate observation.

[0067] Specifically, the final delay observation is determined based on the priori delay model and the residual delay, and the final delay rate observation is determined based on the priori delay model and the residual delay rate.

[0068] Among them, based on the formula: obs =Ψ(t)+τ new , determine the final time-delay observation; Ψ obs is the final time-delay observation, τ new is the residual delay. Based on the formula: Determine the final delay rate observation, τ n ' ew is the residual delay rate, Ψ′ obs is the final delay rate observation.

[0069] It should be noted that the interferometric observation verification was carried out using the measuring station, in which the detector's prior orbit accuracy was low, resulting in large residual delay noise obtained by direct correlation processing, ambiguity in the delay rate, and large interference fringe noise. Figure 2 Schematic diagram of residual delay and delay rate obtained for interference processing, Figure 3 Schematic diagram of the cross power spectrum and interference fringes obtained by interference processing. Figure 4 A schematic diagram showing the residual delay and residual delay rate of the present application is shown, Figure 5 A schematic diagram showing the cross-correlation spectrum and interference fringes of the present application is shown, Figure 6 The phase diagram of the high-order harmonic signal obtained by the present application is shown. It can be seen that the technical solution of this embodiment can improve the observation accuracy when the detector's prior orbit accuracy is limited, and overcome the problem of low observation accuracy caused by large residual delay noise, ambiguity in the delay rate, and large interference fringe noise.

[0070] Figure 7 FIG. 2 is a schematic diagram showing a structure of an embodiment of an interference processing system 200 provided by the present invention, which is applicable to a detector with limited prior orbit accuracy. Figure 7 As shown, the system 200 includes: a data acquisition module 210, a first processing module 220 and a second processing module 230;

[0071] The data acquisition module 210 is used to: obtain a polynomial model of a residual differential delay rate sequence of the target detector's collected signal data in multiple integration periods;

[0072] The first processing module 220 is used to: perform two integration processes on the polynomial model of the residual differential delay rate sequence to obtain and calculate the residual delay and the residual delay rate according to the polynomial model of the residual delay sequence;

[0073] The second processing module 230 is used to determine a final delay observation and a final delay rate observation based on the residual delay and the residual delay rate and in combination with a priori delay model of the target detector.

[0074] In an optional manner, the method further includes: a building module; the building module is used to:

[0075] A residual time delay rate sequence of the target detector's collected signal data in multiple integration periods is constructed, and the residual time delay rate sequence is differentially processed to obtain the residual differential time delay rate sequence, and a polynomial model of the residual differential time delay rate sequence is constructed.

[0076] In an optional manner, the step of performing two integration processes on the polynomial model of the residual differential delay rate sequence in the first processing module 220 to obtain the polynomial model of the residual delay sequence includes:

[0077] The polynomial model of the residual differential delay rate sequence is integrated to obtain the polynomial model of the residual delay rate sequence, and the polynomial model of the residual delay rate sequence is integrated to obtain the polynomial model of the residual delay sequence.

[0078] In an optional manner, the second processing module 230 is specifically used to:

[0079] The final delay observation is determined based on the a priori delay model and the residual delay, and the final delay rate observation is determined based on the a priori delay model and the residual delay rate.

[0080] The technical solution of this embodiment can improve the observation accuracy when the detector's prior orbit accuracy is limited, and overcomes the problem of low observation accuracy caused by large residual delay noise, ambiguity in the delay rate, and large interference fringe noise.

[0081] For the above-mentioned parameters and steps for each module to implement corresponding functions in the interference processing system 200 applicable to the situation where the prior accuracy of the detector is limited, please refer to the parameters and steps in the embodiment of the interference processing method applicable to the situation where the prior accuracy of the detector is limited, and no further details will be given here.

[0082] like Figure 8As shown, an electronic device 300 according to an embodiment of the present invention includes a processor 320, the processor 320 is coupled to a memory 310, and the memory 310 stores at least one computer program 330, and the at least one computer program 330 is loaded and executed by the processor 320, so that the electronic device 300 implements any of the above-mentioned interference processing methods applicable to the detector under the condition of limited prior orbit accuracy, specifically:

[0083] The electronic device 300 may have relatively large differences due to different configurations or performances, and may include one or more processors 320 (Central Processing Units, CPU) and one or more memories 310, wherein the one or more memories 310 store at least one computer program 330, and the at least one computer program 330 is loaded and executed by the one or more processors 320, so that the electronic device 300 implements any of the interference processing methods provided in the above embodiments that are applicable to the detector under the condition of limited prior orbit accuracy. Of course, the electronic device 300 may also have components such as a wired or wireless network interface, a keyboard, and an input and output interface for input and output. The electronic device 300 may also include other components for realizing device functions, which will not be repeated here.

[0084] A computer-readable storage medium according to an embodiment of the present invention stores at least one computer program, and the at least one computer program is loaded and executed by a processor so that a computer implements any one of the above-mentioned interference processing methods applicable to the detector under the condition of limited prior orbit accuracy.

[0085] Optionally, the computer-readable storage medium may be a read-only memory (ROM), a random access memory (RAM), a compact disc (CD-ROM), a magnetic tape, a floppy disk, an optical data storage device, and the like.

[0086] In an exemplary embodiment, a computer program product or a computer program is also provided, the computer program product or the computer program comprising computer instructions, the computer instructions being stored in a computer-readable storage medium. A processor of an electronic device reads the computer instructions from the computer-readable storage medium, and the processor executes the computer instructions, so that the electronic device executes any of the above-mentioned interference processing methods applicable to the detector under the condition of limited a priori orbit accuracy.

[0087] It should be noted that the terms "first", "second", etc. in the specification and claims of the present application are used to distinguish similar objects and represent the definition of a specific order or sequence. The order of use of similar objects can be interchanged where appropriate, so that the embodiments of the present application described herein can be implemented in an order other than the order shown or described.

[0088] Those skilled in the art know that the present invention can be implemented as a system, method or computer program product. Therefore, the present disclosure can be specifically implemented in the following forms, namely: it can be complete hardware, it can be complete software (including firmware, resident software, microcode, etc.), or it can be a combination of hardware and software, which is generally referred to as "circuit", "module" or "system" herein. In addition, in some embodiments, the present invention can also be implemented in the form of a computer program product in one or more computer-readable media, and the computer-readable medium contains computer-readable program code.

[0089] Any combination of one or more computer-readable media can be used. Computer-readable media can be computer-readable signal media or computer-readable storage media. Computer-readable storage media can be, for example, but not limited to, electrical, magnetic, optical, electromagnetic, infrared, or semiconductor systems, devices or devices, or any combination of the above. More specific examples of computer-readable storage media (a non-exhaustive list) include: an electrical connection with one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above. In the present application, a computer-readable storage medium can be any tangible medium containing or storing a program, which can be used by an instruction execution system, device or device or used in combination with it.

[0090] Although the embodiments of the present invention have been shown and described above, it is to be understood that the above embodiments are exemplary and are not to be construed as limitations of the present invention. A person skilled in the art may change, modify, replace and vary the above embodiments within the scope of the present invention.

Claims

1. An interference processing method suitable for detectors with limited prior orbit accuracy, characterized in that: include: Obtaining a polynomial model of a residual differential delay rate sequence of the collected signal data of the target detector in multiple integration periods; Performing two integration processes on the polynomial model of the residual differential delay rate sequence to obtain and calculate the residual delay and the residual delay rate according to the polynomial model of the residual delay sequence; Based on the residual delay and the residual delay rate, and in combination with a priori delay model of the target detector, a final delay observation and a final delay rate observation are determined.

2. The interference processing method according to claim 1, which is suitable for the detector with limited prior orbit accuracy, is characterized in that: Also includes: A residual time delay rate sequence of the target detector's collected signal data in multiple integration periods is constructed, and the residual time delay rate sequence is differentially processed to obtain the residual differential time delay rate sequence, and a polynomial model of the residual differential time delay rate sequence is constructed.

3. The interference processing method applicable to the detector with limited prior orbit accuracy according to claim 1 is characterized in that: The step of performing two integration processes on the polynomial model of the residual differential delay rate sequence to obtain the polynomial model of the residual delay sequence includes: The polynomial model of the residual differential delay rate sequence is integrated to obtain the polynomial model of the residual delay rate sequence, and the polynomial model of the residual delay rate sequence is integrated to obtain the polynomial model of the residual delay sequence.

4. The interference processing method according to claim 1, which is suitable for the detector with limited prior orbit accuracy, is characterized in that: The step of determining a final delay observation and a final delay rate observation based on the residual delay and the residual delay rate and in combination with a priori delay model of the target detector includes: The final delay observation is determined based on the a priori delay model and the residual delay, and the final delay rate observation is determined based on the a priori delay model and the residual delay rate.

5. An interferometric processing system suitable for use with detectors with limited prior orbit accuracy, characterized in that: include: A data acquisition module, a first processing module and a second processing module; The data acquisition module is used to: obtain a polynomial model of a residual differential delay rate sequence of the target detector's collected signal data in multiple integration periods; The first processing module is used to: perform two integration processes on the polynomial model of the residual differential delay rate sequence to obtain and calculate the residual delay and the residual delay rate according to the polynomial model of the residual delay sequence; The second processing module is used to determine a final delay observation and a final delay rate observation based on the residual delay and the residual delay rate and in combination with a priori delay model of the target detector.

6. The interferometric processing system according to claim 5, which is suitable for use with detectors with limited a priori orbit accuracy, is characterized in that: Also includes: Building blocks; the building blocks are used to: A residual time delay rate sequence of the target detector's collected signal data in multiple integration periods is constructed, and the residual time delay rate sequence is differentially processed to obtain the residual differential time delay rate sequence, and a polynomial model of the residual differential time delay rate sequence is constructed.

7. The interferometric processing system suitable for detectors with limited a priori orbit accuracy according to claim 5, characterized in that: The step of performing two integration processes on the polynomial model of the residual differential delay rate sequence in the first processing module to obtain the polynomial model of the residual delay sequence includes: The polynomial model of the residual differential delay rate sequence is integrated to obtain the polynomial model of the residual delay rate sequence, and the polynomial model of the residual delay rate sequence is integrated to obtain the polynomial model of the residual delay sequence.

8. The interferometric processing system suitable for detectors with limited a priori orbit accuracy according to claim 5, characterized in that: The second processing module is specifically used for: The final delay observation is determined based on the a priori delay model and the residual delay, and the final delay rate observation is determined based on the a priori delay model and the residual delay rate.

9. An electronic device, characterized in that: The electronic device includes a processor, the processor is coupled to a memory, at least one computer program is stored in the memory, and the at least one computer program is loaded and executed by the processor so that the electronic device implements the interference processing method as described in any one of claims 1 to 4 under the condition that the prior orbit accuracy of the detector is limited.

10. A computer-readable storage medium, characterized in that: At least one computer program is stored in the computer-readable storage medium, and the at least one computer program is loaded and executed by the processor so that the computer-readable storage medium implements the interference processing method suitable for the detector's prior orbit accuracy being limited as described in any one of claims 1 to 4.

Images