CN108933750B - Distance measurement method and device based on joint estimation and recording medium - Google Patents

Distance measurement method and device based on joint estimation and recording medium Download PDF

Info

Publication number
CN108933750B
CN108933750B CN201810622433.9A CN201810622433A CN108933750B CN 108933750 B CN108933750 B CN 108933750B CN 201810622433 A CN201810622433 A CN 201810622433A CN 108933750 B CN108933750 B CN 108933750B
Authority
CN
China
Prior art keywords
ranging
absolute value
signal
estimation
ranging information
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.)
Active
Application number
CN201810622433.9A
Other languages
Chinese (zh)
Other versions
CN108933750A (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.)
Beijing Research Institute of Telemetry
Aerospace Long March Launch Vehicle Technology Co Ltd
Original Assignee
Beijing Research Institute of Telemetry
Aerospace Long March Launch Vehicle 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 Beijing Research Institute of Telemetry, Aerospace Long March Launch Vehicle Technology Co Ltd filed Critical Beijing Research Institute of Telemetry
Priority to CN201810622433.9A priority Critical patent/CN108933750B/en
Publication of CN108933750A publication Critical patent/CN108933750A/en
Application granted granted Critical
Publication of CN108933750B publication Critical patent/CN108933750B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2689—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation
    • H04L27/2695—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation with channel estimation, e.g. determination of delay spread, derivative or peak tracking
    • 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
    • G01S1/00—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith
    • G01S1/02—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith using radio waves
    • G01S1/08—Systems for determining direction or position line
    • G01S1/20—Systems for determining direction or position line using a comparison of transit time of synchronised signals transmitted from non-directional antennas or antenna systems spaced apart, i.e. path-difference systems
    • G01S1/22—Systems for determining direction or position line using a comparison of transit time of synchronised signals transmitted from non-directional antennas or antenna systems spaced apart, i.e. path-difference systems the synchronised signals being frequency modulations on carrier waves and the transit times being compared by measuring difference of instantaneous frequencies of received carrier waves
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2656—Frame synchronisation, e.g. packet synchronisation, time division duplex [TDD] switching point detection or subframe synchronisation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2657—Carrier synchronisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

A distance measurement method based on joint estimation comprises the following steps: step one, according to a received signal x (kT)s) Obtaining coarse ranging information estimates
Figure DDA0001698363800000011
Calculating the received signal x (kT)s) The sum y (Ln + l) of the real absolute value and the imaginary absolute value of (d); step two, the received signal x (kT) in the step one is compareds) The sum y (Ln + l) of the real part absolute value and the imaginary part absolute value is subjected to weighted accumulation processing, and the accurate ranging information estimation is calculated
Figure DDA0001698363800000012
Step three, estimating according to the coarse ranging information in the step one
Figure DDA0001698363800000013
And, fine ranging information estimation in step two
Figure DDA0001698363800000014
Calculating accurate ranging information estimate between two communication terminals
Figure DDA0001698363800000015
Wherein
Figure DDA0001698363800000016

Description

Distance measurement method and device based on joint estimation and recording medium
Technical Field
The invention relates to a distance measurement method, a distance measurement device and a recording medium based on joint estimation, and belongs to the technical field of digital wireless communication transmission.
Background
In a wireless communication system, because the attitude change between mobile communication terminals is often in a state that navigation information is weak or navigation signals are interfered, position information provided by satellite navigation is difficult to obtain, and distance information and position information between the terminals need to be obtained through a data link and a signal processing means. Especially in the unmanned aerial vehicle data link system, thereby the measurement and control system needs to acquire the position of unmanned aerial vehicle and adjust the parameter configuration and the link strategy of the information transmission system, promotes the survivability of the operation system.
The traditional pseudo-random code ranging technology has the characteristics of interference resistance, multipath fading resistance, good concealment and the like, but the ranging precision of the traditional pseudo-random code ranging technology depends on the code rate; the super-resolution algorithm based on spectrum estimation has higher precision but higher complexity; the distance measurement method adopting the delay locked loop mainly depends on the tracking precision of the delay locked loop; the phase ranging is a ranging method with excellent performance, but certain contradiction exists between the ranging precision and the ranging range, namely the ranging range is smaller when the ranging precision is higher; the ranging method based on received signal energy indication is the simplest to implement, but the least accurate.
PN sequence based methods are studied in "thank you.pn sequence based wireless ranging method study, university of vinpocetine, 2012"; in Yan Star facing the star, the handle of the Yangtze river, the big Dipper, et al. a distance measuring method of a multi-rate OFDM measurement and control system, telemetering and remote control, 2017,38(2):1-6 ", the author provides a distance measuring method of a combined time domain and a frequency domain of an Orthogonal Frequency Division Multiplexing (OFDM) technology;
a two-way ranging method which needs iterative equalization and decoding and performs zero calibration and calibration on a distance value is disclosed in a patent 'unmanned aerial vehicle measurement and control data chain incoherent ranging method' with a publication number of CN105391489A granted by China in 5, month and 4 days of 2018. A frequency domain nonlinear estimation algorithm insensitive to jitter of code loop tracking of spread spectrum codes is designed in ' Wangbuyun ', a frequency domain nonlinear estimation algorithm for pseudo code ranging, and radio engineering 2007 '. However, most of the above documents are ranging methods based on CDMA signals, OFDM signals, PN sequences, or two-way ranging and calibration methods, and the high-precision ranging technology for the linear modulation signals of the common data link system of the measurement and control system is less involved, or the ranging precision is not high.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the method comprises the steps of firstly estimating coarse ranging information by using shorter leader segment data, and then estimating fine ranging information by using data segment oversampling signals, thereby completing accurate ranging between two communication terminals in a data chain system.
The purpose of the invention is realized by the following technical scheme:
a distance measurement method based on joint estimation comprises the following steps:
step one, according to a received signal x (kT)s) Obtaining coarse ranging information estimates
Figure GDA0002861600670000021
Calculating the received signal x (kT)s) The sum y (Ln + l) of the real absolute value and the imaginary absolute value of (d);
step two, the received signal x (kT) in the step one is compareds) The sum y (Ln + l) of the real part absolute value and the imaginary part absolute value is subjected to weighted accumulation processing, and the accurate ranging information estimation is calculated
Figure GDA0002861600670000022
Step three, estimating according to the coarse ranging information in the step one
Figure GDA0002861600670000023
And, fine ranging information estimation in step two
Figure GDA0002861600670000024
Calculating accurate ranging information estimate between two communication terminals
Figure GDA0002861600670000025
Wherein
Figure GDA0002861600670000026
In the above distance measurement method based on joint estimation, the received signal x (kT) in step ones) To carry outTiming synchronization, carrier synchronization and frame synchronization, and then carrying out coarse ranging information estimation according to PN sequence ranging algorithm or CDMA pseudo code ranging algorithm to obtain
Figure GDA0002861600670000027
In the above distance measuring method based on joint estimation, in the first step, the signal x (kT) is receiveds) The sum y (Ln + l) of the absolute value of the real part and the absolute value of the imaginary part is:
when the transmitted signal is a binary phase shift keying BPSK signal
y(Ln+l)＝|xI(Ln+l)|+|xQ(Ln+l)|，1≤l≤L；
When the transmitted signal is a generally linear modulated signal
y(Ln+l)＝|x(Ln+l)|，1≤l≤L；
Wherein x isI(Ln + l) is a received signal x (kT)s) Real part of (x)Q(Ln + l) is a received signal x (kT)s) L is a sampling multiple; l is a first ordinal number, and l is a positive integer.
In the above distance measuring method based on joint estimation, the fine distance measuring information estimation in the second step
Figure GDA0002861600670000031
Comprises the following steps:
Figure GDA0002861600670000032
wherein c is the speed of light, T is the symbol period, arg { } is the phase information of complex signal, and L is the sampling multiple; n is the length of the data symbol of the sending end; l is a first ordinal number, l is a positive integer, and w (l) is a weighting factor.
The distance measurement method based on the joint estimation is applied to a data chain system.
A recording medium having a computer-readable program recorded thereon, the program, when executed by a processor, implementing the steps of:
step one, according to a received signal x (kT)s)x(kTs) Obtaining coarse ranging information estimates
Figure GDA0002861600670000033
Calculating the received signal x (kT)s) The sum y (Ln + l) of the real absolute value and the imaginary absolute value of (d);
step two, the received signal x (kT) in the step one is compareds) The sum y (Ln + l) of the real part absolute value and the imaginary part absolute value is subjected to weighted accumulation processing, and the accurate ranging information estimation is calculated
Figure GDA0002861600670000034
Step three, estimating according to the coarse ranging information in the step one
Figure GDA0002861600670000035
And, fine ranging information estimation in step two
Figure GDA0002861600670000036
Calculating accurate ranging information estimate between two communication terminals
Figure GDA0002861600670000037
Wherein
Figure GDA0002861600670000038
The recording medium, wherein the pair of received signals x (kT) in the first steps) Timing synchronization, carrier synchronization and frame synchronization are carried out, and then coarse ranging information estimation is carried out according to a PN sequence ranging algorithm or a CDMA pseudo code ranging algorithm to obtain
Figure GDA0002861600670000039
The recording medium, wherein in the first step, the signal x (kT) is receiveds) The sum y (Ln + l) of the absolute value of the real part and the absolute value of the imaginary part is:
when the transmitted signal is a binary phase shift keying BPSK signal
y(Ln+l)＝|xI(Ln+l)|+|xQ(Ln+l)|，1≤l≤L；
When the transmitted signal is a generally linear modulated signal
y(Ln+l)＝|x(Ln+l)|，1≤l≤L；
Wherein x isI(Ln + l) is a received signal x (kT)s) Real part of (x)Q(Ln + l) is a received signal x (kT)s) L is a sampling multiple; l is a first ordinal number, and l is a positive integer.
The recording medium, the fine ranging information estimation in the second step
Figure GDA0002861600670000041
Comprises the following steps:
Figure GDA0002861600670000042
wherein c is the speed of light, T is the symbol period, arg { } is the phase information of complex signal, and L is the sampling multiple; n is the length of the data symbol of the sending end; l is a first ordinal number, l is a positive integer, and w (l) is a weighting factor.
A distance measurement device based on joint estimation, which adopts the distance measurement method based on joint estimation as claimed in any one of claims 1-5.
Compared with the prior art, the invention has the following beneficial effects:
(1) the invention relates to a distance measurement method based on joint estimation in a data link system, which is suitable for general linear modulation signals and can more accurately improve the distance measurement accuracy of the data link system compared with the traditional PN sequence or pseudo code distance measurement method;
(2) according to the distance measurement method based on the joint estimation in the data link system, the precise distance measurement is carried out based on the sum of the absolute values of the over-sampling signals, so that the power consumption of the system can be more efficient due to lower complexity;
(3) the invention relates to a distance measurement method based on joint estimation in a data chain system, which is based on the estimation of the accurate distance measurement information of an over-sampling signal and combined with the estimation of the coarse distance measurement information based on a short leading segment, thereby improving the distance measurement action distance and reducing the distance measurement ambiguity.
Drawings
FIG. 1 is a block diagram of a method for joint estimation based ranging in a data link system according to the present invention;
fig. 2 shows the root mean square error performance of the range finding method based on joint estimation in the data link system according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be described in detail with reference to the accompanying drawings.
Fig. 1 is a block diagram of a ranging method based on joint estimation in a data link system according to the present invention. A distance measurement method based on joint estimation comprises the following steps:
step one, according to a received signal x (kT)s) Obtaining coarse ranging information estimates
Figure GDA0002861600670000051
Calculating the received signal x (kT)s) The sum y (Ln + l) of the real absolute value and the imaginary absolute value of (f).
Specifically, the sending end is at T0The transmitting signal is sent out at any moment, and the data frame sent by the sending end is as follows: s ═ a1,a2,...,ap,b1,b2,...,bdContains a leader { a }1,a2,...,apAnd data segment b2,...,bd}. Aircraft to be tested at T0+T1After receiving the signal at time T0+T1+T2Forwarding out all the time, the receiving end is at T0+2T1+T2And receiving the signals at the moment and performing parameter estimation, ranging and other processing on the signals. Received signal x (kT) received by receiving ends) Comprises the following steps:
Figure GDA0002861600670000052
wherein, TsFor sampling clocks, TsT/L; l is a sampling multiple; t is a symbol period; snFor data symbol with length N of sending end, including leading segment and dataA segment; g (t) is a filter; τ T is time delay, z (T) is Gaussian white noise; k is a second ordinal number, and k is a positive integer; n is a third ordinal number, and n is a positive integer.
The received signal may also be expressed as: x (kT)s) X (Ln +1), x (Ln +2),.., x (Ln + L) }. For received signal x (kT)s) Timing synchronization, carrier synchronization and frame synchronization are carried out, and then coarse ranging information estimation is carried out according to a PN sequence ranging algorithm or a CDMA pseudo code ranging algorithm to obtain
Figure GDA0002861600670000053
Receiving a signal x (kT)s) The sum y (Ln + l) of the absolute value of the real part and the absolute value of the imaginary part is:
when the transmitted signal is a binary phase shift keying BPSK signal
y(Ln+l)＝|xI(Ln+l)|+|xQ(Ln+l)|，1≤l≤L；
When the transmitted signal is a generally linear modulated signal
y(Ln+l)＝|x(Ln+l)|，1≤l≤L；
Wherein x isI(Ln + l) is a received signal x (kT)s) Real part of (x)Q(Ln + l) is a received signal x (kT)s) L is a sampling multiple; l is a first ordinal number, and l is a positive integer.
Step two, the received signal x (kT) in the step one is compareds) The sum y (Ln + l) of the real part absolute value and the imaginary part absolute value is subjected to weighted accumulation processing, and the accurate ranging information estimation is calculated
Figure GDA0002861600670000054
Fine ranging information estimation
Figure GDA0002861600670000055
Comprises the following steps:
Figure GDA0002861600670000061
where c is the speed of light, about 3X 108M/s; t is a symbol period, arg { } is phase information of complex signals, and L is a sampling multiple; n is the length of the data symbol of the sending end; l is a first ordinal number, l is a positive integer, w (l) is a weighting factor, w (l) ═ e-j2π(l-1)/L。
For the receiver with L4, the fine ranging information estimation is equivalent to:
Figure GDA0002861600670000062
step three, estimating according to the coarse ranging information in the step one
Figure GDA0002861600670000063
And, fine ranging information estimation in step two
Figure GDA0002861600670000064
Calculating accurate ranging information estimate between two communication terminals
Figure GDA0002861600670000065
Wherein
Figure GDA0002861600670000066
The ranging method based on joint estimation can be applied to a data chain system.
A recording medium having a computer-readable program recorded thereon, the program, when executed by a processor, implementing the steps of:
step one, according to a received signal x (kT)s) Obtaining coarse ranging information estimates
Figure GDA0002861600670000067
Calculating the received signal x (kT)s) The sum y (Ln + l) of the real absolute value and the imaginary absolute value of (f).
Specifically, the sending end is at T0The transmitting signal is sent out at any moment, and the data frame sent by the sending end is as follows: s ═ a1,a2,...,ap,b1,b2,...,bdContains a leader { a }1,a2,...,apAnd data segment b2,...,bd}. Aircraft to be tested at T0+T1After receiving the signal at time T0+T1+T2Forwarding out all the time, the receiving end is at T0+2T1+T2And receiving the signals at the moment and performing parameter estimation, ranging and other processing on the signals. Received signal x (kT) received by receiving ends) Comprises the following steps:
Figure GDA0002861600670000068
wherein, TsFor sampling clocks, TsT/L; l is a sampling multiple; t is a symbol period; snThe data symbol with the length of N at the sending end comprises a leading segment and a data segment; g (t) is a filter; τ T is time delay, z (T) is Gaussian white noise; k is a second ordinal number, and k is a positive integer; n is a third ordinal number, and n is a positive integer.
The received signal may also be expressed as: x (kT)s) X (Ln +1), x (Ln +2),.., x (Ln + L) }. For received signal x (kT)s) Timing synchronization, carrier synchronization and frame synchronization are carried out, and then coarse ranging information estimation is carried out according to a PN sequence ranging algorithm or a CDMA pseudo code ranging algorithm to obtain
Figure GDA0002861600670000071
Receiving a signal x (kT)s) The sum y (Ln + l) of the absolute value of the real part and the absolute value of the imaginary part is:
when the transmitted signal is a binary phase shift keying BPSK signal
y(Ln+l)＝|xI(Ln+l)|+|xQ(Ln+l)|，1≤l≤L；
When the transmitted signal is a generally linear modulated signal
y(Ln+l)＝|x(Ln+l)|，1≤l≤L；
Wherein x isI(Ln + l) is a received signal x (kT)s) Real part of (x)Q(Ln + l) is a received signal x (kT)s) L is a sampling multiple; l is the firstOrdinal number, l is a positive integer.
Step two, the received signal x (kT) in the step one is compareds) The sum y (Ln + l) of the real part absolute value and the imaginary part absolute value is subjected to weighted accumulation processing, and the accurate ranging information estimation is calculated
Figure GDA0002861600670000072
Fine ranging information estimation
Figure GDA0002861600670000073
Comprises the following steps:
Figure GDA0002861600670000074
where c is the speed of light, about 3X 108M/s; t is a symbol period, arg { } is phase information of complex signals, and L is a sampling multiple; n is the length of the data symbol of the sending end; l is a first ordinal number, l is a positive integer, w (l) is a weighting factor, w (l) ═ e-j2π(l-1)/L。
For the receiver with L4, the fine ranging information estimation is equivalent to:
Figure GDA0002861600670000075
step three, estimating according to the coarse ranging information in the step one
Figure GDA0002861600670000076
And, fine ranging information estimation in step two
Figure GDA0002861600670000077
Calculating accurate ranging information estimate between two communication terminals
Figure GDA0002861600670000078
Wherein
Figure GDA0002861600670000079
A distance measuring device based on joint estimation is realized by adopting a distance measuring method based on joint estimation.
Example 1:
comparing the performance simulation result of the ranging method of the present invention with the theoretical performance limit of ranging, fig. 2 is a diagram illustrating the root mean square error performance of the ranging error of the ranging method based on joint estimation in a data link system of the present invention. BPSK signal is adopted, oversampling L is 4 times, roll-off coefficient is alpha 0.1, 0.35, oversampling signal N is 1000, accuracy of coarse ranging method is 0.1 mus x 3 x 10 at 10MHz symbol rate8And m/s is 30 m. As a result of analysis, when α is 0.1, the SNR is 10dB, 20dB, and 30dB, and the mean square error of the ranging error is 0.26m, 0.07m, and 0.02 m; when α is 0.35, the SNR is 10dB, 20dB, and 30dB, and the mean square error of the ranging error is 0.19m, 0.05m, and 0.01 m. SNR at signal-to-noise ratio>When the signal rate is 10dB, the mean square error of the ranging error of the algorithm can approach the theoretical performance limit, and the ranging performance limited by the signal rate of the traditional coarse ranging algorithm is greatly improved.
The above description is only for the best mode of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention.
Those skilled in the art will appreciate that those matters not described in detail in the present specification are well known in the art.

