Interferometric processing method, system, device and medium based on carrier phase reference

Abstract

The application discloses an interference processing method, system, device and medium based on a carrier phase reference, and the method comprises the following steps: constructing a reference phase of each side tone signal based on the phase of a main carrier signal and the design parameter of a detector beacon signal; constructing a reference model of each side tone signal based on the reference phase of each side tone signal, and calculating the cross-correlation residual phase of each side tone signal according to the reference model of each side tone signal; and calculating a group delay observation based on the cross-correlation residual phase of each side tone signal. The application can overcome the problem of precision decline of a traditional algorithm in a low signal-to-noise ratio and a frequency fluctuation scene, improve the accuracy of interference time delay observation, and meet the demand of high-precision measurement of a detector orbit.

Description

Technical Field

[0001] This invention relates to the field of information acquisition and processing technology, and in particular to an interferometric processing method, system, device and medium based on carrier phase reference. Background Technology

[0002] Very Long Baseline Interferometry (VLBI) is a key technique for determining detector orbits, and its accuracy depends on local correlation processing algorithms. However, in practical applications, when the target signal has a weak signal-to-noise ratio or the detector's nominal downlink frequency is not ideally constant, the processing accuracy of traditional local correlation processing algorithms decreases significantly over long integration periods, making it difficult to meet the requirements for high-precision measurement of interferometric time delay observations.

[0003] Therefore, there is an urgent need to provide a technical solution to address the above problems. Summary of the Invention

[0004] To address the aforementioned technical problems, this invention provides an interferometric processing method, system, device, and medium based on carrier phase reference.

[0005] In a first aspect, the present invention provides an interferometric processing method based on carrier phase reference, the technical solution of which is as follows: Based on the phase of the main carrier signal and the design parameters of the detector beacon signal, a reference phase for each sidetone signal is constructed; Based on the reference phase of each sidetone signal, a reference model for each sidetone signal is constructed, and the cross-correlation residual phase of each sidetone signal is calculated according to the reference model of each sidetone signal. The group delay observable is calculated based on the cross-correlation residual phase of each sidetone signal.

[0006] The beneficial effects of the interferometric processing method based on carrier phase reference of the present invention are as follows: The method of this invention can overcome the problem of accuracy degradation of traditional algorithms in low signal-to-noise ratio and frequency fluctuation scenarios, improve the accuracy of interferometric time delay observation, and meet the requirements of high-precision determination of detector orbit.

[0007] Based on the above scheme, the interferometric processing method based on carrier phase reference of the present invention can be further improved as follows.

[0008] In one alternative approach, it also includes: Using a local correlation processing algorithm, the residual phase of the main carrier signal and the theoretical model are calculated; The phase of the main carrier signal is constructed based on the residual phase of the main carrier signal and the theoretical model.

[0009] In one alternative approach, the design parameters include: the frequency ratio of each sidetone signal to the main carrier signal; the step of constructing a reference phase for each sidetone signal based on the phase of the main carrier signal and the design parameters of the detector beacon signal includes: Based on the frequency ratio of any side tone signal to the main carrier signal, and in combination with the phase of the main carrier signal, the reference phase of the side tone signal is determined. Repeat the step of determining the reference phase of the side tone signal based on the frequency ratio of any side tone signal to the main carrier signal and in combination with the phase of the main carrier signal, until the reference phase of each side tone signal is obtained.

[0010] In one alternative approach, the step of constructing a reference model for each sidetone signal based on the reference phase of each sidetone signal includes: The reference phase of any sidetone signal is used as the input parameter of the target exponential function, and a reference model of the sidetone signal is constructed through the target exponential function. Repeat the steps of using the reference phase of any sidetone signal as the input parameter of the target exponential function and constructing a reference model of the sidetone signal through the target exponential function until a reference model for each sidetone signal is constructed.

[0011] In one alternative approach, the step of calculating the cross-correlation residual phase of each sidetone signal based on a reference model for each sidetone signal includes: Based on any side tone signal and its reference model, calculate the cross-correlation residual phase of the side tone signal in each integration period; Repeat the steps of calculating the cross-correlation residual phase of the side tone signal in each integration period based on any side tone signal and the reference model of the side tone signal, until the cross-correlation residual phase of each side tone signal in each integration period is obtained.

