# Patent

## Title
Target curvature estimation considering vehicle dynamics

## Abstract

This document describes target curvature estimation considering vehicle dynamics constraints. Some vehicle dynamics models use several kinematic states to predict the paths of targets identified by an object-tracking system. Rather than rely on predetermined industry measurements as the kinematic inputs to a model, an example system instead estimates some of these states using near real-time data output from a tracking filter. This can include estimating position and velocity states, as well as curvature and centripetal acceleration states. The example system is particularly suited at accurately predicting the curvature and centripetal acceleration states, which greatly improves the accuracy of the object-tracking system, and further improves safety in avoiding a target collision.

## Background
BACKGROUND
Advanced Driver-Assistance Systems (ADAS) have been developed to assist an operator of a vehicle with various driving tasks, including tracking and predicting the paths of target vehicles. Estimating the kinematic states (e.g., position, velocity, acceleration, curvature) of the target vehicle may be necessary to accurately predict a target path and avoid a collision. Without vehicle-to-vehicle communication, these kinematic states may be perceived by an onboard sensor, such as a radar. A radar usually reports multiple detections for a single target vehicle within its field-of-view (FOV). A tracking function is used to fuse multiple detections on the target and estimate other required kinematic states based on these detections. To ensure safety and prevent collisions with the targets, accurate estimates for these states are desirable.
SUMMARY
This document describes target curvature estimation considering vehicle dynamics Many vehicle dynamics models use several kinematic states to predict the paths of targets identified by an object tracking system. Some of these states are estimated using data from a tracking filter. Many of the states are accurately estimated, including position and velocity. One such state is a curvature state, which tends to be less accurate due to its nonlinearity, but it can be improved using the methods and systems described herein.
In one example, a method includes determining, by a processor of a host vehicle, measurement states of a target, the measurement states being derived from radar detections received by an object tracking system of the host vehicle. The method further includes determining, based at least in part on the measurement states, one or more kinematic states of a motion model defined for a tracking filter employed by the object tracking system to track the target. The method further includes predicting, using the one or more kinematic states as inputs to the motion model of the tracking filter, a path of the target. The method further includes controlling, based on the predicted path of the target, the host vehicle to take evasive action and avoid a collision with the target.
In one example, a system includes a processor configured to perform this and other methods. In another example, a system is described including means for performing this and other methods. This document also describes a computer-readable storage medium comprising instructions that, when executed, configure a processor to perform the above-summarized method and other methods set forth herein, in addition to describing systems configured to perform the above-summarized method and the other methods set forth herein.
This summary introduces simplified concepts of target curvature estimation considering vehicle dynamics This summary is not intended to identify essential features of the claimed subject matter, nor is it intended for use in determining the scope of the claimed subject matter. That is, one problem solved by the described techniques is a problem of estimation errors exceeding physical limitations. Therefore, although primarily described in the context of improving radar tracking functionality of an automotive vehicle, improving target curvature estimations can be applied to other applications where accurate path prediction is desired.



BRIEF DESCRIPTION OF THE DRAWINGS
The details of one or more aspects of target curvature estimation considering vehicle dynamics are described in this document with reference to the following figures. The same numbers are often used throughout the drawings to reference like features and components.
FIG. 1 illustrates an example environment in which a vehicle is configured to perform target curvature estimation considering vehicle dynamics, in accordance with techniques of this disclosure.
FIG. 2 illustrates an example of an automotive system configured to perform target curvature estimation considering vehicle dynamics, in accordance with techniques of this disclosure.
FIG. 3 illustrates an example of a vehicle dynamics model that is employed to determine target curvature estimation considering vehicle dynamics, in accordance with techniques of this disclosure.
FIG. 4 illustrates a graph representing the relationship, using the Pacejka Magic Formula, between lateral tire force and tire side-slip angle, in accordance with techniques of this disclosure.
FIG. 5 illustrates the coefficients of lateral tire force versus tire side-slip angle used in the Pacejka Magic Formula, in accordance with techniques of this disclosure.
FIG. 6 illustrates a graph representing a tire friction forces ellipse vehicle dynamics model, in accordance with techniques of this disclosure.
FIG. 7 illustrates ultimate upper-bound values determined from the vehicle dynamics models employed, in accordance with techniques of this disclosure.
FIG. 8 illustrates some kinematic states that may be included in a state vector of a motion model.
FIG. 9 illustrates a flow chart of the Constrained Constant Turn and Constant Acceleration (CCTCA) tracking filter used to determine target curvature estimation considering vehicle dynamics, in accordance with techniques of this disclosure.
FIG. 10 illustrates an example method for target curvature estimation considering vehicle dynamics, in accordance with techniques of this disclosure.
FIGS. 11-1 to 11-3 illustrate example scenarios of target curvature estimation considering vehicle dynamics, in accordance with techniques of this disclosure.



DETAILED DESCRIPTION
Overview
In some advanced driver-assistance systems (ADAS), estimating a target vehicle's kinematic states may be necessary to predict a target vehicle's path. Without vehicle-to-vehicle communication, the target's kinematic states may be perceived by an onboard sensor, a radar being one sensor capable of kinematic-state estimation. The measurement states, derived from radar detections, may include range, range rate, azimuth angle, and elevation angle. With a sufficient resolution, these radars can report multiple detections on a target vehicle inside its field-of-view (FOV). A radar tracking function can fuse multiple radar detections on a single target and report the required kinematic state.
The radar tracking function may include one or more motion models used to characterize the target, a data association between raw sensor detections and targets, and a tracking filter. This techniques, systems, and methods described in this document focus on improving the accuracy of these target kinematic-state estimations, which may be used by the motion model defined for the tracking filter. It is assumed, for simplicity of this disclosure, that the data association between raw sensor detections and targets has been achieved. In some aspects, the motion model is the Constant Turn and Constant Tangential Acceleration (CTCA) motion model and represents the dynamics of a two-axle vehicle. Due to the apparent nonlinearity between the radar measurement states and the CTCA motion state, an Extended Kalman Filter (EKF) may be selected as the tracking filter.
In general, tracking filters estimate the position and velocity states of the target relatively accurately for most driving scenarios. However, the curvature and lateral acceleration reported by the tracking filter may have larger relative errors. For example, a vehicle with a longitudinal velocity of 36.7 m/s and a curvature estimate of 0.0076 rad/m may report an implausible lateral acceleration of 10.2115 m/s2. This may be explained by a higher nonlinearity between the curvature and the radar measurement states. In this example, the target vehicle-heading angle is quite small, and thus the implausible lateral acceleration estimated by the tracker primarily originated from an implausible centripetal acceleration. In the CTCA model, centripetal acceleration, ac, is expressed as Eq. (1):

ac=cs2   (1)

where c is the curvature, and s is the velocity. By performing a simple error analysis, ∂ac/∂c=s2, the curvature estimation is determined to be a major contributor to the centripetal-acceleration error. Since the centripetal-acceleration state is necessary for target-path prediction, a more-accurate centripetal-acceleration state estimation is desirable.

