CN111007551A - Multi-tone ranging ambiguity-resolving method in USB side tone ranging system - Google Patents

Multi-tone ranging ambiguity-resolving method in USB side tone ranging system Download PDF

Info

Publication number
CN111007551A
CN111007551A CN201911353373.6A CN201911353373A CN111007551A CN 111007551 A CN111007551 A CN 111007551A CN 201911353373 A CN201911353373 A CN 201911353373A CN 111007551 A CN111007551 A CN 111007551A
Authority
CN
China
Prior art keywords
sidetone
tone
phi
side tone
phase
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911353373.6A
Other languages
Chinese (zh)
Other versions
CN111007551B (en
Inventor
万扬洋
李志强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Space E Star Communication Technology Co ltd
Original Assignee
Space E Star Communication Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Space E Star Communication Technology Co ltd filed Critical Space E Star Communication Technology Co ltd
Priority to CN201911353373.6A priority Critical patent/CN111007551B/en
Publication of CN111007551A publication Critical patent/CN111007551A/en
Application granted granted Critical
Publication of CN111007551B publication Critical patent/CN111007551B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method

Abstract

The invention discloses a multi-tone distance measurement ambiguity resolving method in a USB side tone distance measurement system. The method comprises a ambiguity resolution stage and a tracking stage, wherein the side tones used for ranging comprise a first side tone to a seventh side tone, in the ambiguity resolution stage, the first side tone to the sixth side tone are respectively sent correspondingly from a first time to a sixth time, the duration of sending each side tone from the first side tone to the sixth side tone is the same, the seventh side tone starts to be sent at the first time and keeps being sent continuously, ambiguity resolution operation is carried out by utilizing the first side tone to the seventh side tone, the accurate phase increment of the seventh side tone is obtained, and the distance from a satellite to the ground is calculated; in the tracking stage, the seventh sidetone is still continuously transmitted, and the phase increment of the seventh sidetone is continuously tracked, so that the continuous tracking measurement of the satellite-ground distance is realized. The fuzzy solution operation method avoids the process of repeated iteration in the prior art, has high operation efficiency and is beneficial to the realization of hardware circuits.

Description

Multi-tone ranging ambiguity-resolving method in USB side tone ranging system
Technical Field
The invention relates to the field of satellite measurement and control, in particular to a multi-tone distance measurement ambiguity resolving method in a USB side tone distance measurement system.
Background
In the satellite USB distance measuring system, a sinusoidal distance measuring sound signal with the phase phi 1 is sent from a ground measurement and control station to a satellite, and the phase of the signal is changed into phi 2 after the signal returns to the ground through a satellite transponder. The distance between the ground station and the satellite can be calculated by measuring the difference phi between the transmitted signal and the received signal phi 2-phi 1. However, due to the periodicity of the phases, the maximum value of the measured phase difference cannot exceed 360 ° because the precise value cannot be determined when the phase difference exceeds 360 °. In practical application, in order to solve the technical problem, when the ranging signal is transmitted, a set of ranging signals with different frequencies is transmitted to perform range ambiguity resolution instead of a single side tone ranging signal with a certain frequency, but when the ranging signals with multiple frequencies are used for performing ambiguity resolution calculation, phase values or virtual phase values of all side tone signals with frequencies higher than the current side tone participate in calculation, so that the calculation process is complex and the efficiency is low.
Disclosure of Invention
The invention mainly solves the technical problem of providing a multi-tone distance measurement ambiguity resolving method in a USB side tone distance measurement system, and solves the problems of complex calculation process and low efficiency of the multi-tone distance measurement ambiguity resolving method in the prior art.
In order to solve the technical problem, one technical scheme adopted by the invention is to provide a multi-tone ranging ambiguity-resolving method in a USB sidetone ranging system, which comprises the following steps: the method comprises the steps that a ground satellite ranging terminal generates and sends sidetone for ranging, the sidetone is modulated by a transmitter of the ground satellite ranging terminal and then is transmitted by an antenna, the sidetone is transmitted by a satellite transponder, and the sidetone is demodulated by a receiver of the ground satellite ranging terminal and then is subjected to ranging calculation; the distance measurement calculation comprises an ambiguity resolution stage and a tracking stage, and the sidetones comprise a first sidetone, a sixth sidetone and a seventh sidetone; in the ambiguity resolution stage, the first sidetone to the sixth sidetone are respectively sent correspondingly from the first time to the sixth time, the duration of sending each sidetone from the first sidetone to the sixth sidetone is the same, the seventh sidetone also starts to be sent at the first time and keeps being sent continuously, ambiguity resolution operation is carried out by using the first sidetone to the seventh sidetone, the accurate phase increment of the seventh sidetone is obtained, and the distance from the satellite to the ground is calculated; in the tracking stage, the seventh sidetone is still continuously transmitted, and the phase increment of the seventh sidetone is continuously tracked, so that the continuous tracking measurement of the satellite-ground distance is realized.
Preferably, in the deblurring phase, the first sidetone f is obtained at a second instant in time1Is phase estimate phi1And seventh lateral soundf7Is phase estimate phi7Then separately calculate
Figure BDA0002335249700000021
Wherein n is1Is a serial number, the value range satisfies n1＝0,1,2,……,(Fb71-1)，Fb71Representing the frequency ratio of the seventh side tone to the first side tone, and further calculating
Figure BDA0002335249700000022
Find the minimum value delta phi among them1,minAnd said minimum value Δ φ1,minCorresponding number value n1,min。
Preferably, in the deblurring phase, the second sidetone f is obtained at a third instant2Is phase estimate phi2And a first sidetone f1Is phase estimate phi1Then separately calculate
Figure BDA0002335249700000023
Wherein n is2Is a serial number, the value range satisfies n2＝0,1,2,……,(Fb12-1)，Fb12Representing the frequency ratio of the first sidetone and the second sidetone, and further calculating
Figure BDA0002335249700000024
Find the minimum value delta phi among them2,minAnd a minimum value of delta phi2,minCorresponding number value n2,min。
Preferably, in the deblurring stage, the third sidetone f is obtained at a fourth moment in time3Is phase estimate phi3And a second sidetone f2Is phase estimate phi2Then separately calculate
Figure BDA0002335249700000025
Wherein n is3Is a serial number, the value range satisfies n3＝0，1，2，……，(Fb23-1)，Fb23Representing the frequency ratio of the second side tone and the third side tone, and further calculating
Figure BDA0002335249700000026
Find the minimum value delta phi among them3,minAnd the minimum value delta phi3,minCorresponding number value n3,min。
Preferably, in the deblurring phase, the fourth side note f is obtained at a fifth instant4Is phase estimate phi4And a third side tone f3Is phase estimate phi3Then separately calculate
Figure BDA0002335249700000031
Wherein n is4Is a serial number, the value range satisfies n4＝0,1,2,……,(Fb34-1)，Fb34Representing the frequency ratio of the third side sound and the fourth side sound, and further calculating
Figure BDA0002335249700000032
Find the minimum value delta phi among them4,minAnd the minimum value delta phi4,minCorresponding number value n4,min。
Preferably, in the deblurring phase, the fifth sidetone f is obtained at a sixth instant5Is phase estimate phi5And a fourth side tone f4Is phase estimate phi4Then separately calculate
Figure BDA0002335249700000033
Wherein n is5Is a serial number, the value range satisfies n5＝0,1,2,……,(Fb45-1)，Fb45Representing the frequency ratio of the fourth side tone and the fifth side tone, and further calculating
Figure BDA0002335249700000034
Find the minimum value delta phi among them5,minAnd the minimum value delta phi5,minCorresponding number value n5,min。
Preferably, in the deblurring phase, a sixth sidetone f is obtained at a seventh instant6Is phase estimate phi6And a fifth sidetone f5Is phase estimate phi5Then separately calculate
Figure BDA0002335249700000035
Wherein n is6Is a serial number, the value range satisfies n6＝0,1,2,……,(Fb56-1)，Fb56Representing the frequency ratio of the fifth side tone and the sixth side tone, and further calculating
Figure BDA0002335249700000036
Find the minimum value delta phi among them6,minAnd a minimum value of delta phi6,minCorresponding number value n6,min。
Preferably, in the stage of deblurring, the phase increment of the seventh sidetone is calculated as:
Δφ＝φ7+(n1,min+n2,min×5+n3,min×52+n4,min×53+n5,min×54+n6,min×55)×2π
then further calculating the distance between the star and the ground according to the delta phi:
Figure BDA0002335249700000037
preferably, the seventh sidetone f7100kHz, first sidetone f120kHz, second sidetone f24kHz, third sidetone f3800Hz, fourth side tone f4160Hz, fifth sidetone f532Hz, sixth sidetone f6The corresponding frequency ratio relationship is as follows: fb71＝5，Fb12＝5，Fb23＝5，Fb34＝5，Fb45＝5，Fb56＝4。
The invention has the beneficial effects that: the invention discloses a multi-tone distance measurement ambiguity resolving method in a USB side tone distance measurement system. The method comprises a ambiguity resolution stage and a tracking stage, wherein the side tones used for ranging comprise a first side tone to a seventh side tone, in the ambiguity resolution stage, the first side tone to the sixth side tone are respectively sent correspondingly from a first time to a sixth time, the duration of sending each side tone from the first side tone to the sixth side tone is the same, the seventh side tone starts to be sent at the first time and keeps being sent continuously, ambiguity resolution operation is carried out by utilizing the first side tone to the seventh side tone, the accurate phase increment of the seventh side tone is obtained, and the distance from a satellite to the ground is calculated; in the tracking stage, the seventh sidetone is still continuously transmitted, and the phase increment of the seventh sidetone is continuously tracked, so that the continuous tracking measurement of the satellite-ground distance is realized. The fuzzy solution operation method avoids the process of repeated iteration in the prior art, has high operation efficiency and is beneficial to the realization of hardware circuits.
Drawings
FIG. 1 is a diagram illustrating an embodiment of a multi-tone ranging disambiguation method in a USB sidetone ranging system according to the present invention;
FIG. 2 is a flowchart illustrating another exemplary method for ambiguity resolution in a multi-tone range finding system according to the USB sidetone range finding system of the present invention.
Detailed Description
In order to facilitate an understanding of the invention, the invention is described in more detail below with reference to the accompanying drawings and specific examples. Preferred embodiments of the present invention are shown in the drawings. This invention may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete.
It is to be noted that, unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
In the embodiment shown in fig. 1, in the USB sidetone ranging system, the first sidetone f1, the second sidetone f2, the second sidetone f … … and the sixth sidetone f6 are respectively transmitted at the first time t1, the second time t2, the … … and the sixth time t6, and the seventh sidetone f7 is continuously transmitted throughout the continuous interval from the first time t1 to the sixth time t 6. The side tones are generated by the ground satellite ranging terminal, are modulated by the transmitter, are transmitted by the ground antenna, are forwarded by the satellite transponder, and enter the ground satellite ranging terminal after being demodulated by the receiver. During each group of side tone continuous transmission period, the matching and de-blurring circuit carries out successive matching on the received side tone signals. It can be seen that the working process of the ranging can be divided into a deblurring stage (referred to as a matching stage) and a tracking stage, which are two stages. In the deblurring phase, the side tones are transmitted two by two, wherein the seventh side tone is continuously transmitted all the time, and the other six tones are transmitted in sequence at different time instants.
Preferably, the seventh side sound f7 and the first side sound f1 are transmitted from the first time t1 to the second time t2, the 1 st matching is performed at the second time t2, and the phase ambiguity of the seventh side sound f7 is resolved by the first side sound f 1.
Further, during the second time t2 to the third time t3, the seventh side tone f7 and the second side tone f2 are transmitted while tracking the phase of the second side tone f2, and the first side tone f1 is virtually tracked by using the tracked phase increment. At time t3, the phase ambiguity of f2 is recovered from the virtually tracked f1 and f2, and then the 2 nd matching is performed. The matching is carried out in 2 steps: step 1 solves the phase ambiguity of f1 for virtual tracking with f 2: step 2 further resolves the phase ambiguity of f7 with a virtual f 1.
Further, during the period from the third time t3 to the third time t4, the seventh side tone f7 and the third side tone f3 are transmitted while continuously tracking the phase of the third side tone f3, and the first side tone f1 and the second side tone f2 are virtually tracked by using the tracked phase increment. At time t4, the phase ambiguity of f3 is recovered from the virtually tracked f1 and f2, and then the 3 rd matching is performed. The matching is carried out in 3 steps: the phase ambiguity of f2 of the virtual tracking is resolved by f3 in step 1, the phase ambiguity of f1 is further resolved by virtual f2 in step 2, and the phase ambiguity of f7 is further resolved by virtual f1 in step 3.
Similarly, the 4 th matching is performed at the time t5, and the phase ambiguities of the virtual phases f4 to f1 and f7 are solved in 4 steps; performing the 5 th matching at the time t6, and resolving the phase ambiguities of the virtual f 5-f 1 and f7 in 5 steps; the 6 th match is made at time t7, and the phase ambiguities for virtual f6 through f1 and f7 are resolved in 6 steps.
After the above-mentioned step-by-step successive matching, the phase of each side tone without ambiguity, that is, the actual distance, can be obtained at time t 7. After time t7, the satellite range can be measured by continuously measuring the phase increment (i.e., the range increment) and adding the range increment to the range at time t 7.
Preferably, the seventh sidetone f7 is 100kHz, the first sidetone f1 is 20kHz, the second sidetone f2 is 4kHz, the third sidetone f3 is 800Hz, the fourth sidetone f4 is 160Hz, the fifth sidetone f5 is 32Hz, and the sixth sidetone f6 is 8 Hz. The frequency ratio relationship between them is: f 7: f1 ═ 5, f 1: f2 ═ 5, f 2: f3 ═ 5, f 3: f4 ═ 5, f 4: f5 ═ 5, f 5: f6 ═ 4.
Specifically, for example, taking the time t1 as an example, if the phase value measured by f7(100kHz) is 184 °, and the phase value measured by f1(20kHz) is 182 °, the matching algorithm can be expressed as follows:
(1)(φ20kHzx frequency ratio)/360 ° (182 ° × 5)/360 ° (2 +190 °), rounding off, i.e. leaving the integer quotient 2, and rounding off the remainder 190 °;
(2) (360 ℃ X. integer quotient 2) + φ100kHz＝904°；
(3)904°÷5＝180.8°。
Through the above 3 steps, the phase precision of the lOOkHz is transferred to the phase of 20kHz, i.e. the phase of 20kHz is determined to be 180.8 °, and the phase ambiguity of the lookkhz is solved, i.e. the phase of the lookkhz is determined to be 904 °. The above process is performed at the following time points, so that the phase ambiguity of each side tone can be resolved step by step.
It can be seen that, although the method adopted in the foregoing embodiment can correctly resolve the fuzzy range of each ranging tone, in the process of engineering implementation, each time the fuzzy calculation is performed, the phase values or virtual phase values of all side tones (f2 and f1) before the current side tone (e.g., f3) need to participate in the calculation, and the calculation process is complex and inefficient.
Further preferably, the invention also provides another embodiment of the deblurring method. Similar to the previous embodiment, as shown in fig. 2, firstly, step S101 is performed, in which a terrestrial satellite ranging terminal generates and transmits a sidetone for ranging, the sidetone is modulated by a transmitter of the terrestrial satellite ranging terminal, transmitted by an antenna, forwarded by a satellite transponder, and demodulated by a receiver of the terrestrial satellite ranging terminal to perform ranging calculation; the distance measurement calculation comprises an ambiguity resolution stage and a tracking stage, and the sidetones comprise a first sidetone, a sixth sidetone and a seventh sidetone; step S102, in the ambiguity resolution stage, the first sidetone to the sixth sidetone are respectively sent correspondingly from the first time to the sixth time, the duration of sending each sidetone from the first sidetone to the sixth sidetone is the same, the seventh sidetone also starts to be sent at the first time and keeps being sent continuously, ambiguity resolution operation is carried out on the first sidetone to the seventh sidetone, accurate phase increment of the seventh sidetone is obtained, and the distance from the satellite to the ground is calculated; and step S103, in the tracking stage, the seventh sidetone is still continuously transmitted, and the phase increment of the seventh sidetone is continuously tracked, so that the continuous tracking measurement of the satellite-ground distance is realized.
Preferably, in the deblurring phase, at a second instant t2Obtaining a first sidetone f at a time1Is phase estimate phi1And seventh sidetone f7Is phase estimate phi7Then separately calculate
Figure BDA0002335249700000071
Wherein n is1Is a serial number, the value range satisfies n1＝0,1,2,……,(Fb71-1)，Fb71Indicating the frequency ratio of the seventh sidetone to the first sidetone, e.g. if Fb71If 5, then there will be phi1，0＝φ1/5，φ1,1＝(φ1+2π)/5，φ1,2＝(φ1+4π)/5，φ1,3＝(φ1+6π)/5，φ1,4＝(φ1+8 π)/5; further separately calculate
Figure BDA0002335249700000072
E.g. having a value of delta phi1,0，Δφ1,1，Δφ1,2，Δφ1,3，Δφ1,4At a value of delta phi1，0～Δφ1，4In which the minimum value delta phi is found1,minE.g. delta phi1,2Is the minimum value among them, the minimum value delta phi1,minThe corresponding serial number value is n1,minE.g. delta phi1,2The number value corresponding to the minimum value is 2, i.e. n1,min＝2。
Further, at a third time t3Obtaining the second sidetone f at a moment2Is phase estimate phi2And a first sidetone f1Is phase estimate phi1Then separately calculate
Figure BDA0002335249700000073
Wherein n is2Is a serial number, the value range satisfies n2＝0,1,2,……,(Fb12-1)，Fb12Indicating the frequency ratio of the first sidetone to the second sidetone, e.g. if Fb12If 5, then there will be phi2,0＝φ2/5，φ2,1＝(φ2+2π)/5，φ2,2＝(φ2+4π)/5，φ2,3＝(φ2+6π)/5，φ2,4＝(φ2+8 π)/5; further separately calculate
Figure BDA0002335249700000074
E.g. having a value of delta phi2,0，Δφ2,1，Δφ2,2，Δφ2,3，Δφ2,4At a value of delta phi2,0～Δφ2,4In which the minimum value delta phi is found2,minE.g. delta phi2,3Is the minimum value among them, the minimum value delta phi2,minThe corresponding serial number value is n2,minE.g. delta phi2,3The number value corresponding to the minimum value is 3, i.e. n2,min＝3。
Further, at a fourth time t4Obtaining a third sidetone f3Is phase estimate phi3And a second sidetone f2Is phase estimate phi2Then separately calculate
Figure BDA0002335249700000075
Wherein n is3Is a serial number and has a value range of fullFoot n3＝0,1,2,……,(Fb23-1)，Fb23Indicating the frequency ratio of the second sidetone to the third sidetone, e.g. if Fb23If 5, then there will be phi3,0＝φ3/5，φ3,1＝(φ3+2π)/5，φ3,2＝(φ3+4π)/5，φ3,3＝(φ3+6π)/5，φ3,4＝(φ3+8 π)/5; further separately calculate
Figure BDA0002335249700000081
E.g. having a value of delta phi3,0，Δφ3,1，Δφ3,2，Δφ3,3，Δφ3,4At a value of delta phi3,0～Δφ3,4In which the minimum value delta phi is found3,minE.g. delta phi3,1Is the minimum value among them, the minimum value delta phi3,minThe corresponding serial number value is n3,minE.g. delta phi3,1The number value corresponding to the minimum value is 1, i.e. n3,min＝1。
Further, at a fifth time t5Obtaining a fourth sidetone f4Is phase estimate phi4And a third side tone f3Is phase estimate phi3Then separately calculate
Figure BDA0002335249700000082
Wherein n is4Is a serial number, the value range satisfies n4＝0,1,2,……,(Fb34-1)，Fb34Indicating the frequency ratio of the third sidetone to the fourth sidetone, e.g. if Fb34If 5, then there will be phi4,0＝φ4/5，φ4,1＝(φ4+2π)/5，φ4,2＝(φ4+4π)/5，φ4,3＝(φ4+6π)/5，φ4,4＝(φ4+8 π)/5; further separately calculate
Figure BDA0002335249700000083
E.g. having a value of delta phi4,0，Δφ4,1，Δφ4,2，Δφ4,3，Δφ4,4At a value of delta phi4,0～Δφ4,4Is found inMinimum value of which delta phi4,minE.g. delta phi4,2Is the minimum value among them, the minimum value delta phi4,minThe corresponding serial number value is n4,minE.g. delta phi4,2Is the minimum value, the corresponding index value is 2, i.e. n4,min＝2。
Further, at a sixth time t6Obtaining a fifth sidetone f5Is phase estimate phi5And a fourth side tone f4Is phase estimate phi4Then separately calculate
Figure BDA0002335249700000084
Wherein n is5Is a serial number, the value range satisfies n5＝0,1,2,……,(Fb45-1)，Fb45Indicating the frequency ratio of the fourth sidetone to the fifth sidetone, e.g. if Fb45If 5, then there will be phi5,0＝φ5/5，φ5,1＝(φ5+2π)/5，φ5,2＝(φ5+4π)/5，φ5,3＝(φ5+6π)/5，φ5,4＝(φ5+8 π)/5; further separately calculate
Figure BDA0002335249700000085
E.g. having a value of delta phi5,0，Δφ5,1，Δφ5,2，Δφ5,3，Δφ5,4At a value of delta phi5,0～Δφ5,4In which the minimum value delta phi is found5,minE.g. delta phi5,0Is the minimum value among them, the minimum value delta phi5,minThe corresponding serial number value is n5,minE.g. delta phi5,0Is the minimum value, the corresponding index value is 0, i.e. n5,min＝0。
Further, at a seventh time t7Obtaining the sixth sidetone f6Is phase estimate phi6And a fifth sidetone f5Is phase estimate phi5Then separately calculate
Figure BDA0002335249700000091
Wherein n is6Is a serial number, the value range satisfies n6＝0,1,2,……,(Fb56-1)，Fb56Indicating the frequency ratio of the fifth sidetone to the sixth sidetone, e.g. if Fb56If 4, then there will be phi6,0＝φ6/5，φ6,1＝(φ6+2π)/5，φ6,2＝(φ6+4π)/5，φ6,3＝(φ6+6 π)/5; further separately calculate
Figure BDA0002335249700000092
E.g. having a value of delta phi6,0，Δφ6,1，Δφ6,2，Δφ6,3At a value of delta phi6，0～Δφ6，3In which the minimum value delta phi is found6,minE.g. delta phi6,3Is the minimum value among them, the minimum value delta phi6,minThe corresponding serial number value is n6,minE.g. delta phi6,3Is the minimum value, the corresponding index value is 3, i.e. n6,min＝3。
From the above calculation, it can be seen that the embodiment only needs to calculate the phase estimation values of the current sidetone and the previous sidetone, and the implementation complexity is greatly simplified. Compared with the current common ambiguity resolution mode, the embodiment not only reduces the calculation complexity of the USB sidetone ranging, but also mainly takes multiplication as a main part, cancels division used in the calculation of the prior art, enables hardware ambiguity resolution to become a simple step in engineering realization, and enables the engineering realization process to be simpler and more flexible.
Based on the above calculation, the seventh side-tone phase difference after the ambiguity resolution is finally obtained as follows:
Δφ＝φ7+(n1,min+n2,min×5+n3,min×52+n4,min×53+n5,min×54+n6,min×55)×2π
then further calculating the distance between the star and the ground according to the delta phi:
Figure BDA0002335249700000093
preferably, the seventh sidetone f7100kHz, first sidetone f120kHz, second sidetone f24kHz, third sidetone f3800Hz, fourth side tone f4160Hz, fifth sidetone f532Hz, sixth sidetone f6The corresponding frequency ratio relationship is as follows: fb71＝5，Fb12＝5，Fb23＝5，Fb34＝5，Fb45＝5，Fb56＝4。
Therefore, the invention discloses a multi-tone distance measurement ambiguity resolution method in a USB side tone distance measurement system. The method comprises a ambiguity resolution stage and a tracking stage, wherein the side tones used for ranging comprise a first side tone to a seventh side tone, in the ambiguity resolution stage, the first side tone to the sixth side tone are respectively sent correspondingly from a first time to a sixth time, the duration of sending each side tone from the first side tone to the sixth side tone is the same, the seventh side tone starts to be sent at the first time and keeps being sent continuously, ambiguity resolution operation is carried out by utilizing the first side tone to the seventh side tone, the accurate phase increment of the seventh side tone is obtained, and the distance from a satellite to the ground is calculated; in the tracking stage, the seventh sidetone is still continuously transmitted, and the phase increment of the seventh sidetone is continuously tracked, so that the continuous tracking measurement of the satellite-ground distance is realized. The fuzzy solution operation method avoids the process of repeated iteration in the prior art, has high operation efficiency and is beneficial to the realization of hardware circuits.
The above description is only an embodiment of the present invention, and not intended to limit the scope of the present invention, and all equivalent structural changes made by using the contents of the present specification and the drawings, or applied directly or indirectly to other related technical fields, are included in the scope of the present invention.

Claims (9)

1. A multi-tone distance measurement ambiguity resolution method in a USB side tone distance measurement system is characterized by comprising the following steps:
the method comprises the steps that a ground satellite ranging terminal generates and sends sidetone for ranging, the sidetone is modulated by a transmitter of the ground satellite ranging terminal and then is transmitted by an antenna, the sidetone is transmitted by a satellite transponder, and the sidetone is demodulated by a receiver of the ground satellite ranging terminal and then is subjected to ranging calculation; the distance measurement calculation comprises an ambiguity resolution stage and a tracking stage, and the sidetones comprise a first sidetone, a sixth sidetone and a seventh sidetone;
in the ambiguity resolution stage, the first sidetone to the sixth sidetone are respectively sent correspondingly from the first time to the sixth time, the duration of sending each sidetone from the first sidetone to the sixth sidetone is the same, the seventh sidetone also starts to be sent at the first time and keeps being sent continuously, ambiguity resolution operation is carried out by using the first sidetone to the seventh sidetone, the accurate phase increment of the seventh sidetone is obtained, and the distance from the satellite to the ground is calculated;
in the tracking stage, the seventh sidetone is still continuously transmitted, and the phase increment of the seventh sidetone is continuously tracked, so that the continuous tracking measurement of the satellite-ground distance is realized.
2. The multi-tone ranging disambiguation method of claim 1 wherein the first sidetone f is obtained at the second time during the disambiguation phase1Is phase estimate phi1And seventh sidetone f7Is phase estimate phi7Then separately calculate
Figure FDA0002335249690000011
Wherein n is1Is a serial number, the value range satisfies n1＝0,1,2,……,(Fb71-1)，Fb71Representing the frequency ratio of the seventh side tone to the first side tone, and further calculating
Figure FDA0002335249690000012
Find the minimum value delta phi among them1,minAnd said minimum value Δ φ1,minCorresponding number value n1,min。
3. The multi-tone ranging disambiguation method of claim 2 wherein the second sidetone f is obtained at a third time during the disambiguation phase2Is phase estimate phi2And a first sidetone f1Phase ofBit estimate phi1Then separately calculate
Figure FDA0002335249690000013
Wherein n is2Is a serial number, the value range satisfies n2＝0,1,2,……,(Fb12-1)，Fb12Representing the frequency ratio of the first sidetone and the second sidetone, and further calculating
Figure FDA0002335249690000021
Find the minimum value delta phi among them2,minAnd a minimum value of delta phi2,minCorresponding number value n2,min。
4. The multi-tone ranging disambiguation method of claim 3 wherein the third sidetone f is obtained at a fourth time during the disambiguation stage3Is phase estimate phi3And a second sidetone f2Is phase estimate phi2Then separately calculate
Figure FDA0002335249690000022
Wherein n is3Is a serial number, the value range satisfies n3＝0，1，2，……，(Fb23-1)，Fb23Representing the frequency ratio of the second side tone and the third side tone, and further calculating
Figure FDA0002335249690000023
Find the minimum value delta phi among them3,minAnd the minimum value delta phi3,minCorresponding number value n3,min。
5. The multi-tone ranging disambiguation method of claim 4 in the USB sidetone ranging regime, wherein the fourth side tone f is obtained at a fifth time during the disambiguation phase4Is phase estimate phi4And a third side tone f3Is phase estimate phi3Then separately calculate
Figure FDA0002335249690000024
Wherein n is4Is a serial number, the value range satisfies n4＝0,1,2,……,(Fb34-1)，Fb34Representing the frequency ratio of the third side sound and the fourth side sound, and further calculating
Figure FDA0002335249690000025
Find the minimum value delta phi among them4,minAnd the minimum value delta phi4,minCorresponding number value n4,min。
6. The multi-tone ranging disambiguation method of claim 5 wherein the fifth sidetone f is obtained at the sixth time during the disambiguation phase5Is phase estimate phi5And a fourth side tone f4Is phase estimate phi4Then separately calculate
Figure FDA0002335249690000026
Wherein n is5Is a serial number, the value range satisfies n5＝0,1,2,……,(Fb45-1)，Fb45Representing the frequency ratio of the fourth side tone and the fifth side tone, and further calculating
Figure FDA0002335249690000027
Find the minimum value delta phi among them5,minAnd the minimum value delta phi5,minCorresponding number value n5,min。
7. The multi-tone ranging disambiguation method of claim 6 wherein the sixth sidetone f is obtained at the seventh time during the disambiguation phase6Is phase estimate phi6And a fifth sidetone f5Is phase estimate phi5Then separately calculate
Figure FDA0002335249690000028
Wherein n is6Is a serial number, the value range satisfies n6＝0,1,2,……,(Fb56-1)，Fb56Representing the frequency ratio of the fifth side tone and the sixth side tone, and further calculating
Figure FDA0002335249690000031
Find the minimum value delta phi among them6,minAnd a minimum value of delta phi6,minCorresponding number value n6,min。
8. The multi-tone ranging ambiguity resolution method in the USB sidetone ranging system according to claim 7, wherein in the ambiguity resolution stage, the phase increment of the seventh sidetone is calculated as:
Δφ＝φ7+(n1,min+n2,min×5+n3,min×52+n4,min×53+n5,min×54+n6,min×55)×2π
then further calculating the distance between the star and the ground according to the delta phi:
Figure FDA0002335249690000032
9. the multi-tone ranging disambiguation method of claim 8 wherein the seventh sidetone f is a multi-tone ranging signal7100kHz, first sidetone f120kHz, second sidetone f24kHz, third sidetone f3800Hz, fourth side tone f4160Hz, fifth sidetone f532Hz, sixth sidetone f6The corresponding frequency ratio relationship is as follows: fb71＝5，Fb12＝5，Fb23＝5，Fb34＝5，Fb45＝5，Fb56＝4。
CN201911353373.6A 2019-12-25 2019-12-25 Multi-tone ranging ambiguity-resolving method in USB side tone ranging system Active CN111007551B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911353373.6A CN111007551B (en) 2019-12-25 2019-12-25 Multi-tone ranging ambiguity-resolving method in USB side tone ranging system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911353373.6A CN111007551B (en) 2019-12-25 2019-12-25 Multi-tone ranging ambiguity-resolving method in USB side tone ranging system

Publications (2)

Publication Number Publication Date
CN111007551A true CN111007551A (en) 2020-04-14
CN111007551B CN111007551B (en) 2021-01-22

Family

ID=70117892

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911353373.6A Active CN111007551B (en) 2019-12-25 2019-12-25 Multi-tone ranging ambiguity-resolving method in USB side tone ranging system

Country Status (1)

Country Link
CN (1) CN111007551B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101726746A (en) * 2009-12-14 2010-06-09 北京航空航天大学 Intermediate frequency direct sequence spread spectrum receiver for satellite ranging
CN101738598A (en) * 2008-11-19 2010-06-16 中国科学院国家天文台 Method for performing real-time measurement on track and position of satellite or aircraft
CN105391489A (en) * 2015-10-18 2016-03-09 中国电子科技集团公司第十研究所 Noncoherent range finding method for data link in measurement and control of unmanned aerial vehicle
CN105785351A (en) * 2016-03-10 2016-07-20 中国人民解放军国防科学技术大学 Method and system based on satellite fuzzy range finding
CN109347539A (en) * 2018-10-10 2019-02-15 北京遥感设备研究所 A kind of absolute time delay high accuracy calibration method of USB TT&C Transponder
CN110138470A (en) * 2019-04-17 2019-08-16 北京遥感设备研究所 A kind of USB TT&C Transponder sidetone power number normalization system and method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101738598A (en) * 2008-11-19 2010-06-16 中国科学院国家天文台 Method for performing real-time measurement on track and position of satellite or aircraft
CN101726746A (en) * 2009-12-14 2010-06-09 北京航空航天大学 Intermediate frequency direct sequence spread spectrum receiver for satellite ranging
CN105391489A (en) * 2015-10-18 2016-03-09 中国电子科技集团公司第十研究所 Noncoherent range finding method for data link in measurement and control of unmanned aerial vehicle
CN105785351A (en) * 2016-03-10 2016-07-20 中国人民解放军国防科学技术大学 Method and system based on satellite fuzzy range finding
CN109347539A (en) * 2018-10-10 2019-02-15 北京遥感设备研究所 A kind of absolute time delay high accuracy calibration method of USB TT&C Transponder
CN110138470A (en) * 2019-04-17 2019-08-16 北京遥感设备研究所 A kind of USB TT&C Transponder sidetone power number normalization system and method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
朱克勤 等: "纯音测距体制中软件解模糊的实现方法", 《上海航天》 *

Also Published As

Publication number Publication date
CN111007551B (en) 2021-01-22

Similar Documents

Publication Publication Date Title
JP7301771B2 (en) PHASE CORRECTOR, RANGING DEVICE AND PHASE VARIATION DETECTION DEVICE
US20090138744A1 (en) Multiplier device with suppression of higher-order distortion
CN101729461B (en) System and method for eliminating single-frequency interference and multi-frequency interference
CN110501729B (en) GNSS signal capturing method based on FPGA step-by-step code phase refinement
US8798125B2 (en) Phase tracking in communications systems
CN111007551B (en) Multi-tone ranging ambiguity-resolving method in USB side tone ranging system
TWI653841B (en) Carrier frequency offset estimation device and carrier frequency offset estimation method
CN101309100B (en) I/Q phase and gain compensating apparatus and method of orthogonal receiver
JP2017183818A (en) Measurement device for synchronous broadcasting and reception device for synchronous broadcasting
JP6457299B2 (en) Radar equipment
JP4169393B2 (en) Angle demodulator
CN113132284B (en) Carrier phase tracking method and device
CN108337202B (en) System and method for frequency estimation using frequency offset bias removal
CN112014864B (en) High dynamic satellite tracking method based on steepest ascent method
JP3945623B2 (en) Frequency synchronization method and OFDM receiver using the same
US20220091253A1 (en) Distance measurement device and distance measurement method
RU2357359C2 (en) Device for synchronising receiver of noise-like signals with minimal frequency-shift keying
CN109729034B (en) Timing method, device and terminal
Guan et al. Intelligent recognition of subcarrier for wireless link of satellite communication
CN111884675B (en) Method and system for tracking multi-system phase hopping spread spectrum modulation signal
RU2423798C1 (en) Clock synchronisation device
CN117590446B (en) Broadband composite navigation signal tracking method and device
JP3088906B2 (en) Synchronizer
JP2005201814A (en) Receiver of spread spectrum system and method of processing receiving signal
JP6269873B1 (en) Receiving apparatus and receiving method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant