CN109946731A - A kind of highly reliable fusion and positioning method of vehicle based on fuzzy self-adaption Unscented kalman filtering - Google Patents

A kind of highly reliable fusion and positioning method of vehicle based on fuzzy self-adaption Unscented kalman filtering Download PDF

Info

Publication number
CN109946731A
CN109946731A CN201910168995.5A CN201910168995A CN109946731A CN 109946731 A CN109946731 A CN 109946731A CN 201910168995 A CN201910168995 A CN 201910168995A CN 109946731 A CN109946731 A CN 109946731A
Authority
CN
China
Prior art keywords
uwb
vehicle
observation
kalman filtering
distance
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
Application number
CN201910168995.5A
Other languages
Chinese (zh)
Other versions
CN109946731B (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.)
Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN201910168995.5A priority Critical patent/CN109946731B/en
Publication of CN109946731A publication Critical patent/CN109946731A/en
Application granted granted Critical
Publication of CN109946731B publication Critical patent/CN109946731B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Navigation (AREA)

Abstract

The invention discloses a kind of highly reliable fusion and positioning methods of the vehicle towards typical urban environment, for under urban environment, Vehicular satellite navigation is easily using limited, the problems such as positioning accuracy and not high reliability, on the basis of traditional Vehicular satellite and inertia combined navigation, introduce ultra wide band (Ultra-Wideband, UWB) location technology, and classified using fuzzy algorithmic approach to UWB accuracy of observation, to the observation noise variance matrix of automatic adjusument UWB, herein on basis, the fusion positioning of vehicle is realized based on Unscented kalman filtering algorithm.Compared to traditional Vehicular satellite and inertia combined navigation, method in the present invention, under the complex environment (such as urban canyons, intersection) that urban environment, especially satellite-signal are seriously blocked, reliability is higher, helps to realize continuous, complete, reliable, the real-time positioning of vehicle.

Description

A kind of highly reliable fusion positioning of the vehicle based on fuzzy self-adaption Unscented kalman filtering Method
Technical field
The present invention relates to automobile navigation positioning field, in particular to the highly reliable fusion of a kind of vehicle towards urban environment is fixed Position method.
Background technique
With the development and progress of economic society, the vehicle guaranteeding organic quantity in China is quicklyd increase, and road traffic is faced with huge Big challenge, in order to solve increasingly serious urban transport problems, intelligent bus or train route cooperative system (Intelligent Vehicle Infrastructure Cooperative Systems, IVICS) it comes into being, and gradually become intelligent transportation system The latest development direction of (Intelligent Transport Systems, ITS) research.The either application of bus or train route collaboration, also It is the realization of intelligent transportation, all be unable to do without high-precision vehicle positioning technology: only in the premise for realizing accurate reliable vehicle positioning Under, comprehensive implementation vehicle vehicle, bus or train route dynamic realtime information exchange, and in full space-time dynamic traffic information collection and the basis merged Upper development vehicle active safety control and road coordinated management are sufficiently realized effective collaboration of people's bus or train route, can effectively be referred to It waves and dispatches buses, improve urban transportation, guarantee vehicle safe driving.Therefore, vehicle positioning technology is bus or train route collaboration or even intelligence One of the basis of the researchs such as traffic and core content.
Automobile navigation location technology common at present includes: dead reckoning (Dead Reckoning, DR), inertial navigation system It unites (Inertial Navigation System, INS), satellite navigation (Global Navigation Satellite System, GNSS) etc..Since the presence of single location technology is respectively insufficient, in order to realize under urban environment it is accurate, reliable, Continuously, it completely positions, mostly uses two kinds and two or more location technologies combines, wherein GNSS/INS fusion positioning is answered With the most extensively: GNSS/INS fusion positioning can solve the orientation problem in the case of satellite-signal short interruptions, and certain The accumulated error that INS is compensated in degree can satisfy the vehicle location demand under opposite free environments, but when satellite-signal is long Between interrupt environment (such as urban canyons, tunnel) in, if GNSS can not work normally, only rely on INS and independently calculate, can still lead Cause biggish position error.
In recent years, based on ultra wide band (Ultra-Wideband, UWB) wireless location technology rise and rapid development be Realize that reliable location of the vehicle under GNSS constrained environment provides new approaches, currently, UWB location technology is primarily used to room Interior positioning field such as personnel, intelligent carriage, robot localization etc..Since the very bandwidth of UWB is wide, pulse signal penetration power is strong, more The well equal technological merits of diameter resolution capability, realize that the vehicle location under outdoor environment also has certain feasibility based on UWB, In sensor layer, if introducing relatively reliable UWB information source, can effectively improve fusion positioning system robustness, but Under complicated urban environment, vehicle, the trees of road both sides and the interference of building that UWB signal is travelled vulnerable to surrounding and go out Existing multipath and non-market value certainly will lead to fusion essence if being used to merge positioning for the poor UWB observed quantity of signal quality Therefore how the decline of degree is directed to relative complex traffic scene (such as right-angled intersection, rotary island crossing), in conjunction with road Means of transportation, rational deployment trackside UWB node are realized people's bus or train route information dynamic interaction under bus or train route collaboration, and then are chosen properly Method, automatic adjusument blending algorithm is to realize the vehicle based on UWB to the degree of dependence of the UWB observation information of different quality Fusion positions critical issue urgently to be resolved.
Summary of the invention
Technical problem: under urban environment, easily using being limited, positioning accuracy and reliability be not high for Vehicular satellite navigation The problems such as, on the basis of traditional Vehicular satellite and inertia combined navigation, introduce ultra wide band (Ultra-Wideband, UWB) location technology, and classified using fuzzy algorithmic approach to UWB accuracy of observation, thus the observation noise variance of automatic adjusument UWB Herein on basis, the fusion positioning of vehicle is realized based on Unscented kalman filtering algorithm for battle array.
Technical solution: to achieve the above object, the present invention adopts the following technical scheme: firstly, arranging UWB in road both sides Stationary nodes, and pass through the position coordinates of high-precision difference GPS acquisition UWB stationary nodes;Secondly, arranging one in vehicle roof UWB mobile node, using reach time difference method, obtain the UWB mobile node to each UWB stationary nodes distance;Then, The rough position coordinate of vehicle is obtained based on ranging localization principle;Then, according to the distance and bearing corner characteristics of UWB signal, base Classify in accuracy of observation of the fuzzy logic classifier algorithm to UWB;And then it is made an uproar according to fuzzy classification result automatic adjusument UWB observation Acoustic matrix;Finally, realize that the fusion positioning of vehicle obtains the exact position of vehicle based on Unscented kalman filtering.
With reference to the accompanying drawing, thinking of the invention is further described:
Process of the invention is as shown in Figure 1.
A kind of highly reliable fusion and positioning method of vehicle based on fuzzy self-adaption Unscented kalman filtering, which is characterized in that In urban environment, on the basis of Vehicular satellite and inertia combined navigation, UWB positioning is introduced, and combine fuzzy algorithmic approach to UWB Accuracy of observation classification, so that the observation noise variance matrix of automatic adjusument UWB, herein on basis, is filtered based on Unscented kalman Wave algorithm realizes the highly reliable fusion positioning of vehicle, and described method includes following steps:
Step 1) arranges UWB stationary nodes in road both sides, and measures the position coordinates of UWB stationary nodes;
Step 2) arranges a UWB mobile node in vehicle roof, using arrival time difference method (Time Difference Of Arrival, TDOA) the UWB mobile node is obtained to the distance of each UWB stationary nodes, and observation of adjusting the distance is counted Data preprocess;
Step 3), apart from observation, establishes equation in coordinates based on ranging localization principle, using minimum according to pretreated Square law calculates the coordinate of UWB mobile node, obtains the rough position coordinate of vehicle;
Step 4) measures the current course angle of vehicle by vehicle electronics compass, and according to the rough position coordinate of vehicle with And the position coordinates of UWB stationary nodes, calculate azimuth of each UWB stationary nodes relative to vehicle current driving direction;
Step 5) is according to two, distance and bearing angle feature, using fuzzy logic classifier algorithm to each UWB stationary nodes Accuracy of observation classify, one is divided into A, B, C, five grades of D, E, the precision apart from observation of corresponding UWB by height to Low, specific step is as follows for fuzzy Classified Algorithms Applied:
51) three subordinating degree functions of distance feature, and adjust the distance classification, respectively N, M and R are determined, distance is represented " close ", " in ", " remote ";
52) two subordinating degree functions of orientation corner characteristics are determined, and azimuthal is classified, respectively L and S represent orientation " big " and " small " at angle;
53) according to Fuzzy classification rule, 5 accuracy classes A, B, C, D, E of fuzzy logic classifier output are determined, are obscured Classifying rules is as follows:
1. when distance is N, when azimuth is L, accuracy class A;
2. when distance is N, when azimuth is S, accuracy class B;
3. when distance is M, when azimuth is L, accuracy class C;
4. when distance is M, when azimuth is S, accuracy class D;
5. when distance is R, when azimuth is L, accuracy class E;
6. when distance is R, when azimuth is S, accuracy class E;
Step 6) combines the observation information of vehicle-mounted GNSS, gyroscope and wheel speed sensors output, utilizes Unscented kalman The accurate positioning of filtering algorithm realization vehicle, the specific steps are as follows:
61) firstly, establishing the Unscented kalman filtering state model of vehicle positioning system
System mode vector is
Wherein,For vehicle latitude information, λ is vehicle longitude information, and h is height of car information;veFor the east orientation speed of vehicle Degree, vnFor the north orientation speed of vehicle, vuFor the vertical velocity of vehicle;P is the pitch angle of vehicle, and r is the angle of heel of vehicle, and A is vehicle Course angle;sfodFor wheel speed dilution of precision error, bωFor Gyroscope Random Drift, bGNSSFor GNSS receiver clock deviation, dGNSS For GNSS receiver clock deviation drift rate, bUWBFor the offset error of UWB;
Exterior input vector is U=[vod aod fx fyωz]
Wherein, vodFor the car speed that wheel speed sensors measure, aodFor longitudinal acceleration of the vehicle, fxFor longitudinal direction of car acceleration Spend measured value, fyFor vehicle lateral acceleration measured value, ωzFor the yaw velocity of gyroscope output;
System state equation are as follows:
Wherein, X (k) is the system mode vector of discrete instants k, and X (k-1), U (k-1), W (k-1), T (k-1) are respectively The system mode vector of discrete instants (k-1), external input vector, system noise vector, external input noise vector;RMFor ground Ball meridian radius, RNFor earth prime vertical radius, Δ t is the sampling interval;γwFor the error of wheel speed sensors and time correlation, σwFor its corresponding noise variance;βzFor Gyroscope Random Drift, σzFor its corresponding noise variance;
62) the Unscented kalman filtering observation model of vehicle positioning system then, is established
Systematic observation vector is
Wherein,The pseudo range observed quantity of respectively m GNSS satellite,It is n UWB stationary nodes apart from observation;
Systematic observation equation are as follows:
Z (k)=H (k) X (k)+V (k)
Wherein, Z (k), H (k) are the systematic observation vector and systematic observation matrix of discrete instants k, when V (k) is discrete Carve the observation noise vector of k;
The corresponding observation noise variance matrix of observation noise vector V (k) is R (k):
R (k)=diag [RGNSS RUWB]
Wherein, RGNSSFor the observation noise variance matrix of GNSS, RUWBFor the observation noise variance matrix of UWB;
63) fuzzy Classified Algorithms Applied according to step 5) classifies to the accuracy of observation of n UWB stationary nodes, And then according to the noise variance matrix R of classification results automatic adjusument UWBUWB, final observation noise variance matrix R (k) is obtained, is had Body method is as follows:
1. calculating the observation noise of i-th of UWB stationary nodes
Wherein, R0Noise variance basic value for UWB apart from observation;μ is adjustment factor, true by the result of fuzzy classification Fixed, five kinds of different corresponding adjustment factors of accuracy class A, B, C, D, E are respectively μa, μb, μc, μd, μe
2. calculating the observation noise variance matrix of UWB
3. updating observation noise variance matrix R (k): R (k)=diag [RGNSS RUWB];
64) System State Model according to step 61), observation model described in step 62), and through step 63) Noise variance matrix after automatic adjusument, the time for carrying out Unscented kalman filtering updates and measurement updaue process, obtains vehicle Accurate location information.
1. constructing the sigma point set ξ of Unscented kalman filteringi(k-1), and its weight coefficient is determined
Wherein, η is scale parameter, and approximation accuracy can be improved by adjusting η,α1It determines The distribution of sigma point is set to a lesser positive number (1e-4≤α1< 1), α3Second scale parameter, is arranged It is 0, α2The dimension that 2, n is system mode vector is set to for Gaussian Profile for state distribution parameter;
2. carrying out Unscented kalman filtering time renewal process, one-step prediction quantity of state is calculated by following equationsWith one-step prediction error covariance matrix P (k, k-1):
ξi(k, k-1)=f (ξi(k-1)) i=0,1 ..., 2n
ζi(k, k-1)=h (ξi(k,k-1))
3. updating observation noise variance matrix R (k) according to step 63) the method;
4. calculating Unscented kalman filtering gain matrix K (k) by following equations:
K (k)=PXZ·PZZ -1
5. carrying out Unscented kalman filtering measurement updaue process, by following equations, state estimation is calculatedAnd estimation Error covariance matrix P (k):
P (k)=P (k, k-1)-K (k) PZZK(k)T
The Unscented kalman filtering state estimationThe vehicle-state vector estimated value of as discrete instants k, thus It can determine the precise position information of vehicle.
The invention has the benefit that
1. the method in the present invention solves under actual traffic environment, UWB location technology is led applied to automobile navigation The highly reliable fusion orientation problem of city vehicle is realized in domain;
2. being classified using fuzzy algorithmic approach to UWB accuracy of observation in the present invention, and then the observation noise of automatic adjusument UWB Variance matrix can effectively exclude the low precision UWB observed quantity as caused by non line of sight or multipath to the shadow of system globe area precision It rings, helps to improve the accuracy and reliability of fusion positioning;
3. realizing the fusion positioning of vehicle in the present invention based on Unscented kalman filtering, traditional nonlinear filtering is overcome For example Extended Kalman filter precision is not high, stability is poor for the method for wave, to target maneuver delay of response the disadvantages of.
Detailed description of the invention
Fig. 1 is flow chart of the invention.
Fig. 2 is the UWB node layout schematic diagram under the right-angled intersection environment of typical urban.
Fig. 3 is the membership function figure of fuzzy logic classifier algorithm.
Fig. 4 is the vehicle location track comparison diagram before and after the algorithm using fuzzy self-adaption.
Fig. 5 is the vehicle location error comparison diagram before and after the algorithm using fuzzy self-adaption.
Specific embodiment
The application is described in further detail with reference to the accompanying drawings and examples.It is understood that this place is retouched The specific embodiment stated is only used for explaining related invention, rather than the restriction to the invention.
With the development and progress of economic society, the vehicle guaranteeding organic quantity in China is quicklyd increase, and urban transportation is faced with day Beneficial severe problem, therefore intelligent transportation system (Intelligent Transport System, ITS) is come into being and gradually The hot spot focused as research.It is fixed all to be unable to do without vehicle for the either realization of ITS or the research in vehicle active safety field Position technology: only under the premise of realizing accurate, reliable location, can effectively command scheduling vehicle, improve urban transportation, Guarantee vehicle safe driving.Therefore, vehicle positioning technology is in the basis and core of the researchs such as ITS and vehicle active safety One of hold.
Automobile navigation location technology common at present includes: dead reckoning (Dead Reckoning, DR), inertial navigation system It unites (Inertial Navigation System, INS), satellite navigation (Global Navigation Satellite System, GNSS) etc..Since the presence of single location technology is respectively insufficient, in order to realize under urban environment it is accurate, reliable, Continuously, it completely positions, mostly uses two kinds and two or more location technologies combines, wherein GNSS/INS fusion positioning is answered With the most extensively: GNSS/INS fusion positioning can solve the orientation problem in the case of satellite-signal short interruptions, and certain The accumulated error that INS is compensated in degree can satisfy the vehicle location demand under opposite free environments, but when satellite-signal is long Between interrupt environment (such as urban canyons, tunnel) in, if GNSS can not work normally, only rely on INS and independently calculate, can still lead Cause biggish position error.
In recent years, the rise and rapid development of the wireless location technology based on UWB are realization vehicle in GNSS constrained environment Under reliable location provide new approaches, currently, UWB location technology be primarily used to indoor positioning field such as personnel, intelligence it is small Vehicle, robot localization etc..The technologies such as since the very bandwidth of UWB is wide, pulse signal penetration power is strong, and multi-path resolved ability is good are excellent Point realizes that the vehicle location under outdoor environment also has certain feasibility based on UWB, in sensor layer, if introducing opposite Reliable UWB information source, then can effectively improve the robustness of fusion positioning system, but under complicated urban environment, UWB signal There is multipath and non-market value vulnerable to vehicle, the trees of road both sides and the interference of building that surrounding travels, If being used to merge positioning for the poor UWB observed quantity of signal quality, the decline of fusion accuracy certainly will be caused, therefore, how to be directed to Relative complex traffic scene (such as right-angled intersection, rotary island crossing), in conjunction with road traffic facility, rational deployment trackside UWB node realizes people's bus or train route information dynamic interaction under bus or train route collaboration, and then chooses suitable method, and automatic adjusument is merged Algorithm is to realize that the Vehicle Fusion based on UWB positions pass urgently to be resolved to the degree of dependence of the UWB observation information of different quality Key problem.
For under urban environment, Vehicular satellite navigation is easily using limited, the problems such as positioning accuracy and not high reliability, On the basis of traditional Vehicular satellite and inertia combined navigation, ultra wide band (Ultra-Wideband, UWB) positioning skill is introduced Art, and classified using fuzzy algorithmic approach to UWB accuracy of observation, thus the observation noise variance matrix of automatic adjusument UWB, in this base On plinth, the fusion positioning of vehicle is realized based on Unscented kalman filtering algorithm.
Referring to FIG. 1, it illustrates the methods according to the present invention to realize that one of vehicle high-precision localization method is implemented The process of example:
Firstly, arranging UWB stationary nodes in road both sides, and UWB stationary nodes are obtained by high-precision difference GPS Position coordinates;Secondly, arranging a UWB mobile node in vehicle roof, using time difference method is reached, the UWB movable joint is obtained Point arrives the distance of each UWB stationary nodes;Then, the rough position coordinate of vehicle is obtained based on ranging localization principle;Then, root According to the distance and bearing corner characteristics of UWB signal, classified based on accuracy of observation of the fuzzy logic classifier algorithm to UWB;And then basis Fuzzy classification result automatic adjusument UWB observation noise battle array;Finally, the fusion positioning of vehicle is realized based on Unscented kalman filtering Obtain the exact position of vehicle.
With reference to the accompanying drawing, thinking of the invention is further described:
Process of the invention is as shown in Figure 1.
A kind of highly reliable fusion and positioning method of vehicle based on fuzzy self-adaption Unscented kalman filtering, which is characterized in that In urban environment, on the basis of Vehicular satellite and inertia combined navigation, UWB positioning is introduced, and combine fuzzy algorithmic approach to UWB Accuracy of observation classification, so that the observation noise variance matrix of automatic adjusument UWB, herein on basis, is filtered based on Unscented kalman Wave algorithm realizes the highly reliable fusion positioning of vehicle, and described method includes following steps:
Step 1) arranges UWB stationary nodes in road both sides, and measures the position coordinates of UWB stationary nodes;
Wherein the number of UWB stationary nodes and installation site will be directed to actual traffic scene, considering cost and fixed Position accuracy requirement, determines reasonable placement scheme, please refers in Fig. 2, and the present invention is what typical urban intersection region proposed A kind of placement scheme: in the road both sides and intersection rotary island region of four different directions of right-angled intersection The heart has laid nine UWB stationary nodes in total, and this layout type can cover entire intersection to greatest extent, realizes Vehicle is in the complete positioning in crossing region, and in the UWB node of rotary island center arrangement, not vulnerable to interference is blocked, it is fixed to facilitate The raising of position precision;
Step 2) arranges a UWB mobile node in vehicle roof, using arrival time difference method (Time Difference Of Arrival, TDOA) the UWB mobile node is obtained to the distance of each UWB stationary nodes, and observation of adjusting the distance is counted Data preprocess;
Step 3), apart from observation, establishes equation in coordinates based on ranging localization principle, using minimum according to pretreated Square law calculates the coordinate of UWB mobile node, obtains the rough position coordinate of vehicle;
Step 4) measures the current course angle of vehicle by vehicle electronics compass, and according to the rough position coordinate of vehicle with And the position coordinates of UWB stationary nodes, calculate azimuth of each UWB stationary nodes relative to vehicle current driving direction;
Step 5) is according to two, distance and bearing angle feature, using fuzzy logic classifier algorithm to each UWB stationary nodes Accuracy of observation classify, one is divided into A, B, C, five grades of D, E, the precision apart from observation of corresponding UWB by height to Low, specific step is as follows for fuzzy Classified Algorithms Applied:
51) three subordinating degree functions of distance feature, and adjust the distance classification, respectively N, M and R are determined, distance is represented " close ", " in ", " remote ";
52) two subordinating degree functions of orientation corner characteristics are determined, and azimuthal is classified, respectively L and S represent orientation " big " and " small " at angle;
53) according to Fuzzy classification rule, 5 accuracy classes A, B, C, D, E of fuzzy logic classifier output are determined, are obscured Classifying rules is as follows:
1. when distance is N, when azimuth is L, accuracy class A;
2. when distance is N, when azimuth is S, accuracy class B;
3. when distance is M, when azimuth is L, accuracy class C;
4. when distance is M, when azimuth is S, accuracy class D;
5. when distance is R, when azimuth is L, accuracy class E;
6. when distance is R, when azimuth is S, accuracy class E;
In the present embodiment, the subordinating degree function of the fuzzy Classified Algorithms Applied please refers to Fig. 3;
Step 6) combines the observation information of vehicle-mounted GNSS, gyroscope and wheel speed sensors output, utilizes Unscented kalman The accurate positioning of filtering algorithm realization vehicle, the specific steps are as follows:
61) firstly, establishing the Unscented kalman filtering state model of vehicle positioning system
System mode vector is
Wherein,For vehicle latitude information, λ is vehicle longitude information, and h is height of car information;veFor the east orientation speed of vehicle Degree, vnFor the north orientation speed of vehicle, vuFor the vertical velocity of vehicle;P is the pitch angle of vehicle, and r is the angle of heel of vehicle, and A is vehicle Course angle;sfodFor wheel speed dilution of precision error, bωFor Gyroscope Random Drift, bGNSSFor GNSS receiver clock deviation, dGNSS For GNSS receiver clock deviation drift rate, bUWBFor the offset error of UWB;
Exterior input vector is U=[vod aod fx fyωz]
Wherein, vodFor the car speed that wheel speed sensors measure, aodFor longitudinal acceleration of the vehicle, fxFor longitudinal direction of car acceleration Spend measured value, fyFor vehicle lateral acceleration measured value, ωzFor the yaw velocity of gyroscope output;
System state equation are as follows:
Wherein, X (k) is the system mode vector of discrete instants k, and X (k-1), U (k-1), W (k-1), T (k-1) are respectively The system mode vector of discrete instants (k-1), external input vector, system noise vector, external input noise vector;RMFor ground Ball meridian radius, RNFor earth prime vertical radius, Δ t is the sampling interval;γwFor the error of wheel speed sensors and time correlation, σwFor its corresponding noise variance;βzFor Gyroscope Random Drift, σzFor its corresponding noise variance;
62) the Unscented kalman filtering observation model of vehicle positioning system then, is established
Systematic observation vector is
Wherein,The pseudo range observed quantity of respectively m GNSS satellite,It is n UWB stationary nodes apart from observation;
Systematic observation equation are as follows:
Z (k)=H (k) X (k)+V (k)
Wherein, Z (k), H (k) are the systematic observation vector and systematic observation matrix of discrete instants k, when V (k) is discrete Carve the observation noise vector of k;
The corresponding observation noise variance matrix of observation noise vector V (k) is R (k):
R (k)=diag [RGNSS RUWB]
Wherein, RGNSSFor the observation noise variance matrix of GNSS, RUWBFor the observation noise variance matrix of UWB;
63) fuzzy Classified Algorithms Applied according to step 5) classifies to the accuracy of observation of n UWB stationary nodes, And then according to the noise variance matrix R of classification results automatic adjusument UWBUWB, final observation noise variance matrix R (k) is obtained, is had Body method is as follows:
1. calculating the observation noise of i-th of UWB stationary nodes
Wherein, R0Noise variance basic value for UWB apart from observation;μ is adjustment factor, true by the result of fuzzy classification Fixed, five kinds of different corresponding adjustment factors of accuracy class A, B, C, D, E are respectively μa, μb, μc, μd, μe
2. calculating the observation noise variance matrix of UWB
3. updating observation noise variance matrix R (k): R (k)=diag [RGNSS RUWB];
64) System State Model according to step 61), observation model described in step 62), and through step 63) Noise variance matrix after automatic adjusument, the time for carrying out Unscented kalman filtering updates and measurement updaue process, obtains vehicle Accurate location information.
1. constructing the sigma point set ξ of Unscented kalman filteringi(k-1), and its weight coefficient is determined
Wherein, η is scale parameter, and approximation accuracy can be improved by adjusting η,α1Determine sigma The distribution of point, is set to a lesser positive number (1e-4≤α1< 1), α3Second scale parameter, is set to 0, α2 The dimension that 2, n is system mode vector is set to for Gaussian Profile for state distribution parameter;
2. carrying out Unscented kalman filtering time renewal process, one-step prediction quantity of state is calculated by following equationsWith one-step prediction error covariance matrix P (k, k-1):
ξi(k, k-1)=f (ξi(k-1)) i=0,1 ..., 2n
ζi(k, k-1)=h (ξi(k,k-1))
3. updating observation noise variance matrix R (k) according to step 63) the method;
4. calculating Unscented kalman filtering gain matrix K (k) by following equations:
K (k)=PXZ·PZZ -1
5. carrying out Unscented kalman filtering measurement updaue process, by following equations, state estimation is calculatedAnd estimation Error covariance matrix P (k):
P (k)=P (k, k-1)-K (k) PZZK(k)T
The Unscented kalman filtering state estimationThe vehicle-state vector estimated value of as discrete instants k, thus It can determine the precise position information of vehicle.
In this embodiment, for the Vehicle Fusion based on fuzzy self-adaption Unscented kalman filtering of inspection institute's proposition The beneficial effect that localization method improves positioning accuracy, has carried out real train test, forward and backward to the algorithm for using fuzzy self-adaption Vehicle location result carried out comparative test.In addition, it is necessary to be pointed out that: herein in test " before fuzzy algorithmic approach processing " Method refer to traditional Unscented kalman filtering fusion method, " using fuzzy algorithmic approach processing after " method refer to this hair Bright middle method.
Fig. 4 is battery of tests comparative result figure, after You Tuzhong's, it is apparent that use fuzzy adaptive algorithm, The track of vehicle location is closer to true value.Fig. 5 is the front and back that one group of data uses fuzzy Classified Algorithms Applied in embodiment Vehicle location error comparison diagram, the position error of You Tuzhong is, it is apparent that use fuzzy Classified Algorithms Applied automatic adjusument After the observation noise matrix of Unscented kalman filtering, Euclidean distance position error reduces, and vehicle location precision significantly improves.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field For art personnel, the invention may be variously modified and varied.All within the spirits and principles of the present invention, made any to repair Change, equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.

Claims (1)

1. a kind of highly reliable fusion and positioning method of vehicle based on fuzzy self-adaption Unscented kalman filtering, which is characterized in that institute The method of stating includes the following steps:
Step 1) arranges UWB stationary nodes in road both sides, and measures the position coordinates of UWB stationary nodes;
Step 2) arranges a UWB mobile node in vehicle roof, obtains the UWB mobile node to respectively using time difference method is reached The distance of a UWB stationary nodes, and observation of adjusting the distance carries out data prediction;
Step 3), apart from observation, establishes equation in coordinates based on ranging localization principle, using least square according to pretreated Method calculates the coordinate of UWB mobile node, obtains the rough position coordinate of vehicle;
Step 4) measures the current course angle of vehicle by vehicle electronics compass, and according to the rough position coordinate of vehicle and The position coordinates of UWB stationary nodes calculate azimuth of each UWB stationary nodes relative to vehicle current driving direction;
Sight of the step 5) according to two, distance and bearing angle feature, using fuzzy logic classifier algorithm to each UWB stationary nodes Precision is surveyed to classify, one is divided into A, B, C, five grades of D, E, the precision apart from observation of corresponding UWB from high to low, mould Pasting sorting algorithm, specific step is as follows:
51) three subordinating degree functions of distance feature, and adjust the distance classification, respectively N, M and R are determined, represent distance " close ", " in ", " remote ";
52) two subordinating degree functions of orientation corner characteristics are determined, and azimuthal is classified, respectively L and S are represented azimuthal " big " and " small ";
53) according to Fuzzy classification rule, 5 accuracy classes A, B, C, D, E of fuzzy logic classifier output, fuzzy classification rule are determined It is then as follows:
<1>is when distance is N, and azimuth is L, accuracy class A;
<2>is when distance is N, and azimuth is S, accuracy class B;
<3>is when distance is M, and azimuth is L, accuracy class C;
<4>is when distance is M, and azimuth is S, accuracy class D;
<5>is when distance is R, and azimuth is L, accuracy class E;
<6>is when distance is R, and azimuth is S, accuracy class E;
Step 6) combines the observation information of vehicle-mounted GNSS, gyroscope and wheel speed sensors output, utilizes Unscented kalman filtering The accurate positioning of algorithm realization vehicle, the specific steps are as follows:
61) firstly, establishing the Unscented kalman filtering state model of vehicle positioning system
System mode vector is
Wherein,For vehicle latitude information, λ is vehicle longitude information, and h is height of car information;veFor the east orientation speed of vehicle, vnFor the north orientation speed of vehicle, vuFor the vertical velocity of vehicle;P is the pitch angle of vehicle, and r is the angle of heel of vehicle, and A is vehicle Course angle;sfodFor wheel speed dilution of precision error, bωFor Gyroscope Random Drift, bGNSSFor GNSS receiver clock deviation, dGNSSFor GNSS receiver clock deviation drift rate, bUWBFor the offset error of UWB;
Exterior input vector is U=[vod aod fx fy ωz]
Wherein, vodFor the car speed that wheel speed sensors measure, aodFor longitudinal acceleration of the vehicle, fxFor longitudinal acceleration of the vehicle survey Magnitude, fyFor vehicle lateral acceleration measured value, ωzFor the yaw velocity of gyroscope output;
System state equation are as follows:
Wherein, X (k) is the system mode vector of discrete instants k, and X (k-1), U (k-1), W (k-1), T (k-1) are respectively discrete The system mode vector at moment (k-1), external input vector, system noise vector, external input noise vector;RMFor earth Noon line radius, RNFor earth prime vertical radius, Δ t is the sampling interval;γwFor the error of wheel speed sensors and time correlation, σwFor Its corresponding noise variance;βzFor Gyroscope Random Drift, σzFor its corresponding noise variance;
62) the Unscented kalman filtering observation model of vehicle positioning system then, is established
Systematic observation vector is
Wherein,The pseudo range observed quantity of respectively m GNSS satellite,For N UWB stationary nodes apart from observation;
Systematic observation equation are as follows:
Z (k)=H (k) X (k)+V (k)
Wherein, Z (k), H (k) are the systematic observation vector and systematic observation matrix of discrete instants k, and V (k) is discrete instants k's Observation noise vector;
The corresponding observation noise variance matrix of observation noise vector V (k) is R (k):
R (k)=diag [RGNSS RUWB]
Wherein, RGNSSFor the observation noise variance matrix of GNSS, RUWBFor the observation noise variance matrix of UWB;
63) fuzzy Classified Algorithms Applied according to step 5) classifies to the accuracy of observation of n UWB stationary nodes, in turn According to the noise variance matrix R of classification results automatic adjusument UWBUWB, final observation noise variance matrix R (k) is obtained, specific side Method is as follows:
<1>calculates the observation noise of i-th of UWB stationary nodes
Wherein, R0Noise variance basic value for UWB apart from observation;μ is adjustment factor, is determined by the result of fuzzy classification, five The different corresponding adjustment factor of accuracy class A, B, C, D, E of kind is respectively μa, μb, μc, μd, μe
<2>calculates the observation noise variance matrix of UWB
<3>updates observation noise variance matrix R (k): R (k)=diag [RGNSS RUWB];
64) System State Model according to step 61), observation model described in step 62), and it is adaptive through step 63) Noise variance matrix after should adjusting, the time for carrying out Unscented kalman filtering updates and measurement updaue process, obtains the essence of vehicle True location information, the specific steps are as follows:
<1>constructs the sigma point set ξ of Unscented kalman filteringi(k-1), and its weight coefficient is determined
Wherein, η is scale parameter, and approximation accuracy can be improved by adjusting η,α1Determine sigma point Distribution, is set to a lesser positive number (1e-4≤α1< 1), α3Second scale parameter, is set to 0, α2For shape State distribution parameter is set to the dimension that 2, n is system mode vector for Gaussian Profile;
<2>carries out Unscented kalman filtering time renewal process, calculates one-step prediction quantity of state by following equations With one-step prediction error covariance matrix P (k, k-1):
ξi(k, k-1)=f (ξi(k-1)) i=0,1 ..., 2n
ζi(k, k-1)=h (ξi(k,k-1))
<3>updates observation noise variance matrix R (k) according to step 63) the method;
<4>calculates Unscented kalman filtering gain matrix K (k) by following equations:
K (k)=PXZ·PZZ -1
<5>carries out Unscented kalman filtering measurement updaue process, by following equations, calculates state estimationIt is missed with estimation Poor covariance matrix P (k):
P (k)=P (k, k-1)-K (k) PZZK(k)T
The Unscented kalman filtering state estimationThe vehicle-state vector estimated value of as discrete instants k, it is possible thereby to Determine the precise position information of vehicle.
CN201910168995.5A 2019-03-06 2019-03-06 Vehicle high-reliability fusion positioning method based on fuzzy self-adaptive unscented Kalman filtering Active CN109946731B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910168995.5A CN109946731B (en) 2019-03-06 2019-03-06 Vehicle high-reliability fusion positioning method based on fuzzy self-adaptive unscented Kalman filtering

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910168995.5A CN109946731B (en) 2019-03-06 2019-03-06 Vehicle high-reliability fusion positioning method based on fuzzy self-adaptive unscented Kalman filtering

Publications (2)

Publication Number Publication Date
CN109946731A true CN109946731A (en) 2019-06-28
CN109946731B CN109946731B (en) 2022-06-10

Family

ID=67009302

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910168995.5A Active CN109946731B (en) 2019-03-06 2019-03-06 Vehicle high-reliability fusion positioning method based on fuzzy self-adaptive unscented Kalman filtering

Country Status (1)

Country Link
CN (1) CN109946731B (en)

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110296706A (en) * 2019-07-11 2019-10-01 北京云迹科技有限公司 Localization method and device, robot for Indoor Robot
CN110940999A (en) * 2019-12-13 2020-03-31 合肥工业大学 Self-adaptive unscented Kalman filtering method based on error model
CN111090113A (en) * 2019-12-31 2020-05-01 南京泰通科技股份有限公司 High-precision train positioning terminal and positioning method based on Beidou, UWB and inertial navigation
CN111551976A (en) * 2020-05-20 2020-08-18 四川万网鑫成信息科技有限公司 Method for automatically completing abnormal positioning by combining various data
CN111811503A (en) * 2020-07-15 2020-10-23 桂林电子科技大学 Unscented Kalman filtering fusion positioning method based on ultra wide band and two-dimensional code
CN111915921A (en) * 2020-08-11 2020-11-10 大陆投资(中国)有限公司 Lane-level precision positioning method using roadside equipment
CN111954153A (en) * 2020-08-12 2020-11-17 南京工程学院 Intelligent taxi calling navigation positioning method based on UWB positioning
CN112406861A (en) * 2019-08-19 2021-02-26 通用汽车环球科技运作有限责任公司 Method and device for Kalman filter parameter selection by using map data
CN112533149A (en) * 2020-11-27 2021-03-19 桂林理工大学 Moving target positioning algorithm based on UWB mobile node
CN112577512A (en) * 2019-09-27 2021-03-30 北京初速度科技有限公司 State quantity error determination method based on wheel speed fusion and vehicle-mounted terminal
CN112577513A (en) * 2019-09-27 2021-03-30 北京初速度科技有限公司 State quantity error determination method and vehicle-mounted terminal
CN113029137A (en) * 2021-04-01 2021-06-25 清华大学 Multi-source information self-adaptive fusion positioning method and system
CN114440876A (en) * 2022-01-21 2022-05-06 北京自动化控制设备研究所 Underground heading machine positioning and guiding method and underground heading machine positioning and guiding system
CN114527487A (en) * 2022-02-24 2022-05-24 东南大学 Active characteristic assisted self-adaptive fusion positioning method in point cloud sparse scene
CN114779307A (en) * 2022-06-17 2022-07-22 武汉理工大学 Port area-oriented UWB/INS/GNSS seamless positioning method
US11718317B2 (en) 2020-10-16 2023-08-08 Ford Global Technologies, Llc Vehicle location correction using roadside devices
CN117148406A (en) * 2023-10-30 2023-12-01 山东大学 Indoor and outdoor seamless elastic fusion positioning method, system, medium and equipment
CN117610601A (en) * 2023-11-08 2024-02-27 深圳市百千诚国际物流有限公司 Logistics equipment intelligent positioning method and system based on Internet of things

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102928816A (en) * 2012-11-07 2013-02-13 东南大学 High-reliably integrated positioning method for vehicles in tunnel environment
CN103323005A (en) * 2013-03-06 2013-09-25 郭雷 Multi-objective optimization anti-interference filtering method for SINS/GPS/polarized light combined navigation system
CN104154913A (en) * 2014-06-28 2014-11-19 南京诺导电子科技有限公司 Autonomous indoor positioning algorithm in all-attitude environment
CN105807301A (en) * 2016-03-03 2016-07-27 东南大学 Enhanced digital map based vehicle optimization oriented satellite selection positioning method
CN108363978A (en) * 2018-02-12 2018-08-03 华南理工大学 Using the emotion perception method based on body language of deep learning and UKF
CN108445480A (en) * 2018-02-02 2018-08-24 重庆邮电大学 Mobile platform based on laser radar adaptively extends Target Tracking System and method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102928816A (en) * 2012-11-07 2013-02-13 东南大学 High-reliably integrated positioning method for vehicles in tunnel environment
CN103323005A (en) * 2013-03-06 2013-09-25 郭雷 Multi-objective optimization anti-interference filtering method for SINS/GPS/polarized light combined navigation system
CN104154913A (en) * 2014-06-28 2014-11-19 南京诺导电子科技有限公司 Autonomous indoor positioning algorithm in all-attitude environment
CN105807301A (en) * 2016-03-03 2016-07-27 东南大学 Enhanced digital map based vehicle optimization oriented satellite selection positioning method
CN108445480A (en) * 2018-02-02 2018-08-24 重庆邮电大学 Mobile platform based on laser radar adaptively extends Target Tracking System and method
CN108363978A (en) * 2018-02-12 2018-08-03 华南理工大学 Using the emotion perception method based on body language of deep learning and UKF

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王宇: "基于车内总线的车辆双模卫星融合定位技术研究", 《中国优秀硕士学位论文全文数据库工程科技II辑》 *

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110296706A (en) * 2019-07-11 2019-10-01 北京云迹科技有限公司 Localization method and device, robot for Indoor Robot
CN112406861A (en) * 2019-08-19 2021-02-26 通用汽车环球科技运作有限责任公司 Method and device for Kalman filter parameter selection by using map data
CN112406861B (en) * 2019-08-19 2024-03-22 通用汽车环球科技运作有限责任公司 Method and device for carrying out Kalman filter parameter selection by using map data
CN112577512A (en) * 2019-09-27 2021-03-30 北京初速度科技有限公司 State quantity error determination method based on wheel speed fusion and vehicle-mounted terminal
CN112577513A (en) * 2019-09-27 2021-03-30 北京初速度科技有限公司 State quantity error determination method and vehicle-mounted terminal
CN110940999A (en) * 2019-12-13 2020-03-31 合肥工业大学 Self-adaptive unscented Kalman filtering method based on error model
CN111090113A (en) * 2019-12-31 2020-05-01 南京泰通科技股份有限公司 High-precision train positioning terminal and positioning method based on Beidou, UWB and inertial navigation
CN111551976A (en) * 2020-05-20 2020-08-18 四川万网鑫成信息科技有限公司 Method for automatically completing abnormal positioning by combining various data
CN111811503A (en) * 2020-07-15 2020-10-23 桂林电子科技大学 Unscented Kalman filtering fusion positioning method based on ultra wide band and two-dimensional code
CN111915921A (en) * 2020-08-11 2020-11-10 大陆投资(中国)有限公司 Lane-level precision positioning method using roadside equipment
CN111954153B (en) * 2020-08-12 2022-11-29 南京工程学院 Intelligent taxi calling navigation positioning method based on UWB positioning
CN111954153A (en) * 2020-08-12 2020-11-17 南京工程学院 Intelligent taxi calling navigation positioning method based on UWB positioning
US11718317B2 (en) 2020-10-16 2023-08-08 Ford Global Technologies, Llc Vehicle location correction using roadside devices
CN112533149B (en) * 2020-11-27 2022-06-07 桂林理工大学 Moving target positioning algorithm based on UWB mobile node
CN112533149A (en) * 2020-11-27 2021-03-19 桂林理工大学 Moving target positioning algorithm based on UWB mobile node
CN113029137B (en) * 2021-04-01 2023-02-28 清华大学 Multi-source information self-adaptive fusion positioning method and system
CN113029137A (en) * 2021-04-01 2021-06-25 清华大学 Multi-source information self-adaptive fusion positioning method and system
CN114440876A (en) * 2022-01-21 2022-05-06 北京自动化控制设备研究所 Underground heading machine positioning and guiding method and underground heading machine positioning and guiding system
CN114440876B (en) * 2022-01-21 2024-04-02 北京自动化控制设备研究所 Positioning and guiding method and system for underground heading machine
CN114527487A (en) * 2022-02-24 2022-05-24 东南大学 Active characteristic assisted self-adaptive fusion positioning method in point cloud sparse scene
CN114779307A (en) * 2022-06-17 2022-07-22 武汉理工大学 Port area-oriented UWB/INS/GNSS seamless positioning method
CN117148406A (en) * 2023-10-30 2023-12-01 山东大学 Indoor and outdoor seamless elastic fusion positioning method, system, medium and equipment
CN117148406B (en) * 2023-10-30 2024-01-30 山东大学 Indoor and outdoor seamless elastic fusion positioning method, system, medium and equipment
CN117610601A (en) * 2023-11-08 2024-02-27 深圳市百千诚国际物流有限公司 Logistics equipment intelligent positioning method and system based on Internet of things

