Nothing Special   »   [go: up one dir, main page]

CN116448111A - Pedestrian indoor navigation method, device and medium based on multi-source information fusion - Google Patents

Pedestrian indoor navigation method, device and medium based on multi-source information fusion Download PDF

Info

Publication number
CN116448111A
CN116448111A CN202310210607.1A CN202310210607A CN116448111A CN 116448111 A CN116448111 A CN 116448111A CN 202310210607 A CN202310210607 A CN 202310210607A CN 116448111 A CN116448111 A CN 116448111A
Authority
CN
China
Prior art keywords
pedestrian
state
matrix
source information
information
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202310210607.1A
Other languages
Chinese (zh)
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.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
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 South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN202310210607.1A priority Critical patent/CN116448111A/en
Publication of CN116448111A publication Critical patent/CN116448111A/en
Pending legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0248Filters characterised by a particular frequency response or filtering method
    • H03H17/0255Filters based on statistics
    • H03H17/0257KALMAN filters
    • 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S11/00Systems for determining distance or velocity not using reflection or reradiation
    • G01S11/02Systems for determining distance or velocity not using reflection or reradiation using radio waves
    • G01S11/06Systems for determining distance or velocity not using reflection or reradiation using radio waves using intensity measurements
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Automation & Control Theory (AREA)
  • Navigation (AREA)

Abstract

The invention discloses a pedestrian indoor navigation method, a device and a medium based on multi-source information fusion, wherein the method comprises the following steps: map information of the current indoor environment is built, and a Wi-Fi fingerprint library is built; acquiring multi-source information; constructing a position model of the pedestrian to obtain position information of the pedestrian; obtaining the movement distance of the pedestrian through gait detection and step estimation; constructing a fusion model based on extended Kalman filtering, and carrying out fusion processing on multi-source information and pedestrian position information to obtain fused position information; and further correcting and updating the position information of the pedestrians according to the map information and the particle filtering algorithm by utilizing the fused position information. According to the invention, a multisource information fusion technology is introduced, an extended Kalman filtering algorithm and a particle filtering algorithm are introduced, wi-Fi fingerprint matching and map information are added, a proper model is created, errors in an inertial navigation system are corrected, and safer, more reliable and more accurate position information is provided for the system.

Description

一种基于多源信息融合的行人室内导航方法、装置及介质A pedestrian indoor navigation method, device and medium based on multi-source information fusion

技术领域technical field

本发明涉及导航技术领域,尤其涉及一种基于多源信息融合的行人室内导航方法、装置及介质。The present invention relates to the technical field of navigation, in particular to a pedestrian indoor navigation method, device and medium based on multi-source information fusion.

背景技术Background technique

随着社会的发展与进步,人们的生活水平日益提高,对于可靠、实时的基于位置服务需求也越来越高。然而,在室内条件下,由于室内环境复杂多变,存在障碍物多、干扰源多等问题,使得卫星导航信号定位性能变差,在室内定位中效果较差。在卫星导航失效时,常见的导航方法有:With the development and progress of society, people's living standards are improving day by day, and the demand for reliable and real-time location-based services is also getting higher and higher. However, under indoor conditions, due to the complex and changeable indoor environment, there are many obstacles and interference sources, etc., which make the positioning performance of satellite navigation signals deteriorate, and the effect in indoor positioning is poor. When satellite navigation fails, common navigation methods are:

(1)惯性导航,基于惯性器件的室内定位技术是以惯性传感器为基础进行的惯性导航,主要包括行人惯性导航系统PINS和行人航位推算PDR两种方式。前者往往要求传感器的精度较高,并且需要通过抑制误差漂移以保证长时间的定位可靠性;后者往往对传感器的精度要求较低,依赖于人行走的特性,通过跨步检测,在行人航向上以步长累积的方式更新行人位置。相比于PINS,PDR往往能够在成本有限的条件下,实现更好的定位效果,获得更高的成本效益。然而,由于低成本的MEMS惯性传感器往往误差较大,并且误差会随着时间而逐渐累积,这对定位结果造成了极为不利的影响。(1) Inertial navigation. Indoor positioning technology based on inertial devices is based on inertial sensors. It mainly includes pedestrian inertial navigation system PINS and pedestrian dead reckoning PDR. The former often requires high sensor accuracy, and needs to ensure long-term positioning reliability by suppressing error drift; the latter often requires lower sensor accuracy and depends on the characteristics of people walking. Updating the pedestrian position in the way of accumulating steps. Compared with PINS, PDR can often achieve better positioning effect and obtain higher cost-effectiveness under the condition of limited cost. However, because low-cost MEMS inertial sensors tend to have large errors, and the errors will gradually accumulate over time, which has an extremely adverse impact on the positioning results.

(2)Wi-Fi定位,近年来由于Wi-Fi技术的普及,大部分室内已经有Wi-Fi信号覆盖,基于Wi-Fi信号实现的定位技术,虽然有大量已有的基础设施提供支持,并且通过智能手机可以直接获取Wi-Fi的信号强度,在很大程度上推动了Wi-Fi定位技术的发展,但由于Wi-Fi信号易受环境干扰,定位精度通常难以保证。(2) Wi-Fi positioning. Due to the popularization of Wi-Fi technology in recent years, most indoors have already been covered by Wi-Fi signals. The positioning technology based on Wi-Fi signals is supported by a large number of existing infrastructures. Moreover, the signal strength of Wi-Fi can be obtained directly through smartphones, which greatly promotes the development of Wi-Fi positioning technology. However, because Wi-Fi signals are susceptible to environmental interference, positioning accuracy is usually difficult to guarantee.

(3)多源信息融合的室内定位技术,多源信息融合技术利用多个传感器获取的信息,按照特定的规则对这些信息进行分类、处理、融合和使用,实现各传感器之间的信息优势互补。与惯性导航优势互补,有效抑制惯性导航累积误差,进一步提高了定位精度。(3) Multi-source information fusion indoor positioning technology, multi-source information fusion technology uses information obtained by multiple sensors, and classifies, processes, fuses and uses the information according to specific rules, so as to realize the complementary advantages of information between sensors . Complementing the advantages of inertial navigation, it can effectively suppress the cumulative error of inertial navigation and further improve the positioning accuracy.

随着人们对室内高精度定位需求的日益增加,定位精度要求也随之提高,但是尚缺少一种室内高精度定位的技术方案。With the increasing demand for indoor high-precision positioning, the requirements for positioning accuracy also increase, but there is still a lack of a technical solution for indoor high-precision positioning.

发明内容Contents of the invention

为至少一定程度上解决现有技术中存在的技术问题之一,本发明的目的在于提供一种基于多源信息融合的行人室内导航方法、装置及介质。In order to solve one of the technical problems in the prior art at least to a certain extent, the object of the present invention is to provide a pedestrian indoor navigation method, device and medium based on multi-source information fusion.

本发明所采用的技术方案是:The technical scheme adopted in the present invention is:

一种基于多源信息融合的行人室内导航方法,包括以下步骤:A pedestrian indoor navigation method based on multi-source information fusion, comprising the following steps:

构建当前室内环境的地图信息,在地图上选取参考点,建立Wi-Fi指纹库;Construct the map information of the current indoor environment, select reference points on the map, and establish a Wi-Fi fingerprint database;

利用智能终端上的传感器获取多源信息,所述多源信息包括Wi-Fi信号强度,角速度,速度,加速度以及姿态信息;Use the sensor on the smart terminal to obtain multi-source information, the multi-source information includes Wi-Fi signal strength, angular velocity, velocity, acceleration and attitude information;

根据所述Wi-Fi信号强度和Wi-Fi指纹库,构建行人的位置模型,得到行人的位置信息;According to described Wi-Fi signal strength and Wi-Fi fingerprint storehouse, build the position model of pedestrian, obtain the position information of pedestrian;

通过所述姿态信息判断智能终端的瞬时姿态,将加速度从载体坐标系转化到导航坐标系,通过步态检测与步长估计获得行人的运动距离;Judging the instantaneous attitude of the intelligent terminal through the attitude information, transforming the acceleration from the carrier coordinate system to the navigation coordinate system, and obtaining the movement distance of the pedestrian through gait detection and step length estimation;

根据多源信息和运动距离,构建行人的运动模型,得到行人的运动信息,所述运动信息包括:位置、速度;根据扩展卡曼滤波算法,构建基于扩展卡尔曼滤波的融合模型,将多源信息以及行人的位置信息进行融合处理,获取融合后的位置信息,完成拓展卡尔曼滤波中状态向量更新及对误差值的最优估计;According to the multi-source information and the movement distance, the pedestrian movement model is constructed to obtain the pedestrian movement information, the movement information includes: position, speed; according to the extended Kalman filter algorithm, the fusion model based on the extended Kalman filter is constructed, and the multi-source information and the location information of pedestrians are fused, and the fused location information is obtained, and the state vector update and the optimal estimation of the error value in the extended Kalman filter are completed;

