A high order harmonic dorn measurement method, system, device and medium

Abstract

The application discloses a high-order overtone DOR measurement method, system, device and medium, relates to the technical field of information acquisition and processing, and comprises the following steps: performing local correlation processing on a carrier channel signal acquired through a satellite antenna to obtain interference phases in multiple first integral periods, and determining a first phase polynomial of a carrier model based on all the interference phases; determining a high-order overtone model based on the phase polynomial of the carrier model; correcting multiple high-order overtone signals acquired through the satellite antenna based on the high-order overtone model; determining a residual phase corresponding to each high-order overtone signal through the multiple corrected high-order overtone signals; and determining a DOR time delay observation corresponding to each high-order overtone signal through a multi-side tone DOR time delay calculation method based on the residual phase corresponding to each high-order overtone signal. The application can realize acquisition of a DOR measurement time delay observation, and has a wide application prospect in the field of deep space interferometric measurement.

Description

Technical Field

[0001] This invention relates to the field of information acquisition and processing technology, and in particular to a method, system, device and medium for measuring high-order harmonic DOR. Background Technology

[0002] Very Long Baseline Interferometry (VLBI) requires downlink multi-tone beacons for spacecraft. For deep space probes without dedicated DOR beacons, interferometry is performed using the higher harmonics of the carrier wave. Compared to dedicated beacon signals, the signal-to-noise ratio of the higher harmonics of the carrier wave is relatively low, making it impossible to effectively perform correlation processing on the interferometric data. Summary of the Invention

[0003] The technical problem to be solved by this invention is to address the shortcomings of existing technologies, specifically the inability to effectively process interferometric measurement data. Specifically, it provides a method, system, device, and medium for measuring high-order harmonic DOR, as detailed below:

[0004] 1) In a first aspect, the present invention provides a method for measuring the DOR of higher-order harmonics, the specific technical solution of which is as follows:

[0005] S1, perform local correlation processing on the carrier channel signal acquired through the satellite antenna to obtain multiple interference phases under the first integration period, and determine the first phase polynomial of the carrier model based on all interference phases;

[0006] S2, Determine the higher-order harmonic model based on the phase polynomial of the carrier model;

[0007] S3, Based on the aforementioned high-order harmonic model, correct multiple high-order harmonic signals acquired through the satellite antenna;

[0008] S4, using multiple corrected higher-order harmonic signals, determines the residual phase corresponding to each higher-order harmonic signal;

[0009] S5. Based on the residual phase corresponding to each higher-order harmonic signal, the DOR delay observation corresponding to each higher-order harmonic signal is determined by the multi-tone DOR delay calculation method.

[0010] The beneficial effects of the high-order harmonic DOR measurement method provided by this invention are as follows:

[0011] Based on the phase model polynomial of the carrier signal, a phase polynomial model of higher-order harmonics is constructed. By optimizing the integration period, phase correction compensation is performed on higher-order harmonics, and residual phase and residual time delay are calculated. This enables the acquisition of DOR measurement time delay observations, which has broad application prospects in the field of deep space interferometry.

[0012] Based on the above solution, the present invention can be further improved as follows.

[0013] Furthermore, the first phase polynomial of the carrier model determined based on all interferometric phases is specifically as follows:

[0014] The first phase polynomial of the carrier model is determined by the first formula, which is:

[0015] k = polyfit(T,φ,N);

