CN106793087A - A kind of array antenna indoor positioning algorithms based on AOA and PDOA - Google Patents

A kind of array antenna indoor positioning algorithms based on AOA and PDOA Download PDF

Info

Publication number
CN106793087A
CN106793087A CN201710156453.7A CN201710156453A CN106793087A CN 106793087 A CN106793087 A CN 106793087A CN 201710156453 A CN201710156453 A CN 201710156453A CN 106793087 A CN106793087 A CN 106793087A
Authority
CN
China
Prior art keywords
theta
pdoa
signal
calculated
aoa
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
CN201710156453.7A
Other languages
Chinese (zh)
Other versions
CN106793087B (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.)
Tianjin University
Original Assignee
Tianjin University
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 Tianjin University filed Critical Tianjin University
Priority to CN201710156453.7A priority Critical patent/CN106793087B/en
Publication of CN106793087A publication Critical patent/CN106793087A/en
Application granted granted Critical
Publication of CN106793087B publication Critical patent/CN106793087B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W64/00Locating users or terminals or network equipment for network management purposes, e.g. mobility management
    • H04W64/006Locating users or terminals or network equipment for network management purposes, e.g. mobility management with additional information processing, e.g. for direction or speed determination
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations
    • G01C21/206Instruments for performing navigational calculations specially adapted for indoor navigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

The present invention relates to a kind of array antenna indoor positioning algorithms based on AOA and PDOA, including:Multipath signal is distinguished using array antenna;Energy size according to the signal for receiving selects two most strong signal paths of energy, and the phase information of two signal paths is obtained from antenna;The angle of arrival of two bars is calculated using AOA methods of estimation;The propagation distance of two bars is calculated using PDOA methods of estimation;According to the position of the total barrier of positioning scene, virtual base station can be set up, non line of sight NLOS paths are converted into sighting distance LOS path;The position coordinates of object to be positioned is calculated using weighted least-squares WLS algorithms;Using residual weighted LS algorithms to obtaining more accurate position coordinates.

Description

