CN112880699B - Vehicle cooperative positioning method based on brain selective attention mechanism - Google Patents

Vehicle cooperative positioning method based on brain selective attention mechanism Download PDF

Info

Publication number
CN112880699B
CN112880699B CN202110071311.7A CN202110071311A CN112880699B CN 112880699 B CN112880699 B CN 112880699B CN 202110071311 A CN202110071311 A CN 202110071311A CN 112880699 B CN112880699 B CN 112880699B
Authority
CN
China
Prior art keywords
vtbp
formula
vehicle
variance
position coordinate
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.)
Active
Application number
CN202110071311.7A
Other languages
Chinese (zh)
Other versions
CN112880699A (en
Inventor
来磊
邹鲲
杨宾锋
李海林
李保中
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Air Force Engineering University of PLA
Original Assignee
Air Force Engineering University of PLA
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 Air Force Engineering University of PLA filed Critical Air Force Engineering University of PLA
Priority to CN202110071311.7A priority Critical patent/CN112880699B/en
Publication of CN112880699A publication Critical patent/CN112880699A/en
Application granted granted Critical
Publication of CN112880699B publication Critical patent/CN112880699B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/26Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
    • G01C21/34Route searching; Route guidance
    • G01C21/3446Details of route searching algorithms, e.g. Dijkstra, A*, arc-flags, using precalculated routes
    • 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/26Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
    • G01C21/34Route searching; Route guidance
    • G01C21/3407Route searching; Route guidance specially adapted for specific applications
    • G01C21/3415Dynamic re-routing, e.g. recalculating the route when the user deviates from calculated route or after detecting real-time traffic data or accidents

Abstract

The invention discloses a vehicle cooperative localization method based on a brain selective attention mechanism, which specifically aims at the influence of a geometric structure, RB position precision and a relative motion state formed between a cooperative reference node and a vehicle to be localized on cooperative localization; adopting a selective attention mechanism of human brain to information processing, taking the three influencing factors as characteristic points of information selection, and carrying out selective filtering processing on information; and then, comprehensively screening out the optimal RB through feature integration so as to further improve the reliability and precision of the cooperative positioning of the vehicle.

Description

Vehicle cooperative positioning method based on brain selective attention mechanism
Technical Field
The invention relates to the technical field of vehicle cooperative positioning, in particular to a vehicle cooperative positioning method based on a brain selective attention mechanism.
Background
Real-time accurate positioning of a moving vehicle is one of essential key technologies for implementation of an Intelligent Transportation System (ITS) and an Automatic Driving Technology (ADT), such as: vehicle information management, traffic control and vehicle information service in the ITS, automatic obstacle avoidance, lane change and other technologies in the ADT all need the self-positioning capability of the vehicle with high reliability and high precision; therefore, the research on a reliable vehicle positioning method has important significance for the realization of an intelligent traffic system;
at present, the mainly applied ground carrier Positioning technology comprises inertial navigation Positioning, visual navigation Positioning and Global Positioning System (GPS) technology, and the inertial navigation Positioning technology is usually found in special ground carriers for military use and the like due to high cost, so that the application of the technology in civil vehicles is restricted by high cost; although the visual navigation positioning technology is widely researched and applied, it can only provide relative position information; the GPS is the most widely used vehicle positioning technology at present due to higher cost performance and maturity; the GPS has positioning accuracy which is different from meter level to centimeter level according to different implementation technologies, for example, the positioning accuracy of the carrier phase differential GPS can reach centimeter level; however, the experiment proves that: the accuracy and reliability of GPS positioning are severely affected by satellite signal blocking and multipath interference phenomena in a densely populated urban or tunnel environment, even with carrier phase differential GPS, and thus are difficult to adapt to the needs of urban ITS.
With the high-speed development of wireless communication technology, a vehicular ad hoc Network (VANET) is responsible for the interaction task of various traffic information and becomes an important component of ITS; meanwhile, a V2V, V2I and V2X communication mode in the VANET also brings a new solution to vehicle navigation positioning and positioning enhancement, namely a vehicle cooperative positioning technology, which is characterized in that a vehicle takes other vehicles or traffic infrastructure as a reference point and obtains information such as the position and relative position of the reference point through a wireless network to position the vehicle; in recent years, the strong demand of the application of artificial intelligence technology prompts experts in a plurality of fields such as neuroscience, cognitive science, computational science and the like to try to disclose the mechanism of brain information perception and processing from different levels and try to simulate the functions of the brain information perception and processing; the researchers found that: the reason why the brain can rapidly process and respond to massive information of the surrounding environment is that the brain can consciously or unconsciously screen the information in a targeted manner, namely, a selective attention mechanism is provided, but in the prior art, the brain selective attention mechanism is not applied to the aspect of vehicle cooperative positioning, and corresponding well-known reports and researches are not provided.
Disclosure of Invention
Aiming at the existing problems, the invention aims to provide a vehicle cooperative localization method based on a brain selective attention mechanism, and the method is characterized in that the influence of a geometric structure, RB position precision and relative motion state formed between a cooperative reference node and a vehicle to be localized on cooperative localization is researched, a selective attention mechanism of human brain on information processing is provided, and the three influencing factors are used as characteristic points of information selection to selectively filter information; and then, the optimal RB is comprehensively screened out through feature integration so as to further improve the reliability and precision of vehicle cooperative positioning, and the method has the characteristics of high positioning precision and no environmental limitation.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows:
a vehicle cooperative positioning method based on brain selective attention mechanism comprises
The method comprises the following steps: in a VANET environment, selecting N RB combinations from VTBP neighbor vehicles or ITS infrastructure nodes, calculating geometric accuracy characteristic values of the selected N RB combinations, and sorting the advantages and the disadvantages according to an ascending order;
step two: introducing a CSM model of a variance adjustment factor as a state equation, taking the position coordinate of the RB as an observed quantity, correcting the position coordinate of the RB, and simultaneously evaluating the precision of the position coordinate of the RB;
step three: comprehensively evaluating the RB relative position characteristic, the position precision characteristic and the relative motion characteristic; establishing a cooperative positioning equation set by using the evaluated optimal RB position coordinates and the relative distance between the optimal RB position coordinates and the VTBP, and solving the VTBP position coordinates;
step four: and establishing a dead reckoning model by taking the solved VTBP position coordinate as an observed quantity, and correcting the VTBP position coordinate, so as to obtain a final VTBP position coordinate value and complete vehicle cooperative positioning.
Preferably, the selection of N RB combinations and the calculation of the geometric precision characteristic value in step one include:
s1.1, setting the position coordinates of VTBP as unknown quantity (x, y), and selecting n neighbor vehicles or ITS infrastructure nodes from the surrounding environment as positioning RB i Wherein i =1,2, \8230;, n; RB (radio B) i Is a known quantity (x) i ,y i );
S1.2 VTBP and RB can be obtained i Relative distance d therebetween i Comprises the following steps:
Figure BDA0002905890320000031
s1.3, calculating the geometric accuracy characteristic values of the selected N RB combinations, and sorting the advantages and the disadvantages according to ascending order.
Preferably, the calculation process of the geometric accuracy characteristic value in step S1.3 includes:
(1) Introducing measurement error e i Then there is an error of (e) in the position of the corresponding VTBP x ,e y ) Then, the formula (1) after adding the error is:
Figure BDA0002905890320000041
converting the formula (2) into a linear equation to obtain
Figure BDA0002905890320000042
(2) Order to
Figure BDA0002905890320000043
a i =(x-x i )/d i ,b i =(y-y i )/d i Equation (3) can be converted to a matrix form:
L=HX+e (4)
in formula (4): l = [ L 1 ,l 2 ,…l n ] T ,X=[e x ,e y ] T ,e=[e 1 ,e 2 ,…e n ] T
Figure BDA0002905890320000044
(3) When H is full rank, H T H is reversible, then the error between the position estimate and the true value is:
Figure BDA0002905890320000045
the magnitude of the error in equation (5) is measured as covariance and yields:
Figure BDA0002905890320000046
in formula (6): sigma 2 The variance of the noise in e when the noise is uncorrelated two by two;
(4) As can be seen from formula (6), (H) T H) -1 Expressed as a magnification factor for the distance measurement error, and therefore defining the geometric position accuracy factor of the cooperative positioning of the vehicle as G, then
Figure BDA0002905890320000051
Wherein, in the formula (7), tr [ ] represents the inverse operation of the matrix;
(5) The operation of the formula (7) is simplified, the matrix inversion is equivalent to the sum of the eigenvalues of the matrix, and the corresponding formula (7) can be written as
Figure BDA0002905890320000052
In formula (8): lambda 1 、λ 2 Is a matrix H T A characteristic value of H;
and finally, obtaining a geometric position precision factor of vehicle cooperative positioning:
Figure BDA0002905890320000053
in the formula (9), det [ ] represents a determinant for matrix calculation.
Preferably, the process of correcting the position coordinates of the RB in real time in step two includes:
s2.1 adopts CSM as an RB carrier motion state equation, and is expressed as follows:
Figure BDA0002905890320000054
in the formula (10), X k+1 Is a state vector, phi k Being a state transition matrix, U k In order to input the control matrix, the control matrix is,
Figure BDA0002905890320000055
mean value of acceleration at time k, W k Mean value of zero and variance of Q k Gaussian distributed noise vector of (2); the above-mentioned
Figure BDA0002905890320000056
Figure BDA0002905890320000057
Figure BDA0002905890320000058
Figure BDA0002905890320000059
In formula (11), x k
Figure BDA0002905890320000061
Position, velocity and acceleration of RB, respectively; in equations (12) and (13), T is a sampling period; in the formula (14), τ is the maneuvering frequency,
Figure BDA0002905890320000062
the variance of the maneuvering acceleration of the carrier is shown, and q is a noise matrix;
s2.2 selection
Figure BDA0002905890320000063
For the mean acceleration at time k, the variance key is introducedSection factor eta k =μ(r k ) And then obtaining:
Figure BDA0002905890320000064
variance of acceleration
Figure BDA0002905890320000065
The update formula of (c) is:
Figure BDA0002905890320000066
in formula (25): a is max And a min Respectively representing the maximum value and the minimum value of the acceleration;
the equation for performing filter correction on the position coordinates of the RB carrier is:
Figure BDA0002905890320000067
preferably, the variance adjustment factor η is described in step S2.2 k The introduction process comprises the following steps:
(1) Let the innovation vector of the filter in KF be defined as:
Figure BDA0002905890320000068
innovation vector d k Is not relevant in the ideal case, and d k Is that the mean is zero and the variance is S k When the RB carrier motion is maneuvered, the maneuver changes the orthogonality of the information so that d k Is changed and is no longer zero, i.e.
Figure BDA0002905890320000069
(2) Normalizing the sequence of the innovation vectors to obtain statistics
Figure BDA0002905890320000071
(3) Establishing window detection statistic, and setting window size as m, window statistic wg of k moment k Is defined as
Figure BDA0002905890320000072
Coefficient of determination r k Is defined as
r k =wg k /wg k-1 (22)
When the RB carrier is not mobile, r k The value is usually close to 1, r when a maneuver occurs k The value increases rapidly and is proportional to the degree of mobility;
(4) In order to convert the mobility of the RB carrier into a variance adjustment factor, a raised-seminormal distribution function mu (u) is introduced,
Figure BDA0002905890320000073
(5) Will r is k As the input variable of the ascending seminormal distribution function mu (u), the variance adjustment factor eta can be obtained k
η k =μ(r k ) (24)。
Preferably, the specific process of step three includes:
s3.1, comprehensively evaluating the relative position characteristic, the position precision characteristic and the relative motion characteristic of the RB combination by adopting a fuzzy evaluation method, and selecting an optimal RB cooperative positioning combination;
and S3.2, establishing a cooperative positioning equation set by utilizing the position coordinate of the optimal RB and the relative distance between the optimal RB and the VTBP, solving the position coordinate of the VTBP by adopting a Taylor series analytic method, and establishing a cooperative positioning equation.
Preferably, the specific process of step four includes:
s4.1, establishing an observation equation by taking the solved VTBP position coordinate as an observed quantity, and establishing a dead reckoning model by the speed and the course angle measured by the VTBP;
s4.2, establishing a measurement model by taking the relative distance and the relative angle between the two vehicles as observed quantities and combining the position coordinates of the reference vehicle transmitted by the reference vehicle;
and S4.3, correcting the VTBP position coordinate by adopting an extended Kalman to obtain a final VTBP position coordinate value, and finishing the vehicle cooperative positioning.
The invention has the beneficial effects that: the invention discloses a vehicle cooperative positioning method based on brain selective attention mechanism, compared with the prior art, the improvement of the invention is as follows:
aiming at the problems in the prior art, the invention provides a Vehicle cooperative positioning method based on a brain selective attention mechanism, which particularly aims at the influence of a geometric structure, RB position precision and a relative motion state formed between a cooperative Reference node (RB) and a Vehicle To Be Positioned (VTBP) on cooperative positioning; adopting a selective attention mechanism of human brain to information processing, taking the three influencing factors as characteristic points of information selection, and carrying out selective filtering processing on the information; and then, the optimal RB is comprehensively screened out through feature integration so as to further improve the reliability and precision of vehicle cooperative positioning, and the method has the advantages of high positioning precision and no environmental limitation.
Drawings
Fig. 1 is a flowchart of a vehicle cooperative localization method based on brain selective attention mechanism according to the present invention.
Fig. 2 is a simplified diagram of a calculation process of the vehicle cooperative localization method based on the brain selective attention mechanism according to the present invention.
Fig. 3 is a schematic diagram of vehicle cooperative positioning in a VANET environment of the present invention.
FIG. 4 is a diagram of an attention model for optimal RB combination selectivity according to the present invention.
Fig. 5 is a schematic diagram of a motion trajectory of a vehicle in experimental verification of embodiment 1 of the present invention.
FIG. 6 is a graph comparing the results of the experiment of example 1 of the method of the present invention.
Detailed Description
In order to make those skilled in the art better understand the technical solution of the present invention, the following description will be made with reference to the accompanying drawings and embodiments.
Referring to fig. 1-4, a method for cooperative localization of a vehicle based on brain selective attention mechanism includes
The method comprises the following steps: setting the position coordinates of a VTBP (vehicle to be positioned) as unknown quantities (x, y) in a VANET environment, selecting N RB (reference vehicle) combinations from neighbor vehicles or ITS infrastructure nodes according to a geometric barycenter method, wherein the selection rule is that the position of the VTBP at the last moment is positioned near the geometric barycenter of a polygon formed by the RB combinations, calculating geometric precision characteristic values of the selected N RB combinations, and performing quality sequencing according to an ascending order, wherein the smaller the geometric precision characteristic value is, the better the RB combination state is; the specific process comprises the following steps:
s1.1 real-time information interaction between vehicles and ITS infrastructure in VANET environment (as shown in FIG. 3), setting the position coordinates of VTBP in VANET environment as unknown quantity (x, y), and selecting n neighbor vehicles or ITS infrastructure nodes from the surrounding environment as positioning RB according to geometric gravity center method i Wherein i =1,2, \8230, N; RB (radio B) i Is a known quantity (x) i ,y i );
S1.2 VTBP and RB can be obtained i Relative distance d between i Comprises the following steps:
Figure BDA0002905890320000091
s1.3, calculating the geometric accuracy characteristic values of the selected N RB combinations, wherein the calculation process is as follows:
(1) The actual measured VTBP and the reference point RB due to the measurement error of the ranging device i There will be an error e in the relative distance between i The position estimate of the corresponding VTBP has an error of (e) x ,e y ) Then, the formula (1) after adding the error is:
Figure BDA0002905890320000101
converting the formula (2) into a linear equation to obtain
Figure BDA0002905890320000102
(2) Order to
Figure BDA0002905890320000103
a i =(x-x i )/d i ,b i =(y-y i )/d i
Equation (3) can be converted to a matrix form:
L=HX+e (4)
in formula (4): l = [ L 1 ,l 2 ,…l n ] T ,X=[e x ,e y ] T ,e=[e 1 ,e 2 ,…e n ] T
Figure BDA0002905890320000104
(3) When H is full rank, H T H is reversible, then the error between the position estimate and the true value is:
Figure BDA0002905890320000105
the magnitude of the error in equation (5) can be measured as covariance, resulting in:
Figure BDA0002905890320000106
in formula (6): sigma 2 The variance of the noise in e when the noise is uncorrelated two by two;
(4) From the formula (6), (H) can be seen T H) -1 Expressed as a magnification of the distance measurement error, thus, the vehicle is drivenThe geometric position precision factor of cooperative positioning is defined as G, then
Figure BDA0002905890320000111
In the formula (7), tr [ ] represents the inverse operation of the matrix;
when the vehicle to be positioned is cooperatively positioned, the combination of the minimum G value is selected as the optimal RB combination, the amplification effect on errors is minimum, but matrix inversion operation is required to be carried out when the formula (7) is adopted for carrying out calculation on G, the real-time performance is reduced due to the large calculation amount, and therefore the calculation of the formula (7) is required to be simplified;
(5) The operation of the formula (7) is simplified, the matrix inversion is equivalent to the sum of the eigenvalues of the matrix, and the corresponding formula (7) can be written as
Figure BDA0002905890320000112
In formula (8): lambda [ alpha ] 1 、λ 2 Is a matrix H T H characteristic value;
can finally obtain
Figure BDA0002905890320000113
In formula (9), det [ ] represents a determinant for solving a matrix;
obtaining a geometric accuracy characteristic value G, sequencing the N RB combinations according to the geometric accuracy characteristic value G, and when the geometric accuracy characteristic value G is smaller, indicating that the RB combination state is better;
step two: as can be seen from equation (1), when the vehicle is co-located, the position coordinate of the RB must be obtained, and usually, the position coordinate of the RB has a certain error, which inevitably affects the accuracy of the co-location; therefore, when the RB is selected, the position coordinate precision of the RB must be evaluated and corrected, and the RB with higher position precision is selected as the optimal RB combination;
through analyzing the motion characteristics of the vehicle, the vehicle motion has a certain motion rule in a relative time period, and the motion rule can be described by adopting a current statistical model of a motor carrier;
therefore, the improved current statistical model is used as a state equation of the vehicle, and the RB position coordinate sent by the RB is used as an observed quantity to carry out error correction and precision evaluation on the RB position coordinate, so that error interference on a cooperative positioning result is eliminated to a certain extent;
the method comprises the steps that a CSM model introducing variance adjustment factors is a state equation, position coordinates of RB are used as observed quantity, the position coordinates of RB are corrected in real time by adopting a formula (17), meanwhile, the precision of the position coordinates of RB is evaluated according to the difference value between the position coordinates of RB and the corrected value, and in addition, the stable stability of the relative motion state of RB is evaluated according to the position coordinates of RB;
the derivation process of utilizing the CSM model introduced with the variance adjustment factor as a state equation and utilizing the position coordinate of the RB as the observed quantity to modify the position coordinate of the RB in real time comprises the following steps:
s2.1 Current statistical motion model
The vehicle has a variable acceleration characteristic, and the current statistical model adopts a time-varying acceleration probability density function and an acceleration non-zero mean value time correlation model, so that the motion characteristic of the moving vehicle can be better described;
and CSM is adopted as an RB carrier motion state equation, and the equation is expressed as follows:
Figure BDA0002905890320000121
in formula (10), X k+1 Is a state vector, phi k Being a state transition matrix, U k In order to input the control matrix, the control matrix is input,
Figure BDA0002905890320000122
mean value of acceleration at time k, W k Mean value of zero and variance of Q k A gaussian distributed noise vector of (a); the described
Figure BDA0002905890320000123
Figure BDA0002905890320000124
Figure BDA0002905890320000131
Figure BDA0002905890320000132
In formula (11), x k
Figure BDA0002905890320000133
Position, velocity and acceleration of RB, respectively; in the formulas (12) and (13), T is a sampling period; in the formula (14), τ is the maneuvering frequency;
Figure BDA0002905890320000134
is the variance of the vehicle maneuvering acceleration; q is a noise matrix;
s2.2 selecting the mean value of the acceleration
Figure BDA0002905890320000135
Is taken as the current acceleration prediction value, i.e.
Figure BDA0002905890320000136
Variance of acceleration
Figure BDA0002905890320000137
Is taken as
Figure BDA0002905890320000138
In formula (16): a is a max And a min Respectively representing the maximum value and the minimum value of the acceleration;
the equation for filter correction of the position coordinates is
Figure BDA0002905890320000139
Correcting the position coordinates of the RB in real time by using the formula (17);
from the analysis of equations (14) and (16) it follows: system variance Q of motion model k Variance with acceleration
Figure BDA00029058903200001310
Is in direct proportion; a is max And a min After the value is determined, when the carrier maneuvers at a small acceleration, the system variance is large, and the filtering precision is low; when the carrier moves at a larger acceleration, the system variance is smaller, and the filtering precision is higher;
therefore, the conventional CSM has high filtering accuracy for a highly mobile carrier and low filtering accuracy for a weakly mobile carrier; CSM cannot describe the acceleration value as the interval [ (4-pi) a) -max /4,(4-π)a max /4]And weak maneuvering, resulting in lower filtering accuracy.
From the analysis of equations (14) and (16) it follows: system variance Q of motion model k Variance with acceleration
Figure BDA0002905890320000141
Proportional ratio of a max And a min After the value is determined, when the carrier maneuvers with a small acceleration, the system variance is large, and the filtering precision is low; when the carrier moves at a larger acceleration, the system variance is smaller, and the filtering precision is higher;
wherein the analysis finds the acceleration variance
Figure BDA0002905890320000142
The incoordination with the carrier mobility is the main reason of poor filtering effect, therefore, the invention designs the variance adjustment factor eta aiming at the problem k Adjusting the variance by a factor eta k Step S2.2 is introduced to the acceleration variance
Figure BDA0002905890320000143
Performing self-adaptive adjustment to improve CSM filtering accuracy, wherein the variance adjustment factor eta k The design process of (2) comprises:
(1) Let the innovation vector of the filter in KF be defined as:
Figure BDA0002905890320000144
innovation vector d k Is not relevant in the ideal case, and d k Is that the mean is zero and the variance is S k When the RB carrier moves, the movement changes the orthogonality of the new information, so that d k The mean value of (a) changes and is no longer zero;
Figure BDA0002905890320000145
(2) Innovation vector d k Obeying a multidimensional Gaussian distribution with a variance S k Subject to a degree of freedom of m (m being the vector dimension) 2 Distribution, it is generally possible to detect the mobility of the support and the divergence of the filters according to this multidimensional distribution property of the innovation, but this straightforward way is complex and therefore it is possible to normalize the sequence of innovation to obtain the statistic g k
Figure BDA0002905890320000146
(3) Establishing window detection statistic, and setting window size as m, window statistic wg of k moment k Is defined as
Figure BDA0002905890320000147
Coefficient of determination r k Is defined as
r k =wg k /wg k-1 (22)
When the RB carrier is not mobile, r k The value is usually close to 1, r when a maneuver occurs k The value increases rapidly and is proportional to the degree of manoeuvre;
(4) In order to convert the mobility of the RB carrier into a variance adjustment factor, a raised-seminormal distribution function mu (u) is introduced,
Figure BDA0002905890320000151
(5) Will r is k The variance adjustment factor eta can be obtained as an input variable of the ascending seminormal distribution function mu (u) k
η k =μ(r k ) (24);
Obtaining variance adjustment factor eta k Adjusting the factor eta according to the variance k Adjusting the maneuvering acceleration variance in real time to obtain an acceleration variance updating formula:
Figure BDA0002905890320000152
in the formula:
Figure BDA0002905890320000153
is the average acceleration at time k; (ii) a
Figure BDA0002905890320000154
Acceleration variance according to equation (25)
Figure BDA0002905890320000155
And system variance Q k Can be self-adaptiveAdjusting, when the mobility of the carrier is weak, the variance adjusting factor eta k Is smaller, so that
Figure BDA0002905890320000156
And Q k The variance adjustment factor is introduced to ensure that the filtering keeps better robustness to model uncertainty.
Step three: comprehensively evaluating RB relative position characteristics, position precision characteristics and relative motion characteristics by adopting a fuzzy evaluation method, and selecting an optimal RB cooperation positioning combination; establishing a cooperative positioning equation set by using the position coordinate of the optimal RB and the relative distance between the optimal RB and the VTBP, and solving the position coordinate of the VTBP by adopting a Taylor series analytic method, wherein the specific process comprises the following steps:
s3.1, comprehensively evaluating the relative position characteristic, the position precision characteristic and the relative motion characteristic of the RB combination by adopting a fuzzy evaluation method, and selecting an optimal RB cooperative positioning combination;
firstly, carrying out normalization processing on relative position characteristics, position precision characteristics and relative motion characteristics of RB combinations; wherein the relative position features are geometric position precision values in the first step; the position precision characteristic is the average value of the difference value between the original value and the corrected value of each RB vehicle position in the step two; the relative motion characteristic is the average value of the difference value between the RB vehicle motion speed and the VTBP vehicle motion speed;
secondly, carrying out weighted summation on the normalized relative position characteristic, position precision characteristic and relative motion characteristic, wherein the summed value is a comprehensive evaluation value, and the smaller the value is, the better the corresponding RB vehicle combination is;
s3.2, establishing a cooperative positioning equation set by utilizing the position coordinates of the optimal RB and the relative distance between the optimal RB and the VTBP, solving the position coordinates of the VTBP by adopting a Taylor series analytic method, and establishing a cooperative positioning equation
Figure BDA0002905890320000161
The positioning coordinate is calculated by adopting a Taylor series method, and the calculation process is as follows:
step1: setting an initial coordinate estimation position (x) 0 ,y 0 ) And an error threshold, the error of the estimated position from the actual value of the coordinates being (e) x ,e y );
Step2: set the localization equations in (x) 0 ,y 0 ) Expanding by Taylor series, neglecting terms of two or more
Figure BDA0002905890320000162
Transforming the above equation into matrix form: a. The b d=b
Figure BDA0002905890320000171
When A is b When full rank, A b Reversible;
it can be derived that:
Figure BDA0002905890320000172
step3: the position of the updated coordinates is x = x 0 +e x ,y=y 0 +e y (ii) a Judgment (e) x ,e y ) If the threshold value is smaller than the set threshold value, the iteration is stopped; if the error is larger than the set threshold value, the Step is switched to Step2 to continue iteration until the error is smaller than the set error threshold value.
Step four: establishing an observation equation by taking the solved VTBP position coordinate as an observed quantity, establishing a dead reckoning model by the speed and course angle measured by the VTBP, and correcting the VTBP position coordinate by adopting an extended Kalman to obtain a final VTBP position coordinate value so as to complete vehicle cooperative positioning, wherein the specific process comprises the following steps:
s4.1, establishing an observation equation by taking the solved VTBP position coordinate as an observed quantity, and establishing a dead reckoning model by the speed and the course angle measured by the VTBP:
Figure BDA0002905890320000173
in formula (30): x is the number of k Is the east position coordinate of the vehicle; y is k Is the north position coordinate of the vehicle; v. of k Is the running speed of the vehicle; theta k Is the driving steering angle.
S4.2 establishing a measurement model
The relative distance and the relative angle between two vehicles are used as observed quantities, and a system measurement equation generated by combining the position coordinates of the reference vehicle transmitted by the reference vehicle is
Figure BDA0002905890320000174
In formula (31):
Figure BDA0002905890320000175
position coordinates of the reference vehicle at the moment k; r is k The relative distance between vehicles at the moment k;
Figure BDA0002905890320000176
is the relative angle between the two workshops;
Figure BDA0002905890320000177
the angle between the reference vehicle and the X-axis direction;
s4.3 correcting the VTBP position coordinate by adopting an extended Kalman
Figure BDA0002905890320000181
Figure BDA0002905890320000182
Figure BDA0002905890320000183
In the filtering iterative process R k In order to measure the variance matrix of the error,in general R k Is determined according to the error of the measured vector, wherein
Figure BDA0002905890320000184
Through the process, the final VTBP position coordinate value is obtained, and vehicle cooperative positioning is completed; by utilizing the Vehicle cooperative positioning method based on the brain selective attention mechanism, the influence of a geometric structure, RB position precision and a relative motion state formed between a cooperative Reference node (RB) and a Vehicle To Be Positioned (VTBP) on cooperative positioning is specifically aimed at, a selective attention mechanism of human brain on information processing is adopted, and the influence factors of the three aspects are used as characteristic points of information selection to selectively filter information; and then, the optimal RB is comprehensively screened out through feature integration, and the reliability and the precision of vehicle cooperative positioning are further improved.
Example 1: data emulation
The coverage area of the Internet of vehicles is assumed to be a road area with the length of 500 meters and the width of 500 meters; the road area has 5 RB vehicles in total, each vehicle moves along different elliptical tracks, and the initial coordinates are (60.0, 15.4) m, (230.0, 400.1) m, (425, 127.9) m, (320.4, 480.0) m, and (480.0, 332.5) m; the movement speeds are respectively 20m/s, 40m/s, 30m/s, 45m/s and 25m/s; the motion track is shown in fig. 5;
the starting position of the VTBP vehicle is (250.0, 90.0), and the movement speed is 50m/s; in the experiment, the method and the traditional Taylor series positioning method are adopted, the number of RB vehicles selected by the two methods is 3, and the traditional Taylor series positioning method randomly selects 3 vehicles from 5 RB vehicles; fig. 6 is a comparison graph of experimental effects, in 150s experimental verification, the average error of the method of the present invention is 1.9m, while the average error of the conventional taylor series positioning method is 4.9m, and the effectiveness of the method of the present invention can be seen from comparison data.
The foregoing shows and describes the general principles, principal features, and advantages of the invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are described in the specification and illustrated only to illustrate the principle of the present invention, but that various changes and modifications may be made therein without departing from the spirit and scope of the present invention, which fall within the scope of the invention as claimed. The scope of the invention is defined by the appended claims and equivalents thereof.