Claims (5)

1. A distance measurement method based on joint estimation is characterized in that: the method comprises the following steps:
step one, according to a received signal x (kT)s) Obtaining coarse ranging information estimates
Figure FDA0002861600660000011
Calculating the received signal x (kT)s) The sum y (Ln + l) of the real absolute value and the imaginary absolute value of (d);
step two, the received signal x (kT) in the step one is compareds) Absolute value of real part ofThe sum y (Ln + l) of the absolute values of the imaginary parts is subjected to weighted accumulation processing to calculate the estimation of the accurate ranging information
Figure FDA0002861600660000012
Step three, estimating according to the coarse ranging information in the step one
Figure FDA0002861600660000013
And, fine ranging information estimation in step two
Figure FDA0002861600660000014
Calculating accurate ranging information estimate between two communication terminals
Figure FDA0002861600660000015
Wherein
Figure FDA0002861600660000016
The pair of received signals x (kT) in the first steps) Timing synchronization, carrier synchronization and frame synchronization are carried out, and then coarse ranging information estimation is carried out according to a PN sequence ranging algorithm or a CDMA pseudo code ranging algorithm to obtain
Figure FDA0002861600660000017
The second step of fine ranging information estimation
Figure FDA0002861600660000018
Comprises the following steps:
Figure FDA0002861600660000019
wherein c is the speed of light, T is the symbol period, arg { } is the phase information of complex signal, and L is the sampling multiple; n is the length of the data symbol of the sending end; l is a first ordinal number, l is a positive integer, and w (l) is a weighting factor.
2. A method for ranging based on joint estimation as claimed in claim 1, wherein: receiving a signal x (kT) in the first steps) The sum y (Ln + l) of the absolute value of the real part and the absolute value of the imaginary part is:
when the transmitted signal is a binary phase shift keying BPSK signal
y(Ln+l)＝|xI(Ln+l)|+|xQ(Ln+l)|，1≤l≤L；
When the transmitted signal is a generally linear modulated signal
y(Ln+l)＝|x(Ln+l)|，1≤l≤L；
Wherein x isI(Ln + l) is a received signal x (kT)s) Real part of (x)Q(Ln + l) is a received signal x (kT)s) The imaginary part of (c).
3. A method for ranging based on joint estimation as claimed in claim 1, wherein: the ranging method is applied to a data chain system.
4. A computer-readable storage medium characterized by: the computer-readable storage medium has a computer-readable program recorded thereon, which when executed by a processor implements the steps of:
step one, according to a received signal x (kT)s)x(kTs) Obtaining coarse ranging information estimates
Figure FDA0002861600660000021
Calculating the received signal x (kT)s) The sum y (Ln + l) of the real absolute value and the imaginary absolute value of (d);
step two, carrying out weighted accumulation processing on the sum y (Ln + l) of the real part absolute value and the imaginary part absolute value of the received signal x (kT) in the step one, and calculating the estimation of the accurate ranging information
Figure FDA0002861600660000022
Step three, estimating according to the coarse ranging information in the step one
Figure FDA0002861600660000023
And, fine ranging information estimation in step two
Figure FDA0002861600660000024
Calculating accurate ranging information estimate between two communication terminals
Figure FDA0002861600660000025
Wherein
Figure FDA0002861600660000026
The pair of received signals x (kT) in the first steps) Timing synchronization, carrier synchronization and frame synchronization are carried out, and then coarse ranging information estimation is carried out according to a PN sequence ranging algorithm or a CDMA pseudo code ranging algorithm to obtain
Figure FDA0002861600660000027
The second step of fine ranging information estimation
Figure FDA0002861600660000028
Comprises the following steps:
Figure FDA0002861600660000029
wherein c is the speed of light, T is the symbol period, arg { } is the phase information of complex signal, and L is the sampling multiple; n is the length of the data symbol of the sending end; l is a first ordinal number, l is a positive integer, and w (l) is a weighting factor.
5. The computer-readable storage medium of claim 4, wherein: receiving a signal x (kT) in the first steps) The sum y (Ln + l) of the absolute value of the real part and the absolute value of the imaginary part is:
when the transmitted signal is a binary phase shift keying BPSK signal
y(Ln+l)＝|xI(Ln+l)|+|xQ(Ln+l)|，1≤l≤L；
When the transmitted signal is a generally linear modulated signal
y(Ln+l)＝|x(Ln+l)|，1≤l≤L；
Wherein x isI(Ln + l) is a received signal x (kT)s) Real part of (x)Q(Ln + l) is a received signal x (kT)s) The imaginary part of (c).
CN201810622433.9A 2018-06-15 2018-06-15 Distance measurement method and device based on joint estimation and recording medium Active CN108933750B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810622433.9A CN108933750B (en) 2018-06-15 2018-06-15 Distance measurement method and device based on joint estimation and recording medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810622433.9A CN108933750B (en) 2018-06-15 2018-06-15 Distance measurement method and device based on joint estimation and recording medium

