Interference fringe detection method and system based on missing clock error information

Abstract

The present invention discloses an interference fringe detection method and system based on missing clock error information, and relates to the field of information acquisition and processing. The method comprises: obtaining the integral periodic complex interference phase of a radio source, performing phase compensation on the integral periodic complex interference phase according to preset clock error values ​​and clock speed values, obtaining a preset complex interference phase through the integral periodic complex interference phase after phase compensation, obtaining an accuracy assessment factor according to the preset complex interference phase, and detecting the corresponding target interference phase according to the accuracy assessment factor. This scheme is based on the integral periodic complex interference phase of a continuous tracking arc segment, and compensates the complex interference phase by setting the clock error and clock speed step length. By evaluating the compliance accuracy within the continuous tracking arc segment interference phase under traversal search, the acquisition of a priori clock error and clock speed information can be realized, and the application prospect in the field of multi-institution joint interference measurement and control support is broad.

Description

Technical Field

[0001] The present invention relates to the field of information acquisition and processing, and in particular to an interference fringe detection method and system based on missing clock error information. Background Art

[0002] Very Long Baseline Interferometry (VLBI) related processing requires a priori models of a certain accuracy as an aid, which generally include geometric delay, medium delay, link delay, and clock error, etc. When conducting multi-institutional joint interferometry data processing, the reliability of clock error information is reduced due to the obvious differences in the time synchronization characteristics between stations of different institutions. This problem is more prominent for stations that are not equipped with dedicated inter-station time synchronization measurement equipment. Due to the lack of clock error, it is impossible to effectively perform interferometric measurement data related processing.

[0003] Usually, the processing is done at the correlation guidance model end, by traversing the clock error and clock speed information to update the iterative correlation guidance model, and then judging the interference fringes through correlation calculation. This processing method requires re-correlation calculation for each iteration, which results in excessive calculation and low efficiency. Summary of the invention

[0004] The technical problem to be solved by the present invention is to provide an interference fringe detection method and system based on the lack of clock error information in view of the deficiencies in the prior art.

[0005] The technical solution of the present invention to solve the above technical problems is as follows:

[0006] An interference fringe detection method based on missing clock error information, comprising:

[0007] Obtain the integrated periodic complex interferometric phase of the radio source;

[0008] Performing phase compensation on the integrated period complex interference phase according to a preset clock error value and a clock speed value;

[0009] Obtaining a preset complex interference phase through the integrated periodic complex interference phase after phase compensation;

[0010] Obtaining an accuracy assessment factor according to the preset complex interference phase;

[0011] The corresponding target interference phase is detected according to the accuracy evaluation factor.

[0012] The beneficial effects of the present invention are as follows: this scheme is based on the integral periodic complex interference phase of the continuous tracking arc segment, and the complex interference phase is compensated by setting the clock error and clock speed step size. By evaluating the compliance accuracy within the continuous tracking arc segment interference phase under traversal search, the prior clock error and clock speed information can be obtained, and the application prospect is broad in the field of multi-mechanism joint interference measurement and control support.

[0013] In the case of missing clock errors, this scheme can be used to effectively perform interferometric measurement data correlation processing, which has the advantages of small computational complexity and high efficiency.

[0014] Furthermore, the step of obtaining the integrated periodic complex interference phase of the radio source specifically includes:

[0015] According to the continuous tracking arc segment of the radio source, the integrated periodic complex interferometric phase is obtained.

[0016] Furthermore, obtaining a preset complex interference phase by using the integrated periodic complex interference phase after phase compensation specifically includes:

[0017] The preset complex interference phase is obtained by the second formula, wherein the second formula is:

[0018]

[0019] Among them, φ scan Represents the preset complex interference phase, Ф icor (m) represents the integral period complex interference phase after phase compensation, m=1, 2, 3...M, M represents the number of data points contained in the complex interference phase, i=1, 2, 3...N, N represents the number of integral periods of continuous tracking arc segments.

[0020] Furthermore, obtaining the accuracy assessment factor according to the preset complex interference phase specifically includes:

[0021] Performing Fourier transform on the preset complex interference phase to obtain an accuracy evaluation factor.

[0022] Furthermore, the method further includes: obtaining the clock difference value and the clock speed value according to a third formula, wherein the third formula is:

[0023]

[0024] Among them, Clock step Represents the clock difference value, ClockR step represents the clock speed value, F s Indicates the acquisition signal bandwidth, T s represents the length of the integration period, l and k represent the scale factors respectively, and f represents the channel local oscillator frequency.

[0025] Another technical solution of the present invention to solve the above technical problems is as follows:

[0026] An interference fringe detection method based on missing clock error information, comprising: a phase acquisition module, a compensation module, a calculation module, a transformation module and an evaluation module;

[0027] The phase acquisition module is used to acquire the integrated periodic complex interference phase of the radio source;

[0028] The compensation module is used to perform phase compensation on the integrated period complex interference phase according to a preset clock error value and a clock speed value;

[0029] The calculation module is used to obtain a preset complex interference phase through the integrated periodic complex interference phase after phase compensation;

[0030] The transformation module is used to obtain an accuracy assessment factor according to the preset complex interference phase;

[0031] The evaluation module is used to detect the corresponding target interference phase according to the accuracy evaluation factor.

[0032] The beneficial effects of the present invention are as follows: this scheme is based on the integral periodic complex interference phase of the continuous tracking arc segment, and the complex interference phase is compensated by setting the clock error and clock speed step size. By evaluating the compliance accuracy within the continuous tracking arc segment interference phase under traversal search, the prior clock error and clock speed information can be obtained, and the application prospect is broad in the field of multi-mechanism joint interference measurement and control support.

[0033] In the case of missing clock errors, this scheme can be used to effectively perform interferometric measurement data correlation processing, which has the advantages of small computational complexity and high efficiency.

[0034] Furthermore, the phase acquisition module is specifically used to acquire the integrated periodic complex interference phase according to the continuous tracking arc segment of the radio source.

[0035] Furthermore, the calculation module is specifically used to obtain a preset complex interference phase through a second formula, wherein the second formula is:

[0036]

[0037] Among them, φ scan Represents the preset complex interference phase, Ф icor (m) represents the integral period complex interference phase after phase compensation, m=1, 2, 3...M, M represents the number of data points contained in the complex interference phase, i=1, 2, 3...N, N represents the number of integral periods of continuous tracking arc segments.

[0038] Furthermore, the transformation module is specifically used to perform Fourier transformation on the preset complex interference phase to obtain an accuracy assessment factor.

[0039] Furthermore, it also includes: a clock difference value and clock speed value acquisition module, which is used to obtain the clock difference value and the clock speed value according to a third formula, wherein the third formula is:

[0040]

[0041] Among them, Clock step Represents the clock difference value, ClockR step represents the clock speed value, F s Indicates the acquisition signal bandwidth, T s represents the length of the integration period, l and k represent the scale factors respectively, and f represents the channel local oscillator frequency.

[0042] Advantages of additional aspects of the present invention will be given in part in the following description, and in part will be obvious from the following description, or will be learned through practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] Figure 1 A schematic diagram of a method for detecting interference fringes based on missing clock error information provided by an embodiment of the present invention;

[0044] Figure 2 An embodiment of the present invention provides a schematic diagram of an interference fringe detection system based on missing clock error information;

[0045] Figure 3 Schematic diagram of interference phase without clock error information of an observed target provided by other embodiments of the present invention

[0046] Figure 4 Schematic diagram of interference fringe detection algorithm provided by other embodiments of the present invention;

[0047] Figure 5 Other embodiments of the present invention provide an internal conformity accuracy graph obtained by traversal search;

[0048] Figure 6 Schematic diagram of stripe detection results provided by other embodiments of the present invention. DETAILED DESCRIPTION

[0049] The principles and features of the present invention are described below in conjunction with the accompanying drawings. The examples given are only used to explain the present invention but not to limit the scope of the present invention.

[0050] like Figure 1 As shown, an interference fringe detection method based on missing clock error information provided by an embodiment of the present invention includes:

[0051] S1, obtain the integrated periodic complex interference phase of the radio source;

[0052] It should be noted that S1 specifically includes: using a priori model without clock error information to perform correlation processing to obtain N integral period complex interference phases Φ i (m); wherein the related processing may include: integer bit compensation, stripe rotation, and fractional bit compensation.