[0012] In one alternative approach, the step of calculating the group delay observable based on the cross-correlation residual phase of each sidetone signal includes: Based on the phase of the main carrier signal and the cross-correlation residual phase of each side tone signal at the first station in any integration period, and the phase of the main carrier signal and the cross-correlation residual phase of each side tone signal at the second station in any integration period, calculate the group delay observation of the first station and the second station in any integration period. Repeat the steps of calculating the group delay observation of the first station and the second station in any integration period based on the phase of the main carrier signal and the cross-correlation residual phase of each side tone signal of the first station in any integration period, and the phase of the main carrier signal and the cross-correlation residual phase of each side tone signal of the second station in any integration period, until the group delay observation of the first station and the second station in each integration period is calculated.

[0013] Secondly, the present invention provides an interferometric processing system based on carrier phase reference, the technical solution of which is as follows: It includes: an acquisition module, a processing module, and a calculation module; The acquisition module is used to: construct a reference phase for each sidetone signal based on the phase of the main carrier signal and the design parameters of the detector beacon signal; The processing module is used to: construct a reference model for each sidetone signal based on the reference phase of each sidetone signal, and calculate the cross-correlation residual phase of each sidetone signal according to the reference model of each sidetone signal; The calculation module is used to calculate the group delay observation based on the cross-correlation residual phase of each sidetone signal.

[0014] The beneficial effects of the interferometric processing system based on carrier phase reference of the present invention are as follows: The system of this invention can overcome the problem of accuracy degradation of traditional algorithms in low signal-to-noise ratio and frequency fluctuation scenarios, improve the accuracy of interferometric time delay observation, and meet the requirements of high-precision determination of detector orbit.

[0015] Based on the above scheme, the interferometric processing system based on carrier phase reference of the present invention can be further improved as follows.

[0016] In an alternative embodiment, it further includes: a building module; the building module is used for: Using a local correlation processing algorithm, the residual phase of the main carrier signal and the theoretical model are calculated; The phase of the main carrier signal is constructed based on the residual phase of the main carrier signal and the theoretical model.

[0017] Thirdly, the technical solution of an electronic device according to the present invention is as follows: It includes a memory, a processor, and a program stored in the memory and running on the processor, wherein the processor executes the program to implement the steps of the carrier phase reference-based interferometric processing method of the present invention.

[0018] Fourthly, the technical solution of a computer-readable storage medium provided by the present invention is as follows: The computer-readable storage medium stores instructions that, when read, cause the computer-readable storage medium to perform the steps of the carrier phase reference-based interferometric processing method of the present invention.

[0019] The above description is merely an overview of the technical solution of the present invention. In order to better understand the technical means of the present invention and to implement it in accordance with the contents of the specification, and in order to make the above and other objects, features and advantages of the present invention more apparent and understandable, specific embodiments of the present invention are described below. Attached Figure Description

[0020] The accompanying drawings are for illustrative purposes only and are not intended to limit the invention. Furthermore, the same reference numerals denote the same parts throughout the drawings. In the drawings: Figure 1 This is a schematic flowchart of an embodiment of the interferometric processing method based on carrier phase reference of the present invention; Figure 2 A schematic diagram of the spectrum of the main carrier signal; Figure 3 A schematic diagram of the residual phase of the main carrier signal; Figure 4 This is a schematic diagram of the residual phase of the sidetone signal; Figure 5 This is a schematic diagram of an embodiment of an interferometric processing system based on carrier phase reference according to the present invention; Figure 6 This is a schematic diagram of an embodiment of an electronic device according to the present invention. Detailed Implementation

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

[0022] Figure 1This diagram illustrates a flowchart of an embodiment of a carrier phase reference-based interferometric processing method provided by the present invention. This method can be executed by an electronic device such as a terminal device or a server. The terminal device can be any fixed or mobile terminal, such as a user equipment (UE), mobile device, user terminal, terminal, cellular phone, cordless phone, personal digital assistant (PDA), handheld device, computing device, in-vehicle device, or wearable device. The server can be a single server or a server cluster consisting of multiple servers. Any electronic device can implement the carrier phase reference-based interferometric processing method by having its processor call computer-readable instructions stored in its memory. Figure 1 As shown, it includes the following steps: S1. Based on the phase of the main carrier signal and the design parameters of the detector beacon signal, construct the reference phase of each sidetone signal.