A kind of array antenna indoor positioning algorithms based on AOA and PDOA
Technical field
The invention belongs to indoor positioning algorithms technical field, particular for multipath and the indoor positioning environment of non line of sight.
Background technology
The research and application of location technology are quite varied, permeate in various aspects, global location such as military, business, lives System (GPS) provides good Global Positioning Service, but it cannot be used under environment indoors.Passive ultra-high frequency radio frequency is known (UHF RFID) indoor positioning technologies are not paid close attention to by increasing people, because passive label has small volume, deployment letter It is single, without power supply, it is cheap, read distant, the advantage such as read or write speed, passive ultra-high frequency RFID location technology faster Context aware field is progressively applied to, such as in terms of logistics management, intelligent transportation, robot, office automatic.
The method of indoor wireless positioning includes being based on arrival time (TOA), reaching time-difference (TDOA), angle of arrival (AOA), arriving signal intensity (RSS) and arrival phase difference (PDOA) etc..However, in the UHF FRID systems of smaller bandwidth, The method of TOA and TDOA is difficult to be implemented;Signal fadeout of the method based on RSS to multipath is very sensitive, so based on simple The alignment system of RSS methods, positioning precision is not very high;Method based on PDOA is also influenceed by multipath so that range finding knot It is really inaccurate;Meanwhile, the method based on AOA is influenceed by non line of sight (NLOS) and multipath.
Common decrease multipath and non line of sight are divided into statistics method and geometry to the method that indoor position accuracy influences Method.No matter this two classes method is in the modeling link of front end, or location algorithm link in rear end suffers from being widely applied. For the method that modeling link suppresses error, statistical modeling is modeled generally be directed to non-market value, counts special with certain Levy to describe non-market value, or scattering object to causing multi-path problem is modeled;And Geometric Modeling, such as ray trace (Ray-tracing), mainly use radio signal propagation characteristic close with the propagation characteristic of light, signal Analysis are in specific environment The path of lower propagation.For the method that location algorithm link suppresses error, statistical method is missed by measuring reception signal Difference cloth is come to determine signal path be sighting distance or non line of sight;And geometric method is that signal road is determined by geometrical relationship Footpath is sighting distance or non line of sight.
The maximum shortcoming of statistical method is to need known signal statistical property under certain circumstances, and this is accomplished by passing through Substantial amounts of reality measures to obtain this statistical property.When the environment changes, statistical property can change, it is necessary to again a large amount of Measurement carry out corrected parameter.And geometric method does not need substantial amounts of measurement, known indoor arrangement figure is only needed.But it is existing Some methods are positioned using direct path or primary scattering path, are not deposited for Multiple Scattering problem, or hypothesis , or filtered by algorithm judgement.
The content of the invention
The present invention provides a kind of influence that can well improve multipath and non line of sight to indoor positioning, effectively improves and determine The indoor positioning algorithms of position precision, technical scheme is as follows:
A kind of array antenna indoor positioning algorithms based on AOA and PDOA, comprise the following steps:
1) multipath signal is distinguished using array antenna;
2) the energy size according to the signal for receiving selects two most strong signal paths of energy, and from antenna The middle phase information for obtaining two signal paths;
3) angle of arrival of two bars is calculated using AOA methods of estimation;
4) propagation distance of two bars is calculated using PDOA methods of estimation;
5) according to the position of the total barrier of positioning scene, virtual base station can be set up, non line of sight NLOS paths is converted into Sighting distance LOS path;
6) according to step 3) in the angle information, the step 4 that obtain) in the range information that obtains, joint step 5) in set up Virtual base station, the position coordinates (x ', y ') of object to be positioned is calculated using weighted least-squares WLS algorithms:
Z1=[x ', y ', x '2+y′2]T=(GTWG)-1GTWH
WhereinDiag { } represents diagonal matrix,
Z=[x, y, x2+y2]T, Xk=[xk,yk]T, k=1,2 ..., N is corresponding with kth signal path virtual Base station, dkIt is the path length of the kth signal path estimated by PDOA, and dk≤dk+1, θkIt is kth signal path Angle of arrival;
7) using residual weighted LS algorithms to step 6) position that obtains optimizes, and obtains more accurate position coordinates (x,y)。
The present invention distinguishes multipath signal using aerial array, and obstructed path is converted into sighting distance using virtual base station Path.The most strong signal path of two energy is chosen, using AOA and PDOA algorithm for estimating, combines weighted least-squares and residual error Weighting algorithm calculates the position of object to be positioned.The influence to indoor positioning of multipath and non line of sight can well be improved, had Improve positioning precision in effect ground.
Brief description of the drawings
Fig. 1 is overall flow block diagram of the invention.
Fig. 2 is positioning scene schematic diagram of the invention.
Obstructed path is converted into Fig. 3 the schematic diagram of los path using virtual base station.
Fig. 4 is the analogous diagram of the positioning precision using the algorithm in the present invention, and the environment of emulation is a room of 10m*10m Son, in the environment of white Gaussian noise, using algorithm of the invention, positioning precision is less than 0.5m with 80% probability.
Specific embodiment
Below with uniform linear array, algorithm of the present invention is described in detail.
1) multipath signal is distinguished using array antenna.
2) the energy size according to the signal for receiving selects two most strong signal paths of energy, and from antenna The middle phase information for obtaining two signal paths.
3) angle of arrival of two bars is calculated using AOA methods of estimation, according to step 2) in the phase information that obtains, The angle of arrival of signal can be just calculated using the MUSIC algorithms based on space smoothing classical in AOA methods of estimation.
4) propagation distance of two bars, the following institute of formula of PDOA methods range finding are calculated using PDOA methods of estimation Show:
Wherein, c is the aerial spread speed of light, and Δ φ represents that signal travels to object to be positioned and produces from array antenna Raw phase difference, π=3.14, Δ f is the difference on the frequency between two signals of different frequency, and d is the biography of the signal being calculated Broadcast distance.
5) according to the position of the total barrier of positioning scene, virtual base station can be set up, by the conversion of non line of sight (NLOS) path It is sighting distance (LOS) path.As shown in figure 3, XbIt is the position of array antenna, XtIt is the position of object to be positioned, XvIt is according to obstacle The virtual base station that the position that thing l sets up is set up.From virtual base station as viewed from object to be positioned, we will obtain a sighting distance Path.
6) according to step 3) in the angle information, the step 4 that obtain) in the range information that obtains, joint step 5) in set up Virtual base station, the position of object to be positioned is calculated using weighted least-squares (WLS) algorithm:
It is provided with being estimated by entity base station apart from d and angle, θ for N signal paths, Xk=[xk,yk]T, k=1, 2 ..., N is the virtual base station corresponding with kth signal path, dkIt is the kth signal path estimated by PDOA Path length, and dk≤dk+1, θkIt is the angle of arrival of kth signal path.With X=[x, y]TRepresent the position of object to be positioned Put, then:
x sin(θi)-y cos(θi)=xi sin(θi)+yi cos(θi)
Write two formula above as matrix form, GZ=H can be obtained, wherein:
Z=[x, y, x2+y2]T
The position of object to be positioned is calculated using weighted least square algorithm, the initial position co-ordinates of object to be positioned are obtained (x′,y′):
Z1=[x ', y ', x '2+y′2]T=(GTWG)-1GTWH,
WhereinDiag { } represents diagonal matrix,
7) setUsing residual weighted (LS) algorithm To step 6) position that obtains optimizes, and obtains the position X=[x, y] of more accurate object to be positionedT, wherein:
X′i=[x 'i,y′i]TI=1,2 ..., N is the object to be positioned derived from the i-th signal path Position.

Claims (1)

1. a kind of array antenna indoor positioning algorithms based on AOA and PDOA, comprise the following steps:
1) multipath signal is distinguished using array antenna.
2) the energy size according to the signal for receiving selects two most strong signal paths of energy, and is obtained from antenna Take the phase information of two signal paths;
3) angle of arrival of two bars is calculated using AOA methods of estimation;
4) propagation distance of two bars is calculated using PDOA methods of estimation;
5) according to the position of the total barrier of positioning scene, virtual base station can be set up, non line of sight NLOS paths is converted into sighting distance LOS path;
6) according to step 3) in the angle information, the step 4 that obtain) in the range information that obtains, joint step 5) in the void set up Intend base station, the position coordinates (x ', y ') of object to be positioned is calculated using weighted least-squares WLS algorithms:
Z1=[x ', y ', x '2+y′2]T=(GTWG)-1GTWH
WhereinDiag { } represents diagonal matrix,
G = - 2 x 1 - 2 y 1 1 . . . . . . . . . - 2 x N - 2 y N 1 sin ( θ 1 ) cos ( θ 1 ) 0 . . . . . . . . . sin ( θ N ) cos ( θ N ) 0 , H = d 1 2 - ( x 1 2 + y 1 2 ) . . . d N 2 - ( x N 2 + y N 2 ) x 1 sin ( θ 1 ) - y 1 cos ( θ 1 ) . . . x N sin ( θ N ) - y N cos ( θ N ) ,
Z=[x, y, x2+y2]T, Xk=[xk,yk]T, k=1,2 ..., N are the virtual base stations corresponding with kth signal path, dkIt is the path length of the kth signal path estimated by PDOA, and dk≤dk+1, θkIt is arriving for kth signal path Up to angle;
7) using residual weighted LS algorithms to step 6) position that obtains optimizes, obtain more accurate position coordinates (x, y)。
CN201710156453.7A 2017-03-16 2017-03-16 Array antenna indoor positioning method based on AOA and PDOA Active CN106793087B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710156453.7A CN106793087B (en) 2017-03-16 2017-03-16 Array antenna indoor positioning method based on AOA and PDOA

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710156453.7A CN106793087B (en) 2017-03-16 2017-03-16 Array antenna indoor positioning method based on AOA and PDOA