[0053] Among them, the integrated period complex interference phase is:

[0054] Ф i (m), m=1, 2...M; i=1, 2...N,

[0055] Wherein, N is the number of integration cycles of the continuous tracking arc segment (called SCAN), and M is the number of data points contained in the complex interferometric phase.

[0056] S2, performing phase compensation on the integrated period complex interference phase according to a preset clock error value and a clock speed value;

[0057] It should be noted that the preset clock difference value and clock speed value can be obtained according to a third formula, wherein the third formula is:

[0058]

[0059] Among them, Clock step Represents the clock difference value, ClockR step represents the clock speed value, F s Indicates the acquisition signal bandwidth, T s represents the length of the integration period, l and k represent the scale factors respectively, and f represents the channel local oscillator frequency.

[0060] S3, obtaining a preset complex interference phase through the integrated periodic complex interference phase after phase compensation; wherein the preset complex interference phase can be a SCAN complex interference phase:

[0061] It should be noted that the compensated integral period complex interference phase is combined with the second formula to obtain a preset complex interference phase, wherein the second formula is:

[0062]

[0063] Among them, φ scan Represents the preset complex interference phase, Ф icor (m) represents the integral period complex interference phase after phase compensation, m = 1, 2, 3...M, M represents the number of data points contained in the complex interference phase, i = 1, 2, 3...N, N represents the number of integral periods of continuous tracking arc segments

[0064] S4, obtaining an accuracy evaluation factor according to the preset complex interference phase; it should be noted that for φ sconPerform Fourier transform to obtain the accuracy assessment factor δ scan =FFT(φ scan ).

[0065] S5, detecting a corresponding target interference phase according to the accuracy evaluation factor.

[0066] It should be noted that the accuracy evaluation factor δ scan When the peak value is the largest and the peak point position is within 1 resolution unit, the prior clock error and clock speed information can be obtained, and the interference phase can be effectively detected; it should be noted that the accuracy evaluation factor δ scan When the above conditions are met, interference fringes can be obtained. Figure 6 As shown, a relatively concentrated straight line trend can be determined.

[0067] Among them, the interferometric phase of a certain observed target without clock error information is as follows Figure 3 As shown, the phase points are completely randomly distributed and no linear trend can be discerned.

[0068] This scheme is based on the integrated periodic complex interferometer phase of the continuous tracking arc segment. By setting the clock error and clock rate step, the complex interferometer phase is compensated. By evaluating the compliance accuracy of the continuous tracking arc segment interferometer phase of the radio source under ergodic search, the prior clock error and clock rate information can be obtained. It has broad application prospects in the field of multi-institution joint interferometer measurement and control support.

[0069] In the case of missing clock errors, this scheme can be used to effectively perform interferometric measurement data correlation processing, which has the advantages of small computational complexity and high efficiency.

[0070] Optionally, in any of the above embodiments, acquiring the integrated periodic complex interferometric phase of the radio source specifically includes:

[0071] According to the continuous tracking arc segment of the radio source, the integrated periodic complex interferometric phase is obtained.

[0072] Optionally, in any of the above embodiments, obtaining a preset complex interference phase by using the integrated periodic complex interference phase after phase compensation specifically includes:

[0073] The preset complex interference phase is obtained by the second formula, wherein the second formula is:

[0074]

[0075] Among them, φ scan Represents the preset complex interference phase, Ф icor (m) represents the integral period complex interference phase after phase compensation, m=1, 2, 3...M, M represents the number of data points contained in the complex interference phase, i=1, 2, 3...N, N represents the number of integral periods of continuous tracking arc segments.

[0076] In one embodiment, the integrated period complex interference phase compensation can be performed according to the fourth formula, wherein the fourth formula is:

[0077]

[0078] Among them, ck and ClockR are the clock difference and clock speed values ​​during a specific traversal.

[0079] Optionally, in any of the above embodiments, obtaining the accuracy assessment factor according to the preset complex interference phase specifically includes:

[0080] Performing Fourier transform on the preset complex interference phase to obtain an accuracy evaluation factor.

[0081] Optionally, in any of the above embodiments, the method further includes: obtaining the clock difference value and the clock speed value according to a third formula, wherein the third formula is:

[0082]