利用融合后的位置信息,根据地图信息和粒子滤波算法对行人的位置信息进行进一步的修正与更新,完成多源信息融合的行人室内导航。Using the fused position information, the pedestrian's position information is further corrected and updated according to the map information and the particle filter algorithm, and the pedestrian indoor navigation with multi-source information fusion is completed.

进一步地,所述行人的位置模型的表达式为:Further, the expression of the position model of the pedestrian is:

其中,表示估计的位置,x(tk)、y(tk)表示前K个位置点对应的位置坐标,ωj(tk)表示相应数据点的权重系数,基于信号强度的权重系数可表示/>RSSI为信号强度。in, Represents the estimated position, x(t k ), y(t k ) represent the position coordinates corresponding to the first K position points, ω j (t k ) represents the weight coefficient of the corresponding data point, and the weight coefficient based on signal strength can be expressed/ > RSSI is signal strength.

进一步地,采用以下方式将加速度从载体坐标系转化到导航坐标系:Further, the acceleration is converted from the carrier coordinate system to the navigation coordinate system in the following way:

其中,为导航系到载体系的方向余弦矩阵,Cr为绕Y轴的旋转矩阵,Cp为绕Z轴的旋转矩阵,Cy为绕X轴的旋转矩阵,γ为绕Y轴的旋转角度,p为绕X轴的旋转角度,y为绕Z轴的旋转角度。in, is the direction cosine matrix from the navigation system to the carrier system, C r is the rotation matrix around the Y axis, C p is the rotation matrix around the Z axis, Cy is the rotation matrix around the X axis, γ is the rotation angle around the Y axis, p is the rotation angle around the X axis, and y is the rotation angle around the Z axis.

进一步地,所述步态检测的步骤包括:Further, the step of described gait detection comprises:

计算整体加速度:Compute the overall acceleration:

式中,ax,ay,az为采集的三轴加速度值,通过计算整体加速度a来降低传感器姿态的影响;In the formula, a x , a y , a z are the collected three-axis acceleration values, and the influence of the sensor attitude is reduced by calculating the overall acceleration a;

根据整体加速度,在滑动窗口中获得潜在峰值,利用预设的加速度阈值进行初次判断;According to the overall acceleration, the potential peak value is obtained in the sliding window, and the initial judgment is made using the preset acceleration threshold;

计算潜在波峰与前一波峰时间差,利用预设的行走一步时间阈值范围进行二次判断;Calculate the time difference between the potential peak and the previous peak, and use the preset walking step time threshold range to make a second judgment;

将潜在波峰处与前后邻域加速度比较,进行三次判断去除伪波峰,若潜在峰值点为最大值则算法记一步,否则不做计步处理。Compare the potential peak with the acceleration of the front and rear neighbors, and make three judgments to remove the false peak. If the potential peak point is the maximum value, the algorithm will count one step, otherwise, no step counting process will be performed.

进一步地,所述步长估计为:Further, the step size is estimated as:

其中SL为步长,amax和amin分别为行人纵向加速度的最大值和最小值。Where SL is the step length, a max and a min are the maximum and minimum values of the pedestrian's longitudinal acceleration, respectively.

进一步地,所述扩展卡曼滤波算法包括:Further, the extended Kalman filter algorithm includes:

Xk+1=f(Xk)+wk X k+1 =f(X k )+w k

其中,Xk+1为第k+1步的状态变量,Xk为第k步的状态变量,wk为过程噪声;f()为扩展卡尔曼滤波中的非线性状态函数;Among them, X k+1 is the state variable of the k+1th step, X k is the state variable of the kth step, w k is the process noise; f() is the nonlinear state function in the extended Kalman filter;

实现系统矩阵线性化的方法包括:Methods to achieve linearization of the system matrix include:

其中,Zk为观测变量,vk为观测噪声,h()为扩展卡尔曼滤波中的量测函数;Among them, Z k is the observation variable, v k is the observation noise, h() is the measurement function in the extended Kalman filter;

实现观测矩阵线性化的方法包括:Methods to achieve linearization of the observation matrix include:

其中,Hk为h(X)对X偏导的雅可比矩阵,为第k步的状态估计。Among them, H k is the Jacobian matrix of the partial derivative of h(X) to X, is the state estimate of the kth step.

进一步地,拓展卡尔曼滤波中状态向量更新的步骤包括:Further, the steps of updating the state vector in the expanded Kalman filter include:

一步状态预测更新:One-step state prediction update:

式中,为上一状态的最优值,/>为当前状态的一步预测值;In the formula, is the optimal value of the previous state, /> is the one-step forecast value of the current state;

一步预测估计误差协方差矩阵更新:One-step forecast estimation error covariance matrix update:

式中,Pk+1|k对应的协方差一步预测值,φk+1|k为状态转移矩阵,Pk|k为/>对应的协方差,/>为φk+1|k的转置,Qk为过程噪声协方差矩阵;In the formula, P k+1|k is Corresponding covariance one-step forecast value, φ k+1|k is the state transition matrix, P k|k is /> Corresponding covariance, /> is the transpose of φ k+1|k , Q k is the process noise covariance matrix;

计算扩展卡尔曼滤波增益:Compute the extended Kalman filter gain:

式中,Kk+1为扩展卡尔曼滤波增益矩阵,Hk+1为观测矩阵,为观测矩阵的转置,Rk+1为测量噪声协方差矩阵;In the formula, K k+1 is the extended Kalman filter gain matrix, H k+1 is the observation matrix, is the transposition of the observation matrix, and R k+1 is the measurement noise covariance matrix;

由观测向量计算新息,并更新状态估计:Compute innovations from observation vectors and update state estimates:

式中,为当前状态的最优估计值,Zk+1为当前观测值,/>为当前观测预测值;更新估计误差协方差:In the formula, is the optimal estimated value of the current state, Z k+1 is the current observed value, /> is the predicted value for the current observation; update the estimated error covariance:

Pk+1=[I-Kk+1Hk+1]Pk+1|k P k+1 =[IK k+1 H k+1 ]P k+1|k

式中,Pk+1对应的协方差,I为单位矩阵。In the formula, P k+1 is The corresponding covariance, I is the identity matrix.

进一步地,所述根据地图信息和粒子滤波算法对行人的位置信息进行进一步的修正与更新,包括:Further, the further correction and updating of the pedestrian's position information according to the map information and the particle filter algorithm includes:

粒子滤波的状态转移方程和观测方程如下:The state transition equation and observation equation of the particle filter are as follows:

x(t)=f(x(t-1),u(t),w(t))x(t)=f(x(t-1),u(t),w(t))

y(t)=f(x(t),e(t))y(t)=f(x(t),e(t))

式中,x(t)为t时刻状态,u(t)为控制量,w(t)和e(t)分别为状态噪音和观测噪音,粒子滤波从观测y(t)和上个时刻状态x(t-1),u(t),w(t)中过滤出t时刻的状态x(t);In the formula, x(t) is the state at time t, u(t) is the control quantity, w(t) and e(t) are the state noise and observation noise respectively, and the particle filter is obtained from the observation y(t) and the state at the last time x(t-1), u(t), w(t) filter out the state x(t) at time t;

粒子滤波中的更新方法主要包含以下步骤:The update method in particle filter mainly includes the following steps:

A1、初始化:将观测值初值xi(k)作为概率密度函数x(0)的初值,由进行N次抽样,初始wi(k)设为/> A1. Initialization: take the initial value x i (k) of the observed value as the initial value of the probability density function x(0), by right Sampling is carried out N times, and the initial w i (k) is set to />

A2、一步预测:对每个粒子xi(k),通过转换公式p(x(k+1)|xi(k))获得一个新的粒子;A2. One-step prediction: For each particle x i (k), obtain a new particle by converting the formula p(x(k+1)| xi (k));

A3、重要性采样:对于任一粒子xi(k+1)求解它们的权值wi(k+1)=p(z(k+1)|xi(k+1));A3. Importance sampling: For any particle x i (k+1), solve their weight w i (k+1)=p(z(k+1)| xi (k+1));

A4、归一化:对权值进行归一化处理:A4. Normalization: Normalize the weights:

A5、重采样:根据权重大小对粒子进行重采样,权重大的粒子重复的概率较大,权重较小的粒子重复概率小;A5. Resampling: Resampling of particles according to the weight size, particles with large weights have a higher probability of repetition, and particles with smaller weights have a lower probability of repetition;