[0016] Where T is the integral period of the carrier channel signal phase; For time vectors, For the phase vector of the carrier model, This represents the carrier signal phase in the m-th integration period, where M is the number of integration periods within the continuous tracking arc, and k = [k N ,k N-1 [,...,k0] is the polynomial coefficient vector, k (n) n = 1, 2, ..., N, k (n) The coefficients of the nth order polynomial are represented by , and polyfit(T,φ,N) represents the least squares fitting calculation of T,φ,N, where N is the order of the polynomial and M>N.

[0017] Furthermore, the determination of the higher-order harmonic model based on the phase polynomial of the carrier model specifically involves:

[0018] The phase polynomial of the higher-order harmonic model is determined by the second formula, which is:

[0019]

[0020] Among them, f c Design the nominal frequency for the carrier wave, f sub Design frequencies for harmonics, This represents a high-order harmonic model of order L, where t represents the time variable. N-1 t represents the N-1 power of the time variable. N The time variable is represented by N raised to the power of L, where L represents the local oscillator frequency of the channel. This represents the initial phase of a higher-order harmonic of order L.

[0021] Furthermore, the specific steps for correcting the multiple high-order harmonic signals acquired through the satellite antenna based on the aforementioned high-order harmonic model are as follows:

[0022] The higher-order harmonic signals are corrected using a third formula, which is:

[0023]

[0024] Where j is the imaginary unit, s sub-L (t) represents the original high-order harmonic signal, s sub-Lcor (t) represents the higher-order harmonic signal after phase model correction.

[0025] Furthermore, determining the residual phase corresponding to each higher-order harmonic signal specifically involves:

[0026] The residual phase is calculated using a fourth formula, which is:

[0027]

[0028] Where J represents the second integration period T sub-L The number of data sampling points contained therein m ind J is the modulation coefficient for higher-order harmonics. L (m ind ) is the independent variable m ind The L-order Bessel function.

[0029] 2) In a second aspect, the present invention also provides a high-order harmonic DOR measurement system, the specific technical solution of which is as follows:

[0030] The acquisition module is used to: perform local correlation processing on the carrier channel signal acquired by the satellite antenna to obtain multiple interference phases under the first integration period, and determine the first phase polynomial of the carrier model based on all interference phases;

[0031] The first determining module is used to: determine the higher-order harmonic model based on the phase polynomial of the carrier model;

[0032] The correction module is used to correct multiple high-order harmonic signals acquired by the satellite antenna based on the high-order harmonic model.

[0033] The second determining module is used to: determine the residual phase corresponding to each higher-order harmonic signal through multiple corrected higher-order harmonic signals;

[0034] The measurement module is used to determine the DOR delay observation corresponding to each higher-order harmonic signal based on the residual phase corresponding to each higher-order harmonic signal and through the multi-tone DOR delay calculation method.

[0035] Based on the above solution, the present invention can be further improved as follows.

[0036] Furthermore, the first phase polynomial of the carrier model determined based on all interferometric phases is specifically as follows:

[0037] The first phase polynomial of the carrier model is determined by the first formula, which is:

[0038] k = polyfit(T,φ,N);

[0039] Where T is the integral period of the carrier channel signal phase; For time vectors, For the phase vector of the carrier model, This represents the carrier signal phase in the m-th integration period, where M is the number of integration periods within the continuous tracking arc, and k = [k N ,k N-1 [,...,k0] is the polynomial coefficient vector, k (n) n = 1, 2, ..., N, k (n) The coefficients of the nth order polynomial are represented by , and polyfit(T,φ,N) represents the least squares fitting calculation of T,φ,N, where N is the order of the polynomial and M>N.

[0040] Furthermore, the determination of the higher-order harmonic model based on the phase polynomial of the carrier model specifically involves:

[0041] The phase polynomial of the higher-order harmonic model is determined by the second formula, which is:

[0042]

[0043] Among them, f c Design the nominal frequency for the carrier wave, f sub Design frequencies for harmonics, This represents a high-order harmonic model of order L, where t represents the time variable. N-1 t represents the N-1 power of the time variable. N The time variable is represented by N raised to the power of L, where L represents the local oscillator frequency of the channel. This represents the initial phase of a higher-order harmonic of order L.

[0044] Furthermore, the specific steps for correcting the multiple high-order harmonic signals acquired through the satellite antenna based on the aforementioned high-order harmonic model are as follows:

[0045] The higher-order harmonic signals are corrected using a third formula, which is:

[0046]

[0047] Where j is the imaginary unit, s sub-L (t) represents the original high-order harmonic signal, s sub-Lcor (t) represents the higher-order harmonic signal after phase model correction.

[0048] Furthermore, determining the residual phase corresponding to each higher-order harmonic signal specifically involves:

[0049] The residual phase is calculated using a fourth formula, which is:

[0050]

[0051] Where J represents the second integration period T sub-L The number of data sampling points contained therein m ind J is the modulation coefficient for higher-order harmonics. L (m ind ) is the independent variable m ind The L-order Bessel function.

[0052] 3) In a third aspect, the present invention also provides a computer device, the computer device including a processor coupled to a memory, the memory storing at least one computer program, the at least one computer program being loaded and executed by the processor to enable the computer device to implement any of the above methods.