[0023] The primary carrier signal refers to the radio frequency signal that serves as the frequency reference in the detector's downlink. The phase of the primary carrier signal refers to the angular position of the primary carrier signal waveform relative to the time reference at a specific moment. The detector beacon signal refers to the composite signal emitted by the detector, which includes the primary carrier and multiple side-tone signals. Design parameters refer to the pre-defined frequency configuration parameters for the detector beacon signal; for example, the frequency ratio between the side-tone signals and the primary carrier signal specified in the design document.

[0024] The sidetone signal refers to a low-frequency sinusoidal signal modulated on the main carrier. The reference phase of the sidetone signal refers to the theoretical phase of the sidetone signal calculated based on the ratio of the main carrier phase to the frequency.

[0025] S2. Based on the reference phase of each sidetone signal, construct a reference model for each sidetone signal, and calculate the cross-correlation residual phase of each sidetone signal according to the reference model of each sidetone signal.

[0026] The reference model for the sidetone signal refers to the mathematical model of the sidetone signal constructed based on the reference phase. The cross-correlation residual phase of the sidetone signal refers to the phase difference obtained after cross-correlation processing between the actual sidetone signal and the reference model of the sidetone signal.

[0027] S3. Calculate the group delay observable based on the cross-correlation residual phase of each sidetone signal.

[0028] Among them, the group time delay observation refers to the observation value that characterizes the time difference of signal propagation between different stations.

[0029] The technical solution of this embodiment can overcome the problem of accuracy degradation of traditional algorithms in low signal-to-noise ratio and frequency fluctuation scenarios, improve the accuracy of interferometric time delay observation, and meet the high-precision measurement requirements of detector orbit.

[0030] In one alternative approach, it also includes: The residual phase of the main carrier signal and the theoretical model are calculated using a local correlation processing algorithm.

[0031] The local correlation processing algorithm is DOR (Differential One-way Ranging) local correlation processing method. Residual phase refers to the difference between the actual signal phase and the theoretical model phase. The theoretical model refers to the ideal signal mathematical model constructed based on the nominal frequency.

[0032] The theoretical model of the main carrier signal is expressed as: ;in, The theoretical model representing the main carrier signal. This represents the natural exponential function. This represents the nominal constant frequency assumed by the local related processing procedures. Represents a time variable. This represents a function for multinomial computation. Represents the coefficients of the constructed polynomial.

[0033] The residual phase representation sequence of the main carrier signal is as follows: ; This represents the residual phase sequence of the main carrier signal. Indicates the number of integration cycles; These represent the residual phase values ​​of the main carrier signal at the 1st to Nth integration period points, respectively.

[0034] The phase of the main carrier signal is constructed based on the residual phase of the main carrier signal and the theoretical model.

[0035] The spectrum of the main carrier signal is as follows: Figure 2 As shown, the residual phase of the main carrier signal is as follows: Figure 3 As shown.

[0036] Specifically, the phase of the main carrier signal is calculated using the following formula, where interpolation is used to obtain the phase at the sampling point from the phase at the integration period point. The formula for calculating the phase of the main carrier signal is expressed as: ;in, Indicates the phase of the main carrier signal; This represents the interpolation function.

[0037] It should be noted that the calculation formula for the radio frequency domain model of the detector's main carrier signal is as follows: ;in, This indicates the received detector main carrier signal. To represent a complex number, This represents the nominal frequency of the main carrier signal. This frequency is not constant in the long integration period calculation under weak signal-to-noise ratio conditions, and the unit is Hz. Indicates the signal propagation delay; This represents Gaussian noise. and Related, This refers to the actual received primary carrier signal from the detector, while It is a theoretical model of the main carrier signal constructed based on the nominal frequency and a polynomial model.

[0038] In the above-mentioned optional methods, the residual phase of the main carrier signal and the theoretical model are further accurately calculated by local correlation processing algorithms. The phase of the main carrier signal is efficiently constructed by combining the residual phase and the theoretical model, thereby improving the accuracy and reliability of signal processing.

[0039] In one alternative approach, the design parameters include: the frequency ratio of each sidetone signal to the main carrier signal; the step of constructing a reference phase for each sidetone signal based on the phase of the main carrier signal and the design parameters of the detector beacon signal includes: Based on the frequency ratio of any side tone signal to the main carrier signal, and in conjunction with the phase of the main carrier signal, the reference phase of the side tone signal is determined.

[0040] The frequency ratio refers to the ratio of the sidetone signal frequency to the main carrier signal frequency. Different sidetone signals have different frequency ratios to the main carrier signal.

[0041] Specifically, the sidetone signal and the main carrier signal share the same frequency origin, and the reference phase is calculated using the frequency ratio determined by the design parameters. The formula for calculating the reference phase of the sidetone signal is as follows: ;in, Indicates the reference phase of the sidetone signal; This represents the frequency ratio of the sidetone signal to the main carrier signal; This indicates the phase of the main carrier signal. The frequency ratio of different sidetone signals. The value of varies and is determined based on the ratio between the frequency of the corresponding side tone signal and the frequency of the main carrier signal.

[0042] Repeat the step of determining the reference phase of the side tone signal based on the frequency ratio of any side tone signal to the main carrier signal and in combination with the phase of the main carrier signal, until the reference phase of each side tone signal is obtained.

[0043] In the above-mentioned optional methods, the reference phase of the side tone signal is further determined one by one based on the frequency ratio of the side tone signal to the main carrier signal in the design parameters, combined with the main carrier phase, to ensure the accuracy and relevance of the reference phase construction and enhance the adaptability of signal processing.

[0044] In one alternative approach, the step of constructing a reference model for each sidetone signal based on the reference phase of each sidetone signal includes: The reference phase of any sidetone signal is used as the input parameter of the target exponential function, and a reference model of the sidetone signal is constructed through the target exponential function.

[0045] The target exponential function refers to the complex exponential function used to construct the signal reference model.

[0046] The reference model for the sidetone signal is constructed using the following formula: ;in, A reference model representing the sidetone signal; To represent a complex number; This indicates the reference phase of the sidetone signal.

[0047] Repeat the steps of using the reference phase of any sidetone signal as the input parameter of the target exponential function and constructing a reference model of the sidetone signal through the target exponential function until a reference model for each sidetone signal is constructed.

[0048] In the above-mentioned optional methods, the reference phase of the sidetone signal is further input into the target exponential function to construct reference models that fit the signal characteristics one by one, providing a high-precision reference for subsequent signal analysis and improving the scientificity and effectiveness of signal processing.

[0049] In one alternative approach, the step of calculating the cross-correlation residual phase of each sidetone signal based on a reference model for each sidetone signal includes: Based on any side tone signal and its reference model, the cross-correlation residual phase of the side tone signal in each integration period is calculated.

[0050] The integration period refers to the fixed time length used in signal processing. The residual phase of the sidetone signal is as follows: Figure 4 As shown.

[0051] The formula for calculating the cross-correlation residual phase is: ;in, The symbol represents the residual phase of the cross-correlation of the m-th sidetone signal at the i-th integration period. This represents the integration operation within the time interval of the i-th integration period; This represents the m-th side note signal. This indicates the conjugate operation; The reference model representing the m-th sidetone signal; symbol Represents the integral variable; , , This represents the total number of integration cycles. This indicates the total number of sidetone signals.

[0052] Repeat the steps of calculating the cross-correlation residual phase of the side tone signal in each integration period based on any side tone signal and the reference model of the side tone signal, until the cross-correlation residual phase of each side tone signal in each integration period is obtained.

[0053] In the above-mentioned optional methods, the cross-correlation residual phase is further calculated periodically for each sidetone signal and its reference model to refine the evaluation of signal deviation and enhance the dynamic adaptability and accuracy of signal processing.

[0054] In one alternative approach, the step of calculating the group delay observable based on the cross-correlation residual phase of each sidetone signal includes: Based on the phase of the main carrier signal and the cross-correlation residual phase of each side tone signal at the first station in any integration period, and the phase of the main carrier signal and the cross-correlation residual phase of each side tone signal at the second station in any integration period, the group delay observation of the first station and the second station in any integration period is calculated.

[0055] In this context, the first station refers to one ground receiving station participating in the interferometry. The second station refers to another ground receiving station participating in the interferometry. For example, ground station A acts as the first station to receive detector signals, and ground station B acts as the second station, forming the interferometric baseline with station A.

[0056] The formula for calculating group delay observations is: Among them, the symbol Represents the group delay observation. Indicates the main delay component. This represents the residual time delay component.

[0057] The formula for calculating the main delay component is: ;in, This indicates the phase of the primary carrier signal at the first station. This indicates the phase of the main carrier signal at the second station; This indicates the nominal frequency of the main carrier signal. The nominal frequencies of the main carrier signals of the first and second stations are the same.

[0058] The formula for calculating the residual time delay component is: ;in, Represents the least squares fitting function; symbol This represents the set of differences in the residual phase of each sidetone signal between the first and second measuring stations. This represents the set of frequencies for each sidetone signal, and 1 indicates the coefficient of the first-order term of the fitted polynomial. Least square fitting is used to obtain the residual time delay, which is taken as the coefficient of the first-order term in the fitting result.

[0059] Repeat the steps of calculating the group delay observation of the first station and the second station in any integration period based on the phase of the main carrier signal and the cross-correlation residual phase of each side tone signal of the first station in any integration period, and the phase of the main carrier signal and the cross-correlation residual phase of each side tone signal of the second station in any integration period, until the group delay observation of the first station and the second station in each integration period is calculated.

[0060] It should be noted that Table 1 compares the performance of the local correlation algorithm and the algorithm after introducing phase reference correction. By introducing phase reference correction, the residual group delay accuracy obtained in this embodiment can be improved by approximately 5%.

[0061] Table 1: Among the above-mentioned optional methods, the main carrier phase and side tone residual phase of the two stations are further integrated, and the group time delay observation is calculated on an integral cycle basis, which effectively improves the observation accuracy of the signal time difference between the stations and meets the requirements of high-precision determination of the detector orbit.

[0062] Figure 5 A schematic diagram of an embodiment of an interferometric processing system 200 based on carrier phase reference provided by the present invention is shown. Figure 5 As shown, the carrier phase reference-based interferometric processing system 200 includes: an acquisition module 201, a processing module 202, and a calculation module 203; The acquisition module 201 is used to: construct a reference phase for each sidetone signal based on the phase of the main carrier signal and the design parameters of the detector beacon signal; The processing module 202 is used to: construct a reference model for each side note signal based on the reference phase of each side note signal, and calculate the cross-correlation residual phase of each side note signal according to the reference model of each side note signal; The calculation module 203 is used to calculate the group delay observation based on the cross-correlation residual phase of each sidetone signal.

[0063] In an alternative embodiment, it further includes: a building module; the building module is used for: Using a local correlation processing algorithm, the residual phase of the main carrier signal and the theoretical model are calculated; The phase of the main carrier signal is constructed based on the residual phase of the main carrier signal and the theoretical model.

[0064] In one alternative approach, the design parameters include: the frequency ratio of each sidetone signal to the main carrier signal; the acquisition module 201 is specifically used for: Based on the frequency ratio of any side tone signal to the main carrier signal, and in combination with the phase of the main carrier signal, the reference phase of the side tone signal is determined. Repeat the step of determining the reference phase of the side tone signal based on the frequency ratio of any side tone signal to the main carrier signal and in combination with the phase of the main carrier signal, until the reference phase of each side tone signal is obtained.

[0065] In an alternative embodiment, the processing module 202 is specifically used for: The reference phase of any sidetone signal is used as the input parameter of the target exponential function, and a reference model of the sidetone signal is constructed through the target exponential function. Repeat the steps of using the reference phase of any sidetone signal as the input parameter of the target exponential function and constructing a reference model of the sidetone signal through the target exponential function until a reference model for each sidetone signal is constructed.

[0066] In an alternative embodiment, the processing module 202 is specifically used for: Based on any side tone signal and its reference model, calculate the cross-correlation residual phase of the side tone signal in each integration period; Repeat the steps of calculating the cross-correlation residual phase of the side tone signal in each integration period based on any side tone signal and the reference model of the side tone signal, until the cross-correlation residual phase of each side tone signal in each integration period is obtained.

[0067] In an alternative embodiment, the computing module 203 is specifically used for: Based on the phase of the main carrier signal and the cross-correlation residual phase of each side tone signal at the first station in any integration period, and the phase of the main carrier signal and the cross-correlation residual phase of each side tone signal at the second station in any integration period, calculate the group delay observation of the first station and the second station in any integration period. Repeat the steps of calculating the group delay observation of the first station and the second station in any integration period based on the phase of the main carrier signal and the cross-correlation residual phase of each side tone signal of the first station in any integration period, and the phase of the main carrier signal and the cross-correlation residual phase of each side tone signal of the second station in any integration period, until the group delay observation of the first station and the second station in each integration period is calculated.

[0068] It should be noted that the beneficial effects of the carrier phase reference-based interferometric processing system 200 provided in the above embodiments are the same as those of the carrier phase reference-based interferometric processing 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.