The techniques disclosed in this document improve the target curvature estimation considering vehicle dynamics. A primary and a secondary upper bound of a centripetal acceleration state are determined from two different vehicle dynamics models for complementary purposes. The primary and secondary upper bounds are compared and the lesser value of the two upper bounds is chosen as an ultimate upper bound. Choosing the ultimate upper bound in this manner ensures that more corner cases are covered. A constrained EKF based on a Constrained Constant Turn and Constant Acceleration (CCTCA) model (e.g., the model using the upper bound) runs on top of the unconstrained CTCA-EKF filter. The CCTCA-EKF checks whether a target state triggers the centripetal acceleration constraints, and if the target state violates the constraint, the CCTCA-EKF projects the unconstrained estimation to the constrained surface.
The highly nonlinear vehicle dynamics models are approximated by using a tire model developed by Pacejka in 1987, which is commonly referred to in the automotive industry as the “Pacejka Magic Formula” and using one or more tire friction force ellipse models. The extreme values of tire side-slip angles and other parameters that determine the upper bound are derived based on reasonable assumptions and publications related to tire and vehicle dynamics. In this way, an example system, through applying these extreme values and approximated vehicle dynamics models, simplifies implementation and reduces computation time.
Example Environment
FIG. 1 illustrates an example environment 100 in which a vehicle 102 is configured for target curvature estimation considering vehicle dynamics, in accordance with techniques of this disclosure. In the depicted environment 100, the vehicle 102 is equipped with an ADAS 104. The ADAS 104 can include sensors (e.g., radar sensor 106) and a radar tracking module 108. Although depicted as a car, the vehicle 102 can represent other types of vehicles and machinery (e.g., a motorcycle, a bus, a tractor, a semi-trailer truck, watercraft, aircraft, or other heavy equipment), including manned and unmanned systems that may be used for a variety of purposes. The vehicle 102 can travel on a road 110. In this example environment, the road 110 is a one-way road with two lanes, a lane 112 and a lane 114. Moving objects may also be in or near the road 110, for example, a moving vehicle 120. Based on radar detections gathered by the radar sensor 106, the radar tracking module may determine a predicted path 122 of the vehicle 120.
The ADAS 104 can track objects in the field-of-view (FOV) based on sensor data obtained from one or more sensors of the vehicle 102. In this example, vehicle 102 is equipped with the radar sensor 106, but the vehicle 102 can be equipped with any number of sensors. The ADAS 104 tracks targets in the FOV based on sensor data received from sensors of the vehicle 102. Predicting the paths (i.e., the predicted path 122) of any targets in the FOV enables the ADAS 104 to reliably track the targets and avoid them while the vehicle 102 navigates the road 110.
In general, manufacturers can mount the ADAS 104 to any moving platform that can travel on the road 110. The ADAS 104 can project its FOV from any exterior surface of the vehicle 102. For example, vehicle manufactures can integrate at least a part of the ADAS 104 into a side mirror, roof, bumper, or any other interior or exterior location where the FOV includes the road 110 and any targets, moving or stationary, that are near the road 110. In the example, a portion of the ADAS 104 is mounted near the front bumper of the vehicle 102.
The ADAS 104 is equipped with the radar tracking module 108 that tracks targets within the FOV through a tracking function. The radar tracking module 108 comprises machine-readable instructions that when executed by a processor or other logic of the ADAS 104, causes the processor or other logic to accept detections from the radar sensor 106 and predict the path of the target. The radar tracking module 108 may include a combination of hardware, software, and/or firmware, for enabling target path predictions.
FIG. 2 illustrates an example of an automotive system 200 configured to perform target curvature estimation considering vehicle dynamics, in accordance with techniques of this disclosure. The automotive system 200 can be integrated within the vehicle 102.
The automotive system 200 includes a controller 202 and an ADAS 104-1. The ADAS 104-1 is an example of the ADAS 104 and can be integrated into an automotive or other vehicular environment. The ADAS 104-1 and the controller 202 communicate over a link 204. The link 204 may be a wired or wireless link and, in some cases, includes a communication bus. The controller 202 performs operations based on information received over the link 204, such as data output from the ADAS 104 as objects in the FOV are identified and tracked.
The controller 202 includes a processor 206 and a computer-readable storage medium (CRM) 208 (e.g., a memory, long-term storage, short-term storage), which stores instructions for an automotive module 210. The ADAS 104-1 includes a radar sensor 106-1. Any number of other sensors, including the radar sensor 106-1, may likewise be used. The ADAS 104-1 may include processing hardware that includes a processor 212 and a computer-readable storage medium (CRM) 214, which stores instructions associated with a radar tracking module 108-1. The radar tracking module 108-1, which is an example of the radar tracking module 108, may include a tracking filter 216 (e.g., EKF) defined by one or more motion models. In this example, the tracking filter 216 is defined by a CTCA model 218 and a CCTCA model 220. The CCTCA model 220 is constrained by an upper-bound parameter 222 of dynamic vehicle states that is determined by techniques disclosed herein.
The processors 206 and 212 can be two separate or a single microprocessor, or a pair of or a single system-on-chip of a computing device, a controller, or a control unit. The processors 206 and 212 execute computer-executable instructions stored within the CRMs 208 and 214. As an example, the processor 206 can execute the automotive module 210 to perform a driving function or other operation of the automotive system 200. Similarly, the processor 212 can execute the radar tracking module 108-1 to infer objects in the FOV based on sensor data obtained from multiple different sensors 106 of the vehicle 102. The automotive module 210, when executing at the processor 206, can receive an indication of one or more objects detected by the radar tracking module 108-1 in response to the radar tracking module 108-1 combining and analyzing sensor data generated at each of the sensors 106 (e.g., the radar sensor 106-1).
The radar tracking module 108-1 can receive target data from the radar sensor 106-1, and in response to receiving the target data, input the target data into the CTCA model 218. The CTCA model 218 can estimate vehicle dynamic states from the target data. For instance, if the CTCA model 218 estimates a vehicle state that exceeds a threshold set by an upper-bound parameter 222, the unconstrained estimation can be projected to a constrained surface that is constrained due to the estimation being replaced by the upper-bound parameter 222, as determined by a CCTCA model 220. By constraining the dynamic vehicle state estimation with a maximum realistic bound (e.g., the upper-bound parameter 222), the dynamic vehicle state estimate may reflect a realistic value. The estimated dynamic vehicle states, determined by either the CTCA model 218 or the CCTCA model 220 (if the one or more states were constrained), can be input to a tracking filter 216. The tracking filter 216 can use the state estimations to predict a path of the target. As the radar tracking module 108-1 tracks the target path, based on the path prediction by the tracking filter 216, the ADAS 104-1 can report instructions for avoiding the target to the automotive module 210.
Example Architecture
The following example architecture includes examples of determining the primary, secondary, and ultimate upper bounds. Once an ultimate upper bound is determined, it may be used by the radar tracking module to constrain the vehicle (e.g., centripetal) acceleration state for path prediction purposes.
FIG. 3 illustrates an example of a vehicle dynamics model (e.g., a bicycle model) that is employed to determine target curvature estimation considering vehicle dynamics constraints, in accordance with techniques of this disclosure. Physical limitations of vehicle dynamics may restrict the achievable centripetal acceleration of a vehicle. The bicycle model and the Pacejka Magic Formula are employed to perform vehicle dynamics analysis and derive a primary upper bound of the achievable centripetal acceleration. The bicycle model combines, for analytical purposes, both front tires into a single front tire and both rear tires into a single rear tire.
The bicycle model in FIG. 3 uses the following variables:

    F: tire force;
    L: distance from Center of Gravity (COG) to front or rear axle;
    v: velocity (note that v without subscripts f or r represents the velocity of vehicle COG);
    ω: yaw rate (heading rate);
    α: tire side-slip angle (steering angle not included);
    φ: steering angle of a front tire;
    subscripts:


x: longitudinal;
y: lateral;
f: front tire;
r: rear tire.
In one example of bicycle model analysis, the target vehicle drives on a flat road with no lateral slope, and the target vehicle has four wheels (two axles) and steers with the front wheels. The vehicle dynamics in the bicycle model can be represented by the Eq. Set (2):




  
    
      
        {
        
          
            
              
                
                  v
                  x
                
                =
                
                  
                    v
                    f
                  
                  ·
                  
                    cos
                    ⁡
                    
                      (
                      
                        φ
                        -
                        
                          α
                          f
                        
                      
                      )
                    
                  
                
              
            
          
          
            
              
                
                  v
                  x
                
                =
                
                  
                    
                      v
                      r
                    
                    ·
                    cos
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    (
                    
                      α
                      r
                    
                    )
                  
                
              
            
          
          
            
              
                
                  
                    v
                    y
                  
                  +
                  
                    ω
                    ·
                    
                      L
                      f
                    
                  
                
                =
                
                  
                    v
                    f
                  
                  ·
                  
                    sin
                    ⁡
                    
                      (
                      
                        φ
                        -
                        
                          α
                          f
                        
                      
                      )
                    
                  
                
              
            
          
          
            
              
                
                  
                    v
                    y
                  
                  -
                  
                    ω
                    *
                    
                      L
                      r
                    
                  
                
                =
                
                  
                    -
                    
                      v
                      r
                    
                  
                  ·
                  
                    sin
                    ⁡
                    
                      (
                      
                        α
                        r
                      
                      )
                    
                  
                
              
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




The tire side-slip angles, αf and αr, can be derived from Eq. Set (2) and represented as Eq. Set (3):





  
    
      
        {
        
          
            
              
                
                  α
                  f
                
                =
                
                  φ
                  -
                  
                    
                      tan
                      
                        -
                        1
                      
                    
                    ⁡
                    
                      (
                      
                        
                          
                            v
                            y
                          
                          +
                          
                            ω
                            ·
                            
                              L
                              f
                            
                          
                        
                        
                          v
                          x
                        
                      
                      )
                    
                  
                
              
            
          
          
            
              
                
                  α
                  r
                
                =
                
                  -
                  
                    
                      tan
                      
                        -
                        1
                      
                    
                    ⁡
                    
                      (
                      
                        
                          
                            v
                            y
                          
                          +
                          
                            ω
                            ·
                            
                              L
                              r
                            
                          
                        
                        
                          v
                          x
                        
                      
                      )
                    
                  
                
              
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




The equations in Eq. Set (3) are nonlinear, and not all the parameters are perceivable from a host sensor. However, general values for most vehicles can be found in industry publications related to tire and vehicle dynamics. Eq. (4) is an example of general values from industry publications used in the techniques described herein:





  
    
      
        
          α
          max
        
        =
        
          {
          
            
              
                
                  
                    3
                    ⁢
                    
                      
                    
                    ⁢
                    degrees
                  
                  -
                  
                    cars
                    ⁢
                    
                      
                    
                    ⁢
                    including
                    ⁢
                    
                      
                    
                    ⁢
                    SUVs
                  
                
              
            
            
              
                
                  
                    2
                    ⁢
                    
                      
                    
                    ⁢
                    degrees
                  
                  -
                  
                    trucks
                    ⁢
                    
                      
                    
                    ⁢
                    and
                    ⁢
                    
                      
                    
                    ⁢
                    buses
                  
                
              
            
          
        
      
    
    
      
        (
        4
        )
      
    
  




where αmax is the maximum tire side-slip angle for most vehicles.

The Pacejka Magic Formula, Eq. (5), is employed to derive the lateral force that each tire can provide, Fy primary.




  
    
      
        
          F
          
            y
            ⁢
            
              
            
            ⁢
            primary
          
        
        =
        
          m
          ·
          g
          ·
          D
          ·
          
            sin
            ⁡
            
              (
              
                C
                ·
                
                  
                    tan
                    
                      -
                      1
                    
                  
                  ⁡
                  
                    (
                    
                      
                        B
                        ·
                        α
                      
                      -
                      
                        E
                        ·
                        
                          (
                          
                            
                              B
                              ·
                              α
                            
                            -
                            
                              
                                tan
                                
                                  -
                                  1
                                
                              
                              ⁡
                              
                                (
                                
                                  B
                                  ·
                                  α
                                
                                )
                              
                            
                          
                          )
                        
                      
                    
                    )
                  
                
              
              )
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




where α is the tire side-slip angle, and the coefficient B represents stiffness factor, coefficient C represents shape factor, m represents vehicle mass loaded on the tire, and g represents gravity. Coefficient D represents peak factor, and coefficient E curvature factor. The subscript y indicates “in the lateral direction,” and primary indicates “for primary bound.”

FIG. 4 illustrates a graph representing the relationship, using the Pacejka Magic Formula, between a lateral tire force and a tire side-slip angle, in accordance with techniques of this disclosure. Fy is the tire lateral force, αD is the tire side-slip angle corresponding to the maximum tire lateral force, and tan−1(Cα) is the slope of the curve at the origin point.
Shown in Eq. (4) is αmax, which is smaller than αD for most vehicles. As a result, Eq. (5) is used to derive the tire lateral force Fy corresponding to αmax. FIG. 5 illustrates the coefficients of lateral tire force versus tire side-slip angle used in the Pacejka Magic Formula, in accordance with techniques of this disclosure. The general coefficient values can be derived from industry publications related to tire and vehicle dynamics The maximum centripetal acceleration, αprimary can be derived from Fy by:




  
    
      
        
          α
          
            p
            ⁢
            rimary
          
        
        =
        
          
            F
            
              y
              ⁢
              
                
              
              ⁢
              primary
            
          
          m
        
      
    
    
      
        (
        6
        )
      
    
  




The applied Pacejka Magic Formula only considers tire lateral forces (which are of the same value but counteracted by tire friction forces from the ground). However, the total tire forces include both lateral and longitudinal force components. Thus, the friction force ellipse model [4] is employed to compute the maximum lateral forces that a tire can provide by considering the contribution of tire longitudinal forces. These forces are used to derive the secondary upper bound.
FIG. 6 illustrates a graph representing a tire friction forces ellipse vehicle dynamics model, in accordance with techniques of this disclosure. In the illustration, the variables have the following meaning:

    Fxf, max: maximum tire longitudinal force;
    Fyf, max: maximum tire lateral force;
    Fxf: longitudinal tire force;
    Fyf, potential: maximum achievable tire lateral force;
    Ff: resultant tire force (combined longitudinal and lateral forces).
    The maximum friction forces that a tire can achieve are represented by the blue ellipse and is given by Eq. Set (7),






  
    
      
        {
        
          
            
              
                
                  F
                  
                    xf
                    ,
                    max
                  
                
                =
                
                  
                    μ
                    x
                  
                  ·
                  
                    F
                    
                      z
                      ⁢
                      f
                    
                  
                
              
            
          
          
            
              
                
                  F
                  
                    yf
                    ,
                    max
                  
                
                =
                
                  
                    μ
                    y
                  
                  ·
                  
                    F
                    
                      z
                      ⁢
                      f
                    
                  
                
              
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




where

    μ: the friction coefficient between the tire and the road;
    Fzf: the vertical force (gravity of Earth) on the tire;
    subscripts


x: longitudinal direction;
y: lateral direction.
When the longitudinal friction forces increase, the maximum achievable lateral friction forces decrease. In other words, the resultant force value of the maximum tire lateral friction forces derived for the primary bound and the tire longitudinal forces never goes beyond the circle. When there are longitudinal tire friction forces, which is true during most of the driving time, the maximum achievable tire lateral friction forces will decrease from Fyf,max to Fyf,potential. Therefore, it may not be necessary to compare the lateral forces derived for the primary bound with the maximum achievable lateral forces and make sure that they do not go beyond the friction force ellipse. The ellipse may be formulated by Eq. (8),




  
    
      
        
          
            
              (
              
                
                  F
                  
                    x
                    ⁢
                    f
                  
                
                
                  F
                  
                    xf
                    ,
                    max
                  
                
              
              )
            
            2
          
          +
          
            
              (
              
                
                  F
                  
                    yf
                    ,
                    potential
                  
                
                
                  F
                  
                    yf
                    ,
                    max
                  
                
              
              )
            
            2
          
        
        =
        1
      
    
    
      
        (
        8
        )
      
    
  




By accounting for tire longitudinal forces, the maximum achievable tire lateral forces may be given by Eqs. (9) and (10),

Fy,secondary=√{square root over ((μmg)2−Fxf2)}  (9)

where

Fxf≅0.125·ρ·CdA·v2+μroll·m·g   (10)

P: air density,

CdA: drag coefficient multiplied by drag area,

μ: friction coefficient between the tire and road,

μroll: rolling friction coefficient between the tire and road,

v: vehicle velocity,

m: lower bound of mass of each certain type of vehicle loaded on the tire,

g: gravity of Earth,

The secondary bound of the centripetal acceleration is derived in Eq (11),




  
    
      
        
          α
          
            s
            ⁢
            e
            ⁢
            condary
          
        
        =
        
          
            F
            
              y
              ⁢
              
                
              
              ⁢
              secondary
            
          
          m
        
      
    
    
      
        (
        11
        )
      
    
  




Note that the longitudinal tire forces Fxf given in Eq. (10) is a conservative estimate as the parameters in Eqs. (9)-(11) are set to the value of a lower bound for certain types of vehicles to cover extreme cases. Fxf may be smaller than real values for most vehicles. As a result, Fy,secondary may be greater than real values. Therefore, when comparing the secondary bound with the primary one, the smaller value is chosen as the value of the ultimate upper bound.
FIG. 7 illustrates ultimate upper-bound values of centripetal acceleration determined from the vehicle dynamics models employed, in accordance with techniques of this disclosure. These ultimate upper-bound values will be applied in a CCTCA-EKF tracking filter, as presented herein.
The vehicle dynamic constraints derived in FIG. 7 may be added to the CTCA tracker filter to provide better curvature estimation. FIG. 8 illustrates some kinematic states that may be included in a state vector of a motion model such as the CTCA model. The CTCA state vector is denoted as X=[x, y, ψ, s, c, a]T, where

    x: longitudinal position,
    y: lateral position,
    ψ: heading angle,
    s: velocity,
    c: curvature,
    a: tangential acceleration.

Taking the centripetal acceleration upper bound into account, the target vehicle follows the CCTCA motion model, Eq. (12):







  
    
      
        
          
            [
            
              
                
                  
                    x
                    .
                  
                
              
              
                
                  
                    y
                    .
                  
                
              
              
                
                  
                    ψ
                    .
                  
                
              
              
                
                  
                    s
                    .
                  
                
              
              
                
                  
                    c
                    .
                  
                
              
              
                
                  
                    a
                    .
                  
                
              
            
            ]
          
          =
          
            
              [
              
                
                  
                    
                      s
                      *
                      
                        cos
                        ⁡
                        
                          (
                          ψ
                          )
                        
                      
                    
                  
                
                
                  
                    
                      s
                      *
                      
                        sin
                        ⁡
                        
                          (
                          ψ
                          )
                        
                      
                    
                  
                
                
                  
                    
                      c
                      ⁢
                      s
                    
                  
                
                
                  
                    a
                  
                
                
                  
                    0
                  
                
                
                  
                    0
                  
                
              
              ]
            
            +
            
              [
              
                
                  
                    0
                  
                
                
                  
                    0
                  
                
                
                  
                    0
                  
                
                
                  
                    0
                  
                
                
                  
                    
                      n
                      1
                    
                  
                
                
                  
                    
                      n
                      2
                    
                  
                
              
              ]
            
          
        
        ,
        
          
            n
            1
          
          ⁢
          
            ∼
          
          ⁢
          
            N
            ⁡
            
              (
              
                0
                ,
                
                  q
                  1
                
              
              )
            
          
        
        ,
        
          
            n
            2
          
          ⁢
          
            ∼
          
          ⁢
          
            N
            ⁡
            
              (
              
                0
                ,
                
                  q
                  2
                
              
              )
            
          
        
      
    
    
      
        (
        12
        )
      
    
  




subject to the centripetal acceleration constraint, Eq. (13):

|ac|=|c*s2|<d   (13)

where the upper bound of the centripetal acceleration is known from FIG. 7. For different types of vehicles at different road conditions, d can take different values. Given the CTCA state vector, both the state transition model and constrain are nonlinear. A constrained EKF is implemented to resolve the problem.

FIG. 9 illustrates a flow chart of the CCTCA tracking filter used to determine target vehicle centripetal acceleration estimation considering vehicle dynamics constraints, in accordance with techniques of this disclosure. It is generally executed on top of the unconstrained CTCA-EKF filter but with two additional steps: 1) it checks whether a target in the application scope triggers the centripetal acceleration constraints, and 2) for any track violating the constraint, it projects the unconstrained estimation to the constrained surface.
Given the CTCA state vector and its state transition function mentioned above, we can use the following EKF filter to estimate its state.




  
    
      
        {
        
          
            
              
                
                  
                    x
                    ^
                  
                  
                    
                      k
                      _
                    
                    +
                    1
                  
                
                =
                
                  
                    F
                    k
                  
                  ⁢
                  
                    
                      x
                      ^
                    
                    k
                  
                
              
            
          
          
            
              
                
                  P
                  
                    
                      k
                      _
                    
                    +
                    1
                  
                
                =
                
                  
                    
                      F
                      k
                    
                    ⁢
                    
                      P
                      k
                    
                    ⁢
                    
                      F
                      k
                      T
                    
                  
                  +
                  Q
                
              
            
          
        
      
    
    
      
        (
        14
        )
      
    
  
  
    
      
        {
        
          
            
              
                
                  K
                  
                    k
                    +
                    1
                  
                
                =
                
                  
                    P
                    
                      
                        k
                        _
                      
                      +
                      1
                    
                  
                  ⁢
                  
                    
                      
                        H
                        T
                      
                      ⁡
                      
                        (
                        
                          
                            H
                            ⁢
                            
                              P
                              
                                
                                  k
                                  _
                                
                                +
                                1
                              
                            
                            ⁢
                            
                              H
                              T
                            
                          
                          +
                          R
                        
                        )
                      
                    
                    
                      -
                      1
                    
                  
                
              
            
          
          
            
              
                
                  
                    x
                    ^
                  
                  
                    k
                    +
                    1
                  
                
                =
                
                  
                    
                      x
                      ^
                    
                    
                      
                        k
                        _
                      
                      +
                      1
                    
                  
                  +
                  
                    
                      K
                      
                        k
                        +
                        1
                      
                    
                    ⁡
                    
                      (
                      
                        
                          z
                          
                            k
                            +
                            1
                          
                        
                        -
                        
                          H
                          ⁢
                          
                            
                              x
                              ^
                            
                            
                              
                                k
                                _
                              
                              +
                              1
                            
                          
                        
                      
                      )
                    
                  
                
              
            
          
          
            
              
                
                  P
                  
                    k
                    +
                    1
                  
                
                =
                
                  
                    (
                    
                      I
                      -
                      k
                      +
                      
                        1
                        ⁢
                        H
                      
                    
                    )
                  
                  ⁢
                  
                    P
                    
                      
                        k
                        _
                      
                      +
                      1
                    
                  
                
              
            
          
        
      
    
    
      
        (
        15
        )
      
    
  




where Eq. (14) is used to predict the state and Eq. (15) is used to update the state, and

    xk: state vector at scan index k,
    Pk: error covariance at scan index k,
    Pk+1: error covariance at scan index k,
    Q: covariance of process noise,
    R: covariance of measurement noise,
    z: measurement vector,
    H: observation matrix, will be Jacobian matrix in the case of nonlinear measurement,
    Fk: state transition matrix (Jacobian matrix).

By solving the Eq. (12) of the CTCA state transition function, the state transition matrix Fk can be expressed as:







  
    
      
        
          F
          k
        
        =
        
          [
          
            
              
                1
              
              
                0
              
              
                
                  F
                  
                    1
                    ⁢
                    3
                  
                
              
              
                
                  F
                  
                    1
                    ⁢
                    4
                  
                
              
              
                
                  F
                  
                    1
                    ⁢
                    5
                  
                
              
              
                
                  F
                  
                    1
                    ⁢
                    6
                  
                
              
            
            
              
                0
              
              
                1
              
              
                
                  F
                  
                    2
                    ⁢
                    3
                  
                
              
              
                
                  F
                  
                    2
                    ⁢
                    4
                  
                
              
              
                
                  F
                  
                    2
                    ⁢
                    5
                  
                
              
              
                
                  F
                  
                    2
                    ⁢
                    6
                  
                
              
            
            
              
                0
              
              
                0
              
              
                1
              
              
                
                  F
                  
                    3
                    ⁢
                    4
                  
                
              
              
                
                  F
                  
                    3
                    ⁢
                    5
                  
                
              
              
                
                  F
                  36
                
              
            
            
              
                0
              
              
                0
              
              
                0
              
              
                1
              
              
                0
              
              
                T
              
            
            
              
                0
              
              
                0
              
              
                0
              
              
                0
              
              
                1
              
              
                0
              
            
            
              
                0
              
              
                0
              
              
                0
              
              
                0
              
              
                0
              
              
                1
              
            
          
          ]
        
      
    
    
      
        (
        16
        )
      
    
  
  
    
      Where
    
    
      
        
      
    
  
  
    
      
        {
        
          
            
              
                
                  F
                  
                    1
                    ⁢
                    3
                  
                
                =
                
                  
                    
                      -
                      
                        sin
                        ⁡
                        
                          (
                          
                            ψ
                            k
                          
                          )
                        
                      
                    
                    ⁢
                    
                      (
                      
                        
                          
                            s
                            k
                          
                          ⁢
                          T
                        
                        +
                        
                          
                            a
                            k
                          
                          ⁢
                          
                            
                              T
                              2
                            
                            2
                          
                        
                      
                      )
                    
                  
                  -
                  
                    
                      cos
                      ⁡
                      
                        (
                        
                          ψ
                          k
                        
                        )
                      
                    
                    ⁢
                    
                      (
                      
                        
                          
                            c
                            k
                          
                          ⁢
                          
                            s
                            k
                            2
                          
                          ⁢
                          
                            
                              T
                              2
                            
                            2
                          
                        
                        +
                        
                          
                            c
                            k
                          
                          ⁢
                          
                            a
                            k
                          
                          ⁢
                          
                            s
                            k
                          
                          ⁢
                          
                            
                              T
                              3
                            
                            2
                          
                        
                        +
                        
                          
                            c
                            k
                          
                          ⁢
                          
                            a
                            k
                            2
                          
                          ⁢
                          
                            
                              T
                              4
                            
                            8
                          
                        
                      
                      )
                    
                  
                
              
            
          
          
            
              
                
                  F
                  
                    1
                    ⁢
                    4
                  
                
                =
                
                  
                    
                      cos
                      ⁡
                      
                        (
                        
                          ψ
                          k
                        
                        )
                      
                    
                    ⁢
                    T
                  
                  -
                  
                    
                      sin
                      ⁡
                      
                        (
                        
                          ψ
                          k
                        
                        )
                      
                    
                    ⁢
                    
                      (
                      
                        
                          
                            c
                            k
                          
                          ⁢
                          
                            s
                            k
                          
                          ⁢
                          
                            T
                            2
                          
                        
                        +
                        
                          
                            c
                            k
                          
                          ⁢
                          
                            a
                            k
                          
                          ⁢
                          
                            
                              T
                              3
                            
                            2
                          
                        
                      
                      )
                    
                  
                
              
            
          
          
            
              
                
                  F
                  
                    1
                    ⁢
                    5
                  
                
                =
                
                  
                    -
                    
                      sin
                      ⁡
                      
                        (
                        
                          ψ
                          k
                        
                        )
                      
                    
                  
                  ⁢
                  
                    (
                    
                      
                        
                          s
                          k
                          2
                        
                        ⁢
                        
                          
                            T
                            2
                          
                          2
                        
                      
                      +
                      
                        
                          a
                          k
                        
                        ⁢
                        
                          s
                          k
                        
                        ⁢
                        
                          
                            T
                            3
                          
                          2
                        
                      
                      +
                      
                        
                          a
                          k
                          2
                        
                        ⁢
                        
                          
                            T
                            4
                          
                          8
                        
                      
                    
                    )
                  
                
              
            
          
          
            
              
                
                  F
                  
                    1
                    ⁢
                    6
                  
                
                =
                
                  
                    
                      cos
                      ⁡
                      
                        (
                        
                          ψ
                          k
                        
                        )
                      
                    
                    ⁢
                    
                      
                        T
                        2
                      
                      2
                    
                  
                  -
                  
                    
                      sin
                      ⁡
                      
                        (
                        
                          ψ
                          k
                        
                        )
                      
                    
                    ⁢
                    
                      (
                      
                        
                          
                            c
                            k
                          
                          ⁢
                          
                            s
                            k
                          
                          ⁢
                          
                            
                              T
                              3
                            
                            2
                          
                        
                        +
                        
                          
                            c
                            k
                          
                          ⁢
                          
                            a
                            k
                          
                          ⁢
                          
                            
                              T
                              4
                            
                            4
                          
                        
                      
                      )
                    
                  
                
              
            
          
          
            
              
                
                  F
                  
                    2
                    ⁢
                    3
                  
                
                =
                
                  
                    
                      cos
                      ⁡
                      
                        (
                        
                          ψ
                          k
                        
                        )
                      
                    
                    ⁢
                    
                      (
                      
                        
                          
                            s
                            k
                          
                          ⁢
                          T
                        
                        +
                        
                          
                            a
                            k
                          
                          ⁢
                          
                            
                              T
                              2
                            
                            2
                          
                        
                      
                      )
                    
                  
                  -
                  
                    
                      sin
                      ⁡
                      
                        (
                        
                          ψ
                          k
                        
                        )
                      
                    
                    ⁢
                    
                      (
                      
                        
                          
                            c
                            k
                          
                          ⁢
                          
                            s
                            k
                            2
                          
                          ⁢
                          
                            
                              T
                              2
                            
                            2
                          
                        
                        +
                        
                          
                            c
                            k
                          
                          ⁢
                          
                            a
                            k
                          
                          ⁢
                          
                            s
                            k
                          
                          ⁢
                          
                            
                              T
                              3
                            
                            2
                          
                        
                        +
                        
                          
                            c
                            k
                          
                          ⁢
                          
                            a
                            k
                            2
                          
                          ⁢
                          
                            
                              T
                              4
                            
                            8
                          
                        
                      
                      )
                    
                  
                
              
            
          
          
            
              
                
                  F
                  
                    2
                    ⁢
                    4
                  
                
                =
                
                  
                    
                      sin
                      ⁡
                      
                        (
                        
                          ψ
                          k
                        
                        )
                      
                    
                    ⁢
                    T
                  
                  -
                  
                    
                      cos
                      ⁡
                      
                        (
                        
                          ψ
                          k
                        
                        )
                      
                    
                    ⁢
                    
                      (
                      
                        
                          
                            c
                            k
                          
                          ⁢
                          
                            s
                            k
                          
                          ⁢
                          
                            T
                            2
                          
                        
                        +
                        
                          
                            c
                            k
                          
                          ⁢
                          
                            a
                            k
                          
                          ⁢
                          
                            
                              T
                              3
                            
                            2
                          
                        
                      
                      )
                    
                  
                
              
            
          
          
            
              
                
                  F
                  
                    2
                    ⁢
                    5
                  
                
                =
                
                  
                    cos
                    ⁡
                    
                      (
                      
                        ψ
                        k
                      
                      )
                    
                  
                  ⁢
                  
                    (
                    
                      
                        
                          s
                          k
                          2
                        
                        ⁢
                        
                          
                            T
                            2
                          
                          2
                        
                      
                      +
                      
                        
                          a
                          k
                        
                        ⁢
                        
                          s
                          k
                        
                        ⁢
                        
                          
                            T
                            3
                          
                          2
                        
                      
                      +
                      
                        
                          a
                          k
                          2
                        
                        ⁢
                        
                          
                            T
                            4
                          
                          8
                        
                      
                    
                    )
                  
                
              
            
          
          
            
              
                
                  F
                  
                    2
                    ⁢
                    6
                  
                
                =
                
                  
                    
                      sin
                      ⁡
                      
                        (
                        
                          ψ
                          k
                        
                        )
                      
                    
                    ⁢
                    
                      
                        T
                        2
                      
                      2
                    
                  
                  +
                  
                    
                      cos
                      ⁡
                      
                        (
                        
                          ψ
                          k
                        
                        )
                      
                    
                    ⁢
                    
                      (
                      
                        
                          
                            c
                            k
                          
                          ⁢
                          
                            s
                            k
                          
                          ⁢
                          
                            
                              T
                              3
                            
                            2
                          
                        
                        +
                        
                          
                            c
                            k
                          
                          ⁢
                          
                            a
                            k
                          
                          ⁢
                          
                            
                              T
                              4
                            
                            4
                          
                        
                      
                      )
                    
                  
                
              
            
          
          
            
              
                
                  F
                  
                    3
                    ⁢
                    4
                  
                
                =
                
                  
                    c
                    k
                  
                  ⁢
                  T
                
              
            
          
          
            
              
                
                  F
                  
                    3
                    ⁢
                    5
                  
                
                =
                
                  
                    
                      s
                      k
                    
                    ⁢
                    T
                  
                  +
                  
                    
                      a
                      k
                    
                    ⁢
                    
                      
                        T
                        2
                      
                      2
                    
                  
                
              
            
          
          
            
              
                
                  F
                  
                    3
                    ⁢
                    6
                  
                
                =
                
                  
                    c
                    k
                  
                  ⁢
                  
                    
                      T
                      2
                    
                    2
                  
                
              
            
          
        
      
    
    
      
        (
        17
        )
      
    
  




Assuming the tracker filter is properly initialized, the process noise Q can be derived by a reasonable noise model, and the measurement vector z and its error covariance R are known from the sensor model, the EKF above may be rolled over to estimate the state at each cycle.
At the end of each tracking cycle, a check can be made as to whether the centripetal acceleration constraint is triggered for each vehicle track in the constraint scope. For the sake of simplicity, in the following analysis it is assumed the curvature is positive, and therefore the sign of centripetal acceleration is also positive. At track index k, for any track that violated the inequality constraint defined in Eq. (13), the track state is projected to the closest surface defined by

f(xk)=ck*sk2=d   (18)

To further simplify the analysis, first linearize the constraint surface defined in the above equation at the estimated state {circumflex over (x)}k is linearized. Ignoring the higher-order terms results in




  
    
      
        
          
            
              
                c
                ^
              
              k
            
            *
            
              
                s
                ^
              
              k
            
          
          +
          
            
              D
              k
            
            ⁡
            
              (
              
                
                  x
                  k
                
                -
                
                  
                    x
                    ^
                  
                  k
                
              
              )
            
          
        
        =
        
          
            d
            ⁢
            
              →
              yields
            
            ⁢
            
              
                D
                k
              
              ⁢
              
                x
                k
              
            
          
          =
          
            d
            l
          
        
      
    
    
      
        (
        19
        )
      
    
  




Where D=[0, 0, 0, 2ŝkĉk, sk2, 0] and dl=d−ĉk*ŝk+Dk{circumflex over (x)}k.

To project the states to the constrained surface, it is preferable the constrained state to have the maximum probability. In the EKF filter, the unconstrained estimation {circumflex over (x)}k follows the normal distribution: N({circumflex over (x)}k, Pk+1). To maximize the probability of the constrained estimation {tilde over (x)}k, the problem turns out to be min[({tilde over (x)}k−{circumflex over (x)}k)T Pk−1({tilde over (x)}k−{circumflex over (x)}k)], subject to Dk{tilde over (x)}k=dl.
The Lagrangian function L=({tilde over (x)}k−{circumflex over (x)}k)T Pk−1({tilde over (x)}k−{circumflex over (x)}k)+2λ(Dk{tilde over (x)}k−dl)




  
    
      
        {
        
          
            
              
                
                  
                    δ
                    ⁢
                    L
                  
                  
                    δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        x
                        ~
                      
                      k
                    
                  
                
                =
                
                  
                    0
                    ⁢
                    
                      →
                      yields
                    
                    ⁢
                    
                      
                        
                          P
                          k
                          
                            -
                            1
                          
                        
                        ⁡
                        
                          (
                          
                            
                              
                                x
                                ~
                              
                              k
                            
                            -
                            
                              
                                x
                                ^
                              
                              k
                            
                          
                          )
                        
                      
                      +
                      
                        
                          D
                          k
                          T
                        
                        ⁢
                        λ
                      
                    
                  
                  =
                  0
                
              
            
          
          
            
              
                
                  
                    δ
                    ⁢
                    
                      
                    
                    ⁢
                    L
                  
                  
                    δ
                    ⁢
                    λ
                  
                
                =
                
                  
                    0
                    ⁢
                    
                      →
                      yields
                    
                    ⁢
                    
                      
                        
                          D
                          k
                        
                        ⁢
                        
                          
                            x
                            ~
                          
                          k
                        
                      
                      -
                      
                        d
                        l
                      
                    
                  
                  =
                  0
                
              
            
          
        
      
    
    
      
        (
        20
        )
      
    
  




solves the optimization problem.

Solving the above joint equation gives the constrained estimation:

{tilde over (x)}k={circumflex over (x)}k+PkDkT(DkPkDkT)−1(dl−Dk{circumflex over (x)}k)   (21)

Example Methods

FIG. 10 illustrates an example method for target curvature estimation considering vehicle dynamics constraints, in accordance with techniques of this disclosure. The method 1000 is shown as a set of operations (or acts) performed in, but not limited to, the order or combinations in which the operations are shown or described. Further, any of the operations may be repeated, combined, or reorganized to provide other methods. In portions of the following discussion, reference may be made to the previously introduced figures in describing some non-limiting examples of method 1000.
At 1002, measurement states of a target are determined by a processor of a host vehicle. The measurement states may be derived from radar detections received by an object tracking system of the host vehicle. The measurement states may include range, range rate, azimuth angle, and elevation angle. At 1004, based at least in part on the measurement states, one or more kinematic states of a motion model are determined. The kinematic states may include longitudinal position, lateral position, heading angle, velocity, curvature and tangential acceleration. The motion model may be defined for a tracking filter employed by the object tracking system to track the target.
At 1006, a path of the target is predicted using the one or more kinematic states as inputs to the motion model of the tracking filter. At 1008, the host vehicle is controlled, based on the predicted path of the target, to take evasive action and avoid a collision with the target.
Large errors in estimated kinematic states may be reduced, in accordance with techniques of this disclosure. In response to reducing the large errors in estimated kinematic states, the object tracking system may be able to predict a target path with accuracy, which may result in less collisions and more safety
Example Scenarios
FIGS. 11-1 to 11-3 illustrate example scenarios of target curvature estimation considering vehicle dynamics, in accordance with techniques of this disclosure. The example scenarios may be based on aspects of the example environment illustrated in FIG. 1 and the example automotive system illustrated in FIG. 2.
In FIG. 11-1, a host vehicle 1102 is equipped with an ADAS 1104, which is an example of the ADAS 104, and a radar sensor 1106-1, which is an example of the radar sensor 106. The host vehicle 1102 is navigating a road 1108 in a lane 1110 on a path 1112. In an adjacent lane 1114 of the road 1110, a target vehicle 1116 is navigating in the same direction as the host vehicle 1102 at a time T-1. Radar detections received by the radar sensor 1106 are used by a radar tracking module (i.e., radar tracking module 108-1) in the ADAS 1104 to predict a path 1118 of the target vehicle 1116.
As time transitions from time T-1 to time T-2, the radar tracking module of the ADAS 1104 may use various measurement states derived from the radar detections received by the radar sensor 1106 to estimate a curvature state of the target vehicle 1116. The curvature state may be used to determine an estimated centripetal acceleration. If the estimated centripetal acceleration state exceeds an upper bound as determined based on the techniques described herein, the upper bound value is used as the centripetal acceleration state value. Based on the predicted path 1118, the host vehicle 1102 may take evasive action and avoid a collision as the target vehicle moves into the lane 1110 at time T-2.
In FIG. 11-2, the host vehicle 1102 is backing into a parking space 1120. At time T-1, the ADAS 1104 may be tracking a pedestrian 1122-1 based on radar detections received by a radar sensor 1106-2, which is an example of the radar sensor 106 and is located in a rear area of the host vehicle 1102. As the host vehicle 1102 follows a path 1124, the ADAS 1104 may predict a path 1126 of the pedestrian 1122-1. By time T-2, the pedestrian 1122 may have walked into the parking space 1120. The host vehicle 1102 may avoid hitting the pedestrian by using the techniques previously described, estimating the centripetal acceleration of the pedestrian 1122, for example, to set the upper bound.
In FIG. 11-3, the host vehicle 1102 is navigating the road 1108 in the lane 1110 on the path 1128. At time T-5, the ADAS 1104 begins tracking a pedestrian 1130 and predicting a path 1132 for the pedestrian 1130. In the same manner as the previous scenarios, the host vehicle 1102 may be able to take evasive action and avoid a collision with the pedestrian 1130 as the pedestrian 1130 walks into the lane 1110 at time T-6
ADDITIONAL EXAMPLES
In the following section, additional examples of target curvature estimation considering vehicle dynamics constraints are provided.
Example 1
A method comprising: determining, by a processor of a host vehicle, measurement states of a target, the measurement states being derived from radar detections received by an object tracking system of the host vehicle; determining, based at least in part on the measurement states, one or more kinematic states of a motion model defined for a tracking filter employed by the object tracking system to track the target; predicting, using the one or more kinematic states as inputs to the motion model of the tracking filter, a path of the target; and controlling, based on the predicted path of the target, the host vehicle to take evasive action and avoid a collision with the target.
Example 2
The method of example 1, further comprising: determining a first upper bound for the estimates of the at least one or more kinematic states of the motion model; determining a second upper bound for the estimates of the at least one or more kinematic states of the motion model; determining a third upper bound for the estimates of the at least one or more kinematic states of the motion model, the third upper bound being equivalent to a lesser of the first upper bound and the second upper bound; and determining a constraint for the at least one or more kinematic states of the motion model, the constraint being equivalent to the third upper bound, wherein determining the one or more kinematic states is further based on the constraint.
Example 3
The method of example 2, wherein: determining the first upper bound comprises determining the first upper bound using the Pacejka Magic Formula.
Example 4
The method of example 2, wherein: determining the second upper bound comprises determining the second upper bound from a friction force ellipse model.
Example 5
The method of example 1, wherein: determining, based on the measurement states, estimates of the at least one or more kinematic states of a motion model, the motion model being a Constant Turn Constant Tangential Acceleration model.
Example 6
The method of example 1, wherein: predicting the path of the target by the tracking filter, based on the motion model, the tracking filter being an extended Kalman Filter.
Example 7
The method of example 2, wherein: determining the first upper bound, the second upper bound, and the third upper bound for the estimates of the at least one or more kinematic states of the motion model, the kinematic state being curvature state.
Example 8
A system comprising: a host vehicle; an object tracking system configured to receive radar detections; and one or more processors configured to: determine measurement states of a target, the measurement states being derived from the radar detections received by the object tracking system; determine one or more kinematic states of a motion model defined for a tracking filter employed by the object tracking system to track the target; predict, using the one or more kinematic states as inputs to the motion model of the tracking filter, a path of the target; and control, based on the predicted path of the target, the host vehicle to take evasive action and avoid a collision with the target.
Example 9
The system of example 8, wherein the one or more processors are further configured to determine a first upper bound for the estimates of the at least one or more kinematic states of the motion model; determine a second upper bound for the estimates of the at least one or more kinematic states of the motion model; determine a third upper bound for the estimates of the at least one or more kinematic states of the motion model, the third upper bound being equivalent to a lesser of the first upper bound and the second upper bound; and determine a constraint for the at least one or more kinematic states of the motion model, the constraint being equivalent to the third upper bound, wherein determining the one or more kinematic states is further based on the constraint.
Example 10
The system of example 9, wherein the one or more processors are further configured to determine the first upper bound using the Pacejka Magic Formula.
Example 11
The system of example 9, wherein the one or more processors are further configured to determine the second upper bound from a friction force ellipse model.
Example 12
The system of example 8, wherein the one or more processors are further configured to determine, based on the measurement states, estimates of the at least one or more kinematic states of a motion model, the motion model being a Constant Turn Constant Tangential Acceleration model.
Example 13
The system of example 8, wherein the one or more processors are further configured to predict the path of the target by the tracking filter, based on the motion model, the tracking filter being an extended Kalman Filter.
Example 14
The system of example 9, wherein the one or more processors are further configured to determine the first upper bound, the second upper bound, and the third upper bound for the estimates of the at least one or more kinematic states of the motion model, the kinematic state being curvature state.
Example 15
The system of example 8, wherein the object tracking system comprises a radar system configured to receive the radar detections.
Example 16
The system of example 8, wherein the target comprises another vehicle in a travel path of the host vehicle.
Example 17
The system of example 16, wherein the travel path comprises a forward trajectory of the host vehicle or a reverse trajectory of the host vehicle.
Example 18
The system of example 8, wherein the target comprises a pedestrian or animal in a travel path of the host vehicle.
Example 19
The system of example 8, wherein the at least one processor is further configured to, after initiating control of the vehicle to take the evasive action and avoid the collision with the target: update the kinematic states of the motion model based on updated measurement states of the target; predict, using the updated kinematic states as the inputs to the motion model, an updated path of the target; and control the host vehicle to take further evasive action and avoid the collision with the target based on the updated path of the target.
Example 20
A system comprising: means for determining, by a processor of a host vehicle, measurement states of a target, the measurement states being derived from radar detections received by an object tracking system of the host vehicle; means for determining, based on the measurement states, one or more kinematic states of a motion model defined for a tracking filter employed by the object tracking system to track the target; means for predicting, using the one or more kinematic states as inputs to the motion model of the tracking filter, a path of the target; and means for controlling, based on the predicted path of the target, the host vehicle to take evasive action and avoid a collision with the target.
CONCLUSION
While various embodiments of the disclosure are described in the foregoing description and shown in the drawings, it is to be understood that this disclosure is not limited thereto but may be variously embodied to practice within the scope of the following claims. From the foregoing description, it will be apparent that various changes may be made without departing from the scope of the disclosure as defined by the following claims. Problems associated with large errors in target curvature estimation can occur in other systems. Therefore, although described as a way to improve target curvature estimation considering vehicle dynamics techniques, the techniques of the foregoing description can be applied to other assignment problems to reduce a total number of assignments between tasks and agents.
The use of “or” and grammatically related terms indicates non-exclusive alternatives without limitation unless the context clearly dictates otherwise. As used herein, a phrase referring to “at least one of” a list of items refers to any combination of those items, including single members. As an example, “at least one of: a, b, or c” is intended to cover a, b, c, a-b, a-c, b-c, and a-b-c, as well as any combination with multiples of the same element (e.g., a-a, a-a-a, a-a-b, a-a-c, a-b-b, a-c-c, b-b, b-b-b, b-b-c, c-c, and c-c-c or any other ordering of a, b, and c).

## Claims

1. A method comprising:
determining, by a processor of a host vehicle, measurement states of a target, the measurement states being derived from radar detections received by an object tracking system of the host vehicle;
determining, based at least in part on the measurement states, estimates of one or more kinematic states of a motion model defined for a tracking filter employed by the object tracking system to track the target, at least one of the kinematic states being curvature;
determining, for a constrained motion model, a respective constrained kinematic state for each respective kinematic state of the one or more kinematic states by:
determining a first upper bound for each respective kinematic state using a first dynamics model;
determining a second upper bound for each respective kinematic state using a second dynamics model;
determining a third upper bound for each respective kinematic state by comparing the first upper bound and the second upper bound and selecting the respective upper bound that has a lesser value to be the third upper bound; and
using the third upper bound for each respective kinematic state as the constrained kinematic state for that respective kinematic state;

comparing, using the estimates of the one or more kinematic states as inputs to the motion model and using the respective constrained kinematic states as inputs to the constrained motion model, the motion model to the constrained motion model;
predicting, using the one or more kinematic states as inputs to the motion model of the tracking filter based on the comparison of the motion model to the constrained motion model, a path of the target; and
controlling, based on the predicted path of the target, the host vehicle to take evasive action and avoid a collision with the target.

2. The method of claim 1, wherein predicting, using the one or more kinematic states as inputs to the motion model of the tracking filter based on the comparison of the motion model to the constrained motion model, a path of the target comprises:
responsive to the one or more kinematic states exceeding the respective constrained kinematic state, projecting an unconstrained estimation to a constrained surface.

3. The method of claim 1, wherein:
determining the first upper bound comprises determining the first upper bound using the Pacejka Magic Formula.

4. The method of claim 1, wherein:
determining the second upper bound comprises determining the second upper bound from a friction force ellipse model.

5. The method of claim 1, wherein:
determining, based on the measurement states, estimates of the one or more kinematic states of a motion model, the motion model being a Constant Turn Constant Tangential Acceleration model.

6. The method of claim 1, wherein:
predicting the path of the target by the tracking filter, based on the motion model, the tracking filter being an extended Kalman Filter.

7. The method of claim 1, wherein:
determining the first upper bound, the second upper bound, and the third upper bound for the estimates of the one or more kinematic states of the motion model, the kinematic state being a centripetal acceleration state.

8. A system comprising:
a host vehicle;
an object tracking system configured to receive radar detections; and
one or more processors configured to:
determine measurement states of a target, the measurement states being derived from the radar detections received by the object tracking system;
determine estimates of one or more kinematic states of a motion model defined for a tracking filter employed by the object tracking system to track the target;
determine, for a constrained motion model, a respective constrained kinematic state for each respective kinematic state of the one or more kinematic states by:
determining a first upper bound for each respective kinematic state using a first dynamics model;
determining a second upper bound for each respective kinematic state using a second dynamics model;
determining a third upper bound for each respective kinematic state by comparing the first upper bound and the second upper bound and selecting the respective upper bound that has a lesser value to be the third upper bound; and
using the third upper bound for each respective kinematic state as the constrained kinematic state for that respective kinematic state;

compare, using the estimates of the one or more kinematic states as inputs to the motion model and using the respective constrained kinematic states as inputs to the constrained motion model, the motion model to the constrained motion model;
predict, using the one or more kinematic states as inputs to the motion model of the tracking filter based on the comparison of the motion model to the constrained motion model, a path of the target; and
control, based on the predicted path of the target, the host vehicle to take evasive action and avoid a collision with the target.

9. The system of claim 8, wherein the one or more processors are configured to predict, using the one or more kinematic states as inputs to the motion model of the tracking filter based on the comparison of the motion model to the constrained motion model, a path of the target by at least:
responsive to the one or more kinematic states exceeding the respective constrained kinematic state, projecting an unconstrained estimation to a constrained surface.

10. The system of claim 8, wherein the one or more processors are further configured to:
determine the first upper bound using the Pacejka Magic Formula.

11. The system of claim 8, wherein the one or more processors are further configured to:
determine the second upper bound from a friction force ellipse model.

12. The system of claim 8, wherein the one or more processors are further configured to:
determine, based on the measurement states, estimates of the one or more kinematic states of a motion model, the motion model being a Constant Turn Constant Tangential Acceleration model.

13. The system of claim 8, wherein the one or more processors are further configured to:
predict the path of the target by the tracking filter, based on the motion model, the tracking filter being an extended Kalman Filter.

14. The system of claim 8, wherein the one or more processors are further configured to:
determine the first upper bound, the second upper bound, and the third upper bound for the estimates of the one or more kinematic states of the motion model, the kinematic state being a centripetal acceleration state.

15. The system of claim 8, wherein the object tracking system comprises a radar system configured to receive the radar detections.

16. The system of claim 8, wherein the target comprises another vehicle in a travel path of the host vehicle.

17. The system of claim 16, wherein the travel path comprises a forward trajectory of the host vehicle or a reverse trajectory of the host vehicle.

18. The system of claim 8, wherein the target comprises a pedestrian or animal in a travel path of the host vehicle.

19. The system of claim 8, wherein the one or more processors are further configured to, after initiating control of the vehicle to take the evasive action and avoid the collision with the target:
update the kinematic states of the motion model based on updated measurement states of the target;
predict, using the updated kinematic states as the inputs to the motion model, an updated path of the target; and
control the host vehicle to take further evasive action and avoid the collision with the target based on the updated path of the target.

20. A non-transitory computer-readable storage media comprising instructions that, when executed, cause a processor to:
determine measurement states of a target, the measurement states being derived from radar detections received by an object tracking system;
determine estimates of one or more kinematic states of a motion model defined for a tracking filter employed by the object tracking system to track the target;
determine, for a constrained motion model, a respective constrained kinematic state for each respective kinematic state of the one or more kinematic states by:
determining a first upper bound for each respective kinematic state using a first dynamics model;
determining a second upper bound for each respective kinematic state using a second dynamics model;
determining a third upper bound for each respective kinematic state by comparing the first upper bound and the second upper bound and selecting the respective upper bound that has a lesser value to be the third upper bound; and
using the third upper bound for each respective kinematic state as the constrained kinematic state for that respective kinematic state;

compare, using the estimates of the one or more kinematic states as inputs to the motion model and using the respective constrained kinematic states as inputs to the constrained motion model, the motion model to the constrained motion model;
predict, using the one or more kinematic states as inputs to the motion model of the tracking filter based on the comparison of the motion model to the constrained motion model, a path of the target; and
control, based on the predicted path of the target, the host vehicle to take evasive action and avoid a collision with the target.