A6、粒子更新同时加入地图信息:A6. The particle update also adds map information:

由于室内行走区域存在限值,根据地图信息进行地图匹配,当一步预测的粒子越过可通行区域时(如越墙),则将它的权重值赋为0:Due to the limitation of the indoor walking area, map matching is performed according to the map information. When the one-step predicted particle crosses the passable area (such as crossing the wall), its weight value is assigned to 0:

A7、采用基本的粒子滤波算法进行计算,根据解算结果完成行人室内定位及导航。A7. The basic particle filter algorithm is used for calculation, and the indoor positioning and navigation of pedestrians are completed according to the calculation results.

本发明所采用的另一技术方案是:Another technical scheme adopted in the present invention is:

一种基于多源信息融合的行人室内导航装置,包括:A pedestrian indoor navigation device based on multi-source information fusion, comprising:

至少一个处理器;at least one processor;

至少一个存储器,用于存储至少一个程序;at least one memory for storing at least one program;

当所述至少一个程序被所述至少一个处理器执行,使得所述至少一个处理器实现上所述方法。When the at least one program is executed by the at least one processor, the at least one processor implements the above method.

本发明所采用的另一技术方案是:Another technical scheme adopted in the present invention is:

一种计算机可读存储介质,其中存储有处理器可执行的程序,所述处理器可执行的程序在由处理器执行时用于执行如上所述方法。A computer-readable storage medium stores a processor-executable program therein, and the processor-executable program is used to perform the above method when executed by a processor.

本发明的有益效果是:本发明引入多源信息融合技术,引入拓展卡尔曼滤波算法和粒子滤波算法,同时加入Wi-Fi指纹匹配与地图信息,创建合适的模型,修正惯性导航系统中的误差,为系统提供更加安全可靠和精准的位置信息。The beneficial effects of the present invention are: the present invention introduces multi-source information fusion technology, introduces extended Kalman filter algorithm and particle filter algorithm, adds Wi-Fi fingerprint matching and map information at the same time, creates a suitable model, and corrects errors in the inertial navigation system , to provide more secure, reliable and accurate location information for the system.

附图说明Description of drawings

为了更清楚地说明本发明实施例或者现有技术中的技术方案,下面对本发明实施例或者现有技术中的相关技术方案附图作以下介绍,应当理解的是,下面介绍中的附图仅仅为了方便清晰表述本发明的技术方案中的部分实施例,对于本领域的技术人员而言,在无需付出创造性劳动的前提下,还可以根据这些附图获取到其他附图。In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following describes the accompanying drawings of the embodiments of the present invention or the related technical solutions in the prior art. It should be understood that the accompanying drawings in the following introduction are only In order to clearly describe some embodiments of the technical solutions of the present invention, those skilled in the art can also obtain other drawings based on these drawings without creative work.

图1是本发明实施例中一种基于多源信息融合的行人室内导航方法的总体流程示意图;FIG. 1 is a schematic diagram of an overall flow of a pedestrian indoor navigation method based on multi-source information fusion in an embodiment of the present invention;

图2是本发明实施例中扩展卡尔曼滤波融合模型构建流程示意图;Fig. 2 is a schematic flow chart of building an extended Kalman filter fusion model in an embodiment of the present invention;

图3是本发明实施例中粒子滤波模型构建流程示意图。Fig. 3 is a schematic diagram of the construction process of the particle filter model in the embodiment of the present invention.

具体实施方式Detailed ways

下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。对于以下实施例中的步骤编号,其仅为了便于阐述说明而设置,对步骤之间的顺序不做任何限定,实施例中的各步骤的执行顺序均可根据本领域技术人员的理解来进行适应性调整。Embodiments of the present invention are described in detail below, examples of which are shown in the drawings, wherein the same or similar reference numerals designate the same or similar elements or elements having the same or similar functions throughout. The embodiments described below by referring to the figures are exemplary only for explaining the present invention and should not be construed as limiting the present invention. For the step numbers in the following embodiments, it is only set for the convenience of illustration and description, and the order between the steps is not limited in any way. The execution order of each step in the embodiments can be adapted according to the understanding of those skilled in the art sexual adjustment.

在本发明的描述中,需要理解的是,涉及到方位描述,例如上、下、前、后、左、右等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。In the description of the present invention, it should be understood that the orientation descriptions, such as up, down, front, back, left, right, etc. indicated orientations or positional relationships are based on the orientations or positional relationships shown in the drawings, and are only In order to facilitate the description of the present invention and simplify the description, it does not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus should not be construed as limiting the present invention.

在本发明的描述中,若干的含义是一个或者多个,多个的含义是两个以上,大于、小于、超过等理解为不包括本数,以上、以下、以内等理解为包括本数。如果有描述到第一、第二只是用于区分技术特征为目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量或者隐含指明所指示的技术特征的先后关系。In the description of the present invention, several means one or more, and multiple means more than two. Greater than, less than, exceeding, etc. are understood as not including the original number, and above, below, within, etc. are understood as including the original number. If the description of the first and second is only for the purpose of distinguishing the technical features, it cannot be understood as indicating or implying the relative importance or implicitly indicating the number of the indicated technical features or implicitly indicating the order of the indicated technical features relation.

本发明的描述中,除非另有明确的限定,设置、安装、连接等词语应做广义理解,所属技术领域技术人员可以结合技术方案的具体内容合理确定上述词语在本发明中的具体含义。In the description of the present invention, unless otherwise clearly defined, words such as setting, installation, and connection should be understood in a broad sense, and those skilled in the art can reasonably determine the specific meanings of the above words in the present invention in combination with the specific content of the technical solution.

本发明引入多源信息融合技术,引入扩展卡曼滤波算法和粒子滤波算法,同时加入Wi-Fi指纹匹配与地图信息,创建合适的模型,修正惯性导航系统中的误差,为系统提供更加安全可靠和精准的位置信息。因此,建立基于多源信息融合的行人室内导航方法,为实现行人室内的精确定位提供有效方法。The present invention introduces multi-source information fusion technology, introduces extended Kalman filter algorithm and particle filter algorithm, adds Wi-Fi fingerprint matching and map information at the same time, creates a suitable model, corrects errors in the inertial navigation system, and provides a safer and more reliable solution for the system. and precise location information. Therefore, an indoor pedestrian navigation method based on multi-source information fusion is established to provide an effective method for precise indoor positioning of pedestrians.

参见图1,本实施例提供一种基于多源信息融合的行人室内导航方法,包括以下步骤:Referring to Figure 1, this embodiment provides a pedestrian indoor navigation method based on multi-source information fusion, including the following steps:

步骤1:构建当前室内环境的地图信息,在地图上选取参考点,建立Wi-Fi指纹库;利用手机内置的传感器获取多源信息,多源信息包括:Wi-Fi信号强度,角速度,速度,磁场强度以及姿态信息。Step 1: Construct the map information of the current indoor environment, select reference points on the map, and establish a Wi-Fi fingerprint library; use the built-in sensor of the mobile phone to obtain multi-source information, and the multi-source information includes: Wi-Fi signal strength, angular velocity, speed, Magnetic field strength and attitude information.

步骤2:处理步骤1采集到的Wi-Fi信号,构建行人的位置模型;得到行人的位置信息。Step 2: Process the Wi-Fi signal collected in step 1 to construct a pedestrian location model; obtain pedestrian location information.

在本实施例中,构建行人位置模型时,对其中的行人位置模型进行改进,用受加权k邻近法模型代替k邻近法模型。In this embodiment, when the pedestrian position model is constructed, the pedestrian position model is improved, and the k-proximity model is replaced by a weighted k-proximity model.

具体地,构建行人位置模型包括:Specifically, building a pedestrian position model includes:

其中表示估计的位置,x(tk)、y(tk)表示前K个位置点对应的位置坐标,ωj(tk)表示相应数据点的权重系数,基于信号强度的权重系数可表示/>RSSI为信号强度。in Represents the estimated position, x(t k ), y(t k ) represent the position coordinates corresponding to the first K position points, ω j (t k ) represents the weight coefficient of the corresponding data point, and the weight coefficient based on signal strength can be expressed/ > RSSI is signal strength.

步骤3:通过步骤1中的方向传感器判断手机瞬时姿态,将加速度从载体坐标系通过一定方法转化为导航坐标系。再通过步态检测与步长估计获得行人的运动距离。Step 3: Determine the instantaneous attitude of the mobile phone through the direction sensor in step 1, and convert the acceleration from the carrier coordinate system to the navigation coordinate system through a certain method. Then the pedestrian's movement distance is obtained through gait detection and step length estimation.