Publications (2)

Publication Number Publication Date
CN108933750A CN108933750A (en) 2018-12-04
CN108933750B true CN108933750B (en) 2021-03-26

Family

ID=64446637

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810622433.9A Active CN108933750B (en) 2018-06-15 2018-06-15 Distance measurement method and device based on joint estimation and recording medium

Country Status (1)

Country Link
CN (1) CN108933750B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111884972A (en) * 2020-06-29 2020-11-03 西南电子技术研究所（中国电子科技集团公司第十研究所） Bidirectional high-precision distance measuring system of OFDM communication system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7511667B2 (en) * 2006-02-22 2009-03-31 Novariant, Inc. Precise local positioning systems using ground-based transmitters
WO2014022239A1 (en) * 2012-07-29 2014-02-06 Qualcomm Incorporated Anatomical gestures detection system using radio signals
US9059828B1 (en) * 2012-05-23 2015-06-16 Jose Marcelo Lima Duarte Full search MIMO detector for recovering single or multiple data stream in a multiple antenna receiver
CN106772511A (en) * 2017-01-19 2017-05-31 湖南北云科技有限公司 A kind of distributed cloud RTK localization methods and system
CN107026810A (en) * 2017-03-17 2017-08-08 西南电子技术研究所（中国电子科技集团公司第十研究所） The PN synchronization method of burst directly-enlarging system and its DS waveform that happens suddenly

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7511667B2 (en) * 2006-02-22 2009-03-31 Novariant, Inc. Precise local positioning systems using ground-based transmitters
US9059828B1 (en) * 2012-05-23 2015-06-16 Jose Marcelo Lima Duarte Full search MIMO detector for recovering single or multiple data stream in a multiple antenna receiver
WO2014022239A1 (en) * 2012-07-29 2014-02-06 Qualcomm Incorporated Anatomical gestures detection system using radio signals
CN106772511A (en) * 2017-01-19 2017-05-31 湖南北云科技有限公司 A kind of distributed cloud RTK localization methods and system
CN107026810A (en) * 2017-03-17 2017-08-08 西南电子技术研究所（中国电子科技集团公司第十研究所） The PN synchronization method of burst directly-enlarging system and its DS waveform that happens suddenly

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"一种无线数据链测距技术的研究与实现";王满 宋鹏齐 建中;《测控技术》;20131231;全文 *