Also Published As

Publication number Publication date
CN109946731B (en) 2022-06-10

Similar Documents

Publication Publication Date Title
CN109946731A (en) A kind of highly reliable fusion and positioning method of vehicle based on fuzzy self-adaption Unscented kalman filtering
CN109946730B (en) Ultra-wideband-based high-reliability fusion positioning method for vehicles under cooperation of vehicle and road
CN104076382B (en) A kind of vehicle seamless positioning method based on Multi-source Information Fusion
CN106840179B (en) Intelligent vehicle positioning method based on multi-sensor information fusion
CN103760904B (en) A kind of voice broadcast type intelligent vehicle path planning apparatus and implementation
CN105792135B (en) A kind of method and device in positioning vehicle place lane
CN108089180A (en) Based on UWB sensors as back indicator to the localization method of GPS and inertial navigation system the suspension type rail vehicle corrected
US11874366B2 (en) High-precision vehicle positioning method based on ultra-wideband in intelligent vehicle infrastructure cooperative systems
CN113706612B (en) Underground coal mine vehicle positioning method fusing UWB and monocular vision SLAM
AU2021106247A4 (en) Vehicle fusion positioning method based on V2X and laser point cloud registration for advanced automatic driving
CN108627864A (en) Localization method and system, pilotless automobile system based on automobile key
CN106469505A (en) A kind of floating wheel paths method for correcting error and device
CN111221020A (en) Indoor and outdoor positioning method, device and system
He et al. Research on GNSS INS & GNSS/INS integrated navigation method for autonomous vehicles: A survey
Cho et al. Enhancing GNSS performance and detection of road crossing in urban area using deep learning
Zhang et al. Using UWB aided GNSS/INS integrated navigation to bridge GNSS outages based on optimal anchor distribution strategy
CN115711618B (en) Mining area high-precision positioning method based on road side guidance
CN204115737U (en) A kind of indoor positioning device based on inertial guidance and radio-frequency (RF) identification
CN115061176B (en) Vehicle GPS enhanced positioning method based on V2V instantaneous data exchange
Wilfinger et al. Indoor position determination using location fingerprinting and vehicle sensor data
CN108696820A (en) Indoor orientation method based on iBeacon ranging technologies
Zhong et al. Collaborative navigation in V2I network based on Chan-Taylor joint iterative algorithm
Ke et al. Integrated positioning method for intelligent vehicle based on GPS and UWB
Sharma et al. Real time location tracking map matching approaches for road navigation applications
CN113156479B (en) Combined positioning method for outdoor country road scene

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