其中,加速度坐标系转化的方法为:Among them, the method of transforming the acceleration coordinate system is:

式中,定义导航系n绕Z轴(正轴向上)旋转y角(航向角、yaw),再绕X轴旋转p角(俯仰角、pitch),最后绕Y轴旋转γ角(横滚角、roll)得到导航系到载体系的方向余弦矩阵而载体系到导航系的方向余弦阵/> In the formula, the navigation system n is defined to rotate around the Z axis (positive axis upward) by an angle of y (yaw angle, yaw), then rotate around the X axis by an angle of p (pitch angle, pitch), and finally rotate around the Y axis by an angle of γ (roll angle, roll) to get the cosine matrix of the direction from the navigation system to the vehicle system And the direction cosine matrix from the carrier system to the navigation system />

步态检测包括:Gait detection includes:

1)计算整体加速度,通过计算整体加速度a来降低传感器姿态的影响。式中,ax,ay,az为采集的三轴加速度值。1) Calculate the overall acceleration, and reduce the influence of the sensor attitude by calculating the overall acceleration a. In the formula, a x , a y , a z are the collected three-axis acceleration values.

2)在滑动窗口中获得潜在峰值,利用加速度阈值[1.2g,3g]进行初次判断。2) Obtain the potential peak in the sliding window, and use the acceleration threshold [1.2g, 3g] for initial judgment.

3)计算潜在波峰与前一波峰时间差,利用行走一步时间阈值范围[0.4s,1s]进行二次判断。3) Calculate the time difference between the potential peak and the previous peak, and use the time threshold range of one step [0.4s, 1s] for secondary judgment.

4)将潜在波峰处与前后邻域加速度比较,进行三次判断去除伪波峰,若潜在峰值点为最大值则算法记一步,否则不做计步处理。4) Compare the potential peak with the acceleration of the front and back neighbors, and perform three judgments to remove the false peak. If the potential peak point is the maximum value, the algorithm will count one step, otherwise, no step counting process will be performed.

步长估计为:The step size is estimated as:

其中SL为步长,amax和amin分别为行人纵向加速度的最大值和最小值。Where SL is the step length, a max and a min are the maximum and minimum values of the pedestrian's longitudinal acceleration, respectively.

步骤4:处理步骤1中采集到的多源信息,构建行人的运动模型;得到行人的运动信息,包括:位置、速度;根据拓展卡尔曼滤波算法,构建基于拓展卡尔曼滤波的融合模型;将步骤1中得到的多源信息以及步骤2中得到的行人的位置信息进行融合处理,获取融合后的位置信息;完成拓展卡尔曼滤波中状态向量更新及对误差值的最优估计。Step 4: Process the multi-source information collected in step 1 to construct a pedestrian motion model; obtain pedestrian motion information, including: position, speed; according to the extended Kalman filter algorithm, construct a fusion model based on the extended Kalman filter; The multi-source information obtained in step 1 and the position information of pedestrians obtained in step 2 are fused to obtain the fused position information; the update of the state vector and the optimal estimation of the error value in the extended Kalman filter are completed.

扩展卡尔曼滤波算法包括:Extended Kalman filter algorithms include:

状态方程为:The state equation is:

X(k+1)=φX(k)+ΓW(k)X(k+1)=φX(k)+ΓW(k)

其中X(k)=[N(k)E(k)vn(k)ve(k)]T,N(k)和E(k)为分别为k时刻的北方向和东方向的位置状态值,vn(k)和ve(k)分别为k时刻的北方向和东方向的速度状态值,W(k)为k时刻的系统噪声,φ和Γ分别为k时刻到k+1时刻的状态转换矩阵和系统噪声系数矩阵。Where X(k)=[N(k)E(k)v n (k)v e (k)] T , N(k) and E(k) are respectively the north and east positions at time k State value, v n (k) and v e (k) are the speed state values in the north direction and east direction at time k respectively, W(k) is the system noise at time k, φ and Γ are respectively State transition matrix and system noise coefficient matrix at time 1.

观测方程为:The observation equation is:

Z(k+1)=H(k+1)X(k+1)+VZ(k+1)=H(k+1)X(k+1)+V

其中Z(k)=[Nwifi Ewifi s]T,Nwifi和Ewifi分别为通过Wi-Fi获取的北方向坐标值和东方向坐标值,s是通过步态检测得到的两秒内平面位移数值。H(k+1)为量测矩阵,可由公式推导并线性化后得到。V为观测噪声,并且假设噪声服从高斯分布。Where Z(k)=[N wifi E wifi s] T , N wifi and E wifi are the north and east coordinates obtained through Wi-Fi, respectively, and s is the plane within two seconds obtained through gait detection displacement value. H(k+1) is the measurement matrix, which can be derived from the formula and obtained after linearization. V is the observation noise, and it is assumed that the noise obeys a Gaussian distribution.

扩展卡尔曼滤波中状态向量更新的方法,包括以下步骤:The method for updating the state vector in the extended Kalman filter includes the following steps:

步骤4-1,一步状态预测更新:Step 4-1, one-step state prediction update:

其中,为上一状态的最优值,/>为当前状态的一步预测值;in, is the optimal value of the previous state, /> is the one-step forecast value of the current state;

步骤4-2,一步预测估计误差协方差矩阵更新:Step 4-2, one-step forecast estimation error covariance matrix update:

其中,Pk+1|k对应的协方差一步预测值,φk+1|k为状态转移矩阵,Pk|k为/>对应的协方差,/>为φk+1|k的转置,Qk为过程噪声协方差矩阵;Among them, P k+1|k is Corresponding covariance one-step forecast value, φ k+1|k is the state transition matrix, P k|k is /> Corresponding covariance, /> is the transpose of φ k+1|k , Q k is the process noise covariance matrix;

步骤4-3,计算扩展卡尔曼滤波增益:Step 4-3, calculate the extended Kalman filter gain:

其中,Kk+1为扩展卡尔曼滤波增益矩阵,Hk+1为观测矩阵,为观测矩阵的转置,Rk+1为测量噪声协方差矩阵;Among them, K k+1 is the extended Kalman filter gain matrix, H k+1 is the observation matrix, is the transposition of the observation matrix, and R k+1 is the measurement noise covariance matrix;

步骤4-4,由观测向量计算新息,并更新状态估计:Step 4-4, calculate the innovation from the observation vector, and update the state estimate:

其中,为当前状态的最优估计值,Zk+1为当前观测值,/>为当前观测预测值。步骤4-5,更新估计误差协方差:in, is the optimal estimated value of the current state, Z k+1 is the current observed value, /> is the predicted value for the current observation. Steps 4-5, update the estimated error covariance:

Pk+1=[I-Kk+1Hk+1]Pk+1|k P k+1 =[IK k+1 H k+1 ]P k+1|k

其中,Pk+1对应的协方差,I为单位矩阵。Among them, P k+1 is The corresponding covariance, I is the identity matrix.

步骤5:利用步骤4获得的基于拓展卡尔曼滤波的位置信息,根据粒子滤波算法,同时加入地图匹配,对行人的位置信息进行进一步的修正与更新,完成多源信息融合的行人室内导航。Step 5: Use the location information based on extended Kalman filter obtained in step 4, according to the particle filter algorithm, and add map matching at the same time, further correct and update the location information of pedestrians, and complete the indoor navigation of pedestrians with multi-source information fusion.

粒子滤波算法包括:Particle filter algorithms include:

粒子滤波器是用随机样本来直接估计近似状态后验概率密度函数:The particle filter uses random samples to directly estimate the approximate state posterior probability density function:

其中,wi(k)为粒子xi k的权值,所有权值之和为1,δ(.)为狄拉克函数。Among them, w i (k) is the weight value of particle x i k , the sum of all values is 1, and δ(.) is the Dirac function.

粒子滤波中的更新方法包含以下步骤:The update method in particle filtering consists of the following steps:

步骤5-1,初始化。Step 5-1, initialization.

将观测值初值xi(k)作为概率密度函数x(0)的初值,由对/>进行N次抽样,初始wi(k)设为/> Taking the initial value x i (k) of the observed value as the initial value of the probability density function x(0), by right /> Sampling is carried out N times, and the initial w i (k) is set to />

步骤5-2,一步预测。Step 5-2, one-step prediction.

对每个粒子xi(k),通过转换公式p(x(k+1)|xi(k))获得一个新的粒子xi(k+1)。For each particle xi (k), a new particle xi (k+1) is obtained by converting the formula p(x(k+1)| xi (k)).

步骤5-3,重要性采样。Step 5-3, importance sampling.

对于任一粒子xi(k+1)求解它们的权值wi(k+1)=p(z(k+1)|xi(k+1))。For any particle x i (k+1), solve their weight value w i (k+1)=p(z(k+1)| xi (k+1)).

步骤5-4,归一化。Step 5-4, normalization.

对权值进行归一化处理。 Normalize the weights.

步骤5-5,重采样。Step 5-5, resampling.

根据权重大小对粒子进行重采样,权重大的粒子重复的概率较大,权重较小的粒子重复概率小。Particles are resampled according to the weight size, the probability of repetition of particles with large weight is higher, and the probability of repetition of particles with smaller weight is small.

步骤5-6,粒子更新。Step 5-6, particle update.

及/> make and />

步骤5-7,采用基本的粒子滤波算法进行计算。In steps 5-7, the basic particle filter algorithm is used for calculation.

使用上述解算结果完成行人室内导航。Use the above calculation results to complete pedestrian indoor navigation.

以下结合附图及具体实施例对上述方法进行详细解释说明。The above method will be explained in detail below in conjunction with the accompanying drawings and specific embodiments.

如图1所示,本实施例提供一种基于多源信息融合的行人室内导航方法,具体方法技术点如下:As shown in Figure 1, this embodiment provides a pedestrian indoor navigation method based on multi-source information fusion. The specific technical points of the method are as follows:

S1、提供Wi-Fi定位位置模型构建功能。S1. Provide a Wi-Fi positioning location model building function.

能够通过手机内部内置传感器获取导航系统内所有的传感器信息,通过选择位置、信号强度信息,构建行人位置模型。It can obtain all the sensor information in the navigation system through the built-in sensor of the mobile phone, and build a pedestrian position model by selecting the position and signal strength information.

多源信息融合技术的关键在于确定多源信息之间的关系。在导航系统中,多源信息主要包括三轴角速度、加速度、Wi-Fi信号强度、磁场强度等。构建行人位置模型具体包括:The key to multi-source information fusion technology is to determine the relationship between multi-source information. In the navigation system, multi-source information mainly includes three-axis angular velocity, acceleration, Wi-Fi signal strength, magnetic field strength, etc. Constructing the pedestrian position model specifically includes:

其中表示估计的位置,x(tk)、y(tk)表示前K个位置点对应的位置坐标,ωj(tk)表示相应数据点的权重系数,基于信号强度的权重系数可表示/>RSSI为信号强度。in Represents the estimated position, x(t k ), y(t k ) represent the position coordinates corresponding to the first K position points, ω j (t k ) represents the weight coefficient of the corresponding data point, and the weight coefficient based on signal strength can be expressed/ > RSSI is signal strength.

S2、扩展卡尔曼滤波融合模型构建功能。S2. Extend the Kalman filter fusion model building function.

利用建立的模型,通过扩展卡尔曼滤波算法融合位置、速度等信息,计算出实时位置结果。Using the established model, the extended Kalman filter algorithm is used to fuse position, velocity and other information to calculate the real-time position result.

传统卡尔曼滤波仅适用于线性系统,而扩展卡尔曼滤波器是卡尔曼滤波算法的一种,是将传统卡尔曼滤波应用到非线性领域中,传统卡尔曼滤波中的系统矩阵F和观测矩阵H将被f(x)和h(x)代替,随后将非线性系统按照一阶泰勒展开进行线性化。The traditional Kalman filter is only suitable for linear systems, and the extended Kalman filter is a kind of Kalman filter algorithm, which applies the traditional Kalman filter to the nonlinear field. The system matrix F and the observation matrix in the traditional Kalman filter H will be replaced by f(x) and h(x), and the nonlinear system is then linearized according to a first-order Taylor expansion.

扩展卡尔曼滤波的系统模型为:The system model of the extended Kalman filter is:

Xk+1=f(Xk)+wk X k+1 =f(X k )+w k

其中,Xk+1为第k+1步的状态变量,Xk为第k步的状态变量,wk为过程噪声;f()为扩展卡尔曼滤波中的非线性状态函数。Among them, X k+1 is the state variable of the k+1th step, X k is the state variable of the kth step, w k is the process noise; f() is the nonlinear state function in the extended Kalman filter.

实现系统矩阵线性化的方法包括:Methods to achieve linearization of the system matrix include:

其中,Fk-1为f对X偏导的雅克比矩阵,为第k-1步的后验状态估计,X为状态变量。Among them, F k-1 is the Jacobian matrix of the partial derivative of f to X, is the posterior state estimation of step k-1, and X is the state variable.

扩展卡尔曼滤波的观测模型为:The observation model of the extended Kalman filter is:

Zk=h(Xk)+vk Z k =h(X k )+v k

其中,Zk为观测变量,vk为观测噪声,h()为扩展卡尔曼滤波中的量测函数;Among them, Z k is the observation variable, v k is the observation noise, h() is the measurement function in the extended Kalman filter;

实现观测矩阵线性化的方法包括:Methods to achieve linearization of the observation matrix include:

其中,Hk为h(X)对X偏导的雅可比矩阵,为第k步的状态估计。Among them, H k is the Jacobian matrix of the partial derivative of h(X) to X, is the state estimate for the kth step.

如图2所示,图2给出了本发明的扩展卡尔曼滤波融合模型构建过程。在本实施例中,利用角速度、加速度、位置等多源信息来构建误差模型,具体的扩展卡尔曼滤波更新主要包括以下五个步骤:As shown in Figure 2, Figure 2 shows the construction process of the extended Kalman filter fusion model of the present invention. In this embodiment, the error model is constructed using multi-source information such as angular velocity, acceleration, and position. The specific extended Kalman filter update mainly includes the following five steps:

步骤1,一步状态预测更新:Step 1, one-step state prediction update:

其中,为上一状态的最优值,/>为当前状态的一步预测值;in, is the optimal value of the previous state, /> is the one-step forecast value of the current state;

步骤2,一步预测估计误差协方差矩阵更新:Step 2, one-step forecast estimation error covariance matrix update:

其中,Pk+1|k对应的协方差一步预测值,φk+1|k为状态转移矩阵,Pk|k为/>对应的协方差,/>为φk+1|k的转置,Qk为过程噪声协方差矩阵;Among them, P k+1|k is Corresponding covariance one-step forecast value, φ k+1|k is the state transition matrix, P k|k is /> Corresponding covariance, /> is the transpose of φ k+1|k , Q k is the process noise covariance matrix;

步骤3,计算扩展卡尔曼滤波增益:Step 3, calculate the extended Kalman filter gain:

其中,Kk+1为扩展卡尔曼滤波增益矩阵,Hk+1为观测矩阵,为观测矩阵的转置,Rk+1为测量噪声协方差矩阵;Among them, K k+1 is the extended Kalman filter gain matrix, H k+1 is the observation matrix, is the transposition of the observation matrix, and R k+1 is the measurement noise covariance matrix;

步骤4,由观测向量计算新息,并更新状态估计:Step 4, calculate the innovation from the observation vector, and update the state estimation:

其中,为当前状态的最优估计值,Zk+1为当前观测值,/>为当前观测预测值。in, is the optimal estimated value of the current state, Z k+1 is the current observed value, /> is the predicted value for the current observation.

步骤5,更新估计误差协方差:Step 5, update the estimated error covariance:

Pk+1=[I-Kk+1Hk+1]Pk+1|k P k+1 =[IK k+1 H k+1 ]P k+1|k

其中,Pk+1对应的协方差,I为单位矩阵。Among them, P k+1 is The corresponding covariance, I is the identity matrix.

基于以上五个步骤,完成了扩展卡尔曼滤波中状态向量的更新及对误差值的最优估计。Based on the above five steps, the update of the state vector and the optimal estimation of the error value in the extended Kalman filter are completed.

S3、基于粒子滤波与地图匹配融合模型构建功能。S3, based on particle filter and map matching fusion model construction function.

基于上述模型的估计结果,在拓展卡尔曼滤波融合模型的基础上进行粒子滤波并且加入地图信息,减少运算量,加快运算速度,降低误差,提高模型精度。Based on the estimation results of the above models, the particle filter is carried out on the basis of the extended Kalman filter fusion model and map information is added to reduce the amount of calculation, speed up the calculation, reduce the error, and improve the accuracy of the model.

