CN112305498A - Heterogeneous TDOA (time difference of arrival) positioning system - Google Patents
Heterogeneous TDOA (time difference of arrival) positioning system Download PDFInfo
- Publication number
- CN112305498A CN112305498A CN202011237117.3A CN202011237117A CN112305498A CN 112305498 A CN112305498 A CN 112305498A CN 202011237117 A CN202011237117 A CN 202011237117A CN 112305498 A CN112305498 A CN 112305498A
- Authority
- CN
- China
- Prior art keywords
- station
- tdoa
- time stamp
- calibration
- mode
- 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.)
- Pending
Links
- 238000005070 sampling Methods 0.000 claims abstract description 79
- 238000005259 measurement Methods 0.000 claims abstract description 31
- 238000004088 simulation Methods 0.000 claims abstract description 25
- 238000012360 testing method Methods 0.000 claims description 32
- 238000006243 chemical reaction Methods 0.000 claims description 22
- 230000000153 supplemental effect Effects 0.000 claims description 18
- 238000000034 method Methods 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims 1
- 238000012952 Resampling Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 101000798707 Homo sapiens Transmembrane protease serine 13 Proteins 0.000 description 1
- 102100032467 Transmembrane protease serine 13 Human genes 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/0205—Details
- G01S5/021—Calibration, monitoring or correction
-
- 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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W64/00—Locating users or terminals or network equipment for network management purposes, e.g. mobility management
- H04W64/003—Locating users or terminals or network equipment for network management purposes, e.g. mobility management locating network equipment
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention provides a heterogeneous TDOA positioning system, which converts sampling rates of all stations into consistent sampling rates, acquires a sampling time stamp calibration value between the heterogeneous stations through a measurement task, further acquires a supplementary time stamp calibration value and verifies the positioning precision through the simulation positioning of equipment in a loop, and corrects the time stamp offset between the heterogeneous stations through the sampling time stamp calibration value and the supplementary sampling time stamp calibration value in an actual positioning task, so that the heterogeneous TDOA positioning system consisting of different types of stations can cooperatively position a radio frequency signal source.
Description
Technical Field
The invention belongs to the field of radio communication, and particularly relates to a heterogeneous TDOA (time difference of arrival) positioning technology.
Background
The TDOA location system consists of a central station, a server, and 3 or more stations with a distance between several kilometers and hundreds of kilometers. The central station is deployed in a computer and mainly used for task control and result display, the server is deployed in a data processing computer and mainly used for calculating sensitive data processing tasks, the central station and the server of some systems are deployed in the same computer, and the measuring station consists of a digital radio frequency receiver, a satellite positioning/timing module, a computer and measuring station application software.
The measuring stations of the TDOA system collect IQ data blocks of measured signals at the same time and upload the IQ data blocks to the server, the server calculates the TDOA between the measuring station pairs, and the signal source coordinates are calculated through a plurality of TDOA measured values and measuring station coordinates. The synchronization precision of the time stamp of the IQ data of the measuring station is the primary factor influencing the TDOA result.
Because different clock frequencies are used by radio frequency receivers of different types of observation stations, the sampling rates which can be realized by the receivers are different; on the other hand, although the 1PPS pulse synchronization precision between the measuring stations can reach within tens of nanoseconds through satellite time service, due to the fact that the digital signal processing module structures and the time stamping mechanism implementation processes at the rear ends of the ADCs of the radio frequency receivers of different models are different, even under the condition that the sampling rates are the same, the time stamps are also different.
In view of this, current TDOA positioning systems require strict consistency in hardware and software specifications for all stations. Existing radio spectrum monitoring systems are equipped with various different models of fixed, vehicle-mounted and airborne stations, and the station resources available to TDOA positioning systems are limited by model consistency due to the problem of timestamp differences. In order to implement the heterogeneous TDOA positioning system, i.e., to enable heterogeneous stations of different models to cooperatively perform a positioning task, problems in the measurement and verification of the timestamp calibration value between the heterogeneous stations, the correction of the sampling timestamp in the actual positioning task, and the like need to be solved.
Disclosure of Invention
In order to solve the above problems, the present invention provides a heterogeneous TDOA positioning system, which converts sampling rates of various stations into uniform sampling rates, obtains a sampling timestamp calibration value between the heterogeneous stations by measurement, further obtains a supplemental timestamp calibration value by simulation positioning of a device in a loop and verifies positioning accuracy, and corrects a timestamp offset between the heterogeneous stations by the sampling timestamp calibration value and the supplemental sampling timestamp calibration value in an actual positioning task, so that the heterogeneous TDOA positioning system composed of different types of stations can cooperatively position a radio frequency signal source location.
In order to achieve the purpose, the invention adopts the following technical scheme:
a heterogeneous TDOA location system, comprising:
s1, a central station task control module: providing a setting interface, setting task type MmodeAnd an operation mode RmodeAnd task parameters for issuing task commands to the survey station and the server, wherein the task type comprises Mmode=1 time stamp calibration value measurement task, Mmode= precision test task, Mmode=3 actual positioning task, operating mode RmodeThe =0 expression sampling rate conversion and delay adjustment module is executed in the observation station, otherwise, the sampling rate conversion and delay adjustment module is executed in the server;
s2, a station testing command receiving module: receiving a central station command, setting working parameters and states of a radio frequency receiver, and starting a station testing task;
s3, the station testing signal sampling receiving module: periodically extracting blocks of samples S from a receiver output IQ streamIQAnd packaged as a timed IQ block STIQWhen R ismodeSample rate conversion and delay adjustment module processing timing IQ block S when =0TIQWill time the IQ block STIQUploading to a server through a network;
s4, a sampling rate conversion and time delay adjustment module: timing IQ block STIQSampling rate conversion and sampling time delay adjustment;
s5, a server TDOA calculation module: receive timestamp alignment NtaskThe measurement stations are timed IQ blocks, and are combined into a measurement station pair when R ismodeWhen not equal to 0, a sampling rate conversion and time delay adjustment module is called to process each timing IQ block STIQComputing TDOA measured value T of the station pairmeasRecording a time stamp calibration value or correcting the TDOA measurement value by using a supplementary time stamp calibration value according to the task type;
s6, a server positioning resolving module: calculating the coordinates of the emission source as O according to the TDOA measured value and the coordinates of the measuring stations of the measuring station pairLOCWhen M ismodeWhen the time difference is not less than 2, recording TDOA errors and positioning errors and counting the supplementary timestamp calibration values;
s7, a central station display module: TDOA hyperbolic measurements and truth, location result measurements and truth are displayed on the map, and TDOA error, timestamp calibration value, supplemental timestamp calibration value, and location error are displayed in a table.
Further, the task parameters of S1 include:
s11, assigning a serial number to each station, establishing a relation table of station serial numbers and station model numbers, and designating the model M0For time-stamped reference stations, model number not M0Namely the heterogeneous station testing, the calibration value of the time stamp of the station testing is stored in the quintupleIn a constructed alignment queue, wherein M1nd、M2ndTwo stations type, FrateIs the sampling rate FS,TcaliFor the time stamp calibration value it is,calibrating values for the supplemental timestamps;
s12, the command parameter includes signal frequency FCSignal bandwidth BWSampling rate FSData acquisition period TmarkThe number of samples L of the timing IQ blockblockExtended sampling number LextTime stamp calibration valueSimulated emission point coordinates S0And simulation station coordinates SkAnd the number of stations N required for positioningtaskWherein K is the number of the stations, K belongs to {1,2, ⋯, K }, and K is the number of stations in the system;
s13, when MmodeTime stamp calibration value for all stations at time of =1When M ismodeWhen not equal to 1, if the calibration queue is empty, reading the record in the time stamp calibration file to the calibration queue, and if the model is MnThe kth station of (1), when Mn=M0Device for placingOtherwise, searching the alignment queue for meeting condition M1nd=M0And M2nd=MnAnd Frate=FSTo obtain a timestamp calibration value。
Further, the step of periodically extracting S3 sample blocks S from the receiver output IQ streamIQAnd packaged as a timed IQ block STIQThe method comprises the following specific steps:
s31, obtaining the actual sampling rate of the receiverFrom a reference time TREFStart with TmarkIs the time point T of the cyclecapThe time for synchronously extracting the sampling blocks of each testing station is used for calculating TcapL=Tcap-Lext/(2×FS) Searching IQ sample stream output from receiver for time stamp TcapLFrom which a sample point of (L) is extractedext+Lblock)×(/FS) Is sampled set SIQWherein T isREFIs a constant;
s32, numbering stations and actual sampling rateTime stamp TcapStation coordinates and sample set SIQPackaged as a timed IQ block STIQWhen M ismodeAnd if not, the coordinate of the measuring station is the simulation measuring station coordinate in the command, otherwise, the coordinate is the coordinate output by the satellite positioning and time service module.
Further, the timing IQ block S of S4TIQThe sampling rate conversion and sampling time delay adjustment method specifically comprises the following steps:
s41 opposite timing IQ block STIQSample execution in (1)To FSSampling rate ofConversion so that each station-timing IQ block has Lext+LblockA sample, at this time, sample STIQ (LextTime stamp of/2) is Tcap;
S42, for the station with the number k, when M ismodeTime P of =3k=0, otherwise, the simulation transmitting point coordinates S are calculated0And simulation survey station coordinate SkA distance D betweenkCalculating DkCorresponding number of samples Pk= ⌊FS×DkC + 0.5 ⌋, wherein C is the electromagnetic wave velocity constant;
s43, timing the samples S in IQ blockTIQ (Lext/2-Pk+) Extracting L as a starting pointblockOne sample as a timed IQ block STIQ。
Further, the specific steps of recording the timestamp calibration value or correcting the TDOA measurement value by using the supplemental timestamp calibration value at S5 are as follows:
s51, when MmodeIf =1, the emission point S is determined according to the simulation0Calculating TDOA truth value T by coordinate distance difference of two simulation test stationstrueCalculating the TDOA error Tdev=Tmeas-TtrueWhen the model of the first station in the station pair is M0Second station non-M0Time, TDOA error TdevI.e. the timestamp calibration value between stations, and the model M of the first station in the pair of stations1The model M of the second testing station2Sampling rate FSAnd timestamp calibration value TdevIs a five-membered groupWill beStoring the data into a calibration queue;
s52, when Mmode=1, at the end of the task, calculating the timestamp calibration of each heterogeneous station at different sampling rates according to the nodes cached in the calibration queueMean value, in 5-tuple QcaliThe format is recorded in a time stamp calibration file, and the time stamp calibration value is supplemented;
S53, when MmodeNot equal to 1, according to the model of the first measuring station, the model of the second measuring station and the sampling rate FSSearching matched nodes from the calibration queue to obtain the supplementary timestamp calibration valueCorrection of TDOA measured value Tmeas= Tmeas-。
Further, the specific steps of recording the TDOA error and the positioning error and counting the supplemental timestamp calibration value in S6 are as follows:
s61, calculating OLOCAnd S0The distance is the positioning error, and the sampling rate F is recordedSThe parameters such as the model, the error and the positioning error of each TDOA measuring station are put into an error analysis file;
s62, when the task is finished, the TDOA error mean value e of each station pair is countedtdoaAnd the mean value of the positioning error and the supplementary time stamp calibration value of each node in the time stamp calibration fileTogether, write the last row of the error analysis file, model and sampling rate F in the time stamp calibration file and the measuring station pairSSupplemental timestamp calibration values for matching nodesUpdating supplemental timestamp calibration values in calibration file nodes=etdoa。
The invention has the beneficial effects that:
1) the time stamp calibration value between the heterogeneous measurement stations is obtained through measurement, the supplementary time stamp calibration value is further obtained through a simulation positioning process, and the time stamp offset between the heterogeneous measurement stations after sampling rate conversion can be fully calibrated through a two-stage time stamp calibration mechanism, so that the positioning precision close to that of a homogeneous TDOA positioning system is realized. If the supplementary timestamp calibration value is not used, the timestamp precision can only be calibrated to the order of magnitude equal to the sampling period, and when the sampling rate is low, the calibrated timestamp synchronization precision is far lower than the PPS pulse synchronization precision of the satellite time service module 1;
2) due to the fact that sampling rate conversion is involved, time stamp calibration is a matter sensitive to the implementation process, and the processes of time stamp calibration value measurement, supplementary time stamp calibration value acquisition, actual TDOA positioning task execution and the like are all completed in the system, consistency of time stamp calibration value acquisition and the application process is guaranteed, and a verifiable and trustable technical method is provided for a heterogeneous TDOA positioning system.
Drawings
FIG. 1 is a block diagram of a heterogeneous TDOA system;
the hardware connection relationship when performing timestamp calibration value measurement and simulation positioning in the embodiment of fig. 2.
Detailed Description
The following describes in detail embodiments of the system according to the present invention with reference to examples.
In the embodiment, the hardware of the heterogeneous TDOA system consists of a central station and 3 measuring stations, a server and the central station are arranged in the same computer, the models of two measuring station radio frequency receivers are BB60C, the model of one measuring station radio frequency receiver is USRP B210, each measuring station further comprises a satellite positioning and time service module and one measuring station computer, the radio frequency receivers receive 1PPS (pulse per second) output by the time service module and serve as a clock tame source, the radio frequency receivers transmit data to the measuring station computers through USB ports, and the heterogeneous TDOA system is used for executing a two-dimensional positioning task and positioning the required measuring stationNumber of stations Ntask=3, set BB60C as time stamp reference station, set running mode Rmode=1, i.e. the sample rate conversion and delay adjustment module, is run in the server;
in order to realize the heterogeneous TDOA positioning, the system firstly executes a time stamp calibration value measurement task, secondly updates the supplementary time stamp calibration value and verifies the positioning precision through simulation positioning, and thirdly executes an actual positioning task;
when time stamp calibration value measurement and simulation positioning are carried out, a central station and three measuring stations are placed at the same place and connected through a local area network, radio frequency input terminals of receivers of the three measuring stations are connected to a digital signal generator through a power divider, the connection relation between the measuring stations and the central station is shown in figure 2, and the symbol rate generated by the digital signal generator is FSThe signal frequency point of the QPSK signal of/2 is fixed to be 1GHz, the output level power is set to be-60 dBm, and 5 sampling rates F are respectively subjected to the sampling rate which can be realized according to BB60CS= (20, 10, 5, 2.5, 1.25) make measurements, sample rate unit MSPS, to obtain timestamp calibration values and supplemental timestamp calibration values at different sample rates;
when an actual positioning task is executed, three measuring stations are deployed in an actual environment, the measuring stations and a central station are interconnected through a public mobile network, a radio frequency receiver of the measuring station receives an air signal through an antenna, a network card of the central station uses an internet static IP address, one of 5 sampling rate measured values is selected when a sampling rate is set, and a timestamp calibration value and a supplementary timestamp calibration value corresponding to the sampling rate can be read from a timestamp calibration file.
A heterogeneous TDOA location system, comprising:
s01, a central station task control module: providing a setting interface, setting task type MmodeAnd an operation mode RmodeAnd task parameters for issuing task commands to the survey station and the server, wherein the task type comprises Mmode=1 time stamp calibration value measurement task, Mmode= precision test task, Mmode=3 actual positioning task, operating mode Rmode=1 expression sample rate conversion and delay adjustment module is executed in server;
s02, a station testing command receiving module: receiving the command of the central station, setting the working parameters and states of the radio frequency receiver, setting the bandwidth of the receiver to be constant F in order to avoid the time stamp change caused by the bandwidth of the digital filter and measure the signal carrier-to-noise ratioSStarting a station testing task;
s03, the station testing signal sampling receiving module: periodically extracting blocks of samples S from a receiver output IQ streamIQAnd packaged as a timed IQ block STIQWill time the IQ block STIQUploading to a server through a network;
s04, a sampling rate conversion and time delay adjustment module: timing IQ block STIQSampling rate conversion and sampling time delay adjustment;
s05, a server TDOA calculation module: receive timestamp alignment NtaskEach test station timing IQ block combines the test stations into a test station pair, preferentially uses a time stamp reference test station as a first test station, and calls a sampling rate conversion and time delay adjustment module to process each timing IQ block STIQUsing a cut-off frequency of BWPerforms filtering on the timed IQ block to calculate the TDOA measurement T of the station pairmeasRecording a time stamp calibration value or correcting the TDOA measurement value by using a supplementary time stamp calibration value according to the task type;
s06, a server positioning resolving module: calculating the coordinates of the emission source as O according to the TDOA measured value and the coordinates of the measuring stations of the measuring station pairLOCWhen M ismodeWhen the time difference is not less than 2, recording TDOA errors and positioning errors and counting the supplementary timestamp calibration values;
s07, a central station display module: TDOA hyperbolic measurements and truth, location result measurements and truth are displayed on the map, and TDOA error, timestamp calibration value, supplemental timestamp calibration value, and location error are displayed in a table.
Further, the task parameters of S01 include:
s011, endowing each survey station with a serial number, establishing a relation table of the serial numbers of the survey stations and the model numbers of the survey stations, and designating the model M0For time-stamped reference stations, model number not M0Namely the heterogeneous station testing, the calibration value of the time stamp of the station testing is stored in the quintupleIn a constructed alignment queue, wherein M1nd、M2ndTwo stations type, FrateIs the sampling rate FS,TcaliFor the time stamp calibration value it is,calibrating values for the supplemental timestamps;
s012, the command parameters include: frequency of signal FC,MmodeTiming signal bandwidth B not equal to 3W=FS/2, sampling rate FSData acquisition period Tmark=500ms, number of timed IQ block samples Lblock=16384, set number of extended samples L based on maximum inter-station distance and potential timestamp offsetext=512, time stamp calibration valueSimulation of the coordinates S of the launch point0And simulation station coordinates SkSetting the number of stations N required for positioningtask=3, where K is the number of stations, K belongs to {1,2, ⋯, K }, and K is the number of stations in the system;
s013, when MmodeTime stamp calibration value for all stations at time of =1When M ismodeWhen not equal to 1, if the calibration queue is empty, reading the record in the time stamp calibration file to the calibration queue, and if the model is MnThe kth station of (1), when Mn=M0Device for placingOtherwise, searching the alignment queue for meeting condition M1nd=M0And M2nd=MnAnd Frate=FSTo obtain a timestamp calibration value。
Further, the step of periodically extracting S03 sample blocks S from the receiver output IQ streamIQAnd packaged as a timed IQ block STIQThe method comprises the following specific steps:
s031, obtain the actual sampling rate of the receiverSetting a reference time constant TREFFor a whole second, starting with every whole second and starting with TmarkIs the time point T of the cyclecapThe time for synchronously extracting the sampling blocks of each testing station is used for calculating TcapL=Tcap-Lext/(2×FS) Searching IQ sample stream output from receiver for time stamp TcapLFrom which a sample point of (L) is extractedext+Lblock)×(/FS) Is sampled set SIQ;
S032, numbering stations and actual sampling rateTime stamp TcapStation coordinates and sample set SIQPackaged as a timed IQ block STIQWhen M ismodeAnd if not, the coordinate of the measuring station is the simulation measuring station coordinate in the command, otherwise, the coordinate is the coordinate output by the satellite positioning and time service module.
Further, the timing IQ block S of S04TIQThe sampling rate conversion and sampling time delay adjustment method specifically comprises the following steps:
s041, pair timing IQ block STIQSample execution in (1)To FSSo that each station-timing IQ block has Lext+LblockA sample, at this time, sample STIQ (LextTime stamp of/2) is TcapWhen is coming into contact withTo FSWhen the transformation can be realized by interpolation and integer extraction which are not more than 20 times, a fractional resampling algorithm is used to reduce the overhead, otherwise, a multi-stage arbitrary resampling algorithm is used;
s042, for the survey station with the number of k, when MmodeTime P of =3k=0, otherwise, the simulation transmitting point coordinates S are calculated0And simulation survey station coordinate SkA distance D betweenkCalculating DkCorresponding number of samples Pk= ⌊FS×DkC + 0.5 ⌋, wherein C is the electromagnetic wave velocity constant;
s043, to time samples S in IQ blocksTIQ (Lext/2-Pk+) Extracting L as a starting pointblockOne sample as a timed IQ block STIQ。
Further, the specific steps of recording the timestamp calibration value or correcting the TDOA measurement value by using the supplemental timestamp calibration value at S05 are as follows:
s051, when MmodeIf =1, the emission point S is determined according to the simulation0Calculating TDOA truth value T by coordinate distance difference of two simulation test stationstrueCalculating the TDOA error Tdev=Tmeas-TtrueWhen the model of the first station in the station pair is M0Second station non-M0Time, TDOA error TdevI.e. the timestamp calibration value between stations, and the model M of the first station in the pair of stations1The model M of the second testing station2Sampling rate FSAnd timestamp calibration value TdevIs a five-membered groupWill beStoring the data into a calibration queue;
s052, Mmode=1, root at task endCalculating the mean value of the time stamp calibration values of each heterogeneous station at different sampling rates according to the cached nodes in the calibration queue, and calculating the mean value by using 5-tuple QcaliThe format is recorded in a time stamp calibration file, and the time stamp calibration value is supplemented;
S053, when MmodeNot equal to 1, according to the model of the first measuring station, the model of the second measuring station and the sampling rate FSSearching matched nodes from the calibration queue to obtain the supplementary timestamp calibration valueCorrection of TDOA measured value Tmeas= Tmeas-。
Further, the specific steps of recording the TDOA error and the positioning error and counting the supplemental timestamp calibration value in S06 are as follows:
s061, calculating OLOCAnd S0The distance is the positioning error, and the sampling rate F is recordedSThe parameters such as the model, the error and the positioning error of each TDOA measuring station are put into an error analysis file;
s062, when the task is finished, counting the TDOA error mean value e of each station measuring pairtdoaAnd the mean value of the positioning error and the supplementary time stamp calibration value of each node in the time stamp calibration fileTogether, write the last row of the error analysis file, model and sampling rate F in the time stamp calibration file and the measuring station pairSSupplemental timestamp calibration values for matching nodesUpdating supplemental timestamp calibration values in calibration file nodes=etdoa。
The above description is only a preferred embodiment of the present application and is not intended to limit the present application, and various modifications and changes may be made by those skilled in the art without departing from the spirit and principle of the present application, and any modifications, equivalents, improvements, etc. made therein are intended to be included within the scope of the present application. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (6)
1. A heterogeneous TDOA location system, said system comprising:
s1, a central station task control module: providing a setting interface, setting task type MmodeAnd an operation mode RmodeAnd task parameters for issuing task commands to the survey station and the server, wherein the task type comprises Mmode=1 time stamp calibration value measurement task, Mmode= precision test task, Mmode=3 actual positioning task, operating mode RmodeThe =0 expression sampling rate conversion and delay adjustment module is executed in the observation station, otherwise, the sampling rate conversion and delay adjustment module is executed in the server;
s2, a station testing command receiving module: receiving a central station command, setting working parameters and states of a radio frequency receiver, and starting a station testing task;
s3, the station testing signal sampling receiving module: periodically extracting blocks of samples S from a receiver output IQ streamIQAnd packaged as a timed IQ block STIQWhen R ismodeSample rate conversion and delay adjustment module processing timing IQ block S when =0TIQWill time the IQ block STIQUploading to a server through a network;
s4, a sampling rate conversion and time delay adjustment module: timing IQ block STIQSampling rate conversion and sampling time delay adjustment;
s5, a server TDOA calculation module: receive timestamp alignment NtaskThe measurement stations are timed IQ blocks, and are combined into a measurement station pair when R ismodeCalling sample rate conversion sum when not equal to 0The delay adjusting module processes each timing IQ block STIQComputing TDOA measured value T of the station pairmeasRecording a time stamp calibration value or correcting the TDOA measurement value by using a supplementary time stamp calibration value according to the task type;
s6, a server positioning resolving module: calculating the coordinates of the emission source as O according to the TDOA measured value and the coordinates of the measuring stations of the measuring station pairLOCWhen M ismodeWhen the time difference is not less than 2, recording TDOA errors and positioning errors and counting the supplementary timestamp calibration values;
s7, a central station display module: TDOA hyperbolic measurements and truth, location result measurements and truth are displayed on the map, and TDOA error, timestamp calibration value, supplemental timestamp calibration value, and location error are displayed in a table.
2. The heterogeneous TDOA location system of claim 1, wherein said task parameters of S1 include:
s11, assigning a serial number to each station, establishing a relation table of station serial numbers and station model numbers, and designating the model M0For time-stamped reference stations, model number not M0Namely the heterogeneous station testing, the calibration value of the time stamp of the station testing is stored in the quintupleIn a constructed alignment queue, wherein M1nd、M2ndTwo stations type, FrateIs the sampling rate FS,TcaliFor the time stamp calibration value it is,calibrating values for the supplemental timestamps;
s12, the command parameter includes signal frequency FCSignal bandwidth BWSampling rate FSData acquisition period TmarkThe number of samples L of the timing IQ blockblockExtended sampling number LextTime stamp calibration valueSimulated emission point coordinates S0And simulation station coordinates SkAnd the number of stations N required for positioningtaskWherein K is the number of the stations, K belongs to {1,2, ⋯, K }, and K is the number of stations in the system;
s13, when MmodeTime stamp calibration value for all stations at time of =1When M ismodeWhen not equal to 1, if the calibration queue is empty, reading the record in the time stamp calibration file to the calibration queue, and if the model is MnThe kth station of (1), when Mn=M0Device for placingOtherwise, searching the alignment queue for meeting condition M1nd=M0And M2nd=MnAnd Frate=FSTo obtain a timestamp calibration value。
3. The heterogeneous TDOA location system of claim 1, wherein said S3 periodically extracts sample blocks S from the receiver output IQ streamIQAnd packaged as a timed IQ block STIQThe method comprises the following specific steps:
s31, obtaining the actual sampling rate of the receiverFrom a reference time TREFStart with TmarkIs the time point T of the cyclecapThe time for synchronously extracting the sampling blocks of each testing station is used for calculating TcapL=Tcap-Lext/(2×FS) Searching IQ sample stream output from receiver for time stamp TcapLFrom which a sample point of (L) is extractedext+Lblock)×(/FS) Is sampled set SIQWherein T isREFIs a constant;
s32, numbering stations and actual sampling rateTime stamp TcapStation coordinates and sample set SIQPackaged as a timed IQ block STIQWhen M ismodeAnd if not, the coordinate of the measuring station is the simulation measuring station coordinate in the command, otherwise, the coordinate is the coordinate output by the satellite positioning and time service module.
4. The heterogeneous TDOA location system as recited in claim 1, wherein said timed IQ block S4 isTIQThe sampling rate conversion and sampling time delay adjustment method specifically comprises the following steps:
s41 opposite timing IQ block STIQSample execution in (1)To FSSo that each station-timing IQ block has Lext+LblockA sample, at this time, sample STIQ (LextTime stamp of/2) is Tcap;
S42, for the station with the number k, when M ismodeTime P of =3k=0, otherwise, the simulation transmitting point coordinates S are calculated0And simulation survey station coordinate SkA distance D betweenkCalculating DkCorresponding number of samples Pk= ⌊FS×DkC + 0.5 ⌋, wherein C is the electromagnetic wave velocity constant;
5. The heterogeneous TDOA positioning system of claim 1, wherein the specific steps of S5 recording the timestamp calibration value or correcting the TDOA measurement value using the supplemental timestamp calibration value are as follows:
s51, when MmodeIf =1, the emission point S is determined according to the simulation0Calculating TDOA truth value T by coordinate distance difference of two simulation test stationstrueCalculating the TDOA error Tdev=Tmeas-TtrueWhen the model of the first station in the station pair is M0Second station non-M0Time, TDOA error TdevI.e. the timestamp calibration value between stations, and the model M of the first station in the pair of stations1The model M of the second testing station2Sampling rate FSAnd timestamp calibration value TdevIs a five-membered groupWill beStoring the data into a calibration queue;
s52, when Mmode=1, at the end of the task, calculating the mean value of the timestamp calibration values of each heterogeneous station at different sampling rates according to the nodes cached in the calibration queue, and taking 5-tuple Q as the reference valuecaliThe format is recorded in a time stamp calibration file, and the time stamp calibration value is supplemented;
S53, when MmodeNot equal to 1, according to the model of the first measuring station, the model of the second measuring station and the sampling rate FSSearching matched nodes from the calibration queue to obtain the supplementary timestamp calibration valueCorrection of TDOA measured value Tmeas= Tmeas-。
6. The heterogeneous TDOA positioning system of claim 1, wherein the steps of S6, recording TDOA errors and positioning errors and counting the supplemental timestamp calibration values, are as follows:
s61, calculating OLOCAnd S0The distance is the positioning error, and the sampling rate F is recordedSThe parameters such as the model, the error and the positioning error of each TDOA measuring station are put into an error analysis file;
s62, when the task is finished, the TDOA error mean value e of each station pair is countedtdoaAnd the mean value of the positioning error and the supplementary time stamp calibration value of each node in the time stamp calibration fileTogether, write the last row of the error analysis file, model and sampling rate F in the time stamp calibration file and the measuring station pairSSupplemental timestamp calibration values for matching nodesUpdating supplemental timestamp calibration values in calibration file nodes=etdoa。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011237117.3A CN112305498A (en) | 2020-11-09 | 2020-11-09 | Heterogeneous TDOA (time difference of arrival) positioning system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011237117.3A CN112305498A (en) | 2020-11-09 | 2020-11-09 | Heterogeneous TDOA (time difference of arrival) positioning system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112305498A true CN112305498A (en) | 2021-02-02 |
Family
ID=74326537
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011237117.3A Pending CN112305498A (en) | 2020-11-09 | 2020-11-09 | Heterogeneous TDOA (time difference of arrival) positioning system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112305498A (en) |
Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020196186A1 (en) * | 2001-06-25 | 2002-12-26 | Harris Corporation | Method and system for calibrating wireless location systems |
US20070053340A1 (en) * | 2005-08-09 | 2007-03-08 | Guilford John H | Time synchronization system and method for synchronizing locating units within a communication system using a known external signal |
JP2009025028A (en) * | 2007-07-17 | 2009-02-05 | Brother Ind Ltd | Positioning system |
CN101472330A (en) * | 2007-12-28 | 2009-07-01 | 三星电子株式会社 | Transmission synchronization method based on time difference orientation |
CN103116083A (en) * | 2013-01-23 | 2013-05-22 | 华中科技大学 | Very low frequency (VLF) / low frequency (LF) frequency range ground lightening detection positioning system |
CN103760522A (en) * | 2014-01-09 | 2014-04-30 | 电子科技大学 | Method and system used for time difference estimation and multi-station clock error correction |
US20160291124A1 (en) * | 2015-04-05 | 2016-10-06 | Nicholaus J. Bauer | Determining a location of a transmitter device |
WO2017063271A1 (en) * | 2015-10-16 | 2017-04-20 | 四川中电昆辰科技有限公司 | Quasi-synchronization structure, positioning device thereof, and positioning method therefor |
CN106646377A (en) * | 2016-12-29 | 2017-05-10 | 西安科技大学 | Positioning method for vibration target based on time sequence similarity search |
CN107422300A (en) * | 2017-09-05 | 2017-12-01 | 芜湖华创光电科技有限公司 | It is a kind of to coordinate wireless communication terminal signal step-out time extracting method to non- |
CN108093474A (en) * | 2017-12-27 | 2018-05-29 | 重庆脉实智能制造有限公司 | The method and system of indoor positioning are synchronously realized using virtual time |
CN109951799A (en) * | 2019-03-22 | 2019-06-28 | 西京学院 | A kind of object impact point positioning method based on sound |
CN110051360A (en) * | 2019-05-10 | 2019-07-26 | 成都信息工程大学 | Herding based on LoRaWAN technology monitors orientation management system |
CN110850364A (en) * | 2019-10-24 | 2020-02-28 | 山西大学 | UWB positioning method for eliminating antenna delay error |
US20200166652A1 (en) * | 2016-06-06 | 2020-05-28 | Brian G. Agee | Multi-Subband Methods for Reduced Complexity, Wideband Blind Resilient Detection and Geo-Observable Estimation of Global Navigation Satellite Signals |
CN111273274A (en) * | 2020-03-12 | 2020-06-12 | 四川九洲电器集团有限责任公司 | Multi-base cooperative positioning method, storage medium, radar and radar positioning system |
CN111343571A (en) * | 2020-03-07 | 2020-06-26 | 郑州联睿电子科技有限公司 | 5G-oriented integrated positioning system and positioning method integrating UWB |
CN111491266A (en) * | 2020-05-16 | 2020-08-04 | 深圳市引力波联科技有限公司 | Time synchronization method for UWB TDOA positioning system base station |
-
2020
- 2020-11-09 CN CN202011237117.3A patent/CN112305498A/en active Pending
Patent Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020196186A1 (en) * | 2001-06-25 | 2002-12-26 | Harris Corporation | Method and system for calibrating wireless location systems |
US20070053340A1 (en) * | 2005-08-09 | 2007-03-08 | Guilford John H | Time synchronization system and method for synchronizing locating units within a communication system using a known external signal |
JP2009025028A (en) * | 2007-07-17 | 2009-02-05 | Brother Ind Ltd | Positioning system |
CN101472330A (en) * | 2007-12-28 | 2009-07-01 | 三星电子株式会社 | Transmission synchronization method based on time difference orientation |
CN103116083A (en) * | 2013-01-23 | 2013-05-22 | 华中科技大学 | Very low frequency (VLF) / low frequency (LF) frequency range ground lightening detection positioning system |
CN103760522A (en) * | 2014-01-09 | 2014-04-30 | 电子科技大学 | Method and system used for time difference estimation and multi-station clock error correction |
US20160291124A1 (en) * | 2015-04-05 | 2016-10-06 | Nicholaus J. Bauer | Determining a location of a transmitter device |
WO2017063271A1 (en) * | 2015-10-16 | 2017-04-20 | 四川中电昆辰科技有限公司 | Quasi-synchronization structure, positioning device thereof, and positioning method therefor |
US20200166652A1 (en) * | 2016-06-06 | 2020-05-28 | Brian G. Agee | Multi-Subband Methods for Reduced Complexity, Wideband Blind Resilient Detection and Geo-Observable Estimation of Global Navigation Satellite Signals |
CN106646377A (en) * | 2016-12-29 | 2017-05-10 | 西安科技大学 | Positioning method for vibration target based on time sequence similarity search |
CN107422300A (en) * | 2017-09-05 | 2017-12-01 | 芜湖华创光电科技有限公司 | It is a kind of to coordinate wireless communication terminal signal step-out time extracting method to non- |
CN108093474A (en) * | 2017-12-27 | 2018-05-29 | 重庆脉实智能制造有限公司 | The method and system of indoor positioning are synchronously realized using virtual time |
CN109951799A (en) * | 2019-03-22 | 2019-06-28 | 西京学院 | A kind of object impact point positioning method based on sound |
CN110051360A (en) * | 2019-05-10 | 2019-07-26 | 成都信息工程大学 | Herding based on LoRaWAN technology monitors orientation management system |
CN110850364A (en) * | 2019-10-24 | 2020-02-28 | 山西大学 | UWB positioning method for eliminating antenna delay error |
CN111343571A (en) * | 2020-03-07 | 2020-06-26 | 郑州联睿电子科技有限公司 | 5G-oriented integrated positioning system and positioning method integrating UWB |
CN111273274A (en) * | 2020-03-12 | 2020-06-12 | 四川九洲电器集团有限责任公司 | Multi-base cooperative positioning method, storage medium, radar and radar positioning system |
CN111491266A (en) * | 2020-05-16 | 2020-08-04 | 深圳市引力波联科技有限公司 | Time synchronization method for UWB TDOA positioning system base station |
Non-Patent Citations (3)
Title |
---|
DAPENG WANG AND HONG DU: "A TDOA Estimation Method Through Phase Difference Square Sums (PDSS)", 《IOP PUBLISHING LTD IOP CONFERENCE SERIES: MATERIALS SCIENCE AND ENGINEERING》, vol. 569, no. 5, pages 052007 * |
罗德与施瓦茨(中国)有限公司: "基于TDOA无线电定位系统的接收机要求——罗德与施瓦茨用于TDOA网络的传感器技术(一)", 《中国无线电》, vol. 224, no. 4, pages 72 - 73 * |
陈金苗: "基于异构天线阵列的多维参数联合估计TDOA测量方法研究", 《中国优秀硕士学位论文全文数据库 信息科学辑》, no. 8, pages 136 - 5 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110351658B (en) | Indoor positioning method based on convolutional neural network | |
CN107113762B (en) | Positioning method, positioning server and positioning system | |
CN110856106B (en) | Indoor high-precision three-dimensional positioning method based on UWB and barometer | |
CN110168396A (en) | Arrival time (TOA) measurement | |
CN104396321B (en) | Method of locating terminal and positioning device | |
US10034262B2 (en) | Beacon sensor based network fingerprint | |
US20110189960A1 (en) | Estimating Whether A Wireless Terminal Is Indoors Using Pattern Classification | |
CN112543078A (en) | Network time server calibration method and device based on satellite common view | |
EP4398629A2 (en) | Fingerprint positioning | |
CN107231636A (en) | A kind of method and apparatus of calibration network coverage evaluating | |
JPWO2021211399A5 (en) | ||
CN102550051A (en) | LTE fingerprinting positioning references for other cellular systems | |
Palamà et al. | From experiments to insights: A journey in 5G new radio localization | |
CN101316151B (en) | System for data traffic delay time test and error cancelling method | |
CN112305498A (en) | Heterogeneous TDOA (time difference of arrival) positioning system | |
CN109669157A (en) | TDOA estimation method based on compression transmission | |
CN111158027A (en) | Multi-terminal synchronization-based network RTK service testing method and device and testing system | |
WO2016106667A1 (en) | Method and apparatus for processing information used for positioning | |
WO2017049914A1 (en) | Terminal positioning method, apparatus, and system | |
CN104459682A (en) | Distance measurement method for automatically eliminating dynamic distance zero value | |
CN101848119A (en) | High-accuracy network round-trip delay measuring method for eliminating system processing time | |
CN111246386A (en) | Terminal positioning method and device | |
CN114089273A (en) | GPS and UWB based motion platform positioning method | |
CN109752594B (en) | Multichannel synchronous data acquisition phase correction method | |
CN111464397A (en) | Method and system for measuring bidirectional distance and clock error |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20210202 |