Claims (4)

1. A vehicle cooperative localization method based on brain selective attention mechanism is characterized in that: comprises the steps of
The method comprises the following steps: in a VANET environment, selecting N RB combinations from VTBP neighbor vehicles or ITS infrastructure nodes, calculating geometric precision characteristic values of the selected N RB combinations, and sorting the advantages and the disadvantages according to an ascending order;
step one, the selection of the N RB combinations and the calculation process of the geometric accuracy characteristic values comprise the following steps:
s1.1, setting the position coordinates of VTBP as unknown quantity (x, y), and selecting n neighbor vehicles or ITS infrastructure nodes from the surrounding environment as positioning RB i Wherein i =1,2, \8230;, n; RB (radio B) i Is a known quantity (x) i ,y i );
S1.2 obtaining VTBP and RB i Relative distance d between i Comprises the following steps:
Figure FDA0003988594070000011
s1.3, calculating geometric accuracy characteristic values of the selected N RB combinations, and performing quality sorting according to an ascending order;
step two: introducing a CSM model of a variance adjustment factor as a state equation, taking the position coordinate of the RB as an observed quantity, correcting the position coordinate of the RB, and simultaneously evaluating the position coordinate precision of the RB;
the process of correcting the position coordinates of the RB in real time comprises the following steps:
s2.1 adopts CSM as an RB carrier motion state equation, and is expressed as follows:
Figure FDA0003988594070000012
in the formula (10), X k+1 Is a state vector, phi k Being a state transition matrix, U k In order to input the control matrix, the control matrix is input,
Figure FDA0003988594070000013
mean value of acceleration at time k, W k Is a mean of zero and a variance of Q k Gaussian distributed noise vector of (2); the above-mentioned
Figure FDA0003988594070000014
Figure FDA0003988594070000021
Figure FDA0003988594070000022
Figure FDA0003988594070000023
In formula (11), x k
Figure FDA0003988594070000024
Position, velocity and acceleration of RB, respectively; in equations (12) and (13), T is a sampling period; in the formula (14), τ is the maneuvering frequency,
Figure FDA0003988594070000025
is the variance of the maneuvering acceleration of the carrier, and q is a noise matrix;
s2.2 selection
Figure FDA0003988594070000026
For the mean acceleration at time k, a variance adjustment factor eta is introduced k =μ(r k ) Obtaining:
Figure FDA0003988594070000027
variance of acceleration
Figure FDA0003988594070000028
The update formula of (2) is:
Figure FDA0003988594070000029
in formula (25): a is max And a min Respectively representing the maximum value and the minimum value of the acceleration;
the equation for performing filtering correction on the position coordinates of the RB carrier is as follows:
Figure FDA00039885940700000210
step S2.2 the variance adjustment factor eta k The introduction process of (2) comprises:
(1) Let the innovation vector of the filter in KF be defined as:
Figure FDA0003988594070000031
innovation vector D k Is not relevant in the ideal case, and D k Is that the mean is zero and the variance is S k When the RB carrier moves, the movement changes the orthogonality of the new information, so that D is k Has changed mean value and is no longer zero, i.e.
Figure FDA0003988594070000032
(2) Normalizing the sequence of the innovation vectors to obtain statistics
Figure FDA0003988594070000033
(3) Establishing window detection statistic, and setting window size as m, window statistic wg of k moment k Is defined as
Figure FDA0003988594070000034
Coefficient of determination r k Is defined as
r k =wg k /wg k-1 (22)
When the RB carrier is not mobile, r k Value close to 1, r when a maneuver occurs k The value increases rapidly and is proportional to the degree of mobility;
(4) In order to convert the mobility of the RB carrier into a variance adjustment factor, a raised-seminormal distribution function mu (u) is introduced,
Figure FDA0003988594070000035
(5) The variance adjusting factor eta can be obtained by taking rk as the input variable of the raised seminormal distribution function mu (u) k
h k =μ(r k ) (24);
Step three: comprehensively evaluating the RB relative position characteristic, the position precision characteristic and the relative motion characteristic; establishing a cooperative positioning equation set by using the evaluated optimal RB position coordinate and the relative distance between the optimal RB position coordinate and the VTBP, and solving the VTBP position coordinate;
step four: and establishing a dead reckoning model by taking the solved VTBP position coordinate as an observed quantity, and correcting the VTBP position coordinate, so as to obtain a final VTBP position coordinate value and complete vehicle cooperative positioning.
2. The brain-selective attention mechanism-based vehicle cooperative localization method according to claim 1, wherein: the calculation process of the geometric accuracy characteristic value in the step S1.3 includes:
(1) Introducing measurement error e i Then there is an error of (e) in the position of the corresponding VTBP x ,e y ) Then, the formula (1) after adding the error is:
Figure FDA0003988594070000041
converting the formula (2) into a linear equation to obtain
Figure FDA0003988594070000042
(2) Order to
Figure FDA0003988594070000043
a i =(x-x i )/d i ,b i =(y-y i )/d i
Equation (3) can be converted to a matrix form:
L=HX+E (4)
in formula (4): l = [ L = 1 ,l 2 ,…l n ] T ,X=[E x ,E y ] T ,E=[E 1 ,E 2 ,…E n ] T
Figure FDA0003988594070000044
(3) When H is full rank, H T H is reversible, then the error between the position estimate and the true value is:
Figure FDA0003988594070000051
the magnitude of the error in equation (5) is measured by covariance and yields:
Figure FDA0003988594070000052
in formula (6): sigma 2 The variance of the noise in e when the noise is uncorrelated two by two;
(4) As can be seen from formula (6), (H) T H) -1 Expressed as a magnification factor for the distance measurement error, and therefore defining the geometric position accuracy factor of the cooperative positioning of the vehicle as G, then
Figure FDA0003988594070000053
Wherein, in the formula (7), tr [ ] represents the inverse operation of the matrix;
(5) The operation of the formula (7) is simplified, matrix inversion is equivalent to the sum of matrix eigenvalues, and the corresponding formula (7) can be written as
Figure FDA0003988594070000054
In formula (8): lambda 1 、λ 2 Is a matrix H T A characteristic value of H;
and finally, obtaining a geometric position precision factor of vehicle cooperative positioning:
Figure FDA0003988594070000055
in the formula (9), det [ ] represents a determinant for matrix calculation.
3. The brain-selective attention mechanism-based vehicle cooperative localization method according to claim 1, wherein: the concrete process of the third step comprises:
s3.1, comprehensively evaluating the relative position characteristic, the position precision characteristic and the relative motion characteristic of the RB combination by adopting a fuzzy evaluation method, and selecting an optimal RB cooperation positioning combination;
and S3.2, establishing a cooperative positioning equation set by utilizing the position coordinate of the optimal RB and the relative distance between the optimal RB and the VTBP, solving the position coordinate of the VTBP by adopting a Taylor series analytic method, and establishing a cooperative positioning equation.
4. The brain-selective attention mechanism-based vehicle cooperative localization method according to claim 1, wherein: the concrete process of the step four comprises:
s4.1, establishing an observation equation by taking the solved VTBP position coordinate as an observed quantity, and establishing a dead reckoning model by the speed and the course angle measured by the VTBP per se;
s4.2, establishing a measurement model by taking the relative distance and the relative angle between the two vehicles as observed quantities and combining the position coordinates of the reference vehicle transmitted by the reference vehicle;
and S4.3, correcting the VTBP position coordinate by adopting an extended Kalman to obtain a final VTBP position coordinate value, and finishing the vehicle cooperative positioning.
CN202110071311.7A 2021-01-19 2021-01-19 Vehicle cooperative positioning method based on brain selective attention mechanism Active CN112880699B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110071311.7A CN112880699B (en) 2021-01-19 2021-01-19 Vehicle cooperative positioning method based on brain selective attention mechanism

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110071311.7A CN112880699B (en) 2021-01-19 2021-01-19 Vehicle cooperative positioning method based on brain selective attention mechanism

