CN112394321A - Multi-base-station real-time positioning method and system based on Bluetooth signals - Google Patents
Multi-base-station real-time positioning method and system based on Bluetooth signals Download PDFInfo
- Publication number
- CN112394321A CN112394321A CN202110078642.3A CN202110078642A CN112394321A CN 112394321 A CN112394321 A CN 112394321A CN 202110078642 A CN202110078642 A CN 202110078642A CN 112394321 A CN112394321 A CN 112394321A
- Authority
- CN
- China
- Prior art keywords
- signal
- base station
- representing
- real
- base
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/06—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
The invention relates to the field of Bluetooth signal positioning, in particular to a multi-base-station real-time positioning method and system based on Bluetooth signals. The multi-base-station real-time positioning method comprises the following steps: step S1, acquiring signal receiving intensity, azimuth angle, pitch angle and spatial spectrum peak value of all base station receiving signals at the current moment; step S2, screening the signal receiving intensity, the azimuth angle, the pitch angle and the spatial spectrum peak value; step S3, acquiring a first predicted position of the signal at the current moment; and step S4, determining the real-time positioning result of the current time according to the first predicted position and the second predicted position. The technical scheme of the invention has the beneficial effects that: the method and the system are applied to a plurality of base stations, azimuth angles and pitch angles which are not accurate enough and can cause misjudgment are screened out, and real-time positioning results of signals are obtained according to the screened signal receiving intensity, the screened azimuth angles, the screened pitch angles and spatial spectrum peak values.
Description
Technical Field
The invention relates to the field of Bluetooth signal positioning, in particular to a multi-base-station real-time positioning method and system based on Bluetooth signals.
Background
The GPS, i.e. global positioning system, is a positioning system of high-precision radio navigation based on air satellite, it can provide accurate geographical position, vehicle speed and accurate time information anywhere in the world and in the near-earth space, the distribution of the satellite can make more than 4 satellites be observed anywhere in the world and anytime, and can maintain the geometry of good positioning resolving precision, thus providing continuous global navigation capability in time. However, the GPS system relies on line-of-sight transmission between the satellites and the receiver, and buildings attenuate GPS satellite signals greatly, so that the results of indoor positioning cannot be obtained accurately.
Disclosure of Invention
In view of the above problems in the prior art, a method and a system for multi-base-station real-time positioning based on bluetooth signals are provided.
The multi-base-station real-time positioning method based on the Bluetooth signals is applied to a plurality of base stations, and all the base stations synchronously receive a signal; the multi-base-station real-time positioning method comprises the following steps:
step S1, acquiring the signal receiving intensity, azimuth angle, pitch angle and spatial spectrum peak value of the signal received by all the base stations at the current moment;
step S2, screening the signal reception intensity, the azimuth angle, the pitch angle, and the spatial spectrum peak;
step S3, performing combined positioning according to the filtered signal receiving intensity, the azimuth angle, the pitch angle and the spatial spectrum peak value to obtain a first predicted position of the signal at the current moment;
step S4, determining a second predicted position of the current time according to the real-time positioning result of the previous time, and determining the real-time positioning result of the current time according to the first predicted position and the second predicted position.
Preferably, the step S1 includes:
step S11, acquiring the signal received by the base station at the current time;
step S12, acquiring the corresponding signal receiving strength according to the energy in the signal;
step S13, obtaining a corresponding spatial spectrum according to orthogonality between the noise subspace in the signal and the array manifold vector, performing peak search on the spatial spectrum to obtain a spatial spectrum peak value and a corresponding angle, and using the corresponding angle as the azimuth angle and the pitch angle at which the base station receives the signal.
Preferably, the step S2 includes:
step S21, obtaining a mean value and a standard deviation of the signal reception intensity at the current time;
step S22, presetting a coefficient;
step S23, sequentially determining whether the difference between the received signal strength and the mean is greater than the product of the coefficient and the standard deviation:
if yes, the signal reception strength is screened out, and then the step S24 is carried out;
if not, the signal receiving strength is retained, and then the step S24 is carried out;
step S24, determining whether all the signal reception intensities have been filtered:
if yes, go to step S25;
if not, returning to the step S23;
step S25, screening the azimuth angle, the pitch angle, and the spatial spectrum peak corresponding to the signal reception intensity according to the screened signal reception intensity.
Preferably, the step S3 includes:
step S31, acquiring the position information of all the base stations;
step S32, sequentially obtaining a first tangent value and a second tangent value according to the position information of each base station and the azimuth angle and the pitch angle at which the base station receives the signal, and respectively constructing a first matrix and a second matrix according to the position information, the first tangent value and the second tangent value;
step S33, obtaining a corresponding weight matrix according to the spatial spectrum peak value of the signal received by each base station;
step S34, obtaining the first predicted position according to the first matrix, the second matrix and the weight matrix.
Preferably, the first tangent value and the second tangent value are expressed by the following formulas:
wherein the content of the first and second substances,
for representing the qth of said first tangent, where q is used to represent the ranking number of said base station among all base stations,for indicating the azimuth angle at which the q-th base station receives the signal;
for representing the second tangent value, whereinSaid pitch angle representing the reception of said signal by the qth of said base station;
x is used to represent the abscissa in a first predicted position of the signal;
y is used to represent the ordinate in a first predicted position of the signal;
Z is used to represent the vertical coordinate in a first predicted location of the signal;
Preferably, the first matrix is expressed by the following formula:
wherein the content of the first and second substances,
for representing the qth of said first tangent, where q is used to represent the ranking number of said base station among all base stations,for indicating the azimuth angle at which the q-th base station receives the signal;
for representing the second tangent value, whereinSaid pitch angle representing the reception of said signal by the qth of said base station;
l is used to indicate the number after screening.
Preferably, the second matrix is expressed by the following formula:
wherein the content of the first and second substances,
for representing the qth of said first tangent, where q is used to represent the ranking number of said base station among all base stations,for indicating the azimuth angle at which the q-th base station receives the signal;
for representing the second tangent value, whereinSaid pitch angle representing the reception of said signal by the qth of said base station;
x is used to represent the abscissa in a first predicted position of the signal;
y is used to represent the ordinate in a first predicted position of the signal;
z is used to represent the vertical coordinate in a first predicted location of the signal;
l is used to indicate the number after screening.
Preferably, the weight is expressed as:
wherein the content of the first and second substances,
is used for representing the q-th weight, wherein q is used for representing the sequence number of the base station in all the base stations;
l is used to indicate the number after screening.
Preferably, the weight matrix is represented as:
wherein the content of the first and second substances,
is used for representing the q-th weight, wherein q is used for representing the sequence number of the base station in all the base stations;
Preferably, the first predicted position is expressed by the following formula:
wherein the content of the first and second substances,
the first predicted position is used for representing the current time, wherein K is used for representing the current time;
m is used to represent the first matrix;
w is used to represent the weight matrix;
n is used to represent the second matrix.
Preferably, the step S4 includes:
preferably, the step S4 includes:
step S41, constructing a prediction state equation, and obtaining a second prediction position of the current time according to the real-time positioning result of the previous time and the prediction state equation;
step S42, determining the real-time positioning result of the current time according to the first predicted position and the second predicted position.
A multi-base station real-time positioning system based on Bluetooth signals is applied to a plurality of base stations, and all the base stations synchronously receive a signal; the multi-base-station real-time positioning system comprises:
the acquisition module is used for acquiring the signal receiving intensity, azimuth angle, pitch angle and spatial spectrum peak value of the signals received by all the base stations at the current moment;
the screening module is connected with the acquisition module and is used for screening the signal receiving intensity, the azimuth angle, the pitch angle and the spatial spectrum peak value;
the prediction module is connected with the screening module and used for carrying out combined positioning according to the screened signal receiving intensity, the azimuth angle, the pitch angle and the spatial spectrum peak value to obtain a first prediction position of the signal at the current moment;
and the filtering module is connected with the prediction module and used for determining a second prediction position of the current moment according to a real-time positioning result of the previous moment and determining a real-time positioning result of the current moment according to the first prediction position and the second prediction position.
The technical scheme has the following advantages or beneficial effects: the method and the system are applied to a plurality of base stations, azimuth angles and pitch angles which are not accurate enough and can cause misjudgment are screened out, and real-time positioning results of signals are obtained according to the screened signal receiving intensity, the screened azimuth angles, the screened pitch angles and spatial spectrum peak values.
Drawings
Fig. 1 is a schematic flowchart of a multi-base-station real-time positioning method based on bluetooth signals according to a preferred embodiment of the present invention;
fig. 2 is a schematic diagram of multiple base stations receiving bluetooth signals in a preferred embodiment of the present invention;
FIG. 3 is a schematic flow chart of step S1 according to the preferred embodiment of the present invention;
FIG. 4 is a schematic flow chart of step S2 according to the preferred embodiment of the present invention;
FIG. 5 is a schematic flow chart of step S3 according to the preferred embodiment of the present invention;
FIG. 6 is a schematic flow chart of step S4 according to the preferred embodiment of the present invention;
fig. 7 is a schematic structural diagram of a multi-base-station real-time positioning system based on bluetooth signals in a preferred embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, are within the scope of the present invention.
In addition, the embodiments of the present invention and the features of the embodiments may be combined with each other without conflict.
The invention is further described with reference to the following drawings and detailed description, but is not intended to be limited thereto.
As shown in fig. 1, a multi-base-station real-time positioning method based on bluetooth signals is applied to a plurality of base stations, and as shown in fig. 2, all the base stations synchronously receive a signal; the multi-base-station real-time positioning method comprises the following steps:
step S1, acquiring signal receiving intensity, azimuth angle, pitch angle and spatial spectrum peak value of all base station receiving signals at the current moment;
step S2, screening the signal receiving intensity, the azimuth angle, the pitch angle and the spatial spectrum peak value;
step S3, performing combined positioning according to the filtered signal receiving intensity, azimuth angle, pitch angle and spatial spectrum peak value to obtain a first predicted position of the signal at the current moment;
and step S4, determining a second predicted position of the current time according to the real-time positioning result of the previous time, and determining the real-time positioning result of the current time according to the first predicted position and the second predicted position.
Specifically, the method for positioning multiple base stations in real time based on the Bluetooth signal is provided in consideration of the fact that indoor real-time positioning cannot be accurately performed in the prior art, and is applied to multiple base stations, a first predicted position of a signal emission source is analyzed, and a real-time positioning result is obtained according to filtering. Moreover, considering that the accuracy of the azimuth angle and the pitch angle of the received signal of the base station obtained through the space spectrum search is reduced along with the increase of the distance between the signal emission source and the base station, if the real-time positioning result is obtained according to all the azimuth angles and the pitch angles, the precision of the real-time positioning result is reduced, and the occupancy rate of the operation resource is increased. Therefore, the method also obtains the signal receiving intensity and the spatial spectrum peak value, the signal receiving intensity is used as the basis for judging the accuracy of the azimuth angle, the pitch angle and the spatial spectrum peak value, the azimuth angle and the pitch angle which are not accurate enough and can cause misjudgment are screened out according to the numerical value of the signal receiving intensity, the resource utilization rate is improved, the spatial spectrum peak value has the corresponding relation with the azimuth angle and the pitch angle, in the process of determining the azimuth angle and the pitch angle, the spatial spectrum search is usually needed, the angle corresponding to the spectrum peak value of the spatial spectrum is used as the azimuth angle and the pitch angle, therefore, the spatial spectrum peak value is also used as the basis for judging the accuracy of different azimuth angles and pitch angles, the spatial spectrum peak value corresponding to the screened azimuth angle and pitch angle is obtained, the corresponding weight is set according to the spatial spectrum peak value, and the weight is weighted to, And the pitch angle further improves the positioning accuracy.
Furthermore, the positioning is carried out based on the Bluetooth signals, so that the deployment is convenient, the positioning accuracy is improved, the power consumption is low, the cost is low, the protocol is simple, and the integration into other systems is convenient.
Wherein, in order to accurately express the signal receiving strength, azimuth angle and pitch angle of the received signals of different base stations, the base stations are numbered as。
In a preferred embodiment of the present invention, as shown in fig. 3, step S1 includes:
step S11, acquiring the signal received by the base station at the current moment;
step S12, acquiring corresponding signal receiving strength according to the energy in the signal;
step S13, obtaining a corresponding spatial spectrum according to the orthogonality of the noise subspace in the signal and the array manifold vector, performing peak search on the spatial spectrum to obtain a spatial spectrum peak value and a corresponding angle, and using the corresponding angle as an azimuth angle and a pitch angle of the signal received by the base station.
Specifically, the antenna arrays on all the base stations receive a signal synchronously through the synchronization header information of the signal, and therefore, in step S11, the signals received by the antennas on the base stations can be specifically expressed as:
wherein the content of the first and second substances,
s (t) for representing signals incident to the base station;
n (t) is a noisy data vector where the antenna noise is white noise and its variance isNoise is statistically independent of signal.
Then, in step S12, the following formula can be used to obtain the signal reception strength of the corresponding base station according to the energy in the signal:
wherein the content of the first and second substances,
Finally, in step S13, the following formula can be used to obtain the corresponding spatial spectrum according to the orthogonality between the noise subspace of the signal and the array manifold vector:
wherein the content of the first and second substances,
for representing the spatial spectrum of the q-th base station, whereinFor the purpose of indicating the azimuth angle,for representing pitch angle;
for representing a feature vector of the noise subspace orthogonal to the array manifold vector at the qth base station.
Spatial spectrum to the q base stationAnd searching to obtain a space spectrum peak value and an angle corresponding to the space spectrum peak value, and taking the angle corresponding to the peak value as an azimuth angle and a pitch angle of the base station for receiving the signal.
Since the specific location information of each base station is different and the signal reception strength, the azimuth angle and the pitch angle of different base stations receiving signals synchronously are also different when the method is applied to a plurality of base stations, the signal reception strength, the azimuth angle, the pitch angle and the spatial spectrum peak corresponding to each base station are obtained by steps S11-S13 for each base station.
In a preferred embodiment of the present invention, as shown in fig. 4, step S2 includes:
step S21, obtaining the mean value and standard deviation of the signal receiving intensity at the current moment;
step S22, presetting a coefficient;
step S23, sequentially determining whether the difference between the received signal strength and the mean is greater than the product of the coefficient and the standard deviation:
if yes, the signal reception strength is screened out, and then the process goes to step S24;
if not, the signal receiving strength is retained, and then the step S24 is proceeded to;
step S24, determining whether all signal reception strengths have been screened:
if yes, go to step S25;
if not, returning to the step S23;
and step S25, screening the azimuth angle, the pitch angle and the spatial spectrum peak value corresponding to the signal receiving intensity according to the screened signal receiving intensity.
Specifically, inaccurate azimuth angles, elevation angles and spatial spectrum peaks can be screened out by using the signal receiving strength, and the screening process can be set correspondingly based on different screening rules.
Here, in step S21-step S25:
step S21, obtaining the signal receiving strength of N base stations, ,…, Mean value ofAnd standard deviation, wherein, mean valueThe following formula can be used to represent:
wherein the content of the first and second substances,
mean value for signal reception strength of all base stations, N forIndicates the number of all base stations,for indicating the signal reception strength of the q-th base station.
wherein the content of the first and second substances,
for indicating the standard deviation of the signal reception strengths of all base stations, N for indicating the number of all base stations,for indicating the signal reception strength of the q-th base station.
Wherein the content of the first and second substances,
n is used to indicate the number of all base stations.
Steps S23-S24, sequentially determining the received signal strength and the mean value of each signalWhether the difference is greater than the coefficientAnd standard deviation ofAnd according to the judgment result, the corresponding signal receiving strength is screened out or reserved, namely the signal receiving strength of the q-th base station meets the requirementIf not, all the signal receiving strength after screening is expressed as,,⋯ 。
Step S25, according to the filtered signal receiving strength, expressing as, ,⋯ Screening and signal receiving strength,,⋯Corresponding azimuth angle,,⋯And a pitch angle,,⋯。
Further, in step S2, different filtering rules may be preset, and different filtering is performed on the spatial spectrum peak, the azimuth angle, and the pitch angle according to the different filtering rules. Specifically, in the screening process, the grubbs criterion construction step S2 may be further applied to obtain the spatial spectrum peak value, the azimuth angle and the pitch angle after the screening.
In a preferred embodiment of the present invention, as shown in fig. 5, step S3 includes:
step S31, acquiring the position information of all base stations;
step S32, respectively acquiring a first tangent value and a second tangent value according to the position information of each base station and the azimuth angle and the pitch angle of the signals received by the base stations, and respectively constructing a first matrix and a second matrix according to the position information, the first tangent value and the second tangent value;
step S33, acquiring a corresponding weight matrix according to the spatial spectrum peak of each base station received signal;
in step S34, a first predicted position is obtained according to the first matrix, the second matrix and the weight matrix.
In a preferred embodiment of the present invention, the first tangent value and the second tangent value are expressed by the following formulas:
wherein the content of the first and second substances,
for representing the qth first tangent value, where q is used to represent the base station's ranking number among all base stations,for indicating the azimuth angle of the q base station receiving signal;
for representing a second tangent value, whereinA pitch angle for representing a signal received by the q-th base station;
x is used to represent the abscissa in a first predicted location of the signal;
y is used to represent the ordinate in a first predicted position of the signal;
z is used to represent the vertical coordinate in a first predicted location of the signal;
Considering the position information, the acquired signal reception strength, the azimuth angle, and the pitch angle for only one base station, a numerical relationship with the first predicted position may be constructed:
and (3) sorting to obtain a first predicted position:
thus, from the position information of the plurality of base stations, the acquired signal reception strength, the azimuth angle, and the pitch angle, a numerical relationship with the first predicted position can be constructed:
wherein the content of the first and second substances,
h is used to represent a first predicted position;
m is used to represent a first matrix;
w is used to represent a weight matrix;
n is used to represent the second matrix.
Specifically, the first matrix is represented by the following formula:
wherein the content of the first and second substances,
for representing the qth first tangent value, where q is used to represent the base station's ranking number among all base stations,for indicating the azimuth angle of the q base station receiving signal;
for representing a second tangent value, whereinA pitch angle for representing a signal received by the q-th base station;
l is used to indicate the number after screening.
Specifically, the second matrix is represented by the following formula:
wherein the content of the first and second substances,
for representing the qth first tangent value, where q is used to represent the base station's ranking number among all base stations,for indicating the azimuth angle of the q base station receiving signal;
for representing a second tangent value, whereinA pitch angle for representing a signal received by the q-th base station;
x is used to represent the abscissa in a first predicted location of the signal;
y is used to represent the ordinate in a first predicted position of the signal;
z is used to represent the vertical coordinate in a first predicted location of the signal;
l is used to indicate the number after screening.
Specifically, the weight matrix is represented as:
wherein the content of the first and second substances,
is used for representing the q weight, wherein q is used for representing the sequence number of the base station in all the base stations;
In a preferred embodiment of the invention, the weights are expressed as:
wherein the content of the first and second substances,
is used for representing the q weight, wherein q is used for representing the sequence number of the base station in all the base stations;
l is used to indicate the number after screening.
In this regard, considering that the accuracy of the azimuth angle and the pitch angle of the received signal of the base station obtained through the spatial spectrum search decreases with the increase of the distance between the signal emission source and the base station, the weight of the base station is determined according to the spatial spectrum peak value of the received signal of each base station and a weight matrix is constructed, so that the ratio of the spatial spectrum peak value, the azimuth angle and the pitch angle corresponding to different base stations is different in the process of obtaining the first predicted position, and the influence of the insufficiently accurate azimuth angle and pitch angle on the first predicted position is avoided.
In a preferred embodiment of the present invention, as shown in fig. 6, step S4 includes:
step S41, constructing a prediction state equation, and obtaining a second prediction position at the current moment according to the real-time positioning result at the previous moment and the prediction state equation;
and step S42, determining the real-time positioning result of the current time according to the first predicted position and the second predicted position.
Specifically, the first predicted position at the current time can be obtained according to the steps S1-S3, and in the real-time obtaining process, the position of the signal emission source changes, and the change speed and the change path are different, and the first predicted positions obtained by the plurality of base stations are likely to be interfered to generate errors, so that step S4 is set, and the real-time positioning result at the current time is obtained according to the real-time positioning result at the previous time and the first predicted position obtained in step S3 through comprehensive analysis.
Step S41, a prediction state equation is constructed, and a second prediction position at the current time is obtained according to the real-time positioning result at the previous time and the prediction state equation:
wherein the content of the first and second substances,
for indicating a second predicted position of the current time, whichInFor representing the current time;
for representingA transformation matrix under a state, namely a basis for predicting a variable with a changed position;
Step S42, determining a real-time positioning result of the current time according to the first predicted position and the second predicted position, which may include:
in step S421, the covariance of the current time is obtained according to the covariance of the previous time, and can be represented by the following formula:
wherein the content of the first and second substances,
a covariance matrix for representing a current time, wherein K is used to represent the current time;
for representing the covariance of the previous time instant, whereinFor representing the current time;
Step S422, acquiring kalman gain according to the covariance of the current time, which may be expressed by the following formula:
wherein the content of the first and second substances,
for representing observation matrices, i.e. different types of position variations of the signal emission sources, such as linear variations, or non-linear variations, along a straight line;
the covariance matrix used for representing the observation noise can be regarded as the observation mean value and can be configured according to the actual situation.
Step S423, obtaining a real-time positioning result of the current time according to the first predicted position, the second predicted position and the kalman gain, and using the following formula to represent:
wherein the content of the first and second substances,the real-time positioning result is used for representing the current moment;
In addition, a process of updating the covariance equation may be included, and the covariance at the current time is obtained according to step S42, which may be expressed by the following formula:
accordingly, a second predicted position at the next time can be obtained according to the covariance at the current time and the real-time positioning result.
A multi-base station real-time positioning system, as shown in fig. 7, is applied to a plurality of base stations 1,2 ⋯ q ⋯, n, all of which receive a signal synchronously; the multi-base station real-time positioning system comprises:
an obtaining module a1, configured to obtain signal receiving intensities, azimuth angles, pitch angles, and spatial spectrum peak values of all base station received signals at the current time;
the screening module A2 is connected with the acquisition module A1 and is used for screening the signal receiving intensity, the azimuth angle, the pitch angle and the spatial spectrum peak value;
the prediction module A3 is connected with the screening module A2 and is used for performing combined positioning according to the screened signal receiving intensity, azimuth angle, pitch angle and spatial spectrum peak value to obtain a first prediction position of the signal at the current moment;
and the filtering module A4 is connected to the predicting module A3 and is used for determining a second predicted position at the current moment according to the real-time positioning result at the previous moment and determining the real-time positioning result at the current moment according to the first predicted position and the second predicted position.
Specifically, in view of the fact that the prior art cannot accurately perform indoor real-time positioning, a multi-base-station real-time positioning system is provided, which is applied to a plurality of base stations, analyzes a first predicted position of a signal emission source, and obtains a real-time positioning result according to filtering, specifically, the real-time positioning method is clear from the above description, and the real-time positioning result is determined and then output through an output module a 5.
The first embodiment is as follows:
when a signal transmitting end is arranged at a fixed position of a space coordinate (20, 30, 25), three base stations are adopted for joint positioning, wherein the distance between one base station and the signal transmitting end is far compared with the distance between the other two base stations, and when the pitch angle and the azimuth angle obtained by the base station are not accurate enough, the final positioning result is as follows:
means for | Positioning result |
Comparative example 1: considering only pitch and azimuth | (21.6,26.8,24.7) |
Comparative example 2: considering only the signal reception strength | (16.2,39.1,30.4) |
The invention | (19.8,29.6,25.0) |
As can be seen from the table, through steps S1-S4 of the present invention, joint positioning is performed according to the filtered signal received strength, azimuth angle, pitch angle and spatial spectrum peak value, and the numerical value of the signal received strength is used as a basis for determining the accuracy of the azimuth angle and pitch angle, so as to screen out the azimuth angle and pitch angle which are not accurate enough and may cause erroneous judgment, and the finally determined real-time positioning result is (19.8, 296, 25.0). In comparative example 1, only the pitch angle and the azimuth angle are considered, or in comparative example 2, only the signal receiving strength is considered, and the positioning result obtained by adopting the method is closest to the correct position. Furthermore, even when the positioning results of comparative examples 1 and 2 are combined, the maximum value, the minimum value or the average value cannot be obtained, so that more accurate positioning results than those of the present invention cannot be obtained.
Example two:
when a signal transmitting end is arranged at a fixed position of a space coordinate (5, 15, 20), three base stations are adopted for joint positioning, wherein the distance between one base station and the signal transmitting end is far compared with the distance between the other two base stations, and when the pitch angle and the azimuth angle obtained by the base station are not accurate enough, the final positioning result is as follows:
means for | Positioning result |
All base stations adopt the same weight | (5.8,13.1,15.9) |
Setting different weights by signal strength | (4.7,14.2,16.4) |
The invention | (4.9,14.3,18.9) |
It can be seen from the table that the present invention determines the weight of the base station according to the spatial spectrum peak of the received signal of each base station, so that the numerical ratios of the received signal strength, the azimuth angle and the pitch angle of different base stations are different, and the finally determined real-time positioning result is (4.9, 14.3, 18.9). The method comprises the steps that the numerical value occupation ratios of the signal receiving strength, the azimuth angle and the pitch angle corresponding to different base stations are the same, the finally determined real-time positioning result is (5.8, 13.1, 15.9), and when the numerical value occupation ratios of the signal receiving strength, the azimuth angle and the pitch angle corresponding to different base stations are determined by the weight values corresponding to the signal receiving strength, the finally determined real-time positioning result is (4.7, 14.2, 16.4).
The technical scheme has the following advantages or beneficial effects: the method and the system are applied to a plurality of base stations, azimuth angles and pitch angles which are not accurate enough and possibly cause misjudgment are screened out, and real-time positioning results of signals are obtained according to the screened signal receiving intensity, the screened azimuth angles and the screened pitch angles.
While the invention has been described with reference to a preferred embodiment, it will be understood by those skilled in the art that various changes in form and detail may be made therein without departing from the spirit and scope of the invention.
Claims (12)
1. A multi-base station real-time positioning method based on Bluetooth signals is applied to a plurality of base stations and is characterized in that all the base stations synchronously receive signals; the multi-base-station real-time positioning method comprises the following steps:
step S1, acquiring the signal receiving intensity, azimuth angle, pitch angle and spatial spectrum peak value of the signal received by all the base stations at the current moment;
step S2, screening the signal reception intensity, the azimuth angle, the pitch angle, and the spatial spectrum peak;
step S3, performing combined positioning according to the filtered signal receiving intensity, the azimuth angle, the pitch angle and the spatial spectrum peak value to obtain a first predicted position of the signal at the current moment;
step S4, determining a second predicted position of the current time according to the real-time positioning result of the previous time, and determining the real-time positioning result of the current time according to the first predicted position and the second predicted position.
2. The multi-base-station real-time positioning method according to claim 1, wherein said step S1 includes:
step S11, acquiring the signal received by the base station at the current time;
step S12, acquiring the corresponding signal receiving strength according to the energy in the signal;
step S13, obtaining a corresponding spatial spectrum according to orthogonality between the noise subspace in the signal and the array manifold vector, performing peak search on the spatial spectrum to obtain a spatial spectrum peak value and a corresponding angle, and using the corresponding angle as the azimuth angle and the pitch angle at which the base station receives the signal.
3. The multi-base-station real-time positioning method according to claim 1, wherein said step S2 includes:
step S21, obtaining a mean value and a standard deviation of the signal reception intensity at the current time;
step S22, presetting a coefficient;
step S23, sequentially determining whether the difference between the received signal strength and the mean is greater than the product of the coefficient and the standard deviation:
if yes, the signal reception strength is screened out, and then the step S24 is carried out;
if not, the signal receiving strength is retained, and then the step S24 is carried out;
step S24, determining whether all the signal reception intensities have been filtered:
if yes, go to step S25;
if not, returning to the step S23;
step S25, screening the azimuth angle, the pitch angle, and the spatial spectrum peak corresponding to the signal reception intensity according to the screened signal reception intensity.
4. The multi-base-station real-time positioning method according to claim 1, wherein said step S3 includes:
step S31, acquiring the position information of all the base stations;
step S32, sequentially obtaining a first tangent value and a second tangent value according to the position information of each base station and the azimuth angle and the pitch angle at which the base station receives the signal, and respectively constructing a first matrix and a second matrix according to the position information, the first tangent value and the second tangent value;
step S33, obtaining a corresponding weight matrix according to the spatial spectrum peak value of the signal received by each base station;
step S34, obtaining the first predicted position according to the first matrix, the second matrix and the weight matrix.
5. A multi-base-station real-time positioning method as claimed in claim 4, wherein said first tangent value and said second tangent value are expressed by the following formulas:
wherein the content of the first and second substances,
for representing the qth of said first tangent, where q is used to represent the ranking number of said base station among all base stations,for indicating the azimuth angle at which the q-th base station receives the signal;
for representing the second tangent value, whereinSaid pitch angle representing the reception of said signal by the qth of said base station;
x is used to represent the abscissa in a first predicted position of the signal;
y is used to represent the ordinate in a first predicted position of the signal;
z is used to represent the vertical coordinate in a first predicted location of the signal;
6. The multi-base-station real-time positioning method according to claim 4, wherein the first matrix is expressed by the following formula:
wherein the content of the first and second substances,
for representing the qth of said first tangent, where q is used to represent the ranking number of said base station among all base stations,for indicating the azimuth angle at which the q-th base station receives the signal;
for representing the second tangent value, whereinSaid pitch angle representing the reception of said signal by the qth of said base station;
l is used to indicate the number after screening.
7. The multi-base-station real-time positioning method according to claim 4, wherein the second matrix is expressed by the following formula:
wherein the content of the first and second substances,
for representing the qth of said first tangent, where q is used to represent the ranking number of said base station among all base stations,for indicating the azimuth angle at which the q-th base station receives the signal;
for representing the second tangent value, whereinSaid pitch angle representing the reception of said signal by the qth of said base station;
x is used to represent the abscissa in a first predicted position of the signal;
y is used to represent the ordinate in a first predicted position of the signal;
z is used to represent the vertical coordinate in a first predicted location of the signal;
l is used to indicate the number after screening.
8. The multi-base-station real-time positioning method according to claim 4, wherein the weight is expressed as:
wherein the content of the first and second substances,
is used for representing the q-th weight, wherein q is used for representing the sequence number of the base station in all the base stations;
l is used to indicate the number after screening.
9. The multi-base-station real-time positioning method according to claim 4, wherein the weight matrix is expressed as:
wherein the content of the first and second substances,
is used for representing the q-th weight, wherein q is used for representing the sequence number of the base station in all the base stations;
10. The multi-base-station real-time positioning method according to claim 4, wherein the first predicted position is expressed by the following formula:
wherein the content of the first and second substances,
the first predicted position is used for representing the current time, wherein K is used for representing the current time;
m is used to represent the first matrix;
w is used to represent the weight matrix;
n is used to represent the second matrix.
11. The multi-base-station real-time positioning method according to claim 1, wherein said step S4 includes:
step S41, constructing a prediction state equation, and obtaining a second prediction position of the current time according to the real-time positioning result of the previous time and the prediction state equation;
step S42, determining the real-time positioning result of the current time according to the first predicted position and the second predicted position.
12. A multi-base station real-time positioning system based on Bluetooth signals is applied to a plurality of base stations and is characterized in that all the base stations synchronously receive a signal; the multi-base-station real-time positioning system comprises:
the acquisition module is used for acquiring the signal receiving intensity, azimuth angle, pitch angle and spatial spectrum peak value of the signals received by all the base stations at the current moment;
the screening module is connected with the acquisition module and is used for screening the signal receiving intensity, the azimuth angle, the pitch angle and the spatial spectrum peak value;
the prediction module is connected with the screening module and used for carrying out combined positioning according to the screened signal receiving intensity, the azimuth angle, the pitch angle and the spatial spectrum peak value to obtain a first prediction position of the signal at the current moment;
and the filtering module is connected with the prediction module and used for determining a second prediction position of the current moment according to a real-time positioning result of the previous moment and determining a real-time positioning result of the current moment according to the first prediction position and the second prediction position.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110078642.3A CN112394321B (en) | 2021-01-21 | 2021-01-21 | Multi-base-station real-time positioning method and system based on Bluetooth signals |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110078642.3A CN112394321B (en) | 2021-01-21 | 2021-01-21 | Multi-base-station real-time positioning method and system based on Bluetooth signals |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112394321A true CN112394321A (en) | 2021-02-23 |
CN112394321B CN112394321B (en) | 2021-04-13 |
Family
ID=74625121
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110078642.3A Active CN112394321B (en) | 2021-01-21 | 2021-01-21 | Multi-base-station real-time positioning method and system based on Bluetooth signals |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112394321B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115128540A (en) * | 2022-09-02 | 2022-09-30 | 珠海翔翼航空技术有限公司 | Positioning and orienting method and device suitable for flight platform |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105933975A (en) * | 2016-04-11 | 2016-09-07 | 南京邮电大学 | WiFi fingerprint-based accuracy improved indoor positioning method |
US20180100915A1 (en) * | 2016-10-07 | 2018-04-12 | Marko Beko | Apparatus and method for RSS/AoA target 3-D localization in wireless networks |
CN110837076A (en) * | 2019-11-18 | 2020-02-25 | 哈尔滨工程大学 | Tensor decomposition-based vector hydrophone array orientation estimation method |
CN111212476A (en) * | 2020-04-21 | 2020-05-29 | 杭州优智联科技有限公司 | Multi-base-station ultra-wideband positioning method based on frequency modulation continuous waves |
-
2021
- 2021-01-21 CN CN202110078642.3A patent/CN112394321B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105933975A (en) * | 2016-04-11 | 2016-09-07 | 南京邮电大学 | WiFi fingerprint-based accuracy improved indoor positioning method |
US20180100915A1 (en) * | 2016-10-07 | 2018-04-12 | Marko Beko | Apparatus and method for RSS/AoA target 3-D localization in wireless networks |
CN110837076A (en) * | 2019-11-18 | 2020-02-25 | 哈尔滨工程大学 | Tensor decomposition-based vector hydrophone array orientation estimation method |
CN111212476A (en) * | 2020-04-21 | 2020-05-29 | 杭州优智联科技有限公司 | Multi-base-station ultra-wideband positioning method based on frequency modulation continuous waves |
Non-Patent Citations (2)
Title |
---|
M.W.KHAN: ""Enhanced Hybrid Positioning in Wireless Network I:AoA-ToA"", 《2014 INTERNATIONAL CONFERENCE ON TECOMMUNIACTIONS AND MULTIMEDIA(TEMU)》 * |
李雪松: ""基于定位的超宽带传感器网络路由算法研究"", 《压电与声光》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115128540A (en) * | 2022-09-02 | 2022-09-30 | 珠海翔翼航空技术有限公司 | Positioning and orienting method and device suitable for flight platform |
CN115128540B (en) * | 2022-09-02 | 2022-12-06 | 珠海翔翼航空技术有限公司 | Positioning and orienting method and device suitable for flight platform |
Also Published As
Publication number | Publication date |
---|---|
CN112394321B (en) | 2021-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6748324B2 (en) | Method for determining location information | |
KR102116824B1 (en) | Positioning system based on deep learnin and construction method thereof | |
CN107613559B (en) | A kind of DOA fingerprint base localization method based on 5G signal | |
EP2672288B1 (en) | Locating apparatus, locating system, locating method and program | |
CN101118280B (en) | Distributed wireless sensor network node self positioning method | |
CN102204373B (en) | Apparatus and method for estimating an orientation of a mobile terminal | |
CN101631349B (en) | Method, device and wireless operation maintenance center for positioning terminal | |
JP4592506B2 (en) | Uplink interference source locating apparatus and method | |
KR20030055412A (en) | Apparatus and method for locating a mobile tranceiver | |
CN112346007B (en) | Direction finding positioning method and system | |
CN1729406A (en) | Methods and systems for determining the position of a mobile terminal using digital television signals | |
CN112394321B (en) | Multi-base-station real-time positioning method and system based on Bluetooth signals | |
JP2000102058A (en) | Method and system for detecting position of mobile terminal | |
CN113747360B (en) | Indoor positioning method based on mixed visible light and BLE | |
US20170336512A1 (en) | Method to increase positioning accuracy of global navigation satellite systems by integration of correction service with a mobile communication network | |
CN106501768B (en) | Localization method in a kind of Dan Zhandong | |
CN112616184B (en) | Mobile equipment position estimation method based on multi-base station channel state information fusion | |
CN103781095A (en) | Calibration method for TDOA (Time Difference Of Arrival) measurement errors, transmit points and system | |
CN106872939A (en) | A kind of alignment system and localization method | |
CN108226912B (en) | Sparse network-based non-contact object perception positioning method and system | |
WO2019138376A1 (en) | Method of using a multi-input and multi-output antenna (mimo) array | |
CN113068127B (en) | Positioning method of NB-IoT terminal | |
CN107113690A (en) | A kind of method and device for obtaining positional information | |
CN109459016A (en) | A kind of micro-nano satellite cluster relative positioning method based on location fingerprint | |
CN107948916B (en) | Multi-beam positioning method, device and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20210223 Assignee: SUZHOU PANCHIP MICROELECTRONICS Co.,Ltd. Assignor: SHANGHAI PANCHIP MICROELECTRONICS Co.,Ltd. Contract record no.: X2022990000576 Denomination of invention: A method and system for real-time positioning of multiple base stations based on bluetooth signal Granted publication date: 20210413 License type: Common License Record date: 20220829 |