与卡尔曼滤波相比较粒子滤波的思想基于蒙特卡洛方法,它是利用粒子集来表示概率,可以用在任何形式的状态空间模型上。其核心思想是通过从后验概率中抽取的随机状态粒子来表达其分布,是一种顺序重要性采样法。简单来说,粒子滤波法是指通过寻找一组在状态空间传播的随机样本对概率密度函数进行近似,以样本均值代替积分运算,从而获得状态最小方差分布的过程。这里的样本即指粒子,当样本数量N→∞时可以逼近任何形式的概率密度分布。Compared with Kalman filtering, the idea of particle filtering is based on the Monte Carlo method, which uses particle sets to represent probability and can be used in any form of state space model. Its core idea is to express its distribution through random state particles drawn from the posterior probability, which is a sequential importance sampling method. In simple terms, the particle filter method refers to the process of approximating the probability density function by finding a group of random samples propagated in the state space, and replacing the integral operation with the sample mean value, so as to obtain the minimum variance distribution of the state. The samples here refer to particles, which can approach any form of probability density distribution when the number of samples is N→∞.

粒子滤波的状态转移方程和观测方程如下:The state transition equation and observation equation of the particle filter are as follows:

x(t)=f(x(t-1),u(t),w(t))x(t)=f(x(t-1),u(t),w(t))

y(t)=f(x(t),e(t))y(t)=f(x(t),e(t))

其中的x(t)为t时刻状态,u(t)为控制量,w(t)和e(t)分别为状态噪音和观测噪音,粒子滤波从观测y(t)和上个时刻状态x(t-1),u(t),w(t)中过滤出t时刻的状态x(t)。Among them, x(t) is the state at time t, u(t) is the control quantity, w(t) and e(t) are the state noise and observation noise respectively, and the particle filter is obtained from the observation y(t) and the state x at the last time (t-1), u(t), w(t) filter out the state x(t) at time t.

参见图3,粒子滤波中的更新方法主要包含以下步骤:Referring to Figure 3, the update method in the particle filter mainly includes the following steps:

步骤1,初始化。Step 1, initialization.

将观测值初值xi(k)作为概率密度函数x(0)的初值,由对/>进行N次抽样,初始wi(k)设为/> Taking the initial value x i (k) of the observed value as the initial value of the probability density function x(0), by right /> Sampling is carried out N times, and the initial w i (k) is set to />

步骤2,一步预测。Step 2, one-step prediction.

对每个粒子xi(k),通过转换公式p(x(k+1)|xi(k))获得一个新的粒子xi(k+1)。For each particle xi (k), a new particle xi (k+1) is obtained by converting the formula p(x(k+1)| xi (k)).

步骤3,重要性采样。Step 3, importance sampling.

对于任一粒子xi(k+1)求解它们的权值wi(k+1)=p(z(k+1)|xi(k+1))。For any particle x i (k+1), solve their weight value w i (k+1)=p(z(k+1)| xi (k+1)).

步骤4,归一化。Step 4, normalization.

对权值进行归一化处理。 Normalize the weights.

步骤5,重采样。Step 5, resampling.

根据权重大小对粒子进行重采样,权重大的粒子重复的概率较大,权重较小的粒子重复概率小。Particles are resampled according to the weight size, the probability of repetition of particles with large weight is higher, and the probability of repetition of particles with smaller weight is small.

步骤6,粒子更新同时加入地图信息。Step 6, update the particles and add map information at the same time.

及/> make and />

由于室内行走区域存在限值,根据地图信息进行地图匹配,当一步预测的粒子越过可通行区域时(如越墙),则将它的权重值赋为0:Due to the limitation of the indoor walking area, map matching is performed according to the map information. When the one-step predicted particle crosses the passable area (such as crossing the wall), its weight value is assigned to 0:

步骤7,采用基本的粒子滤波算法进行计算。In step 7, the basic particle filter algorithm is used for calculation.

S4、基于多源信息融合的行人室内导航方法。S4. A pedestrian indoor navigation method based on multi-source information fusion.

在上述三个模型构建完成后,基于多源信息融合的行人室内导航方法基本生成。在卫星导航失效时引入该导航模型,预测误差,实时修正各种因素对惯性导航系统造成的定位误差,提高室内行人的定位精度。After the above three models are constructed, the pedestrian indoor navigation method based on multi-source information fusion is basically generated. When the satellite navigation fails, the navigation model is introduced to predict the error, correct the positioning error caused by various factors to the inertial navigation system in real time, and improve the positioning accuracy of indoor pedestrians.

在Wi-Fi定位位置模型、扩展卡尔曼滤波融合模型及粒子滤波与地图匹配融合模型构建完成之后,本发明提出的基于多源信息融合的行人室内导航方法基本完成,具体的流程图如图1所示。在行人室内导航系统中,利用多源信息构建好扩展卡尔曼滤波融合模型及粒子滤波与地图匹配融合模型,在常规卫星导航失效时,利用手机中的内置传感器获得导航系统所需要信息,进而通过本发明的导航方法进行行人室内导航获得行人当前位置信息。After the Wi-Fi positioning model, the extended Kalman filter fusion model and the particle filter and map matching fusion model are constructed, the pedestrian indoor navigation method based on multi-source information fusion proposed by the present invention is basically completed. The specific flow chart is shown in Figure 1 shown. In the pedestrian indoor navigation system, the multi-source information is used to construct the extended Kalman filter fusion model and the particle filter and map matching fusion model. The navigation method of the present invention performs indoor navigation of pedestrians to obtain the current position information of pedestrians.

综上所述,本实施例方法相对于现有技术,至少具有如下优点及有益效果:In summary, compared with the prior art, the method of this embodiment has at least the following advantages and beneficial effects:

(1)本实施例方法通过引入拓展卡尔曼滤波与粒子滤波,改进传统的PDR定位,并且加入地图匹配,减少运算量。(1) The method of this embodiment improves the traditional PDR positioning by introducing extended Kalman filter and particle filter, and adds map matching to reduce the amount of computation.

(2)本实施例方法基于多源信息融合的行人室内导航方法能够提高行人导航定位效率和精度。(2) The method of this embodiment is based on the pedestrian indoor navigation method based on multi-source information fusion, which can improve the efficiency and accuracy of pedestrian navigation and positioning.

(3)本实施例方法中的模型和方法能够适用于大部分行人室内导航系统。(3) The model and method in the method of this embodiment can be applied to most indoor navigation systems for pedestrians.

本实施例还提供一种基于多源信息融合的行人室内导航装置,包括:This embodiment also provides a pedestrian indoor navigation device based on multi-source information fusion, including:

至少一个处理器;at least one processor;

至少一个存储器,用于存储至少一个程序;at least one memory for storing at least one program;

当所述至少一个程序被所述至少一个处理器执行,使得所述至少一个处理器实现上所述方法。When the at least one program is executed by the at least one processor, the at least one processor implements the above method.

本实施例的一种基于多源信息融合的行人室内导航装置,可执行本发明方法实施例所提供的一种基于多源信息融合的行人室内导航方法,可执行方法实施例的任意组合实施步骤,具备该方法相应的功能和有益效果。A pedestrian indoor navigation device based on multi-source information fusion in this embodiment can execute a pedestrian indoor navigation method based on multi-source information fusion provided by the method embodiment of the present invention, and can execute any combination of implementation steps of the method embodiment , have the corresponding functions and beneficial effects of the method.

本申请实施例还公开了一种计算机程序产品或计算机程序,该计算机程序产品或计算机程序包括计算机指令,该计算机指令存储在计算机可读存介质中。计算机设备的处理器可以从计算机可读存储介质读取该计算机指令,处理器执行该计算机指令,使得该计算机设备执行图1所示的方法。The embodiment of the present application also discloses a computer program product or computer program, where the computer program product or computer program includes computer instructions, and the computer instructions are stored in a computer-readable storage medium. The processor of the computer device can read the computer instruction from the computer-readable storage medium, and the processor executes the computer instruction, so that the computer device executes the method shown in FIG. 1 .

本实施例还提供了一种存储介质,存储有可执行本发明方法实施例所提供的一种基于多源信息融合的行人室内导航方法的指令或程序,当运行该指令或程序时,可执行方法实施例的任意组合实施步骤,具备该方法相应的功能和有益效果。This embodiment also provides a storage medium, which stores an instruction or program that can execute a pedestrian indoor navigation method based on multi-source information fusion provided by the method embodiment of the present invention. When the instruction or program is executed, it can execute Any combination of implementation steps of the method embodiments has the corresponding functions and beneficial effects of the method.