Publications (2)

Publication Number Publication Date
CN106793087A true CN106793087A (en) 2017-05-31
CN106793087B CN106793087B (en) 2020-01-21

Family

ID=58966896

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710156453.7A Active CN106793087B (en) 2017-03-16 2017-03-16 Array antenna indoor positioning method based on AOA and PDOA

Country Status (1)

Country Link
CN (1) CN106793087B (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107831517A (en) * 2017-10-18 2018-03-23 国家新闻出版广电总局广播科学研究院 Combined positioning method to failing to understand radio and television intrusion signal multi-point cooperative direction finding
CN107883959A (en) * 2017-10-31 2018-04-06 西安交通大学 More people's localization methods in a kind of WiFi rooms based on phased array principle
CN108811082A (en) * 2018-06-26 2018-11-13 中国人民解放军陆军工程大学 Wireless positioning method, system and positioning terminal
GB2563834A (en) * 2017-06-23 2019-01-02 Decawave Ltd Wideband antenna array
CN109669176A (en) * 2019-01-17 2019-04-23 普联技术有限公司 Method, apparatus and storage medium based on wireless multi-path measurement house type
CN109698999A (en) * 2017-10-23 2019-04-30 深圳市优必选科技有限公司 Positioning method and terminal equipment
CN109991572A (en) * 2019-04-25 2019-07-09 哈尔滨工业大学 A kind of two-shipper passive location method based on azimuth and pitch angle information
CN111615050A (en) * 2020-04-24 2020-09-01 郑州联睿电子科技有限公司 Emergency positioning searching device, method and system
CN112505624A (en) * 2020-12-09 2021-03-16 电子科技大学 Positioning system based on environment backscattering
CN116106823A (en) * 2023-02-18 2023-05-12 中国人民解放军海军工程大学 TDOA-PDOA combined positioning method based on particle swarm optimization
CN117939630A (en) * 2024-01-23 2024-04-26 电子科技大学 Indoor complex environment passive positioning system and method based on ultrahigh frequency RFID technology

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104812063A (en) * 2015-04-08 2015-07-29 天津大学 TOA (time of arrival) positioning method based on virtual sensors in indoor environment
CN105572635A (en) * 2016-01-25 2016-05-11 西安电子科技大学 Single-station passive quick positioning method based on least square method
CN105960013A (en) * 2016-05-06 2016-09-21 华东交通大学 AOA-based cooperative localization method under non line-of-sight environment
KR101658465B1 (en) * 2015-03-05 2016-09-22 국방과학연구소 Location estimation method based on direction finding using tdoa of signal and direction finder thereof
CN106291540A (en) * 2016-09-14 2017-01-04 河北省电力勘测设计研究院 A kind of multiple-input and multiple-output GPR backwards projection target imaging method estimated based on DOA

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101658465B1 (en) * 2015-03-05 2016-09-22 국방과학연구소 Location estimation method based on direction finding using tdoa of signal and direction finder thereof
CN104812063A (en) * 2015-04-08 2015-07-29 天津大学 TOA (time of arrival) positioning method based on virtual sensors in indoor environment
CN105572635A (en) * 2016-01-25 2016-05-11 西安电子科技大学 Single-station passive quick positioning method based on least square method
CN105960013A (en) * 2016-05-06 2016-09-21 华东交通大学 AOA-based cooperative localization method under non line-of-sight environment
CN106291540A (en) * 2016-09-14 2017-01-04 河北省电力勘测设计研究院 A kind of multiple-input and multiple-output GPR backwards projection target imaging method estimated based on DOA

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
SAMEH SOROUR ; YVES LOSTANLEN ; SHAHROKH VALAEE ; KHAQAN MAJEED: "Joint Indoor Localization and Radio Map Construction with Limited Deployment Load", 《IEEE》 *
YISHUANG GENG ; KAVEH PAHLAVAN: "Design, Implementation, and Fundamental Limits of Image and RF Based Wireless Capsule Endoscopy Hybrid Localization", 《IEEE》 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2563834A (en) * 2017-06-23 2019-01-02 Decawave Ltd Wideband antenna array
US11128058B2 (en) 2017-06-23 2021-09-21 DecaWave, Ltd. Wideband antenna array
CN107831517A (en) * 2017-10-18 2018-03-23 国家新闻出版广电总局广播科学研究院 Combined positioning method to failing to understand radio and television intrusion signal multi-point cooperative direction finding
CN109698999B (en) * 2017-10-23 2021-04-16 深圳市优必选科技有限公司 Positioning method and terminal equipment
CN109698999A (en) * 2017-10-23 2019-04-30 深圳市优必选科技有限公司 Positioning method and terminal equipment
CN107883959A (en) * 2017-10-31 2018-04-06 西安交通大学 More people's localization methods in a kind of WiFi rooms based on phased array principle
CN108811082A (en) * 2018-06-26 2018-11-13 中国人民解放军陆军工程大学 Wireless positioning method, system and positioning terminal
CN108811082B (en) * 2018-06-26 2020-07-21 中国人民解放军陆军工程大学 Wireless positioning method, system and positioning terminal
CN109669176A (en) * 2019-01-17 2019-04-23 普联技术有限公司 Method, apparatus and storage medium based on wireless multi-path measurement house type
CN109991572A (en) * 2019-04-25 2019-07-09 哈尔滨工业大学 A kind of two-shipper passive location method based on azimuth and pitch angle information
CN111615050A (en) * 2020-04-24 2020-09-01 郑州联睿电子科技有限公司 Emergency positioning searching device, method and system
CN112505624A (en) * 2020-12-09 2021-03-16 电子科技大学 Positioning system based on environment backscattering
CN112505624B (en) * 2020-12-09 2023-01-24 电子科技大学 Positioning system based on environment backscattering
CN116106823A (en) * 2023-02-18 2023-05-12 中国人民解放军海军工程大学 TDOA-PDOA combined positioning method based on particle swarm optimization
CN117939630A (en) * 2024-01-23 2024-04-26 电子科技大学 Indoor complex environment passive positioning system and method based on ultrahigh frequency RFID technology

Also Published As

Publication number Publication date
CN106793087B (en) 2020-01-21

Similar Documents

Publication Publication Date Title
CN106793087A (en) A kind of array antenna indoor positioning algorithms based on AOA and PDOA
Poulose et al. Localization error analysis of indoor positioning system based on UWB measurements
CN109275095B (en) Bluetooth-based indoor positioning system, positioning equipment and positioning method
Thomas et al. Performance of a TDOA-AOA hybrid mobile location system
CN110045324B (en) Indoor positioning fusion method based on UWB and Bluetooth technology
CN103501538B (en) Based on the indoor orientation method of multipath energy fingerprint
CN113271540B (en) Bluetooth signal positioning method, multi-signal fusion positioning method and system
CN109341679A (en) A kind of smart machine air navigation aid and navigation system
CN103995250B (en) Radio-frequency (RF) tag trajectory track method
CN102928813A (en) RSSI (Received Signal Strength Indicator) weighted centroid algorithm-based passive RFID (Radio Frequency Identification Device) label locating method
Sung et al. Accurate indoor positioning for UWB-based personal devices using deep learning
CN107071732B (en) RSSI (received Signal Strength indicator) -based MLE-PSO (Multi-layer programmable Gate array) indoor positioning method
CN204166130U (en) Radio frequency locating device and system
CN103813448A (en) Indoor positioning method based on RSSI
CN112929826B (en) Indoor hybrid positioning method based on UWB triangulation positioning and fingerprint information
CN104330771A (en) Indoor RFID precise positioning method and device
CN104181500A (en) Real-time locating method based on inertia information and chance wireless signal characteristics
CN113794991B (en) Single-base-station wireless positioning system based on UWB and LoRa
CN104869636A (en) Indoor positioning method based on distance measurement information fusion
CN107356903A (en) Passive RFID localization method and device based on phase difference measurement
EP4137835A1 (en) Ultra-wide band distance determination with an angle-of-arrival based disturbance compensation
CN110187333A (en) A kind of RFID label tag localization method based on synthetic aperture radar technique
Wang et al. Single base station positioning based on multipath parameter clustering in NLOS environment
CN108344411A (en) Indoor high-accuracy position system
CN109254264A (en) Improved UHF-RFID positioning 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