[0083] Among them, Clock step Represents the clock difference value, ClockR step represents the clock speed value, F s Indicates the acquisition signal bandwidth, T s represents the length of the integration period, l and k represent the scale factors respectively, and f represents the channel local oscillator frequency.

[0084] In one embodiment, if Figure 4 As shown, the interference fringe detection algorithm may include:

[0085] 1) Using the prior model without clock error information to perform correlation processing to obtain the integrated periodic complex interference phase;

[0086] 2) Set the clock error and clock speed step;

[0087] 3) Compensate for the complex interference phase of each integration period;

[0088] 4) Calculate the SCAN complex interferometer phase;

[0089] 5) FFT transformation;

[0090] 6) Calculate the interference phase internal coincidence accuracy.

[0091] Among them, the internal accuracy graph obtained by traversal search is as follows: Figure 5As shown in the figure, it is a three-dimensional graph, where x represents the clock error index and y represents the clock speed index. Each traversal search corresponds to an x / y coordinate, and the z coordinate value is obtained at the same time. The z in the left figure is the index of the peak value of the precision assessment factor, and the z in the right figure is the peak value of the precision assessment factor; among them, the clock error traversal range is -10 / Fs---10 / Fs; the clock speed traversal feedback is -1 / Tsf---1 / Tsf, and the fringe detection result is as follows Figure 6 Show.

[0092] In one embodiment, if Figure 2 As shown, an interference fringe detection method based on missing clock error information includes: a phase acquisition module, a compensation module, a calculation module, a transformation module and an evaluation module;

[0093] The phase acquisition module is used to acquire the integrated periodic complex interference phase of the radio source;

[0094] The compensation module is used to perform phase compensation on the integrated period complex interference phase according to a preset clock error value and a clock speed value;

[0095] The calculation module is used to obtain a preset complex interference phase through the integrated periodic complex interference phase after phase compensation;

[0096] The transformation module is used to obtain an accuracy assessment factor according to the preset complex interference phase;

[0097] The evaluation module is used to detect the corresponding target interference phase according to the accuracy evaluation factor.

[0098] This scheme is based on the integrated periodic complex interference phase of the continuous tracking arc segment. By setting the clock error and clock speed step, the complex interference phase is compensated. By evaluating the conformity accuracy of the continuous tracking arc segment interference phase under traversal search, the prior clock error and clock speed information can be obtained. It has broad application prospects in the field of multi-institution joint interference measurement and control support.

[0099] In the case of missing clock errors, this scheme can be used to effectively perform interferometric measurement data correlation processing, which has the advantages of small computational complexity and high efficiency.

[0100] Optionally, in any of the above embodiments, the phase acquisition module is specifically used to acquire the integrated periodic complex interference phase according to the continuous tracking arc segment of the radio source.

[0101] Optionally, in any of the above embodiments, the calculation module is specifically used to obtain a preset complex interference phase through a second formula, wherein the second formula is:

[0102]

[0103] Among them, φ scan Represents the preset complex interference phase, Фicor (m) represents the integral period complex interference phase after phase compensation, m=1, 2, 3...M, M represents the number of data points contained in the complex interference phase, i=1, 2, 3...N, N represents the number of integral periods of continuous tracking arc segments.

[0104] Optionally, in any of the above embodiments, the transformation module is specifically used to perform Fourier transform on the preset complex interference phase to obtain an accuracy assessment factor.

[0105] Optionally, in any of the above embodiments, the method further comprises: a clock difference value and clock speed value obtaining module, configured to obtain the clock difference value and the clock speed value according to a third formula, wherein the third formula is:

[0106]

[0107] Among them, Clock step Represents the clock difference value, ClockR step represents the clock speed value, F s Indicates the acquisition signal bandwidth, T s represents the length of the integration period, l and k represent the scale factors respectively, and f represents the channel local oscillator frequency.

[0108] It can be understood that in some embodiments, some or all of the optional implementation modes in the above embodiments may be included.

[0109] It should be noted that the above embodiments are product embodiments corresponding to the previous method embodiments. For the description of the optional implementation methods in the product embodiments, reference can be made to the corresponding description in the above method embodiments, which will not be repeated here.