在一些可选择的实施例中,在方框图中提到的功能/操作可以不按照操作示图提到的顺序发生。例如,取决于所涉及的功能/操作,连续示出的两个方框实际上可以被大体上同时地执行或所述方框有时能以相反顺序被执行。此外,在本发明的流程图中所呈现和描述的实施例以示例的方式被提供,目的在于提供对技术更全面的理解。所公开的方法不限于本文所呈现的操作和逻辑流程。可选择的实施例是可预期的,其中各种操作的顺序被改变以及其中被描述为较大操作的一部分的子操作被独立地执行。In some alternative implementations, the functions/operations noted in the block diagrams may occur out of the order noted in the operational diagrams. For example, two blocks shown in succession may, in fact, be executed substantially concurrently or the blocks may sometimes be executed in the reverse order, depending upon the functionality/operations involved. Furthermore, the embodiments presented and described in the flowcharts of the present invention are provided by way of example in order to provide a more comprehensive understanding of the technology. The disclosed methods are not limited to the operations and logical flow presented herein. Alternative embodiments are contemplated in which the order of various operations is changed and in which sub-operations described as part of larger operations are performed independently.

此外,虽然在功能性模块的背景下描述了本发明,但应当理解的是,除非另有相反说明,所述的功能和/或特征中的一个或多个可以被集成在单个物理装置和/或软件模块中,或者一个或多个功能和/或特征可以在单独的物理装置或软件模块中被实现。还可以理解的是,有关每个模块的实际实现的详细讨论对于理解本发明是不必要的。更确切地说,考虑到在本文中公开的装置中各种功能模块的属性、功能和内部关系的情况下,在工程师的常规技术内将会了解该模块的实际实现。因此,本领域技术人员运用普通技术就能够在无需过度试验的情况下实现在权利要求书中所阐明的本发明。还可以理解的是,所公开的特定概念仅仅是说明性的,并不意在限制本发明的范围,本发明的范围由所附权利要求书及其等同方案的全部范围来决定。Furthermore, although the invention has been described in the context of functional modules, it should be understood that one or more of the described functions and/or features may be integrated into a single physical device and/or unless stated to the contrary. or software modules, or one or more functions and/or features may be implemented in separate physical devices or software modules. It will also be appreciated that a detailed discussion of the actual implementation of each module is not necessary to understand the present invention. Rather, given the attributes, functions and internal relationships of the various functional blocks in the devices disclosed herein, the actual implementation of the blocks will be within the ordinary skill of the engineer. Accordingly, those skilled in the art can implement the present invention set forth in the claims without undue experimentation using ordinary techniques. It is also to be understood that the particular concepts disclosed are illustrative only and are not intended to limit the scope of the invention which is to be determined by the appended claims and their full scope of equivalents.

所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。If the functions described above are realized in the form of software function units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the essence of the technical solution of the present invention or the part that contributes to the prior art or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including Several instructions are used to make a computer device (which may be a personal computer, a server, or a network device, etc.) execute all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disk or optical disk and other media that can store program codes. .

在流程图中表示或在此以其他方式描述的逻辑和/或步骤,例如,可以被认为是用于实现逻辑功能的可执行指令的定序列表,可以具体实现在任何计算机可读介质中,以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用,或结合这些指令执行系统、装置或设备而使用。就本说明书而言,“计算机可读介质”可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。The logic and/or steps represented in the flowcharts or otherwise described herein, for example, can be considered as a sequenced listing of executable instructions for implementing logical functions, can be embodied in any computer-readable medium, For use with instruction execution systems, devices, or devices (such as computer-based systems, systems including processors, or other systems that can fetch instructions from instruction execution systems, devices, or devices and execute instructions), or in conjunction with these instruction execution systems, devices or equipment used. For the purposes of this specification, a "computer-readable medium" may be any device that can contain, store, communicate, propagate or transmit a program for use in or in conjunction with an instruction execution system, device or device.

计算机可读介质的更具体的示例(非穷尽性列表)包括以下:具有一个或多个布线的电连接部(电子装置),便携式计算机盘盒(磁装置),随机存取存储器(RAM),只读存储器(ROM),可擦除可编辑只读存储器(EPROM或闪速存储器),光纤装置,以及便携式光盘只读存储器(CDROM)。另外,计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质,因为可以例如通过对纸或其他介质进行光学扫描,接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序,然后将其存储在计算机存储器中。More specific examples (non-exhaustive list) of computer-readable media include the following: electrical connection with one or more wires (electronic device), portable computer disk case (magnetic device), random access memory (RAM), Read Only Memory (ROM), Erasable and Editable Read Only Memory (EPROM or Flash Memory), Fiber Optic Devices, and Portable Compact Disc Read Only Memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program can be printed, as it may be possible, for example, by optically scanning the paper or other medium, followed by editing, interpreting, or other suitable processing if necessary. The program is processed electronically and stored in computer memory.

应当理解,本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中,多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。例如,如果用硬件来实现,和在另一实施方式中一样,可用本领域公知的下列技术中的任一项或他们的组合来实现:具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路,具有合适的组合逻辑门电路的专用集成电路,可编程门阵列(PGA),现场可编程门阵列(FPGA)等。It should be understood that various parts of the present invention can be realized by hardware, software, firmware or their combination. In the embodiments described above, various steps or methods may be implemented by software or firmware stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, it can be implemented by any one or combination of the following techniques known in the art: Discrete logic circuits, ASICs with suitable combinational logic gates, programmable gate arrays (PGAs), field programmable gate arrays (FPGAs), etc.

在本说明书的上述描述中,参考术语“一个实施方式/实施例”、“另一实施方式/实施例”或“某些实施方式/实施例”等的描述意指结合实施方式或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施方式或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施方式或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施方式或示例中以合适的方式结合。In the above description of this specification, the description with reference to the terms "one embodiment/example", "another embodiment/example" or "some embodiments/example" means that the description is described in conjunction with the embodiment or example. A particular feature, structure, material, or characteristic is included in at least one embodiment or example of the invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the described specific features, structures, materials or characteristics may be combined in any suitable manner in any one or more embodiments or examples.

尽管已经示出和描述了本发明的实施方式,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施方式进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同物限定。Although the embodiments of the present invention have been shown and described, those skilled in the art can understand that various changes, modifications, substitutions and variations can be made to these embodiments without departing from the principle and spirit of the present invention. The scope of the invention is defined by the claims and their equivalents.

以上是对本发明的较佳实施进行了具体说明,但本发明并不限于上述实施例,熟悉本领域的技术人员在不违背本发明精神的前提下还可做作出种种的等同变形或替换,这些等同的变形或替换均包含在本申请权利要求所限定的范围内。The above is a specific description of the preferred implementation of the present invention, but the present invention is not limited to the above-mentioned embodiments, and those skilled in the art can also make various equivalent deformations or replacements without violating the spirit of the present invention. Equivalent modifications or replacements are all within the scope defined by the claims of the present application.

Claims (10)