[0069] The carrier phase reference-based interferometric processing system 200 of the present invention can be a computer program (including program code) running on a computer device. For example, the carrier phase reference-based interferometric processing system 200 of the present invention is an application software that can be used to execute the corresponding steps in the carrier phase reference-based interferometric processing method of the present invention.

[0070] In some embodiments, the carrier phase reference-based interferometric processing system 200 of the present invention can be implemented in a combination of hardware and software. As an example, the carrier phase reference-based interferometric processing system 200 of the present invention can be a processor in the form of a hardware decoding processor, which is programmed to execute the carrier phase reference-based interferometric processing method of the present invention. For example, the processor in the form of a hardware decoding processor can be one or more application specific integrated circuits (ASICs), DSPs, programmable logic devices (PLDs), complex programmable logic devices (CPLDs), field-programmable gate arrays (FPGAs), or other electronic components.

[0071] The modules described in the embodiments of this invention can be implemented in software or hardware. The names of the modules are not, in some cases, limiting the scope of the module itself.

[0072] An electronic device according to an embodiment of the present invention includes a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements any of the above-mentioned interferometric processing methods based on carrier phase reference. That is, an electronic device according to an embodiment of the present invention may include, but is not limited to: a processor and a memory; the memory is used to store the computer program; the processor is used to execute the interferometric processing method based on carrier phase reference shown in any embodiment of the present invention by calling the computer program.

[0073] In one alternative embodiment, an electronic device is provided, such as Figure 6 As shown, Figure 6 The illustrated electronic device 4000 includes a processor 4001 and a memory 4003. The processor 4001 and the memory 4003 are connected, for example, via a bus 4002. Optionally, the electronic device 4000 may further include a transceiver 4004, which can be used for data interaction between the electronic device and other electronic devices, such as sending and / or receiving data. It should be noted that in practical applications, the transceiver 4004 is not limited to one type, and the structure of the electronic device 4000 does not constitute a limitation on the embodiments of the present invention.

[0074] Processor 4001 may be a CPU (Central Processing Unit), a general-purpose processor, a DSP (Digital Signal Processor), an ASIC (Application Specific Integrated Circuit), an FPGA (Field Programmable Gate Array), or other programmable logic devices, transistor logic devices, hardware components, or any combination thereof. It can implement or execute the various exemplary logic blocks, modules, and circuits described in conjunction with the disclosure of this invention. Processor 4001 may also be a combination that implements computational functions, such as including one or more microprocessor combinations, a combination of a DSP and a microprocessor, etc.

[0075] Bus 4002 may include a path for transmitting information between the aforementioned components. Bus 4002 may be a PCI (Peripheral Component Interconnect) bus or an EISA (Extended Industry Standard Architecture) bus, etc. Bus 4002 can be divided into address bus, data bus, control bus, etc. For ease of representation, Figure 6The bus 4002 is represented by only one thick line, but this does not mean that there is only one bus or one type of bus.

[0076] The memory 4003 may be ROM (Read Only Memory) or other types of static storage devices capable of storing static information and instructions, RAM (Random Access Memory) or other types of dynamic storage devices capable of storing information and instructions, or EEPROM (Electrically Erasable Programmable Read Only Memory), CD-ROM (Compact Disc Read Only Memory) or other optical disc storage, optical disc storage (including compressed optical discs, laser discs, optical discs, digital universal optical discs, Blu-ray discs, etc.), magnetic disk storage media or other magnetic storage devices, or any other medium capable of carrying or storing desired program code in the form of instructions or data structures and accessible by a computer, but not limited thereto.

[0077] The memory 4003 stores application code (computer program) for executing the present invention, and its execution is controlled by the processor 4001. The processor 4001 executes the application code stored in the memory 4003 to implement the content shown in the foregoing method embodiments.

[0078] Among them, electronic devices can also be terminal devices. A terminal device can be any terminal device that can install applications and access web pages through applications, including at least one of smartphones, tablets, laptops, desktop computers, smart speakers, smartwatches, smart TVs, and smart in-vehicle devices.

[0079] It should be noted that, Figure 6 The electronic device shown is merely an example and should not be construed as limiting the functionality and scope of use of the embodiments of the present invention.

[0080] An embodiment of the present invention provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements any of the above-described interferometric processing methods based on carrier phase reference.