[0053] 4) In a fourth aspect, the present invention also provides a computer-readable storage medium storing at least one computer program, which is loaded and executed by a processor to enable a computer to implement any of the above methods.

[0054] It should be noted that the beneficial effects of the technical solutions of the second to fourth aspects of the present invention and their corresponding possible implementations can be found in the above description of the technical effects of the first aspect and its corresponding possible implementations, and will not be repeated here. Attached Figure Description

[0055] Other features, objects, and advantages of the invention will become more apparent from the following detailed description of non-limiting embodiments with reference to the accompanying drawings:

[0056] Figure 1 This is a flowchart illustrating a method for measuring high-order harmonic DOR according to an embodiment of the present invention.

[0057] Figure 2 This is a structural framework diagram of a high-order harmonic DOR measurement system according to an embodiment of the present invention;

[0058] Figure 3 This is a schematic diagram of a computer device structure for a high-order harmonic DOR measurement method according to an embodiment of the present invention;

[0059] Figure 4This is one of the schematic diagrams of the carrier and harmonic spectrum of the observation target provided in the embodiments of the present invention;

[0060] Figure 5 This is the second schematic diagram of the carrier and harmonic spectrum of the observation target provided in the embodiments of the present invention;

[0061] Figure 6 This is the third schematic diagram of the carrier and harmonic spectrum of the observation target provided in the embodiments of the present invention;

[0062] Figure 7 This is the fourth schematic diagram of the carrier and harmonic spectrum of the observation target provided in the embodiments of the present invention;

[0063] Figure 8 This is a flowchart provided for an embodiment of the present invention;

[0064] Figure 9 A schematic diagram of the carrier signal phase provided in an embodiment of the present invention;

[0065] Figure 10 This is a schematic diagram of the phase of a high-order harmonic signal provided in an embodiment of the present invention. Detailed Implementation

[0066] To make the objectives, technical solutions, and advantages of the present invention clearer, the embodiments of the present invention will be described in further detail below with reference to the accompanying drawings.

[0067] like Figure 1 As shown, an embodiment of the present invention provides a method for measuring the DOR of higher-order harmonics, comprising the following steps:

[0068] S1, perform local correlation processing on the carrier channel signal acquired through the satellite antenna to obtain multiple interference phases under the first integration period, and determine the first phase polynomial of the carrier model based on all interference phases;

[0069] S2, Determine the higher-order harmonic model based on the phase polynomial of the carrier model;

[0070] S3, Based on the aforementioned high-order harmonic model, correct multiple high-order harmonic signals acquired through the satellite antenna;

[0071] S4, using multiple corrected higher-order harmonic signals, determines the residual phase corresponding to each higher-order harmonic signal;

[0072] S5. Based on the residual phase corresponding to each higher-order harmonic signal, the DOR delay observation corresponding to each higher-order harmonic signal is determined by the multi-tone DOR delay calculation method.

[0073] The beneficial effects of the high-order harmonic DOR measurement method provided by this invention are as follows:

[0074] Based on the phase model polynomial of the carrier signal, a phase polynomial model of higher-order harmonics is constructed. By optimizing the integration period, phase correction compensation is performed on higher-order harmonics, and residual phase and residual time delay are calculated. This enables the acquisition of DOR measurement time delay observations, which has broad application prospects in the field of deep space interferometry.

[0075] like Figure 8 As shown, S1 performs local correlation processing on the carrier channel signal acquired through the satellite antenna to obtain multiple interferometric phases under the first integration period, and determines the first phase polynomial of the carrier model based on all interferometric phases. Wherein:

[0076] Carrier channel signal refers to the signal transmitted by the satellite antenna, which includes both carrier signal and higher-order harmonics. The carrier signal mainly supports satellite telemetry and control communication. The carrier signal is usually in the form of a sine wave, which can be a pure single carrier or modulated telemetry information.