[0110] The reader should understand that in the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "examples", "specific examples", or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine different embodiments or examples described in this specification and the features of different embodiments or examples without contradiction.

[0111] In the several embodiments provided in this application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the method embodiments described above are only illustrative, for example, the division of steps is only a logical function division, and there may be other division methods in actual implementation, such as multiple steps can be combined or integrated into another step, or some features can be ignored or not executed.

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

[0113] The above are only specific embodiments of the present invention, but the protection scope of the present invention is not limited thereto. Any technician familiar with the technical field can easily think of various equivalent modifications or replacements within the technical scope disclosed by the present invention, and these modifications or replacements should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention shall be based on the protection scope of the claims.

Claims

1. A method for detecting interference fringes based on missing clock error information, characterized in that: include: Obtain the integrated periodic complex interferometric phase of the radio source; Performing phase compensation on the integrated period complex interference phase according to a preset clock error value and a clock speed value; Obtaining a preset complex interference phase through the integrated periodic complex interference phase after phase compensation; Obtaining an accuracy assessment factor according to the preset complex interference phase; The corresponding target interference phase is detected according to the accuracy assessment factor; The step of obtaining the integrated periodic complex interference phase of the radio source specifically includes: The integrated periodic complex interference phase is obtained according to the continuous tracking arc segments of the radio source.

2. The interference fringe detection method based on missing clock error information according to claim 1, characterized in that: The step of obtaining a preset complex interference phase through the integrated periodic complex interference phase after phase compensation specifically includes: The preset complex interference phase is obtained by the second formula, wherein the second formula is: ， in, represents the preset complex interference phase, It represents the integral period complex interference phase after phase compensation, m=1,2,3…M, M represents the number of data points contained in the complex interference phase, i=1,2,3…N, N represents the number of integral periods of continuous tracking arc segments.

3. The interference fringe detection method based on missing clock error information according to any one of claims 1 to 2, characterized in that: The obtaining of the accuracy evaluation factor according to the preset complex interference phase specifically includes: Performing Fourier transform on the preset complex interference phase to obtain an accuracy evaluation factor.

4. The interference fringe detection method based on missing clock error information according to any one of claims 1 to 2, characterized in that: Also includes: According to the third formula, the clock difference value and the clock speed value are obtained, wherein the third formula is: ， in, represents the clock difference value, represents the clock speed value, represents the acquisition signal bandwidth, represents the length of the integration period, , represent the scale factors, Indicates the channel local oscillator frequency.

5. A method for detecting interference fringes based on missing clock error information, characterized in that: include: Phase acquisition module, compensation module, calculation module, transformation module and evaluation module; The phase acquisition module is used to acquire the integrated periodic complex interference phase of the radio source; The compensation module is used to perform phase compensation on the integrated period complex interference phase according to a preset clock error value and a clock speed value; The calculation module is used to obtain a preset complex interference phase through the integrated periodic complex interference phase after phase compensation; The transformation module is used to obtain an accuracy assessment factor according to the preset complex interference phase; The evaluation module is used to detect the corresponding target interference phase according to the accuracy evaluation factor; The phase acquisition module is specifically used to acquire the integrated periodic complex interference phase according to the continuous tracking arc segment of the radio source.

6. The interference fringe detection method based on missing clock error information according to claim 5, characterized in that: The calculation module is specifically used to obtain a preset complex interference phase through a second formula, wherein the second formula is: ， in, represents the preset complex interference phase, It represents the integral period complex interference phase after phase compensation, m=1,2,3…M, M represents the number of data points contained in the complex interference phase, i=1,2,3…N, N represents the number of integral periods of continuous tracking arc segments.

7. The interference fringe detection method based on missing clock error information according to any one of claims 5 to 6, characterized in that: The transformation module is specifically used to perform Fourier transformation on the preset complex interference phase to obtain an accuracy assessment factor.

8. The interference fringe detection method based on missing clock error information according to any one of claims 5 to 6, characterized in that: Also includes: The clock difference value and clock speed value acquisition module is used to obtain the clock difference value and the clock speed value according to a third formula, wherein the third formula is: ， in, represents the clock difference value, represents the clock speed value, represents the acquisition signal bandwidth, represents the length of the integration period, , represent the scale factors, Indicates the channel local oscillator frequency.

Images