CN103379619B  A kind of localization method and system  Google Patents
A kind of localization method and system Download PDFInfo
 Publication number
 CN103379619B CN103379619B CN201210111075.8A CN201210111075A CN103379619B CN 103379619 B CN103379619 B CN 103379619B CN 201210111075 A CN201210111075 A CN 201210111075A CN 103379619 B CN103379619 B CN 103379619B
 Authority
 CN
 China
 Prior art keywords
 mrow
 mtd
 msub
 course angle
 mtr
 Prior art date
Links
 230000004807 localization Effects 0.000 title claims abstract description 13
 239000002245 particles Substances 0.000 claims description 110
 238000001914 filtration Methods 0.000 claims description 47
 280000638271 Reference Point companies 0.000 claims description 41
 239000011159 matrix materials Substances 0.000 claims description 25
 230000000875 corresponding Effects 0.000 claims description 14
 230000001133 acceleration Effects 0.000 claims description 10
 238000000034 methods Methods 0.000 claims description 9
 230000004927 fusion Effects 0.000 claims description 8
 238000010606 normalization Methods 0.000 claims description 7
 238000005096 rolling process Methods 0.000 claims description 2
 238000005259 measurements Methods 0.000 claims 1
 238000007667 floating Methods 0.000 abstract description 3
 239000000969 carriers Substances 0.000 description 16
 230000000295 complement Effects 0.000 description 7
 238000009826 distribution Methods 0.000 description 7
 241001269238 Data Species 0.000 description 5
 244000171263 Ribes grossularia Species 0.000 description 5
 238000004422 calculation algorithm Methods 0.000 description 5
 238000004364 calculation methods Methods 0.000 description 4
 238000003860 storage Methods 0.000 description 4
 238000005516 engineering processes Methods 0.000 description 3
 241000209094 Oryza Species 0.000 description 2
 235000007164 Oryza sativa Nutrition 0.000 description 2
 238000010586 diagrams Methods 0.000 description 2
 230000005484 gravity Effects 0.000 description 2
 238000002156 mixing Methods 0.000 description 2
 235000009566 rice Nutrition 0.000 description 2
 240000004678 Panax pseudoginseng Species 0.000 description 1
 235000003140 Panax quinquefolius Nutrition 0.000 description 1
 208000008425 Protein Deficiency Diseases 0.000 description 1
 230000003044 adaptive Effects 0.000 description 1
 235000005035 ginseng Nutrition 0.000 description 1
 235000008434 ginseng Nutrition 0.000 description 1
 235000009808 lpulo Nutrition 0.000 description 1
 238000004088 simulation Methods 0.000 description 1
 238000003892 spreading Methods 0.000 description 1
 238000004642 transportation engineering Methods 0.000 description 1
 239000011901 water Substances 0.000 description 1
Classifications

 H—ELECTRICITY
 H04—ELECTRIC COMMUNICATION TECHNIQUE
 H04W—WIRELESS COMMUNICATION NETWORKS
 H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
 H04W4/02—Services making use of location information
 H04W4/025—Services making use of location information using location based information parameters
 H04W4/027—Services making use of location information using location based information parameters using movement velocity, acceleration information