[0077] Local correlation processing refers to the process where various signals transmitted by a satellite antenna are received and recorded by a ground-based VLBI antenna. The basic idea of ​​the local correlation processing algorithm is to construct a reference signal at the data processing end using auxiliary information such as prior ephemeris, perform cross-correlation calculations on the reference signal and the received signal, and extract the phase information of the target signal.

[0078] It should be noted that satellites transmit signals via antennas on the satellite, and the signals are received, collected, and recorded by ground-based VLBI antennas before being processed.

[0079] The interference phases under multiple first integration periods are:

[0080]

[0081] Where M is the number of integration cycles within the continuous tracking arc segment.

[0082] The first phase polynomial of the carrier model, based on all interferometric phases, is determined by the following equation:

[0083] k = polyfit(T,φ,N);

[0084] Where T is the integral period of the carrier channel signal phase; For time vectors, For the phase vector of the carrier model, This represents the carrier signal phase in the m-th integration period, where M is the number of integration periods within the continuous tracking arc, and k = [k N ,k N-1 [,...,k0] is the polynomial coefficient vector, k (n) n = 1, 2, ..., N, k (n)The coefficients of the nth order polynomial are represented by , and polyfit(T,φ,N) represents the least squares fitting calculation of T,φ,N, where N is the order of the polynomial and M>N.

[0085] S2, determine the higher-order harmonic model based on the phase polynomial of the carrier model. Wherein:

[0086] The higher-order harmonic model is determined by the following formula:

[0087]

[0088] Among them, f c Indicates the nominal design frequency of the carrier wave, in Hz, f. sub Indicates the harmonic design frequency. This represents a high-order harmonic model of order L, where t represents the time variable. N-1 t represents the N-1 power of the time variable. N The time variable is represented by N raised to the power of L, where L represents the local oscillator frequency of the channel. This represents the initial phase of a higher-order harmonic of order L.

[0089] S3, Based on the aforementioned high-order harmonic model, correct multiple high-order harmonic signals acquired through the satellite antenna. Wherein:

[0090] The correction is made using the following formula:

[0091]

[0092] Where j is the imaginary unit, s sub-L (t) represents the original high-order harmonic signal, s sub-Lcor (t) represents the signal after phase model correction.

[0093] S4 determines the residual phase corresponding to each higher-order harmonic signal using multiple corrected higher-order harmonic signals. Among them,

[0094] The residual phase is calculated using the following formula:

[0095]

[0096] Where J represents the second integration period T sub-L The number of data sampling points contained within, and the higher-order harmonic integral period T. sub-L for:

[0097]

[0098] Where m ind J is the modulation coefficient for higher-order harmonics, which affects the energy distribution of higher-order harmonics. L (mind ) is the independent variable m ind The L-order Bessel function.

[0099] S5. Based on the residual phase corresponding to each higher-order harmonic signal, the DOR delay observation corresponding to each higher-order harmonic signal is determined by the multi-tone DOR delay calculation method.

[0100] The multi-tone DOR delay calculation method refers to:

[0101]

[0102] Where f sub-L1 f sub-L2 For the frequency of higher-order harmonics, That is M represents the phase ambiguity, which is obtained based on the satellite's prior orbital information.

[0103] Example 1: An observation verification was carried out using a certain measuring station. The four channel signals collected by the station were: carrier wave, +2nd harmonic, -14th harmonic, and +20th harmonic. It can be observed from the image that the signal-to-noise ratio of the higher-order harmonic signals is extremely low, and the signals are submerged in noise and cannot be effectively identified. Figure 4-7 The image shows the signal spectrum for each channel. The carrier phase results are as follows: Figure 9 As shown, the phase results for higher-order harmonics are as follows: Figure 10 As shown.

[0104] Furthermore, the first phase polynomial of the carrier model determined based on all interferometric phases is specifically as follows:

[0105] The first phase polynomial of the carrier model is determined by the first formula, which is:

[0106] k = polyfit(T,φ,N);

[0107] Where T is the integral period of the carrier channel signal phase; For time vectors, For the phase vector of the carrier model, This represents the carrier signal phase in the m-th integration period, where M is the number of integration periods within the continuous tracking arc, and k = [k N ,k N-1 [,...,k0] is the polynomial coefficient vector, k (n) n = 1, 2, ..., N, k (n) The coefficients of the nth order polynomial are represented by , and polyfit(T,φ,N) represents the least squares fitting calculation of T,φ,N, where N is the order of the polynomial and M>N.

[0108] Furthermore, the determination of the higher-order harmonic model based on the phase polynomial of the carrier model specifically involves:

[0109] The phase polynomial of the higher-order harmonic model is determined by the second formula, which is:

[0110]

[0111] Among them, f c Design the nominal frequency for the carrier wave, f sub Design frequencies for harmonics, This represents a high-order harmonic model of order L, where t represents the time variable. N-1 t represents the N-1 power of the time variable. N The time variable is represented by N raised to the power of L, where L represents the local oscillator frequency of the channel. This represents the initial phase of a higher-order harmonic of order L.

[0112] Furthermore, the specific steps for correcting the multiple high-order harmonic signals acquired via the satellite antenna based on the aforementioned high-order harmonic model are as follows:

[0113] The higher-order harmonic signals are corrected using a third formula, which is:

[0114]

[0115] Where j is the imaginary unit, s sub-L (t) represents the original higher-order harmonic signal, s sub-Lcor (t) represents the higher-order harmonic signal after phase model correction.

[0116] Furthermore, determining the residual phase corresponding to each higher-order harmonic signal specifically involves:

[0117] The residual phase is calculated using a fourth formula, which is:

[0118]

[0119] Where J represents the second integration period T sub-L The number of data sampling points contained therein m ind J is the modulation coefficient for higher-order harmonics. L (m ind ) is the independent variable m ind The L-order Bessel function.

[0120] like Figure 2 As shown, the present invention also provides a high-order harmonic DOR measurement system 200, the specific technical solution of which is as follows:

[0121] The acquisition module 210 is used to: perform local correlation processing on the carrier channel signal acquired by the satellite antenna to obtain multiple interference phases under the first integration period, and determine the first phase polynomial of the carrier model based on all interference phases;

[0122] The first determining module 220 is used to: determine the higher-order harmonic model based on the phase polynomial of the carrier model;

[0123] The correction module 230 is used to: correct multiple high-order harmonic signals acquired by the satellite antenna based on the high-order harmonic model;

[0124] The second determining module 240 is used to: determine the residual phase corresponding to each higher-order harmonic signal through multiple corrected higher-order harmonic signals;

[0125] The measurement module 250 is used to: determine the DOR delay observation corresponding to each higher-order harmonic signal based on the residual phase corresponding to each higher-order harmonic signal, using a multi-tone DOR delay calculation method.

[0126] Furthermore, the first phase polynomial of the carrier model determined based on all interferometric phases is specifically as follows:

[0127] The first phase polynomial of the carrier model is determined by the first formula, which is:

[0128] k = polyfit(T,φ,N);

[0129] Where T is the integral period of the carrier channel signal phase; For time vectors, For the phase vector of the carrier model, This represents the carrier signal phase in the m-th integration period, where M is the number of integration periods within the continuous tracking arc, and k = [k N ,k N-1 [,...,k0] is the polynomial coefficient vector, k (n) (,n=1,2,...N,k (n) The coefficients of the nth order polynomial are represented by , and polyfit(T,φ,N) represents the least squares fitting calculation of T,φ,N, where N is the order of the polynomial and M>N.

[0130] Furthermore, the determination of the higher-order harmonic model based on the phase polynomial of the carrier model specifically involves:

[0131] The phase polynomial of the higher-order harmonic model is determined by the second formula, which is:

[0132]

[0133] Among them, f c Design the nominal frequency for the carrier wave, f sub Design frequencies for harmonics, This represents a high-order harmonic model of order L, where t represents the time variable. N-1 t represents the N-1 power of the time variable. N The time variable is represented by N raised to the power of L, where L represents the local oscillator frequency of the channel. This represents the initial phase of a higher-order harmonic of order L.

[0134] Furthermore, the specific steps for correcting the multiple high-order harmonic signals acquired via the satellite antenna based on the aforementioned high-order harmonic model are as follows:

[0135] The higher-order harmonic signals are corrected using a third formula, which is:

[0136]

[0137] Where j is the imaginary unit, s sub-L (t) represents the original higher-order harmonic signal, s sub-Lcor (t) represents the higher-order harmonic signal after phase model correction.

[0138] Furthermore, determining the residual phase corresponding to each higher-order harmonic signal specifically involves:

[0139] The residual phase is calculated using a fourth formula, which is:

[0140]

[0141] Where J represents the second integration period T sub-L The number of data sampling points contained therein m ind J is the modulation coefficient for higher-order harmonics. L (m ind ) is the independent variable m ind The L-order Bessel function.

[0142] In the above embodiments, although the steps are numbered S1, S2, etc., they are only specific embodiments given by the present invention. Those skilled in the art can adjust the execution order of S1, S2, etc. according to the actual situation, which is also within the protection scope of the present invention. It can be understood that in some embodiments, some or all of the above embodiments may be included.

[0143] It should be noted that the beneficial effects of the high-order harmonic DOR measurement system provided in the above embodiments are the same as those of the high-order harmonic DOR measurement method described above, and will not be repeated here. Furthermore, the system provided in the above embodiments is only illustrated by the division of the above functional modules. In practical applications, the above functions can be assigned to different functional modules as needed, that is, the system can be divided into different functional modules according to the actual situation to complete all or part of the functions described above. In addition, the system and method embodiments provided in the above embodiments belong to the same concept, and their specific implementation process is detailed in the method embodiments, and will not be repeated here.

[0144] like Figure 3 As shown, an embodiment of the present invention provides a computer device 300, which includes a processor 320 coupled to a memory 310. The memory 310 stores at least one computer program 330, which is loaded and executed by the processor 320 to enable the computer device 300 to implement any of the above-described methods. Specifically:

[0145] The computer device 300 can vary considerably due to differences in configuration or performance. It may include one or more processors 320 (Central Processing Units, CPUs) and one or more memories 310. The one or more memories 310 store at least one computer program 330, which is loaded and executed by the one or more processors 320 to enable the computer device 300 to implement the high-order harmonic DOR measurement method provided in the above embodiments. Of course, the computer device 300 may also have wired or wireless network interfaces, a keyboard, and input / output interfaces for input and output. The computer device 300 may also include other components for implementing device functions, which will not be elaborated here.

[0146] An embodiment of the present invention provides a computer-readable storage medium storing at least one computer program, which is loaded and executed by a processor to enable a computer to implement any of the above-described methods.

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

[0148] In an exemplary embodiment, a computer program product or computer program is also provided, the computer program product or computer program including computer instructions stored in a computer-readable storage medium. A processor of a computer device reads the computer instructions from the computer-readable storage medium, and the processor executes the computer instructions, causing the computer device to perform any of the methods described above.

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

[0150] Those skilled in the art will recognize that this invention can be implemented as a system, method, or computer program product. Therefore, this disclosure can be specifically implemented in the following forms: it can be entirely hardware, entirely software (including firmware, resident software, microcode, etc.), or a combination of hardware and software, generally referred to herein as a "circuit," "module," or "system." Furthermore, in some embodiments, this invention can also be implemented as a computer program product in one or more computer-readable media containing computer-readable program code.

[0151] Any combination of one or more computer-readable media may be used. A computer-readable medium can be a computer-readable signal medium or a computer-readable storage medium. A computer-readable storage medium can be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples (a non-exhaustive list) of computer-readable storage media include: an electrical connection having one or more wires, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage device, magnetic storage device, or any suitable combination thereof. In this document, a computer-readable storage medium can be any tangible medium that contains or stores a program that can be used by or in connection with an instruction execution system, apparatus, or device.

[0152] Although embodiments of the present invention have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present invention. Those skilled in the art can make changes, modifications, substitutions and variations to the above embodiments within the scope of the present invention.

Claims

1. A method for measuring the DOR of higher-order harmonics, characterized in that, include: S1, perform local correlation processing on the carrier channel signal acquired through the satellite antenna to obtain multiple interference phases under the first integration period, and determine the first phase polynomial of the carrier model based on all interference phases; S2, Determine the higher-order harmonic model based on the phase polynomial of the carrier model; S3, Based on the aforementioned high-order harmonic model, correct multiple high-order harmonic signals acquired through the satellite antenna; S4, using multiple corrected higher-order harmonic signals, determines the residual phase corresponding to each higher-order harmonic signal; S5. Based on the residual phase corresponding to each higher-order harmonic signal, the DOR delay observation corresponding to each higher-order harmonic signal is determined by the multi-tone DOR delay calculation method. The first phase polynomial of the carrier model determined based on all interferometric phases is specifically: The first phase polynomial of the carrier model is determined by the first formula, which is: ； in, For time vectors, The integral period of the carrier channel signal phase; For the phase vector of the carrier model, Where M is the number of integration cycles within the continuously tracking arc segment. Let be the order of the polynomial. Let be the polynomial coefficient vector. (T, ) indicates that for T, Perform least squares fitting calculations; M > N. The determination of the higher-order harmonic model based on the phase polynomial of the carrier model specifically involves: The phase polynomial of the higher-order harmonic model is determined by the second formula, which is: ； in, Design the nominal frequency for the carrier wave. Design frequencies for harmonics, This represents a high-order harmonic model of order L, where t represents the time variable. N-1 t represents the N-1 power of the time variable. N The time variable is represented by N raised to the power of L, where L represents the local oscillator frequency of the channel. This represents the initial phase of a higher-order harmonic of order L.

2. The method for measuring high-order harmonic DOR according to claim 1, characterized in that, The specific steps for correcting multiple high-order harmonic signals acquired via satellite antenna based on the aforementioned high-order harmonic model are as follows: The higher-order harmonic signals are corrected using a third formula, which is: ； Where j is the imaginary unit, The original high-order harmonic signal, This refers to the higher-order harmonic signal after phase model correction.

3. The method for measuring high-order harmonic DOR according to claim 2, characterized in that, The determination of the residual phase corresponding to each higher-order harmonic signal specifically involves: The residual phase is calculated using a fourth formula, which is: ； Where J represents the second integration period. The number of data sampling points contained therein ;m ind The modulation coefficients for higher-order harmonics, It is The order Bessel function.

4. A high-order harmonic DOR measurement system, characterized in that, include: The acquisition module is used to: perform local correlation processing on the carrier channel signal acquired by the satellite antenna to obtain multiple interference phases under the first integration period, and determine the first phase polynomial of the carrier model based on all interference phases; The first determining module is used to: determine the higher-order harmonic model based on the phase polynomial of the carrier model; The correction module is used to correct multiple high-order harmonic signals acquired by the satellite antenna based on the high-order harmonic model. The second determining module is used to: determine the residual phase corresponding to each higher-order harmonic signal through multiple corrected higher-order harmonic signals; The measurement module is used to: determine the DOR delay observation corresponding to each higher-order harmonic signal based on the residual phase corresponding to each higher-order harmonic signal, using the multi-tone DOR delay calculation method; The first phase polynomial of the carrier model determined based on all interferometric phases is specifically: The first phase polynomial of the carrier model is determined by the first formula, which is: ； in, For time vectors, The integral period of the carrier channel signal phase; For the phase vector of the carrier model, Where M is the number of integration cycles within the continuously tracking arc segment. Let be the order of the polynomial. Let be the polynomial coefficient vector. (T, ) indicates that for T, Perform least squares fitting calculations; M > N. The determination of the higher-order harmonic model based on the phase polynomial of the carrier model specifically involves: The phase polynomial of the higher-order harmonic model is determined by the second formula, which is: ； in, Design the nominal frequency for the carrier wave. Design frequencies for harmonics, This represents a high-order harmonic model of order L, where t represents the time variable. N-1 t represents the N-1 power of the time variable. N The time variable is represented by N raised to the power of L, where L represents the local oscillator frequency of the channel. This represents the initial phase of a higher-order harmonic of order L.

5. A computer device, characterized in that, The computer device includes a processor coupled to a memory storing at least one computer program, which is loaded and executed by the processor to enable the computer device to perform the method as claimed in any one of claims 1 to 3.

6. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores at least one computer program, which is loaded and executed by a processor to enable the computer to perform the method as claimed in any one of claims 1 to 3.

Images