Field of the Invention
This invention relates generally to a system and method for estimating vehicle speed and, more particularly, to a system and method for a secondary or redundant estimation of vehicle speed in the event that a primary speed estimation module fails, where a secondary speed estimation module employs sensor inputs from a 6 degree-of-freedom (DOF) inertial measurement unit (IMU), a GPS receiver and range sensors.
Discussion of the Related Art
Vehicles are becoming more autonomous or cognitive with the goal being a completely autonomously driven vehicle, i.e., vehicles that are able to provide driving control with minimal or no driver intervention. Adaptive cruise control systems have been available for a number of years where not only does the system maintain a set speed, but also will automatically slow the vehicle down in the event that a slower moving vehicle is detected in front of the subject vehicle. Vehicle control systems currently exist that include autonomous parking where the vehicle will automatically provide the steering control for parking the vehicle. Also, control systems exist that may intervene if the driver makes harsh steering changes that may affect vehicle stability and lane centering capabilities, where the vehicle system attempts to maintain the vehicle near the center of the travel lane. Future vehicles will likely employ autonomous systems for lane changing, passing, turns away from traffic, turns into traffic, merging into traffic, passing through or turning at intersections, etc.
Various active safety control systems, driver assist systems and autonomous driving operations on vehicles, such as electronic stability control (ECS), adaptive cruise control (ACC), lane keeping (LK), lane changing (LC), etc., require highly robust and precise modules for estimating various vehicle dynamics. Such modules are necessary to provide knowledge of the vehicle position and velocity to control the vehicle.
Active safety control for the systems discussed above, and others, rely on accurate vehicle speed estimation for proper performance. Currently, these types of proposed systems rely on wheel speed sensors and other vehicle kinematic inputs to provide the vehicle speed estimation. However, sometimes sensors and control modules that determine vehicle speed fail or operate incorrectly, where loss of vehicle speed could be serious. Certain automotive safety requirements, such as ASIL-D, require redundant vehicle speed estimation processes in the event of failure of the primary estimation processor. For example, for those systems that require active control, the control systems are required to provide accurate speed estimation for five seconds after the failure event so as to give the driver time to take control of the vehicle. It would be desirable to provide such redundant speed estimation processes using existing hardware on the vehicle to reduce vehicle cost.
The present disclosure describes a system and method for providing redundant vehicle speed estimation. The method includes providing sensor output signals from a plurality of primary sensors and providing inertial measurement signals from an inertial measurement unit. The method also includes estimating the vehicle speed in a primary module using the primary sensor signals, and buffering the estimated vehicle speed values from the primary module for a predetermined period of time. The method further includes determining that one or more of the primary sensors or the primary module has failed, and if so, estimating the vehicle speed in a secondary module using the buffered vehicle speed values and the inertial measurement signals. The method can use GPS signal data and/or velocity data provided by range sensors, such as radar, lidar and vision systems, from static objects to improve the estimated vehicle speed in the secondary module if they are available.
Additional features of the present invention will become apparent from the following description and appended claims, taken in conjunction with the accompanying drawings.
The following discussion of the embodiments of the invention directed to a secondary technique for estimating vehicle speed in the event of primary vehicle speed estimation failure is merely exemplary in nature, and is in no way intended to limit the invention or its applications or uses. For example, as discussed, the system and method has particular application for estimating vehicle speed. However, as will be appreciated by those skilled in the art, the system and method may have application for other mobile platforms, such as on trains, machines, tractors, boats, recreation vehicles, etc.
As will be discussed in detail below, the present invention proposes a secondary system and method for estimating vehicle speed in the event that the primary system for estimating vehicle speed fails. In one embodiment, the secondary speed estimation system has application for providing the vehicle speed estimation for a certain period of time, such as five seconds, after the primary technique has failed as required by automotive standards for vehicle autonomous or semi-autonomous driving control systems to allow time for the vehicle driver to take control of the vehicle. The secondary speed estimation technique employs data from existing vehicle sensors that are not used to determine the vehicle speed in the primary speed estimation system. In one embodiment, the secondary speed estimation system employs signals from a 6-DOF IMU, a GPS receiver and range sensors, such as radar, lidar, vision systems, etc., that can be fused together to determine vehicle speed. During the time that the primary speed estimation system is operational, the calculated vehicle speed is buffered for a number of sample points, where the buffered speed information is continually updated for each sample period. The 6-DOF IMU sensor data is used to extrapolate vehicle speed for some period of time after failure of the primary speed estimation system by using the buffered vehicle speed data calculated by the primary system. This data from the 6-DOF IMU can be used in combination with the GPS signals if they are available. Further, velocity measurement data from stationary objects around the vehicle can be used in combination with one or both of the extrapolated 6-DOF IUM data and the GPS signals.
The secondary module 34 also receives the speed signals from the primary function processor 38 and buffers those speed signals for at least some rolling period of time, such as five seconds. If a fault in the primary module 32 is detected, then the secondary module 34 reads the buffered data and extrapolates that data using angular rate and acceleration signals from the 6-DOF IMU 50 to provide the vehicle speed estimation. The 6-DOF IMU 50 is a well known sensor available on certain vehicles that provides six rate of change measurements, particularly, angular rate measurements to of vehicle roll, pitch and yaw, and acceleration measurements f of longitudinal acceleration, lateral acceleration and up/down acceleration. The extrapolation of the vehicle speed using the stored speed values and the 6-DOF IMU data in the secondary module 34 can be corrected or augmented using the GPS data from the GPS receiver 52 and/or the data from the range sensors 54 that detect stationary objects if available, where the range rate of change of the distance between the vehicle and the detected stationary object can be used to help determine the vehicle speed. If the primary module 32 recovers from the failure or fault, operation of the system 30 transitions back to the primary module 32, where the state variables of the primary module 32, such as vehicle longitudinal speed, vehicle lateral speed, roll angle and pitch angle, are initialized with the latest state variables from the secondary module 34 over the communications bus 48.
The following discussion describes how the rate change data from the 6-DOF IMU 50 can be used to predict the vehicle speed for some period of time after the primary module 32 has failed using the buffered vehicle speed data. The vehicle kinematics of the ordinary differential equation (ODE) for vehicle speed in the vehicle body frame b is defined as:
{dot over (v)}b=−Ωvb+Rabga+fb, (1)
where vb is the vehicle velocity at the center of the IMU 50 in the vehicle body frame b, which corresponds to vehicle longitudinal, lateral and up/down speed, i.e., vb=[vlon,vlat,vup]T, where it is assumed that the wheels of the vehicle always touch the ground, particularly a “virtual” measurement vup=0 in Kalman filtering; Ω is a skew symmetric matrix of the angular velocity ω in the b-frame; ga is the Earth gravitational vector in the inertia frame a, specifically (0,0,−9.80665)T m/s2; Rab is the rotation matrix from the b-frame to the inertia frame a; fb is the vehicle acceleration measured by accelerometers in the b-frame (IMU data), and ω is the angular rates measured by gyros in the b-frame (IMU data), and where the skew matrix Ω is:
and the angular rates ω are:
ω=({dot over (φ)},{dot over (θ)},{dot over (ψ)})T, (3)
where φ is the roll angle, θ is the pitch angle and ψ is the yaw angle.
If the value va is the velocity of the vehicle's center of gravity (CG) in the inertial frame (a-frame), by Newton's laws:
{dot over (v)}a=ga+fa (4)
since:
va=Rbavb, (5)
{dot over (v)}a=Rba{dot over (v)}b=Rba{dot over (v)}b+RbaΩvb=ga+fa, (6)
Note that the property of the rotation matrix Rab and the acceleration vector observed in the vehicle body frame. Multiplying equation (6) by the rotation matrix Rab derives the ODE of the vehicle speed in vehicle body frame b.
v0b=[{circumflex over (v)}lon,{circumflex over (v)}lat,0]T (7)
The vehicle velocity in the body frame is integrated at box 64 to obtain the vehicle speed value vb. The speed value vb is used to obtain the up-speed error −vup at box 66, which is multiplied by a Kalman filtering gain K provided at box 68 in a multiplier 70 to obtain −Kvup. The filtered up-speed error −Kvup is then subtracted from the vehicle speed value vb in a subtractor 72 to obtain the estimated longitudinal and lateral vehicle speed signal on line 74 as the output from the secondary function processor 42.
A discussion of how the speed value {dot over (v)}b is obtained is provided as follows. The pitch, roll and yaw angle rates ω from the 6-DOF IMU 50 are provided at box 76, and are used to calculate the skew matrix Ω at box 78. The skew matrix Ω is multiplied by the vehicle speed value vb in multiplier 80 and is inverted by inverter 82 to obtain the value −Ωvb as the first component of equation (1).
The prediction of the vehicle speed by the secondary function processor 42 in the secondary module 34 requires an initial value of the vehicle attitude, i.e., roll, pitch and yaw, which determines the rotation matrix Rab from the inertial coordinate system (a-frame) to the vehicle coordinate system (b-frame). In this example, the value Rabga is invariant to azimuth rotation, and thus only the roll and pitch angles of the vehicle are needed, where the yaw angle is set to zero. For a stationary vehicle, the roll and pitch attitude of the vehicle can be computed from two-dimensional angles between the gravitational vector and the IMU measured acceleration. However, for a moving vehicle, the attitude for roll and pitch can be determined using the vehicle speed vector value vb from the primary module 32. The skew matrix Ω from the box 78 is inverted by an inverter 58 to obtain the inverse of the skew matrix −Ω, which is multiplied by an estimate of the rotation matrix Rab provided on line 56 in a multiplier 98 to obtain the value −ΩRab. The last buffered rotation matrix Rab provided at box 84 before the detected fault of the primary module 32 and the value −ΩRab from the multiplier 98 are provided to an integral operator 86 to generate a new rotation matrix {circumflex over (R)}ab.
The space formed by rotation matrices is non-Euclidean, and thus special care is needed when handling vehicle attitude estimations. For example, the value {circumflex over (R)}ab from the integral operator 86 may not be a rotation matrix. A singular value decomposition (SVD) projection is needed to provide normalization of the value {circumflex over (R)}ab to identify it as a rotation matrix. SVD matrices U, S and V are computed such that:
USVT={circumflex over (R)}ab, (8)
The normalized rotation matrix Rab is obtained at box 88 by computing:
Rab=USVT. (10)
The vehicle attitude including the roll φ, pitch θ, and yaw ψ angles is obtained from the rotation matrix Rab. The vehicle coordinate system (b-frame) is attached to the body of the vehicle, and the a-frame is the inertia coordinate system.
The rotation matrix Rba can be expressed by three Euler angles as:
θ=(φ,θ,ψ)T. (11)
The rotation matrix Rab from the a-frame to the b-frame is provided as:
where cθ=cos(θ) and sθ=sin(θ).
The derivative of the rotation matrix Rba from the b-frame to the a-frame is determined as follows. Considering infinitesimal angles δθ=(δφ, δθ, δψ)T for the roll φ, pitch θ and yaw ψ motion of the vehicle, the corresponding rotation matrix R can be approximated by Rba≈I3+ΔΘ, where ΔΘ=[δθx] is the skew symmetric matrix representation of the rotation angles θba as:
where x represents a skew symmetric matrix operator, i.e., axb≡[ax]b.
At time t, a vector p in the b-frame can be expressed in the a-frame as q(t)=Rba(t)p. Now consider at time t+Δt that:
q(t+Δt)=Rba(t+Δt)p=(I3+ΔΘ)Rba(t)p. (14)
The time derivative of the rotation matrix Rba is defined as:
is the skew symmetric matrix angular rate ω=({dot over (φ)},{dot over (θ)},{dot over (ψ)})T, i.e., Ω=[ωx]. Note that equation (15) is equivalent to the vehicle attitude:
{dot over (θ)}(t)=ω. (16)
The rotation matrix Rab from the box 88 and the gravitational vector ga from box 90 are multiplied by multiplier 92 to obtain the second component of equation (1), specifically Rabga. The three acceleration values fb from the 6-DOF IMU 50 are provided at box 94 for the third component of equation (1) and all of the components −Ωvb, Rab, ga and fb are added together by adder 96 to obtain the velocity value {dot over (v)}b that is sent to the integral box 64.
Rabga=β, (17)
where:
β={dot over (v)}b+Ωvb−fb. (18)
The SVD matrices U, S and V are computed as:
USVT=β(ga)T, (19)
where:
The estimate for the rotation matrix Rab is:
Re=UCVT. (21)
The vehicle attitude θ for the roll angle φ and the pitch angle θ is computed as:
φ=−sin−1(R(1,3)), (22)
As mentioned above, if the GPS data from the GPS receiver 52 is available it can be used to improve the predicted vehicle speed in combination with the data provided by the IMU 50 because the prediction performance degrades as errors are accumulated. The GPS signal is used as the measurement of the vehicle's longitudinal speed as vGPS=√{square root over (vlon2+vlat2)}+ε, where vlon is the vehicle longitudinal speed, vlat is the vehicle lateral speed, and E is the zero-mean Gaussian distributed error. A Kalman-Bucy filter is applied to correct vehicle velocity estimation when the GPS corrections occur when new GPS data arrives and signal reception is above a predetermined signal-to-noise ratio.
Also, as mentioned above, velocity data from range sensors can be used to improve the vehicle speed estimator. As an exemplar embodiment, the correction using the radar data is discussed below with reference to
where εi and εj are zero-mean Gaussian distributed error terms introduced by radar Doppler measurements. Therefore, a Kalman-Bucy filter can be designed to correct vehicle speed prediction based on the measurement equations.
As will be well understood by those skilled in the art, the several and various steps and processes discussed herein to describe the invention may be referring to operations performed by a computer, a processor or other electronic calculating device that manipulate and/or transform data using electrical phenomenon. Those computers and electronic devices may employ various volatile and/or non-volatile memories including non-transitory computer-readable medium with an executable program stored thereon including various code or executable instructions able to be performed by the computer or processor, where the memory and/or computer-readable medium may include all forms and types of memory and other computer-readable media.
The foregoing discussion disclosed and describes merely exemplary embodiments of the present invention. One skilled in the art will readily recognize from such discussion and from the accompanying drawings and claims that various changes, modifications and variations can be made therein without departing from the spirit and scope of the invention as defined in the following claims.
Number | Name | Date | Kind |
---|---|---|---|
6718248 | Lu | Apr 2004 | B2 |
8855848 | Zeng | Oct 2014 | B2 |
9228836 | Girod et al. | Jan 2016 | B2 |
20030236606 | Lu | Dec 2003 | A1 |
20100017128 | Zeng | Jan 2010 | A1 |
20140278206 | Girod | Sep 2014 | A1 |
Number | Date | Country |
---|---|---|
101320089 | Dec 2008 | CN |
102008026397 | Jan 2009 | DE |
WO 2014145409 | Sep 2014 | WO |
Entry |
---|
Plane-based calibration for multibeam sonar mounting angles; Chau-Chang Wang; Kuang-Che Yang; Sen-Wen Shyue; Liu-Huang Chon-Lon; Oceans 2003. Proceedings; Year: 2003, vol. 2; p. 987 vol. 2, DOI: 10.1109/OCEANS.2003.178465. |
A strategy for selecting optimal vergence pose and calibration patterns for stereo-vision calibration; K S Chidanand Kumar 2015 IEEE International Conference on Computer Graphics, Vision and Information Security (CGVIS); Year: 2015 pp. 44-48, DOI: 10.1109/CGVIS.2015.7449890. |
Estimating desired yaw rate and control strategy analysis on developed air ESC system for performance evaluation; Seunghwan Chung; Hyeongcheol Lee; Control, Automation and Systems (ICCAS), 2015 15th International Conference on; Year: 2015 pp. 2005-2010, DOI: 10.1109/ICCAS.2015.7364697. |
Design and simulation of an innovative micro-magnetic module with dual-function of actuator and sensor; Nan-Chyuan Tsai; Bing-Hong Liou; Chih-Che Lin; 2009 IEEE/ASME International Conference on Advanced Intelligent Mechatronics; Year: 2009 pp. 1903-1908, DOI: 10.1109/AIM.2009.5229767. |
Real time vision for intelligent vehicles; D. M. Gavrila; U. Franke; C. Wohler; S. Gorzig; IEEE Instrumentation & Measurement Magazine; Year: 2001, vol. 4, Issue: 2; pp. 22-27, DOI: 10.1109/5289.930982. |
Experimental investigation of axisymmetric turbulent boundary layers; Kimberly Cipolla; William Keith; Alia Foley; Oceans'11 MTS/IEEE KONA; Year: 2011; pp. 1-6. |
SPIV measurements of axisymmetric turbulent boundary layers; Kimberly Cipolla; William Keith; Alia Foley; Oceans 2010 MTS/IEEE Seattle; Year: 2010; pp. 1-7, DOI: 10.1109/OCEANS.2010.5663848. |
Song, Xiaofeng F. et al. U.S. Appl. No. 14/614,164, filed Feb. 4, 2015, entitled “Fail Operational Vehicle Speed Estimation Through Data Fusion of 6-DOF IMU, GPS, and Radar”. |
Number | Date | Country | |
---|---|---|---|
20160299234 A1 | Oct 2016 | US |