1. The pedestrian indoor navigation method based on multi-source information fusion is characterized by comprising the following steps of:
constructing map information of the current indoor environment, selecting a reference point on a map, and constructing a Wi-Fi fingerprint library;
acquiring multi-source information by using a sensor on an intelligent terminal, wherein the multi-source information comprises Wi-Fi signal strength, angular velocity, speed, acceleration and gesture information;
constructing a pedestrian position model according to the Wi-Fi signal intensity and the Wi-Fi fingerprint library to obtain pedestrian position information; judging the instantaneous gesture of the intelligent terminal through the gesture information, converting the acceleration from a carrier coordinate system to a navigation coordinate system, and obtaining the movement distance of the pedestrian through gait detection and step estimation;
constructing a motion model of the pedestrian according to the multi-source information and the motion distance to obtain motion information of the pedestrian; according to an extended Kalman filtering algorithm, a fusion model based on the extended Kalman filtering is constructed, multi-source information and pedestrian position information are fused, fused position information is obtained, and state vector updating and optimal estimation of error values in the extended Kalman filtering are completed;
and further correcting and updating the pedestrian position information according to the map information and the particle filter algorithm by utilizing the fused position information, so as to complete the pedestrian indoor navigation with multi-source information fusion.
2. The pedestrian indoor navigation method based on multi-source information fusion according to claim 1, wherein the expression of the pedestrian position model is:
wherein,,represents the estimated position, x (t k )、y(t k ) Representing the position coordinates, ω, corresponding to the first K position points j (t k ) Weight coefficients representing the corresponding data points, the weight coefficients based on signal strength may represent +.>RSSI is the signal strength.
3. The pedestrian indoor navigation method based on multi-source information fusion according to claim 1, wherein the acceleration is converted from the carrier coordinate system to the navigation coordinate system by:
wherein,,for navigating a directional cosine matrix from the navigation system to the carrier system, C r C for rotating matrix around Y axis p For a rotation matrix about the Z-axis, C y For a rotation matrix about the X-axis, γ is the rotation angle about the Y-axis, p is the rotation angle about the X-axis, and Y is the rotation angle about the Z-axis.
4. The pedestrian indoor navigation method based on multi-source information fusion of claim 1, wherein the step of gait detection comprises:
calculating the overall acceleration:
wherein a is x ,a y ,a z For the collected triaxial acceleration value, the influence of the sensor posture is reduced by calculating the overall acceleration a;
according to the overall acceleration, a potential peak value is obtained in the sliding window, and primary judgment is carried out by utilizing a preset acceleration threshold value;
calculating the time difference between the potential wave crest and the previous wave crest, and performing secondary judgment by using a preset walking one-step time threshold range;
and comparing the potential peak with the front and rear neighborhood acceleration, judging for three times to remove the pseudo peak, if the potential peak point is the maximum value, recording one step by the algorithm, otherwise, not performing step counting processing.
5. The pedestrian indoor navigation method based on multi-source information fusion of claim 1, wherein the step size estimation is:
where SL is the step size, a max And a min The maximum and minimum of the longitudinal acceleration of the pedestrian, respectively.
6. The pedestrian indoor navigation method based on multi-source information fusion of claim 1, wherein the extended kalman filtering algorithm comprises:
X k+1 =f(X k )+w k
wherein X is k+1 X is the state variable of the (k+1) th step k Is the state variable of the kth step, w k Is process noise; f () is a nonlinear state function in extended kalman filtering;
the method for realizing the linearization of the system matrix comprises the following steps:
wherein Z is k To observe the variable, v k H () is a measurement function in extended kalman filtering for observing noise;
the method for realizing the linearization of the observation matrix comprises the following steps:
wherein H is k Is a jacobian matrix of h (X) to X bias,and (5) estimating the state of the kth step.
7. The pedestrian indoor navigation method based on multi-source information fusion of claim 1, wherein the step of expanding the state vector update in the kalman filter comprises:
one-step state prediction update:
in the method, in the process of the invention,is the optimal value of the last state, +.>A one-step predicted value for the current state;
one-step predictive estimation error covariance matrix update:
wherein P is k+1|k Is thatCorresponding covariance one-step predictor, phi k+1|k For state transition matrix, P k|k Is->Corresponding covariance,/>Is phi k+1|k Transpose of Q k A process noise covariance matrix;
calculating an extended Kalman filtering gain:
wherein K is k+1 To expand the Kalman filter gain matrix, H k+1 In order to observe the matrix,for transpose of the observation matrix, R k+1 Measuring a noise covariance matrix;
calculating information from the observation vectors and updating the state estimate:
in the method, in the process of the invention,z is the optimal estimated value of the current state k+1 For the current observation +.>The current observation predicted value;
updating the estimated error covariance:
P k+1 =[I-K k+1 H k+1 ]P k+1|k
wherein P is k+1 Is thatCorresponding covariance, I is the identity matrix.
8. The pedestrian indoor navigation method based on multi-source information fusion according to claim 1, wherein the further correcting and updating of the pedestrian position information according to the map information and the particle filter algorithm comprises:
the state transfer equation and the observation equation for particle filtering are as follows:
x(t)=f(x(t-1),u(t),w(t))
y(t)=f(x(t),e(t))
wherein x (t) is a state at time t, u (t) is a control quantity, w (t) and e (t) are state noise and observation noise respectively, and particle filtering filters the state at time t x (t) from observation y (t) and the state at the last time x (t-1), u (t) and w (t);
the updating method in particle filtering mainly comprises the following steps:
a1, initializing: initial value x of observed value i (k) As an initial value of the probability density function x (0), byFor->Sampling N times, initial w i (k) Set to->
A2, one-step prediction: for each particle x i (k) By converting the formula p (x (k+1) |x i (k) Obtaining a new particle;
a3, importance sampling: for any particle x i (k+1) solving for their weights w i (k+1)=p(z(k+1)|x i (k+1));
A4, normalization: and (3) carrying out normalization processing on the weight values:
a5, resampling: resampling the particles according to the weight, wherein the probability of repeating the particles with the weight is high, and the probability of repeating the particles with the weight is low;
a6, updating particles and simultaneously adding map information:
order theIs->
Because of the limit value of the indoor walking area, map matching is carried out according to map information, and when the particles predicted in one step cross the passable area, the weight value of the particles is given as 0:
and A7, calculating by adopting a basic particle filtering algorithm, and completing the indoor positioning and navigation of the pedestrians according to the calculation result.
9. A pedestrian indoor navigation device based on multi-source information fusion, comprising:
at least one of a processor;
at least one memory for storing at least one program;
the at least one program, when executed by the at least one processor, causes the at least one processor to implement the method of any one of claims 1-8.
10. A computer readable storage medium, in which a processor executable program is stored, characterized in that the processor executable program is for performing the method according to any of claims 1-8 when being executed by a processor.
CN202310210607.1A 2023-03-06 2023-03-06 Pedestrian indoor navigation method, device and medium based on multi-source information fusion Pending CN116448111A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310210607.1A CN116448111A (en) 2023-03-06 2023-03-06 Pedestrian indoor navigation method, device and medium based on multi-source information fusion

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310210607.1A CN116448111A (en) 2023-03-06 2023-03-06 Pedestrian indoor navigation method, device and medium based on multi-source information fusion

Publications (1)

Publication Number Publication Date
CN116448111A true CN116448111A (en) 2023-07-18

Family

ID=87129134

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310210607.1A Pending CN116448111A (en) 2023-03-06 2023-03-06 Pedestrian indoor navigation method, device and medium based on multi-source information fusion

Country Status (1)

Country Link
CN (1) CN116448111A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118408552A (en) * 2024-07-01 2024-07-30 山东科技大学 Multi-sensor positioning method for unmanned ships based on dual-threshold event triggering mechanism
CN118444647A (en) * 2024-07-08 2024-08-06 微晶数实(山东)装备科技有限公司 Intelligent device control method, system and medium based on multi-source information fusion

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118408552A (en) * 2024-07-01 2024-07-30 山东科技大学 Multi-sensor positioning method for unmanned ships based on dual-threshold event triggering mechanism
CN118444647A (en) * 2024-07-08 2024-08-06 微晶数实(山东)装备科技有限公司 Intelligent device control method, system and medium based on multi-source information fusion
CN118444647B (en) * 2024-07-08 2024-09-27 微晶数实(山东)装备科技有限公司 Intelligent device control method, system and medium based on multi-source information fusion

Similar Documents

Publication Publication Date Title
CN112639502B (en) Robot pose estimation
CN111402339B (en) Real-time positioning method, device, system and storage medium
CN108759833B (en) Intelligent vehicle positioning method based on prior map
WO2018018994A1 (en) Method and system for indoor positioning
CN112965063B (en) Robot mapping and positioning method
CN114323033B (en) Positioning method and equipment based on lane lines and feature points and automatic driving vehicle
CN107036594A (en) The positioning of intelligent Power Station inspection intelligent body and many granularity environment perception technologies
CN116448111A (en) Pedestrian indoor navigation method, device and medium based on multi-source information fusion
CN110763239B (en) Filter combination laser SLAM mapping method and device
CN112444246B (en) Laser fusion positioning method in high-precision digital twin scene
CN113739795B (en) Underwater synchronous positioning and mapping method based on polarized light/inertia/vision integrated navigation
CN109141413B (en) EFIR filtering algorithm and system with data missing UWB pedestrian positioning
CN111964667B (en) geomagnetic-INS (inertial navigation System) integrated navigation method based on particle filter algorithm
CN108387236B (en) A Polarized Light SLAM Method Based on Extended Kalman Filtering
CN112797985A (en) Indoor positioning method and indoor positioning system based on weighted extended Kalman filter
CN114061611A (en) Target object positioning method, apparatus, storage medium and computer program product
CN115435775A (en) Multi-sensor Fusion SLAM Method Based on Extended Kalman Filter
CN113379915A (en) Driving scene construction method based on point cloud fusion
CN115711616B (en) Smooth positioning method and device for indoor and outdoor traversing unmanned aerial vehicle
CN115435782A (en) Anti-interference position estimation method and device under multi-source information constraint
CN114047766B (en) Mobile robot data collection system and method for long-term application in indoor and outdoor scenes
CN117687006A (en) External parameter calibration method from laser radar to inertial measurement unit and device thereof
CN115900697B (en) Object motion trail information processing method, electronic equipment and automatic driving vehicle
CN117075158A (en) Position and orientation estimation method and system of unmanned deformable motion platform based on lidar
CN116929343A (en) Pose estimation method, related equipment and storage medium

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