[0081] 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.

[0082] In an exemplary embodiment, a computer program product or computer program is also provided, which includes computer instructions stored in a computer-readable storage medium. A processor of an electronic device reads the computer instructions from the computer-readable storage medium and executes the computer instructions, causing the electronic device to perform the aforementioned carrier phase reference-based interferometry processing method.

[0083] Computer program code for performing the operations of this invention can be written in one or more programming languages ​​or a combination thereof, including object-oriented programming languages ​​such as Java, Smalltalk, and C++, and conventional procedural programming languages ​​such as C or similar languages. The program code can be executed entirely on the user's computer, partially on the user's computer, as a standalone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving remote computers, the remote computer can be connected to the user's computer via any type of network—including a local area network (LAN) or a wide area network (WAN)—or can be connected to an external computer (e.g., via the Internet using an Internet service provider).

[0084] It should be understood that the flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of methods and computer program products according to various embodiments of the present invention. In this regard, each block in a flowchart or block diagram may represent a module, segment, or portion of code containing one or more executable instructions for implementing the specified logical function. It should also be noted that in some alternative implementations, the functions indicated in the blocks may occur in a different order than those indicated in the drawings. For example, two consecutively indicated blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in the block diagrams and / or flowcharts, and combinations of blocks in the block diagrams and / or flowcharts, may be implemented using a dedicated hardware-based system that performs the specified function or operation, or using a combination of dedicated hardware and computer instructions.

[0085] The computer-readable storage medium provided in this invention can be, but is not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples of a computer-readable storage medium may include, but are not limited to: 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 invention, a computer-readable storage medium can be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, apparatus, or device.

[0086] The aforementioned computer-readable storage medium carries one or more programs, which, when executed by the electronic device, cause the electronic device to perform the method shown in the above embodiments.

[0087] The above description is merely a preferred embodiment of the present invention and an explanation of the technical principles employed. Those skilled in the art should understand that the scope of disclosure in this invention is not limited to technical solutions formed by specific combinations of the above-described technical features, but should also cover other technical solutions formed by arbitrary combinations of the above-described technical features or their equivalents without departing from the above-disclosed concept. For example, technical solutions formed by substituting the above features with (but not limited to) technical features with similar functions disclosed in this invention.

[0088] 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.

[0089] Those skilled in the art will recognize that this invention can be implemented as a system, method, or computer program product. Therefore, this invention 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 contained in one or more computer-readable media, which includes computer-readable program code.

[0090] 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. An interferometric processing method based on carrier phase reference, characterized in that, include: Based on the phase of the main carrier signal and the design parameters of the detector beacon signal, a reference phase for each sidetone signal is constructed; Based on the reference phase of each sidetone signal, a reference model for each sidetone signal is constructed, and the cross-correlation residual phase of each sidetone signal is calculated according to the reference model of each sidetone signal. The group delay observable is calculated based on the cross-correlation residual phase of each sidetone signal; The design parameters include: the frequency ratio of each sidetone signal to the main carrier signal; the steps for constructing a reference phase for each sidetone signal based on the phase of the main carrier signal and the design parameters of the detector beacon signal include: Based on the frequency ratio of any side tone signal to the main carrier signal, and in combination with the phase of the main carrier signal, the reference phase of the side tone signal is determined. Repeat the step of determining the reference phase of the side tone signal based on the frequency ratio of any side tone signal to the main carrier signal and in combination with the phase of the main carrier signal, until the reference phase of each side tone signal is obtained; The formula for calculating group delay observations is: ;in, Represents the group delay observation. Indicates the main delay component. Indicates the residual delay component; The formula for calculating the main delay component is: ;in, This indicates the phase of the primary carrier signal at the first station. This indicates the phase of the main carrier signal at the second station; This indicates the nominal frequency of the main carrier signal. The nominal frequencies of the main carrier signals of the first and second stations are the same. The formula for calculating the residual time delay component is: ;in, This represents the least squares fitting function; This represents the set of differences in the residual phase of each sidetone signal between the first and second measuring stations. This represents the set of frequencies for each sidetone signal, and 1 represents the coefficient of the first-order term of the fitted polynomial. Least square fitting is used to obtain the residual time delay, which is taken as the coefficient of the first-order term of the fitting result.

