Systems and methods for doppler-enhanced radar tracking

1. A method for Doppler-enhanced radar tracking, comprising:

transmitting a probe signal;

in response to reflection of the probe signal by a tracking target, at a radar array comprising a first plurality of radar elements positioned substantially along a first radar axis, receiving a reflected probe signal, wherein the tracking target and radar array cooperatively define a target vector from the radar array to the tracking target;

based on the reflected probe signal, determining a target range between the radar array and the tracking target;

based on the reflected probe signal, determining a first target angle between a first reference vector and a first projected target vector; wherein the first projected target vector is a projection of the target vector onto a first reference plane, the first reference plane containing both of the first radar axis and the first reference vector;

based on the reflected probe signal, determining Doppler shift information, comprising:

determining a set of predicted target composite angles; and

for each predicted target composite angle of the set, determining a respective set of frequency domain outputs;

based on the reflected probe signal and the Doppler shift information, determining a target composite angle; wherein the target composite angle is an angle between the target vector and a composite reference vector; and

based on the target range, the first target angle, and the target composite angle, determining a three-dimensional position of the tracking target relative to the radar array.

2. The method of claim 1 , wherein determining the three-dimensional position comprises implementing a Kalman filter.

3. The method of claim 2 , wherein implementing the Kalman filter comprises:

receiving a prior position estimate associated with the tracking target;

predicting a predicted position based on the prior position estimate;

determining error bounds associated with the target composite angle; and

determining an updated position estimate based on the prior position estimate, the target range, the first target angle, the target composite angle, and the error bounds.

4. The method of claim 1 , further comprising receiving egovelocity data, wherein determining the target composite angle is performed based further on the egovelocity data.

5. The method of claim 1 , wherein the first target angle is determined based on phase information of the reflected probe signal.

6. The method of claim 1 , wherein, for each predicted target composite angle of the set, determining the respective set of frequency domain outputs comprises determining a Fourier transform based on a respective scaled time variable, wherein the respective scaled time variable is determined based on the predicted target composite angle.

7. The method of claim 6 , wherein, for each predicted target composite angle of the set, the respective scaled time variable is scaled by a cosine of the predicted target composite angle.

8. The method of claim 6 , wherein, for each predicted target composite angle of the set, determining the Fourier transform based on the respective scaled time variable comprises pre-warping data indicative of the reflected probe signal to account for change in the target composite angle over time.

9. The method of claim 8 , wherein pre-warping the data comprises selecting a set of samples from the data, wherein the set of samples defines a non-uniform sampling rate.

10. The method of claim 8 , wherein pre-warping the data comprises resampling the data at a non-uniform sampling rate.

11. The method of claim 8 , wherein pre-warping the data is performed based on an assumption of a substantially linear change in the target composite angle over time.

12. The method of claim 1 , wherein:

the radar array further comprises a second plurality of radar elements positioned substantially along a second radar axis, the second radar axis not substantially parallel to the first radar axis;

the method further comprises, based on the reflected probe signal, determining a second target angle between a second reference vector and a second projected target vector; wherein the second projected target vector is a projection of the target vector onto a second reference plane, the second reference plane containing both of the second radar axis and the second reference vector; and

the three-dimensional position is determined based further on the second target angle.

13. The method of claim 12 , wherein the second radar axis is substantially perpendicular to the first radar axis.

14. The method of claim 13 , wherein:

the first target angle is determined based on phase differences of the reflected probe signal between the radar elements of the first plurality; and

the second target angle is determined based on phase differences of the reflected probe signal between the radar elements of the second plurality.

15. The method of claim 12 , wherein the first and second target angles are determined based on phase information of the reflected probe signal.

16. The method of claim 15 , further comprising, based on the target composite angle, filtering information associated with the reflected probe signal, wherein the first and second target angles are determined based on the filtered information.

17. The method of claim 12 , wherein, for each predicted target composite angle of the set, determining the respective set of frequency domain outputs comprises determining a Fourier transform based on a respective scaled time variable, wherein the respective scaled time variable is determined based on the predicted target composite angle.

18. The method of claim 17 , wherein, for each predicted target composite angle of the set, the respective scaled time variable is scaled by a cosine of the predicted target composite angle.

19. The method of claim 12 , wherein:

the first target angle is an elevation angle; and

the second target angle is an azimuth angle.

20. The method of claim 12 , wherein the first and second radar axes are substantially coplanar.

21. The method of claim 1 , wherein determining the three-dimensional position comprises determining error bounds associated with the target composite angle and, based on the error bounds, implementing a Kalman filter.