Abstract
Description
Technical field
The present invention relates to wireless network positioning field, more particularly to a kind of localization method and system.
Background technology
Alignment system based on WLAN has that cost is low, and precision is high, have a wide range of application (indoor and outdoors) the advantages that, in base As emergency relief, intelligent transportation and indoor positioning navigation etc. achieve very big success in the service of position.But still Following two problem urgent need to resolve so be present：(1) received signal strength (Received caused by the factor such as multipath jamming Signal Strength, RSS) floating severe exacerbation WLAN positioning precisions；(2) in wireless access node (Access Point, AP) region that is not covered with, because AP missings cause WLAN positioning failures.
In order to solve the above problems, there has been proposed a variety of methods, can be divided into following three class：
1st, the WLAN alignment systems based on time diversity and probability Distribution Model
The basic thought of WLAN alignment systems based on time diversity and probability Distribution Model is：It is fixed in localization region Position obtains multiple samples of received signal strength using time diversity, and it is strong to establish reception signal according to multiple sample informations The probability Distribution Model of degree, the probability Distribution Model of received signal strength is stored into property data base；In positioning stage, move Movingtarget obtains multiple samples of received signal strength using time diversity, passes through the reception for asking for sample average to obtain stable Signal intensity is positioned.Because time diversity needs to consume the substantial amounts of time, positioning delay is added, can not be realized in real time Positioning, can not be used in running fix.
2. the WLAN alignment systems based on Kalman filtering
The basic thought of WLAN alignment systems based on Kalman filtering is：Moved first with WLAN location algorithms The location estimation of target, then using the speed of mobile target trajectory continuity or the mobile target of hypothesis within the specific limits, The state equation and observational equation of Kalman filter are constructed, processing is filtered to the location estimation of user.Though this method The positioning precision of WLAN alignment systems is so improved, but due to setting mobile target velocity in advance, therefore can not realize certainly Adaptive filtering, limit application in practice.Can not solve the WLAN positioning failure caused by AP is lacked simultaneously.
3.WLAN/GPS integrated positionings
Because GPS is in outdoor open space, good positioning precision can be obtained, therefore in the region of AP missings, can be with Accurate positional information is obtained using GPS；Simultaneously in builtup urban district, can be positioned using WLAN to make up GPS location Deficiency, it may be said that WLAN/GPS integrated positionings achieve good positioning performance in an outdoor environment, but environment indoors Under, because gps signal is blocked, WLAN indoor position accuracies can not be improved.
The content of the invention
The technical problem to be solved in the present invention is to provide a kind of localization method and system, solves existing alignment system positioning accurate Spend the problem of inadequate.
In order to solve the above problems, the invention provides a kind of localization method, including：
Access point based on WLAN is positioned, and obtains the initial estimated location of user equipment；
Obtain the course angle and velocity information of the user equipment；
The initial estimated location is modified according to the course angle and velocity information, obtains final position information.
Further, the above method can also be had the characteristics that, the access point based on WLAN is positioned Including：
Reference point is chosen, the signal intensity from each access point is measured in each reference point, by the position of the reference point, institute Mark/the position for stating signal intensity and corresponding access point is stored into database；Respectively connect around user equipment measurement to be positioned The signal intensity of access point, searching data storehouse obtain corresponding to refer to point set, with it is described with reference to point set match determine use The initial estimated location of family equipment.
Further, the above method can also have the characteristics that, described to determine the use with reference to point set match The initial estimated location of family equipment includes：
The reference point of the Euclidean distance minimum of m received signal strength of selection, uses the line of the position of the m reference point Property initial estimated location of the weighted sum as the user equipment, the m is more than or equal to 1.
Further, the above method can also have the characteristics that, the course angle and speed for obtaining the user equipment Information includes：
The course angle and velocity information are obtained according to the metrical information of MARG sensors.
Further, the above method can also have the characteristics that, described to obtain institute according to the metrical information of MARG sensors Stating course angle includes：
First course angle φ is obtained according to the metrical information of the magnetometer of the MARG sensors_{mag}, passed according to the MARG The metrical information of the gyroscope of sensor obtains the second course angleAccording to obtaining first course angle and the second course angle The course angle φ of user equipment：
Wherein, the W is default weighted value, 0≤W≤1.
Further, the above method can also have the characteristics that methods described also includes, according to the MARG sensors Metrical information obtains the roll angle of the user equipment, the angle of pitch；
It is described according to the course angle and velocity information the initial estimated location is modified including：
The roll angle obtained according to the metrical information of the MARG sensors, the angle of pitch, course angle and the speed are believed The input as Kalman filter is ceased, Kalman filtering is carried out, exports new course angle and velocity information；
The course angle and velocity information and the initial estimated location that the Kalman filter is exported are as particle The input of wave filter, carry out particle filter, output position information, course angle and velocity information, using the positional information of output as The final position information of the user equipment.
Further, the above method can also have the characteristics that, the progress Kalman filtering includes：
The state onestep prediction of Kalman filtering is carried out,
Calculate prediction varivance matrix
Calculate filtering gain matrix
Carry out state estimation
Calculate estimation error variance
Wherein, it is describedRepresent the angular speed of the gyroscope output of the MARG sensors, φ_{k1}Represent root during moment k1 The course angle obtained according to the MARG sensors, Δ T represent the time of measuring interval of the MARG sensors, and Q and R are represented respectively The covariance matrix of process noise and measurement noise, K_{k}For Kalman filter gain,And P_{k}Represent varivance matrix, φ_{PF} The course angle of particle filter output, first φ during Kalman filtering when being positioned for the last time_{PF}For designated value.
Further, the above method can also have the characteristics that, the progress particle filter includes：
, it is necessary to initialize particle when carrying out particle filter first, the probability density of particle is initialized using Gaussian Profile Function；
According to the course angle and velocity information, and the initial estimated location, predict the user equipment in next step Status information：
The weight of each particle and normalization are calculated, it is as follows：
Carry out particle resampling, the particle as particle filter next time；
Wherein, described [x_{k}, y_{k}]^{T}For the state vector of each particle, T_{s}Represent the last access based on WLAN The time interval of the positioning and the positioning of this access point based on WLAN of point, φ_{k}Represent the course angle, v_{k}Represent The velocity information, [η_{x}, η_{y}]^{T}Acceleration is represented, is simulated with the Gaussian noise of zeromean, variance is sensed by the MARG The metrical information estimation of device,To input the state value of the particle filter,Represent state of ith of particle in moment k Value, σ represent the noise variance of signal strength measurement.
The present invention also provides a kind of alignment system, including：
WLAN locating modules, positioned for the access point based on WLAN, obtain the initial of user equipment and estimate Count position；
Sensor locating module, for obtaining the course angle and velocity information of the user equipment；
Fusion Module, for being modified according to the course angle and velocity information to the initial estimated location, obtain Final position information.
Further, said system can also have the characteristics that, the WLAN locating modules include：
Database, for being stored in the signal intensity from each access point of each reference point measurement, the position of the reference point Put and correspond to mark/position of access point；
RSS measuring units, for measuring the signal intensity of each access point around user equipment measurement to be positioned；
Positioning unit, for the signal intensity measured according to the RSS measuring units, searching data storehouse obtains corresponding join Examination point set, with the initial estimated location for reference to point set match determination user equipment.
Further, said system can also have the characteristics that, the positioning unit with reference to point set match really The initial estimated location of the fixed user equipment includes：
The positioning unit selects the minimum reference point of the Euclidean distance of m received signal strength, uses described m reference The linear weighted function of the position of point and the initial estimated location as the user equipment, the m are more than or equal to 1.
Further, said system can also have the characteristics that, the sensor locating module includes：MARG sensors and Data processing unit, wherein：
The MARG sensors are used for, and the user equipment is measured, and obtain metrical information；
The data processing unit is used for, and obtaining the course angle and speed according to the metrical information of MARG sensors believes Breath.
Further, said system can also have the characteristics that, the data processing unit is according to the surveys of MARG sensors Course angle described in amount acquisition of information includes：
The data processing unit obtains the first course angle according to the metrical information of the magnetometer of the MARG sensors φ_{mag}, the second course angle is obtained according to the metrical information of the gyroscope of the MARG sensorsAccording to first course angle The course angle φ of the user equipment is obtained with the second course angle：
Wherein, the W is default weighted value, 0≤W≤1.
Further, said system can also have the characteristics that, the Fusion Module includes：Kalman filter and particle Wave filter, wherein：
The data processing unit is additionally operable to, and the user equipment is obtained according to the metrical information of the MARG sensors Roll angle, the angle of pitch；
The Kalman filter is used for, by it is described according to the metrical information of the MARG sensors obtain roll angle, The angle of pitch, course angle and velocity information input as the state value of Kalman filter, carry out Kalman filtering, export new boat To angle and velocity information；
The particle filter is used for, the course angle and velocity information that the Kalman filter is exported, and described Initial estimated location inputs as state value, carries out particle filter, output position information, course angle and velocity information, will export Final position information of the positional information as the user equipment.
Further, said system can also be had the characteristics that, the Kalman filter, which carries out Kalman filtering, to be included：
The state onestep prediction of Kalman filtering is carried out,
Calculate prediction varivance matrix
Calculate filtering gain matrix
Carry out state estimation
Calculate estimation error variance
Wherein, it is describedRepresent the angular speed of the gyroscope output of the MARG sensors, φ_{k1}Represent root during moment k1 The course angle obtained according to the MARG sensors, Δ T represent the time of measuring interval of the MARG sensors, and Q and R are represented respectively The covariance matrix of process noise and measurement noise, K_{k}For Kalman filter gain,And P_{k}Represent varivance matrix, φ_{PF} The course angle of particle filter output, first φ during Kalman filtering when being positioned for the last time_{PF}For designated value.
Further, said system can also have the characteristics that, the progress particle filter includes：
Particle is initialized, the probability density function of particle is initialized using Gaussian Profile；
According to the course angle and velocity information, and the initial estimated location, predict the user equipment in next step Status information：
The weight of each particle and normalization are calculated, it is as follows：
Carry out particle resampling, the particle as particle filter next time；
Wherein, described [x_{k}, y_{k}]^{T}For the state vector of each particle, T_{s}Represent the last access based on WLAN The time interval of the positioning and the positioning of this access point based on WLAN of point, φ_{k}Represent the course angle, v_{k}Represent The velocity information, [η_{x}, η_{y}]^{T}Acceleration is represented, is simulated with the Gaussian noise of zeromean, variance is sensed by the MARG The metrical information estimation of device,Currently to input the state value of the particle filter,Represent shape of ith of particle in moment k State value, σ represent the noise variance of signal strength measurement.
The present invention utilizes MARG (Magnetic, Angular Rate, and Gravity, magnetometer, gyroscope and acceleration Degree meter) sensor auxiliary WLAN (WirelessLAN) alignment system, devise one and be based on particle filter and Kalman filtering Data anastomosing algorithm, the blending algorithm makes full use of WLAN (WirelessLAN) and the complementary characteristic of MARG location technologies, Position error and the accumulated error as caused by sensor noise as caused by floating receive information intensity are effectively corrected, is realized The WLAN/MARG integrated positioning systems of one lowcost and highprecision.
Brief description of the drawings
Fig. 1 is alignment system block diagram of the embodiment of the present invention；
Fig. 2 is localization method flow chart of the embodiment of the present invention.
Embodiment
For the object, technical solutions and advantages of the present invention are more clearly understood, below in conjunction with accompanying drawing to the present invention Embodiment be described in detail.It should be noted that in the case where not conflicting, in the embodiment and embodiment in the application Feature can mutually be combined.
The embodiment of the present invention provides a kind of localization method, including：
Access point based on WLAN is positioned, and obtains the initial estimated location of user equipment；
Obtain the course angle and velocity information of the user equipment；
The initial estimated location is modified according to the course angle and velocity information, obtains final position information.
Wherein, the access point based on WLAN, which carries out positioning, includes：
Reference point is chosen, the signal intensity from each access point is measured in each reference point, by the position of the reference point, institute Mark/the position for stating signal intensity and corresponding access point is stored into database；Respectively connect around user equipment measurement to be positioned The signal intensity of access point, searching data storehouse obtain corresponding to refer to point set, with it is described with reference to point set match determine use The initial estimated location of family equipment.
Wherein, it is described to determine that the initial estimated location of the user equipment includes with reference to point set match：
The reference point of the Euclidean distance minimum of m received signal strength of selection, uses the line of the position of the m reference point Property initial estimated location of the weighted sum as the user equipment, the m is more than or equal to 1.
Wherein, the course angle for obtaining the user equipment and velocity information include：
The course angle and velocity information are obtained according to the metrical information of MARG sensors.
Wherein, it is described to be included according to the metrical information of the MARG sensors acquisition course angle：
First course angle φ is obtained according to the metrical information of the magnetometer of the MARG sensors_{mag}, passed according to the MARG The metrical information of the gyroscope of sensor obtains the second course angleAccording to obtaining first course angle and the second course angle The course angle φ of user equipment：
Wherein, the W is default weighted value, 0≤W≤1.
Wherein methods described also includes, and the rolling of the user equipment is obtained according to the metrical information of the MARG sensors Angle, the angle of pitch；
It is described according to the course angle and velocity information the initial estimated location is modified including：
The roll angle obtained according to the metrical information of the MARG sensors, the angle of pitch, course angle and the speed are believed The input as Kalman filter is ceased, Kalman filtering is carried out, exports new course angle and velocity information；
The course angle and velocity information and the initial estimated location that the Kalman filter is exported are as particle The input of wave filter, carry out particle filter, output position information, course angle and velocity information, using the positional information of output as The final position information of the user equipment.
It is of course also possible to without Kalman filtering, directly the course angle and speed that are obtained according to MARG sensors are believed Breath carries out particle filter, obtains final positional information.
In addition to Kalman filtering and particle filter, Bayesian filter, complementary filter, spreading kalman filter also can be used The blending algorithms such as ripple, Federated Kalman Filtering.
Wherein, the progress Kalman filtering includes：
The state onestep prediction of Kalman filtering is carried out,
Calculate prediction varivance matrix
Calculate filtering gain matrix
Carry out state estimation
Calculate estimation error variance
Wherein, it is describedRepresent the angular speed of the gyroscope output of the MARG sensors, φ_{k1}Represent root during moment k1 The course angle obtained according to MARG sensors, Δ T represent the time of measuring interval of the MARG sensors, and Q and R represent process respectively The covariance matrix of noise and measurement noise, K_{k}For Kalman filter gain,And P_{k}Represent varivance matrix, φ_{PF}To be upper The course angle that particle filter exports during onetime positioning, first φ during Kalman filtering_{PF}For designated value.
Wherein, the progress particle filter includes：
, it is necessary to initialize particle when carrying out particle filter first, the probability density of particle is initialized using Gaussian Profile Function；
According to the course angle and velocity information, and the initial estimated location, predict the user equipment in next step Status information：
The weight of each particle and normalization are calculated, it is as follows：
Carry out particle resampling, the particle as particle filter next time；
Wherein, described [x_{k}, y_{k}]^{T}For the state vector of each particle, T_{s}Represent the last access based on WLAN The time interval of the positioning and the positioning of this access point based on WLAN of point, φ_{k}Represent the course angle, v_{k}Represent The velocity information, [η_{x}, η_{y}]^{T}Acceleration is represented, is simulated with the Gaussian noise of zeromean, variance is sensed by the MARG The metrical information estimation of device,Currently to input the state value of the particle filter,Represent shape of ith of particle in moment k State value, σ represent the noise variance of signal strength measurement.
Alignment system provided in an embodiment of the present invention realized shown in block diagram as accompanying drawing 1, including：WLAN locating modules, sensing Device locating module and Fusion Module, wherein：
WLAN locating modules, positioned for the access point based on WLAN, obtain the initial of user equipment and estimate Count position；
Sensor locating module, for obtaining the course angle and velocity information of the user equipment；
Fusion Module, for being modified according to the course angle and velocity information to the initial estimated location, obtain Final position information.
The WLAN locating modules include：
Database, for being stored in the signal intensity from each access point of each reference point measurement, the position of the reference point Put and correspond to mark/position of access point；
RSS measuring units, for measuring the signal intensity of each access point around user equipment measurement to be positioned；
Positioning unit, for the signal intensity measured according to the RSS measuring units, searching data storehouse obtains corresponding join Examination point set, with the initial estimated location for reference to point set match determination user equipment.
The positioning unit determines that the initial estimated location of the user equipment includes with reference to point set match：
The positioning unit selects the minimum reference point of the Euclidean distance of m received signal strength, uses described m reference The linear weighted function of the position of point and the initial estimated location as the user equipment, the m are more than or equal to 1.
The sensor locating module includes：MARG sensors and data processing unit, wherein：
The MARG sensors are used for, and the user equipment is measured, and obtain metrical information；
The data processing unit is used for, and the course angle and speed are obtained according to the metrical information of the MARG sensors Information.
The data processing unit obtains the course angle according to the metrical information of MARG sensors to be included：
The data processing unit obtains the first course angle according to the metrical information of the magnetometer of the MARG sensors φ_{mag}, the second course angle is obtained according to the metrical information of the gyroscope of the MARG sensorsAccording to first course angle The course angle φ of the user equipment is obtained with the second course angle：
Wherein, the W is default weighted value, 0≤W≤1.
The Fusion Module includes：Kalman filter and particle filter, wherein：
The data processing unit is additionally operable to, and the user equipment is obtained according to the metrical information of the MARG sensors Roll angle, the angle of pitch；
The Kalman filter is used for, by it is described according to the metrical information of the MARG sensors obtain roll angle, The angle of pitch, course angle and velocity information input as the state value of Kalman filter, carry out Kalman filtering, export new boat To angle and velocity information；
The particle filter is used for, the course angle and velocity information that the Kalman filter is exported, and described Initial estimated location inputs as state value, carries out particle filter, output position information, course angle and velocity information, will export Final position information of the positional information as the user equipment.
Particle filter and the specific method of Kalman filtering are referring to embodiment of the method.
The present invention is further illustrated below by concrete application example.
1) WLAN is positioned
Fingerprinting (fingerprint characteristic) position fixing process based on WLAN is broadly divided into training and two stages of positioning As shown in Figure 2.
(1) training stage
It, which is aimed at, establishes a location fingerprint identification database.Groundwork is respectively referred in collection areaofinterest The fingerprint feature information of point (Reference Point, RP) position  RSS (Received Signal Strength, is received Signal intensity).As shown in figure 1, measurement comes from different AP to mobile subscriber (Mobile User, MU) in each reference point successively RSS values, and by corresponding MAC Address and latitude and longitude coordinates information storage into database, until traversal areaofinterest in institute Some reference points, this process complete RSS measurement and the foundation of RP databases.
Specifically：
23 rice, outdoor selection 811 rice are chosen in localization region selection reference point (RP), the spacing room of reference；
The signal intensity samples of WLAN access points, AP addresses or mark are measured in reference point by space diversity technology The positional information of (such as MAC Address) and corresponding reference point (can be represented) with latitude and longitude coordinates；
Filtering is weighted by signal intensity samples, obtains reference point RP finger print information, finger print information number is arrived in storage According in storehouse.
(2) positioning stage
(1) mobile subscriber measures RSS (Received Signal Strength, the reception signal of surrounding WLAN access points Intensity), and data are filtered with processing；
(2) set of fingerprint information corresponding to being searched according to MAC Address in fingerprint database；
(3) matching primitives are carried out with the reference point in set, to determine the position of mobile subscriber.Matching primitives principle is to connect The Euclidean distance of signal intensity is received, the minimum reference point of m Euclidean distance is selected, with the linear weighted function of the coordinate of these reference points With the position for representing mobile subscriber, referred to as initial estimated location.
Except abovementioned matching primitives principle, also the matching based on nerual network technique, based on histogram probabilistic synchronization algorithm Matching principle, matching principle based on SVMs etc..
One example is as follows：
Mobile subscriber measures surrounding AP RSS, and it is carried out into matching meter with the RSS vectors being stored in advance in database Calculate, matching principle is the Euclidean distance of received signal strength, as shown in formula (1)：
Wherein D_{j}Represent the Euclidean distance or similarity of signal intensity between reference point j and mobile subscriber, D_{j}It is smaller to show two Apart nearer between person；Rss=(rss_{i}, rss_{2}..., rss_{n}) vector representation mobile subscriber currently measure to n AP's RSS；RSS=(RSS_{j1}, RSS_{j2}..., RSS_{jn}) vector representation in jth reference point storage to the RSS in database, represents ginseng Examination point j finger print information.
Minimum reference point (the x of m Euclidean distance of selection_{1}, x_{2}..., x_{m}), with these with reference to point coordinates linear weighted function and Represent the current position coordinates x of mobile subscriber_{0}=(x_{0}, y_{0}), calculation formula is as follows：
Wherein w_{k}Expression reference point k weight, calculation formula are as follows：
Abovementioned WLAN localization methods are merely illustrative, and other WLAN localization methods are also applicable in invention.
2) MARG sensing datas are handled
The processing of MARG sensing datas is the information such as the acceleration and angular speed provided using sensor, to obtain carrier The information such as attitude angle and relative position.The intelligence of threedimensional gyroscope, threedimensional accelerometer and threedimensional magnetometer will be integrated with first For energy terminal definitions to an xyz coordinate system, commonly referred to as carrier coordinate system, the center of gravity for taking carrier is carrier coordinate system origin, Three axles coincide with the longitudinal axis, transverse axis and vertical pivot of carrier respectively.Corresponding absolute coordinate system is commonly referred to as XYZ navigation East, north, day are pointed in coordinate system, X, Y, Z axis distribution, it then follows the righthand rule.
Such as：Roll anglePitching angle theta and course angle φ represent what carrier coordinate system rotated around xaxis, yaxis and zaxis respectively Corner, for representing orientation of the carrier coordinate system relative to navigational coordinate system, the also referred to as attitude angle of carrier.Roll angle and pitching The calculation formula at angle is as follows：
Wherein WithRepresent accelerometer under carrier coordinate system along the output valve of x, y and z axes.The output of course angle It can be obtained by the output of magnetometer or the output of gyroscope.With magnetometer m=[m_{x}, m_{y}, m_{z}]^{T}When asking for course angle, Need to make the zaxis of carrier coordinate system and navigational coordinate system Z axis align by spin matrix, then ask for navigating by water angle, formula is as follows：
M '=R ' m (7)
Wherein m '_{x}With m '_{y}Represent earth magnetic field intensity component along the x ' axles and y ' axis components after alignment.
Accelerometer measures obtain carrier in x, y, the acceleration in zaxis, by acceleration obtain carrier roll angle and The angle of pitch, and the velocity information of carrier.
Magnetometer survey obtains magnetic field of the earth in x, y, the magnetic field strength component in zaxis, obtains the course angle of carrier, is referred to as First course angle.
Gyroscope measures carrier angular velocity information, calculates roll angle, the angle of pitch and the second course angle for obtaining carrier.
Navigated by complementary filter to merge the second course angle that gyroscope calculates and first that magnetometer calculates To angle, course angle is obtained.It is of course also possible to without complementary filter, the first course angle or the second course angle are directly used.When So, course angle can also be obtained using GPS technology.
The course angle calculated using gyroscope corrects the output of magnetometer, and formula is as follows：
WhereinRepresent the course angle calculated by gyroscope, ω_{k}Angular speed of the gyroscope in the k periods is represented, W is institute The complementary weight of the complementary filter of design, 0≤W≤1.
3) the particle filter design based on WLAN and MARG
Particle filter step：
(1) course angle of Kalman filter output and velocity information are inputted as the state value of particle filter；
(2) particle is initialized：The probability density function of particle is initialized using Gaussian Profile, average is target initial shape State, the original state include course angle and velocity information；
(3) predict：Using course angle and velocity information, and WLAN positioning results, particle filter prediction target is next The status information of step；
(4) granular Weights Computing and normalization：The power of each particle is asked for by measurement model and present measured value Weight, when particle position estimated state current closer to target, the weight that particle obtains is bigger；
(5) resampling：New particle is produced according to posterior probability density function, solves the problems, such as sample degeneracy.
Particle filter output position information, course angle information, velocity information.Position of the positional information as mobile subscriber Estimation.
Specific filtering method is as follows：
Particle filter be particle collection that is being randomly selected using one group from probability density function and attaching related weights come Approach posterior probability density function：
Wherein x_{k}Represent target in moment k state vector, z_{0:k}The measurement value sequence before moment k+1 is represented,Table Show ith of particle or sample point,For its weight, N is population.The particle filter used in the system is divided into following four Step：
1) initialize
According to initial probability density function Pr (x_{0}) produce N number of particleWherein Pr (x_{0}) use Gauss Distribution, average is target original state, and the value can be set as needed.Carrying out particle filter first needs to be initialized, The particle obtained before followup use after a filtering resampling.
2) predict
With reference to MARG sensing data results, particle filter predicts the status information (x of target next step_{k+1}, y_{k+1}), Formula is as follows：
Wherein [x_{k}, y_{k}]^{T}For the state vector of each particle, T_{s}Represent that 1 WLAN positioning of kth positions with kth time WLAN Time interval, φ_{k}Represent target of the MARG sensing datas by Kalman filtering acquisition around the anglec of rotation (the i.e. course of zaxis Angle), v_{k}Represent the target velocity obtained after the processing of MARG sensing datas, [η_{x}, η_{y}]^{T}The acceleration of target is represented, with zero The Gaussian noise simulation of average, variance can be estimated by MARG sensing datas.
3) weight calculation and normalization
The weight of particle is asked for by measurement model and Current observation value：
Wherein z_{k}The RSS that target currently measures is represented,The current status information (position, course angle) of target is represented,Table Show status information of ith of particle in moment k, σ represents measurement noise variance, is selected according to the variance of the floatings of RSS in practice. Formula (13) represents that when particle position is closer to target current estimated location the weight that particle obtains is bigger, accurate so as to obtain Posterior probability distribution.
4) resampling
Resampling is the key of particle filter, according to probability density function Pr (x_{k}z_{k}) produce N number of new particleTo solve the problems, such as sample degeneracy, a kind of method for resampling is as follows：
4) the Kalman filtering design based on WLAN and MARG
The fine or not heavy dependence course angle φ of the particle filter degree of accuracy, tire out to further eliminate because gyroscope is present Product error and magnetometer are easily disturbed by local magnetic field around, cause the course angle φ errors obtained after MARG data processings, this In embodiment, a Kalman filter is may also provide, (including particle filter is defeated for the attitude information obtained using particle filter The positional information gone out, bearer rate information, course angle information etc.) course angle φ is corrected, so as to obtaining reliable and stable course Angle information.
The step of Kalman filtering：
(1) roll angle, the angle of pitch and the course angle that MARG sensors obtain, and bearer rate information, as Kalman The input of filter state value.
(2) onestep prediction in the state of progress Kalman filtering；
(3) prediction varivance matrix is calculated；
(4) filtering gain matrix are calculated；
(5) state estimation；
(6) estimation error variance calculates.
Kalman filtering algorithm such as following formula：
WhereinRepresent the angular speed of gyroscope output, φ_{k1}The course angle predicted during moment k1 is represented, Δ T represents MARG The time of measuring interval of sensor, Q and R represent process noise and the covariance matrix of measurement noise, K respectively_{k}Filtered for Kalman Ripple device gain,And P_{k}Represent varivance matrix, φ_{PF}For the course angle of particle filter estimation.
The WLAN/MARG integrated positioning systems based on data fusion that the embodiment of the present invention proposes, the system utilize MARG Sensor obtains the information such as the speed of mobile target, posture, is melted by data such as complementary filter, Kalman filtering and particle filters Hop algorithm improves WLAN positioning precisions, realizes a lowcost and highprecision WLAN/MARG integrated positioning system.
One of ordinary skill in the art will appreciate that all or part of step in the above method can be instructed by program Related hardware is completed, and described program can be stored in computerreadable recording medium, such as readonly storage, disk or CD Deng.Alternatively, all or part of step of abovedescribed embodiment can also be realized using one or more integrated circuits.Accordingly Ground, each module/unit in abovedescribed embodiment can be realized in the form of hardware, can also use the shape of software function module Formula is realized.The present invention is not restricted to the combination of the hardware and software of any particular form.
Claims (16)
 A kind of 1. localization method, it is characterised in that including：Access point based on WLAN is positioned, and obtains the initial estimated location of user equipment；Obtain the course angle and velocity information of the user equipment；The initial estimated location is modified according to the course angle and velocity information, obtains final position information；It is described according to the course angle and velocity information the initial estimated location is modified including：Kalman filtering is carried out to the course angle and velocity information, exports new course angle and velocity information；Particle filter is carried out to the new course angle and velocity information and the initial estimated location, by the position of output Information is as final position information.
 2. the method as described in claim 1, it is characterised in that the access point based on WLAN carries out positioning bag Include：Reference point is chosen, the signal intensity from each access point is measured in each reference point, by the position of the reference point, the letter Mark/position of number intensity and corresponding access point is stored into database；Each access point around user equipment measurement to be positioned Signal intensity, searching data storehouse obtain corresponding to refer to point set, match determination user with the reference point set sets Standby initial estimated location.
 3. method as claimed in claim 2, it is characterised in that described to determine that the user sets with reference to point set match Standby initial estimated location includes：The reference point of the Euclidean distance minimum of m received signal strength of selection, added using the linear of position of the m reference point Power and the initial estimated location as the user equipment, the m are more than or equal to 1.
 4. the method as described in claim 1, it is characterised in that the course angle and velocity information for obtaining the user equipment Including：The course angle and velocity information are obtained according to the metrical information of MARG sensors.
 5. method as claimed in claim 4, it is characterised in that described that the boat is obtained according to the metrical information of MARG sensors Include to angle：First course angle φ is obtained according to the metrical information of the magnetometer of the MARG sensors_{mag}, according to the MARG sensors Gyroscope metrical information obtain the second course angleThe use is obtained according to first course angle and the second course angle The course angle φ of family equipment：<mrow> <mi>&phi;</mi> <mo>=</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo></mo> <mi>W</mi> <mo>)</mo> </mrow> <msubsup> <mi>&phi;</mi> <mi>k</mi> <mrow> <mi>g</mi> <mi>y</mi> <mi>r</mi> </mrow> </msubsup> <mo>+</mo> <msub> <mi>W&phi;</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>g</mi> </mrow> </msub> </mrow>Wherein, the W is default weighted value, 0≤W≤1.
 6. the method as described in claim 4 or 5, it is characterised in thatMethods described also includes, and roll angle, the pitching of the user equipment are obtained according to the metrical information of the MARG sensors Angle；It is described according to the course angle and velocity information the initial estimated location is modified including：The roll angle obtained according to the metrical information of the MARG sensors, the angle of pitch, course angle and the velocity information are made For the input of Kalman filter, Kalman filtering is carried out, exports new course angle and velocity information；The course angle and velocity information and the initial estimated location that the Kalman filter is exported are as particle filter The input of device, carry out particle filter, output position information, course angle and velocity information, using the positional information of output as described in The final position information of user equipment.
 7. method as claimed in claim 6, it is characterised in that the progress Kalman filtering includes：The state onestep prediction of Kalman filtering is carried out,Calculate prediction varivance matrixCalculate filtering gain matrixCarry out state estimationCalculate estimation error varianceWherein, it is describedRepresent the angular speed of the gyroscope output of the MARG sensors, φ_{k1}According to institute during expression moment k1 The course angle of MARG sensors acquisition is stated, Δ T represents the time of measuring interval of the MARG sensors, and Q and R represent process respectively The covariance matrix of noise and measurement noise, K_{k}For Kalman filter gain,And P_{k}Represent varivance matrix, φ_{PF}For The course angle of particle filter output, first φ during Kalman filtering during last time positioning_{PF}For designated value.
 8. method as claimed in claim 6, it is characterised in that the progress particle filter includes：, it is necessary to initialize particle when carrying out particle filter first, the probability density function of particle is initialized using Gaussian Profile；According to the course angle and velocity information, and the initial estimated location, the shape of the user equipment next step is predicted State information：<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>&CenterDot;</mo> <mi>cos</mi> <mrow> <mo>(</mo> <msub> <mi>&phi;</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>&CenterDot;</mo> <mi>sin</mi> <mrow> <mo>(</mo> <msub> <mi>&phi;</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mi>k</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mi>k</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>z</mi> <mi>k</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mfrac> <msubsup> <mi>T</mi> <mi>s</mi> <mn>2</mn> </msubsup> <mn>2</mn> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mfrac> <msubsup> <mi>T</mi> <mi>s</mi> <mn>2</mn> </msubsup> <mn>2</mn> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>&eta;</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&eta;</mi> <mi>y</mi> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow>The weight of each particle and normalization are calculated, it is as follows：<mrow> <msubsup> <mi>&omega;</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>i</mi> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msqrt> <mrow> <mn>2</mn> <mi>&pi;</mi> </mrow> </msqrt> <mi>&sigma;</mi> </mrow> </mfrac> <mi>exp</mi> <mo>&lsqb;</mo> <mo></mo> <mfrac> <mrow> <mo></mo> <mo></mo> <msubsup> <mi>x</mi> <mi>k</mi> <mi>z</mi> </msubsup> <mo></mo> <msubsup> <mi>x</mi> <mi>k</mi> <mi>i</mi> </msubsup> <mo></mo> <mo></mo> </mrow> <mrow> <mn>2</mn> <msup> <mi>&sigma;</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>&rsqb;</mo> </mrow><mrow> <msubsup> <mi>&omega;</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>i</mi> </msubsup> <mo>=</mo> <mfrac> <msubsup> <mi>&omega;</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>i</mi> </msubsup> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mi>&omega;</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>j</mi> </msubsup> </mrow> </mfrac> </mrow>Carry out particle resampling, the particle as particle filter next time；Wherein, described [x_{k},y_{k}]^{T}For the state vector of each particle, T_{s}Represent the last access point based on WLAN The time interval of positioning and the positioning of this access point based on WLAN, φ_{k}Represent the course angle, v_{k}Described in expression Velocity information, [η_{x},η_{y}]^{T}Acceleration is represented, is simulated with the Gaussian noise of zeromean, variance is by the MARG sensors Metrical information estimation,To input the state value of the particle filter,Represent state value of ith of particle in moment k, σ Represent the noise variance of signal strength measurement.
 A kind of 9. alignment system, it is characterised in that including：WLAN locating modules, positioned for the access point based on WLAN, obtain the initial estimation position of user equipment Put；Sensor locating module, for obtaining the course angle and velocity information of the user equipment；Fusion Module, for being modified according to the course angle and velocity information to the initial estimated location, obtain final Positional information；The Fusion Module includes：Kalman filter and particle filter, wherein：The Kalman filter is used for, and the course angle and velocity information to the acquisition carry out Kalman filtering, export new Course angle and velocity information；The particle filter is used for, and grain is carried out to the new course angle and velocity information, and the initial estimated location Son filtering, using the positional information of output as final position information.
 10. system as claimed in claim 9, it is characterised in that the WLAN locating modules include：Database, for being stored in the signal intensity from each access point of each reference point measurement, the position of the reference point and Mark/position of corresponding access point；RSS measuring units, for measuring the signal intensity of each access point around user equipment measurement to be positioned；Positioning unit, for the signal intensity measured according to the RSS measuring units, reference point corresponding to the acquisition of searching data storehouse Set, with the initial estimated location for reference to point set match determination user equipment.
 11. system as claimed in claim 10, it is characterised in that the positioning unit with reference to point set with carrying out matching determination The initial estimated location of the user equipment includes：The positioning unit selects the minimum reference point of the Euclidean distance of m received signal strength, uses the m reference point The linear weighted function of position and the initial estimated location as the user equipment, the m are more than or equal to 1.
 12. system as claimed in claim 9, it is characterised in that the sensor locating module includes：MARG sensor sums According to processing unit, wherein：The MARG sensors are used for, and the user equipment is measured, and obtain metrical information；The data processing unit is used for, and the course angle and velocity information are obtained according to the metrical information of MARG sensors.
 13. system as claimed in claim 12, it is characterised in that the data processing unit is according to the measurements of MARG sensors Course angle includes described in acquisition of information：The data processing unit obtains the first course angle φ according to the metrical information of the magnetometer of the MARG sensors_{mag}, root The second course angle is obtained according to the metrical information of the gyroscope of the MARG sensorsAccording to first course angle and second Course angle obtains the course angle φ of the user equipment：<mrow> <mi>&phi;</mi> <mo>=</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo></mo> <mi>W</mi> <mo>)</mo> </mrow> <msubsup> <mi>&phi;</mi> <mi>k</mi> <mrow> <mi>g</mi> <mi>y</mi> <mi>r</mi> </mrow> </msubsup> <mo>+</mo> <msub> <mi>W&phi;</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>g</mi> </mrow> </msub> </mrow>Wherein, the W is default weighted value, 0≤W≤1.
 14. the system as described in claim 12 or 13, it is characterised in thatThe data processing unit is additionally operable to, and the rolling of the user equipment is obtained according to the metrical information of the MARG sensors Angle, the angle of pitch；The Kalman filter is used for, by the roll angle obtained according to the metrical information of the MARG sensors, pitching Angle, course angle and velocity information input as the state value of Kalman filter, carry out Kalman filtering, export new course angle And velocity information；The particle filter is used for, the course angle and velocity information that the Kalman filter is exported, and described initial Estimated location inputs as state value, carries out particle filter, output position information, course angle and velocity information, by the position of output Confidence ceases the final position information as the user equipment.
 15. system as claimed in claim 14, it is characterised in that the Kalman filter, which carries out Kalman filtering, to be included：The state onestep prediction of Kalman filtering is carried out,Calculate prediction varivance matrixCalculate filtering gain matrixCarry out state estimationCalculate estimation error varianceWherein, it is describedRepresent the angular speed of the gyroscope output of the MARG sensors, φ_{k1}According to institute during expression moment k1 The course angle of MARG sensors acquisition is stated, Δ T represents the time of measuring interval of the MARG sensors, and Q and R represent process respectively The covariance matrix of noise and measurement noise, K_{k}For Kalman filter gain,And P_{k}Represent varivance matrix, φ_{PF}For The course angle of particle filter output, first φ during Kalman filtering during last time positioning_{PF}For designated value.
 16. system as claimed in claim 14, it is characterised in that the progress particle filter includes：Particle is initialized, the probability density function of particle is initialized using Gaussian Profile；According to the course angle and velocity information, and the initial estimated location, the shape of the user equipment next step is predicted State information：<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>&CenterDot;</mo> <mi>cos</mi> <mrow> <mo>(</mo> <msub> <mi>&phi;</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>&CenterDot;</mo> <mi>sin</mi> <mrow> <mo>(</mo> <msub> <mi>&phi;</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mi>k</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mi>k</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>z</mi> <mi>k</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mfrac> <msubsup> <mi>T</mi> <mi>s</mi> <mn>2</mn> </msubsup> <mn>2</mn> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mfrac> <msubsup> <mi>T</mi> <mi>s</mi> <mn>2</mn> </msubsup> <mn>2</mn> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>&eta;</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&eta;</mi> <mi>y</mi> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow>The weight of each particle and normalization are calculated, it is as follows：<mrow> <msubsup> <mi>&omega;</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>i</mi> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msqrt> <mrow> <mn>2</mn> <mi>&pi;</mi> </mrow> </msqrt> <mi>&sigma;</mi> </mrow> </mfrac> <mi>exp</mi> <mo>&lsqb;</mo> <mo></mo> <mfrac> <mrow> <mo></mo> <mo></mo> <msubsup> <mi>x</mi> <mi>k</mi> <mi>z</mi> </msubsup> <mo></mo> <msubsup> <mi>x</mi> <mi>k</mi> <mi>i</mi> </msubsup> <mo></mo> <mo></mo> </mrow> <mrow> <mn>2</mn> <msup> <mi>&sigma;</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>&rsqb;</mo> </mrow><mrow> <msubsup> <mi>&omega;</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>i</mi> </msubsup> <mo>=</mo> <mfrac> <msubsup> <mi>&omega;</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>i</mi> </msubsup> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mi>&omega;</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>j</mi> </msubsup> </mrow> </mfrac> </mrow>Carry out particle resampling, the particle as particle filter next time；Wherein, described [x_{k},y_{k}]^{T}For the state vector of each particle, T_{s}Represent the last access point based on WLAN The time interval of positioning and the positioning of this access point based on WLAN, φ_{k}Represent the course angle, v_{k}Described in expression Velocity information, [η_{x},η_{y}]^{T}Acceleration is represented, is simulated with the Gaussian noise of zeromean, variance is by the MARG sensors Metrical information estimation,Currently to input the state value of the particle filter,Represent state of ith of particle in moment k Value, σ represent the noise variance of signal strength measurement.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201210111075.8A CN103379619B (en)  20120416  20120416  A kind of localization method and system 
Applications Claiming Priority (2)
Application Number  Priority Date  Filing Date  Title 

CN201210111075.8A CN103379619B (en)  20120416  20120416  A kind of localization method and system 
PCT/CN2013/073085 WO2013155919A1 (en)  20120416  20130322  Positioning method and system 
Publications (2)
Publication Number  Publication Date 

CN103379619A CN103379619A (en)  20131030 
CN103379619B true CN103379619B (en)  20171128 
Family
ID=49382898
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201210111075.8A CN103379619B (en)  20120416  20120416  A kind of localization method and system 
Country Status (2)
Country  Link 

CN (1)  CN103379619B (en) 
WO (1)  WO2013155919A1 (en) 
Families Citing this family (17)
Publication number  Priority date  Publication date  Assignee  Title 

CN103644911A (en) *  20131127  20140319  南京城际在线信息技术有限公司  Gyroscope assisted positioning method 
CN103674785B (en) *  20131204  20170201  杭州腾海科技有限公司  Floated online tracking monitoring device and monitoring method for flow velocity and flow direction of superficial water flow 
CN104853431A (en) *  20140218  20150819  中兴通讯股份有限公司  Geographical location information positioning method and device 
CN104105060B (en) *  20140718  20180213  南京猫酷科技股份有限公司  Localization method and equipment 
CN104482933B (en) *  20141203  20171229  北京航空航天大学  A kind of method based on particle filter reckoning and WLAN integrated positioning 
CN104501796B (en) *  20141216  20170308  重庆邮电大学  A kind of interior WLAN/MEMS merges across floor 3dimensional localization method 
CN104635203B (en) *  20150212  20170419  国家无线电监测中心  Radio interference source directionfinding and positioning method based on particle filter algorithm 
CN104796866B (en) *  20150506  20180911  北京我联科技有限公司  Indoor orientation method and device 
CN106325287A (en) *  20150614  20170111  复旦大学  Intelligent mower straight line walking control system based on inertial/magnetic sensor MARG attitude detection 
CN105357753B (en) *  20151016  20180918  杭州电子科技大学  A kind of indoor orientation method based on multimodality fusion recursive iteration 
CN105510870B (en) *  20151219  20180518  长安大学  A kind of smart machine indoor locating system and method 
CN105589064B (en) *  20160108  20180323  重庆邮电大学  WLAN location fingerprint database is quickly established and dynamic update system and method 
CN105768480A (en) *  20160225  20160720  北京石油化工学院  Positioning method, device and system for suitcase 
CN105872981B (en) *  20160330  20181225  河海大学常州校区  Indoor orientation method based on signal synthesis and artificial neural network 
CN106535133A (en) *  20161111  20170322  天津大学  Indoor telephone traffic accurate location method based on machine learning in cellular network 
CN107576330A (en) *  20170907  20180112  西北大学  A kind of localization method of the indoor dynamic sensing strategy based on WLAN fingerprints 
CN110099443B (en) *  20190513  20200407  电子科技大学  Load balancing method for node tracking in wireless sensor network 
Citations (3)
Publication number  Priority date  Publication date  Assignee  Title 

CN101208880A (en) *  20050331  20080625  日本电气株式会社  Mobile radio terminal and its moving speed detection method 
CN101285686A (en) *  20080529  20081015  中国农业大学  Agricultural machines navigation hierarchical positioning process and system 
WO2010059934A2 (en) *  20081121  20100527  Qualcomm Incorporated  Wireless position determination using adjusted round trip time measurements 
Family Cites Families (1)
Publication number  Priority date  Publication date  Assignee  Title 

KR100689517B1 (en) *  20041213  20070302  삼성전자주식회사  System of assisted gps using network time protocol server and method measuring location of mobile phone thereof 

2012
 20120416 CN CN201210111075.8A patent/CN103379619B/en active IP Right Grant

2013
 20130322 WO PCT/CN2013/073085 patent/WO2013155919A1/en active Application Filing
Patent Citations (3)
Publication number  Priority date  Publication date  Assignee  Title 

CN101208880A (en) *  20050331  20080625  日本电气株式会社  Mobile radio terminal and its moving speed detection method 
CN101285686A (en) *  20080529  20081015  中国农业大学  Agricultural machines navigation hierarchical positioning process and system 
WO2010059934A2 (en) *  20081121  20100527  Qualcomm Incorporated  Wireless position determination using adjusted round trip time measurements 
NonPatent Citations (1)
Title 

大区域宽带自组织网络中定位信息的作用（上）;马正新;《数字通信世界》;20110615;参见期刊第68页 * 
Also Published As
Publication number  Publication date 

CN103379619A (en)  20131030 
WO2013155919A1 (en)  20131024 
Similar Documents
Publication  Publication Date  Title 

Zhuang et al.  Evaluation of two WiFi positioning systems based on autonomous crowdsourcing of handheld devices for indoor navigation  
US9551561B2 (en)  Determining location using magnetic fields from AC power lines  
Liu et al.  Improving positioning accuracy using GPS pseudorange measurements for cooperative vehicular localization  
Li et al.  An improved inertial/wifi/magnetic fusion structure for indoor navigation  
Li et al.  How feasible is the use of magnetic field alone for indoor positioning?  
CN105588566B (en)  A kind of indoor locating system merged based on bluetooth with MEMS and method  
Jekabsons et al.  An analysis of WiFi based indoor positioning accuracy  
CN102435188B (en)  Monocular vision/inertia autonomous navigation method for indoor environment  
CN104121905B (en)  Course angle obtaining method based on inertial sensor  
CN102426019B (en)  Unmanned aerial vehicle scene matching auxiliary navigation method and system  
Wang et al.  Urban positioning on a smartphone: Realtime shadow matching using GNSS and 3D city models  
CA2701142C (en)  Apparatus and method for determining a coincidence of a position with a reference position  
CN103561462B (en)  Indoor positioning system and method totally based on smart mobile terminal platform  
CN103363988B (en)  A kind of method utilizing intelligent mobile phone sensor to realize the navigation of earth magnetism indoor positioning  
CN104819716A (en)  Indoor and outdoor personal navigation algorithm based on INS/GPS (inertial navigation system/global position system) integration of MEMS (microelectromechanical system)  
CN105547305B (en)  A kind of pose calculation method based on wireless location and laser map match  
Mohssen et al.  It's the human that matters: accurate user orientation estimation for mobile computing applications  
CN104469942B (en)  A kind of indoor orientation method based on HMM  
CN104849740A (en)  Indoor and outdoor seamless positioning system integrated with satellite navigation and bluetooth technology, and method thereof  
CN103282741B (en)  Position determination using horizontal angles  
CN106556854B (en)  A kind of indoor and outdoor navigation system and method  
CN104869639B (en)  A kind of indoor orientation method and device  
CN107194969A (en)  The pick up calibration determined based on end point and location estimation  
CN103744052B (en)  A kind of double star being applied to aerial target location surveys timeofarrival direction finding method and device  
CN106574975A (en)  Trajectory matching using peripheral signal 
Legal Events
Date  Code  Title  Description 

PB01  Publication  
C06  Publication  
SE01  Entry into force of request for substantive examination  
C10  Entry into substantive examination  
GR01  Patent grant  
GR01  Patent grant 