2. The interferometric processing method based on carrier phase reference according to claim 1, characterized in that, Also includes: Using a local correlation processing algorithm, the residual phase of the main carrier signal and the theoretical model are calculated; The phase of the main carrier signal is constructed based on the residual phase of the main carrier signal and the theoretical model.

3. The interferometric processing method based on carrier phase reference according to claim 1, characterized in that, The step of constructing a reference model for each sidetone signal based on the reference phase of each sidetone signal includes: The reference phase of any sidetone signal is used as the input parameter of the target exponential function, and a reference model of the sidetone signal is constructed through the target exponential function. Repeat the steps of using the reference phase of any sidetone signal as the input parameter of the target exponential function and constructing a reference model of the sidetone signal through the target exponential function until a reference model for each sidetone signal is constructed.

4. The interferometric processing method based on carrier phase reference according to any one of claims 1 to 3, characterized in that, The step of calculating the cross-correlation residual phase of each sidetone signal based on a reference model for each sidetone signal includes: Based on any side tone signal and its reference model, calculate the cross-correlation residual phase of the side tone signal in each integration period; Repeat the steps of calculating the cross-correlation residual phase of the side tone signal in each integration period based on any side tone signal and the reference model of the side tone signal, until the cross-correlation residual phase of each side tone signal in each integration period is obtained.

5. The interferometric processing method based on carrier phase reference according to claim 4, characterized in that, The step of calculating the group delay observable based on the cross-correlation residual phase of each sidetone signal includes: Based on the phase of the main carrier signal and the cross-correlation residual phase of each side tone signal at the first station in any integration period, and the phase of the main carrier signal and the cross-correlation residual phase of each side tone signal at the second station in any integration period, calculate the group delay observation of the first station and the second station in any integration period. Repeat the steps of calculating the group delay observation of the first station and the second station in any integration period based on the phase of the main carrier signal and the cross-correlation residual phase of each side tone signal of the first station in any integration period, and the phase of the main carrier signal and the cross-correlation residual phase of each side tone signal of the second station in any integration period, until the group delay observation of the first station and the second station in each integration period is calculated.

6. An interferometric processing system based on carrier phase reference, characterized in that, include: Acquisition module, processing module, and calculation module; The acquisition module is used to: construct a reference phase for each sidetone signal based on the phase of the main carrier signal and the design parameters of the detector beacon signal; The processing module is used to: construct a reference model for each sidetone signal based on the reference phase of each sidetone signal, and calculate the cross-correlation residual phase of each sidetone signal according to the reference model of each sidetone signal; The calculation module is used to: calculate the group delay observation based on the cross-correlation residual phase of each sidetone signal; The design parameters include: the frequency ratio of each sidetone signal to the main carrier signal; the acquisition module is specifically used for: Based on the frequency ratio of any side tone signal to the main carrier signal, and in combination with the phase of the main carrier signal, the reference phase of the side tone signal is determined. Repeat the step of determining the reference phase of the side tone signal based on the frequency ratio of any side tone signal to the main carrier signal and in combination with the phase of the main carrier signal, until the reference phase of each side tone signal is obtained; The formula for calculating group delay observations is: ;in, Represents the group delay observation. Indicates the main delay component. Indicates the residual delay component; The formula for calculating the main delay component is: ;in, This indicates the phase of the primary carrier signal at the first station. This indicates the phase of the main carrier signal at the second station; This indicates the nominal frequency of the main carrier signal. The nominal frequencies of the main carrier signals of the first and second stations are the same. The formula for calculating the residual time delay component is: ;in, This represents the least squares fitting function; This represents the set of differences in the residual phase of each sidetone signal between the first and second measuring stations. This represents the set of frequencies for each sidetone signal, and 1 represents the coefficient of the first-order term of the fitted polynomial. Least square fitting is used to obtain the residual time delay, which is taken as the coefficient of the first-order term of the fitting result.

7. The interferometric processing system based on carrier phase reference according to claim 6, characterized in that, Also includes: Build module; the build module is used for: Using a local correlation processing algorithm, the residual phase of the main carrier signal and the theoretical model are calculated; The phase of the main carrier signal is constructed based on the residual phase of the main carrier signal and the theoretical model.

8. An electronic device, characterized in that, The electronic 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 electronic device to implement the carrier phase reference-based interferometric processing method as described in any one of claims 1 to 5.

9. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores at least one computer program, which, when executed by a processor, implements the interferometric processing method based on carrier phase reference as described in any one of claims 1 to 5.

Images