Also Published As

Publication number Publication date
CN108933750A (en) 2018-12-04

Similar Documents

Publication Publication Date Title
CN105162493A (en) Doppler domain and delay domain two-dimension acquiring method and device
Wang et al. Doppler estimation and timing synchronization of underwater acoustic communication based on hyperbolic frequency modulation signal
CN103428153B (en) GMSK signal acceptance method in satellite mobile communication
CN107493117B (en) The two-dimentional joint acquisition method of direct expansion msk signal under a kind of high dynamic
CN101141425A (en) Time-division pilot based channel estimation method of mobile communication system
US6505053B1 (en) Method for sinusoidal modeling and prediction of fast fading processes
CN103023831B (en) A kind of carrier frequency bias estimation being applicable to burst waveform
CN102035770B (en) Method for estimating channel by means of correlation
CN104219761A (en) Ultra-wideband wireless positioning method based on maximum slope
CN108933750B (en) Distance measurement method and device based on joint estimation and recording medium
US20040179631A1 (en) UWB receiver architecture
CN1972137A (en) A receiving signal frequency offset estimation and compensation method for WCDMA
CN112953593A (en) LoRa advanced receiver
CN117040564A (en) Communication and distance measurement integrated method based on template reconstruction
CN111654308A (en) Precision frequency offset estimation method for burst spread spectrum weak signal
CN115242367B (en) Data error correction method for industrial wireless channel impulse response
CN108401582B (en) A kind of GEO satellite mobile communication system initial frequency synchronization method adapting to TDMA/FDMA systems
CN104779993A (en) Deep space measurement and control system and method on basis of frequency-domain equalization
CN112630761B (en) SC-FDE signal joint ranging method in unmanned aerial vehicle bee colony network
CN113702901B (en) Cross-correlation time delay estimation method based on signal correlation enhancement
CN109547129B (en) Stepping weighting Tong detection method
CN108924073B (en) A kind of quick self-adapted Doppler estimation synchronous based on pseudo-random sequence
CN114500192B (en) Channel estimation method and terminal equipment
CN104219762A (en) Ultra-wideband wireless positioning method based on mean square error
CN115250125B (en) RAKE receiver and improved receiving method for satellite communication system

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