Disclosure of Invention
The invention aims to provide an IMU/Wi-Fi signal ultra-tight combination indoor navigation method based on CKF, and aims to solve the problems of poor stability and poor filtering precision of the existing indoor pedestrian combination navigation method.
The implementation of the invention comprises the following steps:
the method comprises the following steps: selecting a certain density indoorsMeasure and record the signal strength received from the access point APs and the corresponding location information P for each reference pointDBEstablishing an RSS Fingerprint database;
step two: the IMU is fixed on the foot of the pedestrian, the IMU is used for measuring the acceleration a and the rotating speed omega of the motion information of the pedestrian during the advancing process, and the velocity V of the pedestrian is obtained through inertial navigation calculationIMUPosition PIMUAnd attitude AIMU;
Step three: utilizing the acceleration a obtained in the step two to carry out zero-speed correction on the accelerometer, carrying out inertial navigation solution on the corrected motion information, and obtaining the corrected position velocity V through filtering of the extended filterJIMUAnd the velocity V obtained in the step twoIMUObtaining a speed error delta V by difference;
step four: using the position information P obtained in step twoIMUPosition information P of reference point in RSS finger rprint databaseDBComparing to screen out n points near the measuring pointAnd then respectively with the IMU measured position PIMUMaking difference to obtain the distance between pedestrian and reference point
To which is relatedThe expression is as follows:
wherein, XIMUAnd YIMUIs measured by IMU to obtain position coordinates PIMUAnd has a component ofIMU={XIMU,YIMU};Andis the coordinates of the selected reference pointAnd has a component of
Step five: the method comprises the steps of receiving Wi-Fi signals by using a Wi-Fi antenna carried by a pedestrian, measuring the intensity of the received signals, calculating the distance d between a measuring point and an access point APs according to a mathematical model of the signal intensity and the distance, and calculating the position P of the measuring point by using a trilateration methodWi-FiAnd making a difference with the position measured by the IMU to obtain a position error delta P;
the mathematical model between the referred Wi-Fi signal strength RSS and distance d is:
in the formula (d)0For a known reference distance, RSS0To be at a reference distance d0The average signal strength, p is the signal attenuation exponent;
step six: using the position information P from step fiveWi-FiPosition information P of reference point in RSS finger rprint databaseDBComparing to screen out n points near the measuring pointThen respectively locating the obtained positions P with Wi-FiWi-FiMaking difference to obtain the distance between pedestrian and reference pointNamely:
to which is relatedThe expression is as follows:
wherein, XWi-FiAnd YWi-FiIs the measured position coordinate PWi-FiAnd has a component ofWi-Fi={XWi-Fi,YWi-Fi};Andis the selected reference point coordinatesAnd has a component of
Step seven: the distances obtained in the fourth step and the sixth stepMaking a difference to obtain a distance difference delta di;
Step eight: using the velocity error Δ V obtained in step three and the position error Δ P obtained in step five as the state quantity X of CKF ═ Δ V, Δ P]Using the difference Δ d obtained in step seveniMeasuring Z ═ d as CKF quantityi]Performing filtering calculation by using CKF;
step nine: correcting the speed V and the position P measured by the IMU by using the filtered speed error delta V and the filtered position error delta P obtained in the step eight, and comparing the filtered position error delta P obtained in the step eight with the position P of the measuring point in the step fiveWi-FiCorrecting to obtain final indoor navigation information of the pedestrian
In the first step, when a reference point is selected, one reference point is taken every 5 meters, one reference point is taken every 10 meters, 5 groups of data are measured in each direction at each reference point, 40 groups of data are measured for processing, and the Wi-Fi signal strength of the reference point is determined.
The invention has the beneficial effects that:
the method has the advantages of high convergence speed, high filtering precision and strong robustness, so that the indoor navigation precision of the pedestrian can be effectively improved.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail with reference to the following embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
The invention discloses an IMU/Wi-Fi signal ultra-tight combination indoor navigation method based on CKF, which comprises the following steps: establishing an RSS finger printing database of the Wi-Fi signals; utilizing IMU to carry out positioning to obtain the position P of the pedestrianIMUAnd information such as speed, acceleration, and attitude; correcting the IMU accelerometer drift through zero speed correction, and obtaining corrected speed; positioning by utilizing Wi-Fi signals to obtain the position P of the pedestrianWi-Fi(ii) a Constructed numberRespectively finding P in databaseIMUCloser n reference points and PWi-FiThe n closer reference points respectively obtain the distance between the measuring point and the reference point according to the positions of the known reference pointsThereby obtaining a distance difference Deltadi(ii) a And performing information fusion by using a Cubature Kalman Filter (CKF) to finally obtain the position, acceleration, speed and attitude information of the pedestrian. According to the invention, information fusion is carried out on the pedestrian navigation system which is formed by tightly combining the indoor Wi-Fi signal and the IMU by using the CKF, the Wi-Fi is widely applied indoors and the Wi-Fi signal is relatively stable, and the pedestrian navigation system can be effectively assisted to navigate pedestrians by the inertial navigation system. Compared with EKF, UKF and other filtering algorithms, the CKF algorithm has the advantages of high filtering precision, high convergence speed and strong robustness, so that the method can effectively improve the navigation precision.
The IMU/Wi-Fi signal tight combination indoor navigation method based on the CKF comprises the following steps:
the method comprises the following steps: selecting reference Points with certain density indoors, measuring and recording the Received Signal Strength (Received Signal Strength RSS) and corresponding position information P of each reference point from access Points (Access Points aps)DBEstablishing an RSS Fingerprint database;
step two: the IMU is fixed on one foot of the pedestrian, the IMU is used for measuring the acceleration a and the rotating speed omega of the movement information of the pedestrian during the advancing process, and the velocity V of the pedestrian is obtained through inertial navigation calculationIMUPosition PIMUAnd attitude AIMU;
Step three: utilizing the acceleration a obtained in the step two to carry out zero-speed correction on the accelerometer, carrying out inertial navigation solution on the corrected motion information, and obtaining the corrected position velocity V through filtering of the extended filterJIMUAnd the velocity V obtained in the step twoIMUObtaining a speed error delta V by difference;
step four: using the position information P obtained in step twoIMUAnd RSS Fingerposition information P of reference point in print databaseDBComparing to screen out n points near the measuring pointAnd then respectively with the IMU measured position PIMUMaking difference to obtain the distance between pedestrian and reference pointNamely:
wherein, XIMUAnd YIMUIs measured by IMU to obtain position coordinates PIMUAnd has a component ofIMU={XIMU,YIMU};Andis the coordinates of the selected reference pointAnd has a component of
Step five: the method comprises the steps of receiving Wi-Fi signals by using a Wi-Fi antenna carried by a pedestrian, measuring the intensity of the received signals, calculating the distance d between a measuring point and an access point APs according to a mathematical model of the signal intensity and the distance, and calculating the position P of the measuring point by using a trilateration methodWi-FiAnd making a difference with the position measured by the IMU to obtain a position error delta P;
by the formula:
establishing a mathematical model between Wi-Fi signal strength RSS and distance d, where d0For a known reference distance, RSS0To be at a reference distance d0Where p is the signal attenuation exponent, v is a gaussian distributed random variable with a mean of 0 and a mean square error σRSSIndoor obstacles. The maximum likelihood estimation is carried out on the formula to obtain a mathematical model between the Wi-Fi signal strength RSS and the distance d as follows:
the distance d between the measurement point and the APs can be obtained by the above formula, so that the position P of the measurement point is solved by using the trilateration methodWi-FiAnd obtaining a position error delta P by subtracting the position measured by the IMU;
step six: using the position information P from step fiveWi-FiPosition information P of reference point in RSS finger rprint databaseDBComparing to screen out n points near the measuring pointThen respectively locating the obtained positions P with Wi-FiWi-FiMaking difference to obtain the distance between pedestrian and reference pointNamely:
wherein, XWi-FiAnd YWi-FiIs the measured position coordinate PWi-FiAnd has a component ofWi-Fi={XWi-Fi,YWi-Fi};Andis the selected reference point coordinatesAnd has a component of
Step seven: the distances obtained in the fourth step and the sixth stepMaking a difference to obtain a distance difference delta di;
Step eight: using the velocity error Δ V obtained in step three and the position error Δ P obtained in step five as the state quantity X of CKF ═ Δ V, Δ P]Using the difference Δ d obtained in step seveniMeasuring Z ═ d as CKF quantityi]Filtering by using CKF;
step nine: correcting the speed V and the position P measured by the IMU by using the filtered speed error delta V and the filtered position error delta P obtained in the step eight, and comparing the filtered position error delta P obtained in the step eight with the position P of the measuring point in the step fiveWi-FiAnd correcting to obtain final indoor navigation information of the pedestrian.
In the first step, when a reference point is selected, one reference point is taken every 5 meters in a room, one reference point is taken every 10 meters in a corridor, 5 groups of data are measured in each direction at each reference point, 40 groups of data are measured for processing, and the Wi-Fi signal strength of the reference point is determined.
In step eight, a volume kalman filter (CKF) is used for filtering, and the velocity error Δ V and the position error Δ P obtained in step three and step five are used as the state quantity X of the volume kalman filter, [ Δ V, Δ P ═ V]Distance difference Δ d obtained in step seveniMeasurement of Z ═ d as a volumetric Kalman filteri]The processing flow of the CKF algorithm is as follows:
1) time updating
Selecting state volume points:
wherein S isk-1=chol{Pk-1Is a matrix Pk-1Cholesky decomposition of (i.e. error covariance matrix)
State volume point propagated through state equation:
and (3) state prediction:
error covariance prediction:
2) measurement update
State volume point:
wherein,
measuring volume points:
Zj,k=h(Xj,k)
measurement and prediction:
innovation variance:
covariance matrix:
kalman gain:
and (3) updating the state:
error covariance:
FIG. 1 shows a flow of an IMU/Wi-Fi signal ultra-tight combination indoor navigation method based on CKF provided by the invention. For convenience of explanation, only portions relevant to the present invention are shown.
The IMU/Wi-Fi signal ultra-tight combination indoor navigation method based on the CKF comprises the following steps:
the method comprises the following steps: selecting reference Points with certain density indoors, measuring and recording the Received Signal Strength (Received Signal Strength RSS) and corresponding position information P of each reference point from access Points (Access Points aps)DBEstablishing an RSS Fingerprint database;
step two: the IMU is fixed on one foot of the pedestrian, the IMU is used for measuring the acceleration a and the rotating speed omega of the movement information of the pedestrian during the advancing process, and the velocity V of the pedestrian is obtained through inertial navigation calculationIMUPosition PIMUAnd attitude AIMU;
Step three: utilizing the acceleration a obtained in the step two to carry out zero-speed correction on the accelerometer, carrying out inertial navigation solution on the corrected motion information, and obtaining the corrected position velocity V through filtering of the extended filterJIMUAnd the velocity V obtained in the step twoIMUObtaining a speed error delta V by difference;
step four: using the position information P obtained in step twoIMUPosition information P of reference point in RSS finger rprint databaseDBComparing to screen out n points near the measuring pointAnd then respectively with the IMU measured position PIMUMaking difference to obtain the distance between pedestrian and reference pointNamely:
wherein, XIMUAnd YIMUIs measured by IMU to obtain position coordinates PIMUAnd has a component ofIMU={XIMU,YIMU};Andis the coordinates of the selected reference pointAnd has a component of
Step five: the method comprises the steps of receiving Wi-Fi signals by using a Wi-Fi antenna carried by a pedestrian, measuring the intensity of the received signals, calculating the distance d between a measuring point and an access point APs according to a mathematical model of the signal intensity and the distance, and calculating the position P of the measuring point by using a trilateration methodWi-FiAnd making a difference with the position measured by the IMU to obtain a position error delta P;
by the formula:
establishing a mathematical model between Wi-Fi signal strength RSS and distance d, where d0For a known reference distance, RSS0To be at a reference distance d0Where p is the signal attenuation exponent, v is a gaussian distributed random variable with a mean of 0 and a mean square error σRSSIndoor obstacles. The maximum likelihood estimation is carried out on the formula to obtain a mathematical model between the Wi-Fi signal strength RSS and the distance d as follows:
the distance d between the measurement point and the APs can be obtained by the above formula, so that the position P of the measurement point is solved by using the trilateration methodWi-FiAnd obtaining a position error delta P by subtracting the position measured by the IMU;
step six: by usingPosition information P obtained in step fiveWi-FiPosition information P of reference point in RSS finger rprint databaseDBComparing to screen out n points near the measuring pointThen respectively locating the obtained positions P with Wi-FiWi-FiMaking difference to obtain the distance between pedestrian and reference pointNamely:
wherein, XWi-FiAnd YWi-FiIs the measured position coordinate PWi-FiAnd has a component ofWi-Fi={XWi-Fi,YWi-Fi};Andis the selected reference point coordinatesAnd has a component of
Step seven: the distances obtained in the fourth step and the sixth stepMaking a difference to obtain a distance difference delta di;
Step eight: using the velocity error Δ V obtained in step three and the position error Δ P obtained in step five as the state quantity X of CKF ═ Δ V, Δ P]Using the difference Δ d obtained in step seveniAs a measure of CKFAmount Z ═ di]Filtering by using CKF;
step nine: correcting the speed V and the position P measured by the IMU by using the filtered speed error delta V and the filtered position error delta P obtained in the step eight, and comparing the filtered position error delta P obtained in the step eight with the position P of the measuring point in the step fiveWi-FiAnd correcting to obtain final indoor navigation information of the pedestrian.
As an optimization scheme of the embodiment of the invention, in the step one, when the reference point is selected, one reference point is taken every 5 meters in a room, one reference point is taken every 10 meters in a corridor, 5 groups of data are measured in each direction at each reference point, 40 groups of data are measured for processing, and the Wi-Fi signal strength of the reference point is determined. In order to make the RSS Fingerprint database more accurate, the reference point is measured twice, once when someone walks indoors and once when no one exists indoors.
As an optimized solution of the embodiment of the present invention, in step eight, filtering is performed by using a volume kalman filter (CKF), and the speed error Δ V and the position error Δ P obtained in step three and step five are used as the state quantity X ═ Δ V, Δ P of the volume kalman filter]Distance difference Δ d obtained in step seveniMeasurement of Z ═ d as a volumetric Kalman filteri]The processing flow of the CKF algorithm is as follows:
1) time updating
Selecting state volume points:
wherein S isk-1=chol{Pk-1Is a matrix Pk-1Cholesky decomposition of (i.e. error covariance matrix)
State volume point propagated through state equation:
and (3) state prediction:
error covariance prediction:
2) measurement update
State volume point:
wherein,
measuring volume points:
Zj,k=h(Xj,k)
measurement and prediction:
innovation variance:
covariance matrix:
kalman gain:
and (3) updating the state:
error covariance:
as shown in fig. 1, the IMU/Wi-Fi signal ultra-tight combination indoor navigation method based on CKF of the embodiment of the present invention includes the following steps:
s101: establishing an RSS finger printing database of the Wi-Fi signals;
s102: the IMU is fixed on one foot of the pedestrian, the IMU is used for measuring the acceleration a and the rotating speed omega of the movement information of the pedestrian during the advancing process, and the velocity V of the pedestrian is obtained through inertial navigation calculationIMUPosition PIMUAnd attitude AIMU;
S103: the obtained acceleration a is used for carrying out zero-speed correction on the accelerometer, inertial navigation calculation is carried out on the corrected motion information, and the corrected position velocity V is obtained through filtering of the expansion filterJIMUAnd the velocity V obtained in the step twoIMUObtaining a speed error delta V by difference;
s104: using the position information P obtained in step twoIMUParticipating in RSS Fingerprint databasePosition information P of examination pointDBComparing, screening out the point near to the measurement point, and comparing with the position P measured by IMUIMUMaking a difference to obtain the distance between the pedestrian and a reference point;
s105: the method comprises the steps of receiving Wi-Fi signals by using a Wi-Fi antenna carried by a pedestrian, measuring the intensity of the received signals, calculating the distance d between a measuring point and an access point APs according to a mathematical model of the signal intensity and the distance, and calculating the position P of the measuring point by using a trilateration methodWi-FiAnd making a difference with the position measured by the IMU to obtain a position error delta P;
s106: using the position information P from step fiveWi-FiPosition information P of reference point in RSS finger rprint databaseDBComparing, screening out the point close to the measuring point, and locating with Wi-Fi to obtain position PWi-FiMaking a difference to obtain the distance between the pedestrian and a reference point;
s107: the distances obtained in the fourth step and the sixth stepMaking difference to obtain distance difference;
s108: using the obtained speed error Δ V and the obtained position error Δ P as the state quantity X of CKF ═ Δ V, Δ P]Using the resulting difference in distance Δ diMeasuring Z ═ d as CKF quantityi]Filtering by using CKF;
s109: correcting the speed V and the position P measured by the IMU by using the filtered speed error delta V and the filtered position error delta P, and correcting the filtered position error delta P to the position P of the measuring pointWi-FiAnd correcting to obtain final indoor navigation information of the pedestrian.
The method comprises the following specific steps:
the method comprises the following steps: selecting reference Points with certain density indoors, measuring and recording the Received Signal Strength (Received Signal Strength RSS) and corresponding position of each reference point from access Points (Access Points aps)Letter PDBEstablishing an RSS Fingerprint database;
when the reference point is selected, one reference point is taken at every 5 steps in a room, one reference point is taken at every 10 steps in a hall, 5 groups of data are measured in each direction at each reference point, 40 groups of data are measured for processing, and finally the Wi-Fi signal strength of the reference point is determined. In order to make the RSS Fingerprint database more accurate, the reference point is measured twice, once when someone walks indoors and once when no one exists indoors.
Step two: the IMU is fixed on one foot of the pedestrian, the IMU is used for measuring the acceleration a and the rotating speed omega of the movement information of the pedestrian during the advancing process, and the velocity V of the pedestrian is obtained through inertial navigation calculationIMUPosition PIMUAnd attitude AIMU;
Step three: utilizing the acceleration a obtained in the step two to carry out zero-speed correction on the accelerometer, carrying out inertial navigation solution on the corrected motion information, and obtaining the corrected position velocity V through filtering of the extended filterJIMUAnd the velocity V obtained in the step twoIMUObtaining a speed error delta V by difference;
step four: using the position information P obtained in step twoIMUPosition information P of reference point in RSS finger rprint databaseDBComparing to screen out n points near the measuring pointAnd then respectively with the IMU measured position PIMUMaking difference to obtain the distance between pedestrian and reference pointNamely:
wherein, XIMUAnd YIMUIs measured by IMU to obtain position coordinates PIMUAnd has a component ofIMU={XIMU,YIMU};Andis the coordinates of the selected reference pointAnd has a component of
Step five: the method comprises the steps of receiving Wi-Fi signals by using a Wi-Fi antenna carried by a pedestrian, measuring the intensity of the received signals, calculating the distance d between a measuring point and an access point APs according to a mathematical model of the signal intensity and the distance, and calculating the position P of the measuring point by using a trilateration methodWi-FiAnd making a difference with the position measured by the IMU to obtain a position error delta P;
by the formula:
establishing a mathematical model between Wi-Fi signal strength RSS and distance d, where d0For a known reference distance, RSS0To be at a reference distance d0Where p is the signal attenuation exponent, v is a gaussian distributed random variable with a mean of 0 and a mean square error σRSSIndoor obstacles. The maximum likelihood estimation is carried out on the formula to obtain a mathematical model between the Wi-Fi signal strength RSS and the distance d as follows:
the distance d between the measurement point and the APs can be obtained by the above formula, so that the measurement point is solved by using the trilateration methodPosition P ofWi-FiAnd obtaining a position error delta P by subtracting the position measured by the IMU;
step six: using the position information P from step fiveWi-FiPosition information P of reference point in RSS finger rprint databaseDBComparing to screen out n points near the measuring pointThen respectively locating the obtained positions P with Wi-FiWi-FiMaking difference to obtain the distance between pedestrian and reference pointNamely:
wherein, XWi-FiAnd YWi-FiIs the measured position coordinate PWi-FiAnd has a component ofWi-Fi={XWi-Fi,YWi-Fi};Andis the selected reference point coordinatesAnd has a component of
Step seven: the distances obtained in the fourth step and the sixth stepMaking a difference to obtain a distance difference delta di;
Step eight: using the velocity obtained in step threeThe error Δ V and the position error Δ P obtained in step five are used as the state quantity X of CKF ═ Δ V, Δ P]Using the difference Δ d obtained in step seveniMeasuring Z ═ d as CKF quantityi]Filtering by using CKF, wherein the CKF algorithm flow is as follows;
1) time updating
Selecting state volume points:
wherein S isk-1=chol{Pk-1Is a matrix Pk-1Cholesky decomposition of (i.e. error covariance matrix)
State volume point propagated through state equation:
and (3) state prediction:
error covariance prediction:
2) measurement update
State volume point:
wherein,
measuring volume points:
Zj,k=h(Xj,k)
measurement and prediction:
innovation variance:
covariance matrix:
kalman gain:
and (3) updating the state:
error covariance:
step nine: using the filtered speed error Δ V obtained in step eight,Correcting the speed V and the position P measured by the IMU by the position error delta P, and comparing the position error delta P obtained in the step eight after filtering with the position P of the measuring point in the step fiveWi-FiAnd correcting to obtain final indoor navigation information of the pedestrian.