Publications (2)

Publication Number Publication Date
CN112880699A CN112880699A (en) 2021-06-01
CN112880699B true CN112880699B (en) 2023-03-10

Family

ID=76049965

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110071311.7A Active CN112880699B (en) 2021-01-19 2021-01-19 Vehicle cooperative positioning method based on brain selective attention mechanism

Country Status (1)

Country Link
CN (1) CN112880699B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116482716B (en) * 2023-06-26 2023-08-29 北京航空航天大学 Node fault detection method for space-based navigation enhanced ad hoc network

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007303841A (en) * 2006-05-08 2007-11-22 Toyota Central Res & Dev Lab Inc Vehicle position estimation device
US8457891B1 (en) * 2012-06-19 2013-06-04 Honeywell International Inc. Systems and methods for compensating nonlinearities in a navigational model
CN104990554A (en) * 2015-05-04 2015-10-21 南京邮电大学 Inertial navigation positioning method in GNSS blind area based on cooperation between VANET vehicles
CN107728138A (en) * 2017-09-15 2018-02-23 电子科技大学 A kind of maneuvering target tracking method based on current statistical model
CN108519738A (en) * 2018-04-13 2018-09-11 中国科学院微电子研究所 A kind of state of motion of vehicle information delay compensation method and device
CN110986957A (en) * 2019-12-24 2020-04-10 中国人民解放军空军工程大学 Three-dimensional flight path planning method and device for unmanned aerial vehicle
CN111163419A (en) * 2020-02-07 2020-05-15 北京大学 Malicious user detection method based on state mean value in vehicle cooperation dynamic tracking
CN111586632A (en) * 2020-05-06 2020-08-25 浙江大学 Cooperative neighbor vehicle positioning method based on communication sensing asynchronous data fusion

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6697736B2 (en) * 2002-02-06 2004-02-24 American Gnc Corporation Positioning and navigation method and system thereof
US7894512B2 (en) * 2007-07-31 2011-02-22 Harris Corporation System and method for automatic recovery and covariance adjustment in linear filters
CN103383261B (en) * 2013-07-02 2015-11-18 河海大学 A kind of modified can't harm the indoor moving targets location method of Kalman filtering
CN103914985B (en) * 2014-04-25 2015-10-28 大连理工大学 A kind of hybrid power passenger car following speed of a motor vehicle trajectory predictions method
US11713967B2 (en) * 2017-10-13 2023-08-01 JVC Kenwood Corporation Angular speed derivation device and angular speed derivation method for deriving angular speed based on output value of triaxial gyro sensor
US11636375B2 (en) * 2018-02-27 2023-04-25 Toyota Research Institute, Inc. Adversarial learning of driving behavior

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007303841A (en) * 2006-05-08 2007-11-22 Toyota Central Res & Dev Lab Inc Vehicle position estimation device
US8457891B1 (en) * 2012-06-19 2013-06-04 Honeywell International Inc. Systems and methods for compensating nonlinearities in a navigational model
CN104990554A (en) * 2015-05-04 2015-10-21 南京邮电大学 Inertial navigation positioning method in GNSS blind area based on cooperation between VANET vehicles
CN107728138A (en) * 2017-09-15 2018-02-23 电子科技大学 A kind of maneuvering target tracking method based on current statistical model
CN108519738A (en) * 2018-04-13 2018-09-11 中国科学院微电子研究所 A kind of state of motion of vehicle information delay compensation method and device
CN110986957A (en) * 2019-12-24 2020-04-10 中国人民解放军空军工程大学 Three-dimensional flight path planning method and device for unmanned aerial vehicle
CN111163419A (en) * 2020-02-07 2020-05-15 北京大学 Malicious user detection method based on state mean value in vehicle cooperation dynamic tracking
CN111586632A (en) * 2020-05-06 2020-08-25 浙江大学 Cooperative neighbor vehicle positioning method based on communication sensing asynchronous data fusion

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Sang-Hyuk Park ; Young-Joong Kim ; Hoo-Cheol Lee ; Myo-Taeg Lim.Improved adaptive particle filter using adjusted variance and gradient data.《2008 IEEE International Conference on Multisensor Fusion and Integration for Intelligent Systems》.2008,全文. *
基于速度约束与模糊自适应滤波的车载组合导航;胡杰,严勇杰,王子卉;《兵工学报》;20200229;第41卷(第2期);全文 *

Also Published As

Publication number Publication date
CN112880699A (en) 2021-06-01

Similar Documents

Publication Publication Date Title
Balico et al. Localization prediction in vehicular ad hoc networks
Chelouah et al. Localization protocols for mobile wireless sensor networks: A survey
CN109946731B (en) Vehicle high-reliability fusion positioning method based on fuzzy self-adaptive unscented Kalman filtering
Liu et al. Mercury: An infrastructure-free system for network localization and navigation
Schubert et al. Comparison and evaluation of advanced motion models for vehicle tracking
CN108534779B (en) Indoor positioning map construction method based on track correction and fingerprint improvement
CN110264721B (en) Urban intersection surrounding vehicle track prediction method
Xu et al. Spatial-temporal constrained particle filter for cooperative target tracking
CN111324848B (en) Vehicle-mounted track data optimization method of mobile laser radar measurement system
Hao et al. Modal activity-based stochastic model for estimating vehicle trajectories from sparse mobile sensor data
CN112367614B (en) LSTM-based Wi-Fi and geomagnetic field fusion indoor positioning algorithm
CN110849355B (en) Bionic navigation method for geomagnetic multi-parameter multi-target rapid convergence
Baek et al. Accurate vehicle position estimation using a Kalman filter and neural network-based approach
CN112880699B (en) Vehicle cooperative positioning method based on brain selective attention mechanism
CN111307143A (en) Bionic navigation algorithm for multi-target evolution search based on geomagnetic gradient assistance
Wu et al. CLSTERS: A general system for reducing errors of trajectories under challenging localization situations
KR102119196B1 (en) Method and system for indoor positioning based on machine learning
Ten Kathen et al. A comparison of pso-based informative path planners for detecting pollution peaks of the ypacarai lake with autonomous surface vehicles
Magnano et al. Movement prediction in vehicular networks
Shi et al. Indoor localization scheme using magnetic map for smartphones
Zhang et al. An adaptive road-constrained IMM estimator for ground target tracking in GSM networks
Xiao et al. A GPR-PSO incremental regression framework on GPS/INS integration for vehicle localization under urban environment
Li et al. Sliding mode control for vehicular platoon based on v2v communication
Goli et al. An accurate multi-sensor multi-target localization method for cooperating vehicles
CN113852911B (en) Fusion positioning method based on fingerprint library and PDR calculation and fingerprint library updating method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant