Radar system and method for identifying multiple targets in a beam response spectrum
11313948 · 2022-04-26
Assignee
Inventors
Cpc classification
G01S13/42
PHYSICS
International classification
G01S7/41
PHYSICS
Abstract
A radar system may include a transmitter, a receiver, and a controller. The controller may calculate a received beam response spectrum based on the received reflected radar signal, detect a first maximum value of the received beam response spectrum, identify an angle corresponding to the first maximum value as a first target angle, obtain a threshold envelope based on the first maximum value and the first target angle, detect a second maximum value in a portion of the received beam response spectrum being greater than the threshold envelope, identify an angle corresponding to the second maximum value as a second target angle, and output the first target angle as the angle of arrival of the reflected radar signal from the first target and the second target angle as the angle of arrival of the reflected radar signal from the second target.
Claims
1. A radar system comprising: a transmitter structured to transmit a radar signal; a receiver structured to receive a reflected radar signal; and a controller configured to identify an angle of arrival of a reflected radar signal from a first target and an angle of arrival of a reflected radar signal from a second target by performing: calculating a received beam response spectrum based on the received reflected radar signal, the received beam response spectrum comprising a set of power values as a function of angle of arrival; detecting a first maximum value of the received beam response spectrum; identifying an angle corresponding to the first maximum value as a first target angle; obtaining a threshold envelope based on the first maximum value and the first target angle; detecting a second maximum value in a portion of the received beam response spectrum being greater than the threshold envelope; identifying an angle corresponding to the second maximum value as a second target angle; and outputting the first target angle as the angle of arrival of the reflected radar signal from the first target and the second target angle as the angle of arrival of the reflected radar signal from the second target.
2. The radar system of claim 1, wherein the received beam response spectrum comprises a target main lobe including the first maximum value; the threshold envelope is a function having a value equal to the first maximum value for angles corresponding to a width of the target main lobe, and, for all other angles, having a value equal to a highest sidelobe peak level among sidelobes immediately adjacent to the target main lobe in the received beam response spectrum.
3. The radar system of claim 1, further comprising: a storage device storing a plurality of stored threshold envelopes, each stored threshold envelope of the plurality of stored threshold envelopes being associated with a different angle of arrival; wherein controller is configured such that the obtaining a threshold envelope based on the first maximum value and the first target angle comprises: identifying a first stored threshold envelope among the plurality of stored threshold envelopes having an associated angle of arrival corresponding to the first target angle; setting the first stored threshold envelope as the threshold envelope.
4. The radar system of claim 1, wherein the controller is further configured to generate a cumulative threshold envelope initially equal to the threshold envelope; and wherein the controller is further configured to recursively perform the following until no subsequent maximum values are detected: obtaining an additional threshold envelope based on an additional maximum value and an additional target angle; modifying the cumulative threshold envelope based on the additional threshold envelope; and detecting a subsequent maximum value in a portion of the received beam response spectrum being greater than the cumulative threshold envelope; wherein the second maximum value is used as an initial value of the additional maximum value and the second target angle is used as an initial value of the additional target area; identifying an angle corresponding to the second maximum value as a second target angle.
5. The radar system of claim 1, wherein the threshold envelope is calculated based on an arbitrary beam response of the radar system at the first target angle and is a function that, for all angles, has a value greater than a value of the arbitrary beam response.
6. The radar system of claim 1, wherein the threshold envelope is calculated based on an arbitrary beam response of the radar system at the first target angle and is calculated by performing smoothing with a Gaussian window over the arbitrary beam response and adding a bias value.
7. The radar system of claim 1, wherein the threshold envelope is calculated based on an arbitrary beam response of the radar system at the first target angle, the arbitrary beam response comprising a plurality of lobes, each lobe comprising a maximum value and a width; and wherein the threshold envelope is a function such that, for each lobe of the arbitrary beam response, the threshold envelope has a value equal to the maximum value of the lobe plus a bias value for all angles corresponding to the width of the lobe.
8. The radar system of claim 1, wherein the threshold envelope monotonically increases for angles less than the first target angle and monotonically decreases for angles greater than the first target angle.
9. The radar system of claim 1, further comprising: a storage device storing a plurality of stored threshold envelopes, each stored threshold envelope of the plurality of stored threshold envelopes being associated with a different angle of arrival, based on a confusion matrix given by the following equation:
ConfusionMat=CalibrationMat.Math.CalibrationMat.sup.H; wherein ConfusionMat is the confusion matrix, CalibrationMat is a calibration matrix of the radar system, and CalibrationMat.sup.H is a Hermitian transpose of the calibration matrix; wherein the plurality of stored threshold envelopes comprises, for each angle of arrival θ, a stored threshold envelope calculated by: generating an arbitrary beam response B.sub.θ given by the equation B.sub.θ=|A.sub.θ|.sup.2, where A.sub.θ=ConfusionMat[θ,:]; generating an array E equal to the arbitrary beam response B.sub.θ; identifying an index M of a maximum value of the arbitrary beam response B.sub.θ; for each element of array E from E.sub.1 to E.sub.M−1, setting E.sub.i+1 equal to E.sub.i if E.sub.i+1 is less than E.sub.i, i being an index of array E; for each element of array E from E.sub.I to E.sub.M+1, I being a length of array E, setting E.sub.i equal to E.sub.i+1 if E.sub.i+1 is greater than E.sub.i; for each element of array E from E.sub.N to E.sub.M−1, N being a distance between nulls in the arbitrary beam response B.sub.θ, setting elements of array E from E.sub.i−N+1 to E.sub.i equal to E.sub.i+1 if E.sub.i+1 is equal to E.sub.i; for each element of array E from E.sub.I−N to E.sub.M+1, setting elements of array E from E.sub.i+N to E.sub.i equal to E.sub.i+1 if E.sub.i+1 is equal to E.sub.i; setting element of array E from E.sub.M−2N to E.sub.M+2N equal to E.sub.M; adding a bias value to each element of array E, setting array E as the stored threshold envelope for angle θ; wherein the controller is configured such that the obtaining a threshold envelope based on the first maximum value and the first target angle comprises: identifying a first stored threshold envelope among the plurality of stored threshold envelopes having an associated angle of arrival corresponding to the first target angle; setting the first stored threshold envelope as the threshold envelope.
10. A method for identifying targets in a beam response of a radar system, the method comprising: receiving a reflected radar signal with a radar receiver; calculating a received beam response spectrum based on the received reflected radar signal, the received beam response spectrum comprising a set of power values as a function of angle of arrival; detecting a first maximum value of the received beam response spectrum; identifying an angle corresponding to the first maximum value as a first target angle; obtaining a threshold envelope based on the first maximum value and the first target angle; detecting a second maximum value in a portion of the received beam response spectrum being greater than the threshold envelope; identifying an angle corresponding to the second maximum value as a second target angle; and outputting the first target angle as an angle of arrival of the reflected radar signal from a first target and the second target angle as an angle of arrival of the reflected radar signal from a second target.
11. The method of claim 10, wherein the received beam response spectrum comprises a target main lobe including the first maximum value; the threshold envelope is a function having a value equal to the first maximum value for angles corresponding to a width of the target main lobe, and, for all other angles, having a value equal to a highest sidelobe peak level among sidelobes immediately adjacent to the target main lobe in the received beam response spectrum.
12. The method of claim 10, further comprising: previously storing a plurality of stored threshold envelopes, each stored threshold envelope of the plurality of stored threshold envelopes being associated with a different angle of arrival; wherein the obtaining a threshold envelope based on the first maximum value and the first target angle comprises: identifying a first stored threshold envelope among the plurality of stored threshold envelopes having an associated angle of arrival corresponding to the first target angle; setting the first stored threshold envelope as the threshold envelope.
13. The method of claim 10, further comprising: generating a cumulative threshold envelope initially equal to the threshold envelope; and recursively performing the following until no subsequent maximum values are detected: obtaining an additional threshold envelope based on an additional maximum value and an additional target angle; modifying the cumulative threshold envelope based on the additional threshold envelope; and detecting a subsequent maximum value in a portion of the received beam response spectrum being greater than the cumulative threshold envelope; wherein the second maximum value is used as an initial value of the additional maximum value and the second target angle is used as an initial value of the additional target area.
14. The method of claim 10, wherein the threshold envelope is calculated based on an arbitrary beam response of the radar system at the first target angle and is a function that, for all angles, has a value greater than a value of the arbitrary beam response.
15. The method of claim 10, wherein the threshold envelope is calculated based on an arbitrary beam response of the radar system at the first target angle and is calculated by performing smoothing with a Gaussian window over the arbitrary beam response and adding a bias value.
16. The method of claim 10, wherein the threshold envelope is calculated based on an arbitrary beam response of the radar system at the first target angle, the arbitrary beam response comprising a plurality of lobes, each lobe comprising a maximum value and a width; and wherein the threshold envelope is a function such that, for each lobe of the arbitrary beam response, the threshold envelope has a value equal to the maximum value of the lobe plus a bias value for all angles corresponding to the width of the lobe.
17. The method of claim 10, wherein the threshold envelope monotonically increases for angles less than the first target angle and monotonically decreases for angles greater than the first target angle.
18. The method of claim 10, further comprising: calibrating the radar system to generate a calibration matrix; calculating a confusion matrix based on the calibration matrix according to the following equation:
ConfusionMat=CalibrationMat.Math.CalibrationMat.sup.H; wherein ConfusionMat is the confusion matrix, CalibrationMat is the calibration matrix, and CalibrationMat.sup.H is a Hermitian transpose of the calibration matrix; calculating a plurality of stored threshold envelopes based on the confusion matrix, each stored threshold envelope of the plurality of stored threshold envelopes being associated with a different angle of arrival; wherein the calculating a plurality of stored threshold envelopes comprises, for each angle of arrival θ, calculating a stored threshold envelope, wherein the calculating the stored threshold envelope comprises: generating an arbitrary beam response B.sub.θ given by the equation B.sub.θ=|A.sub.θ|.sup.2, where A.sub.θ=ConfusionMat[θ,:]; generating an array E equal to the arbitrary beam response B.sub.θ; identifying an index M of a maximum value of the arbitrary beam response B.sub.θ; for each element of array E from E.sub.1 to E.sub.M−1, setting E.sub.i+1 equal to E.sub.i if E.sub.i+1 is less than E.sub.i, i being an index of array E; for each element of array E from E.sub.I to E.sub.M+1, I being a length of array E, setting E.sub.i equal to E.sub.i+1 if E.sub.i+1 is greater than E.sub.i; for each element of array E from E.sub.N to E.sub.M−1, N being a distance between nulls in the arbitrary beam response B.sub.θ, setting elements of array E from E.sub.i−N+1 to E.sub.i equal to E.sub.i+1 if E.sub.i+1 is equal to E.sub.i; for each element of array E from E.sub.I−N to E.sub.M+1, setting elements of array E from E.sub.i+N to E.sub.i equal to E.sub.i+1 if E.sub.i+1 is equal to E.sub.i; setting element of array E from E.sub.M−2N to E.sub.M+2N equal to E.sub.M; adding a bias value to each element of array E, setting the stored threshold envelope to array E; wherein the obtaining a threshold envelope based on the first maximum value and the first target angle comprises: identifying a first stored threshold envelope among the plurality of stored threshold envelopes having an associated angle of arrival corresponding to the first target angle; setting the first stored threshold envelope as the threshold envelope.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) Other features, advantages and details appear, by way of example only, in the following detailed description, the detailed description referring to the drawings in which:
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
DETAILED DESCRIPTION
(14) The following description is merely exemplary in nature and is not intended to limit the present disclosure, its application or uses. It should be understood that throughout the drawings, corresponding reference numerals indicate like or corresponding parts and features.
(15) In radar systems, a CFAR threshold may be selected so as to exclude all points in the detected beam response spectrum below a predetermined value when identifying detection targets as points in the detected beam response that are above the predetermined value. In other words, a CFAR threshold may be a constant value for all angles of arrival. However, applying an angle-invariant threshold such as a CFAR threshold may result in the false identification of side lobes of a beam response as a valid target, or may exclude secondary detected targets if the power of the beam response for the secondary detected targets is below the CFAR threshold.
(16) Accordingly, as an improvement over using an invariant CFAR threshold, it may be desirable to provide a radar system and method that uses a threshold envelope that varies as a function of angle of arrival according to one or more embodiments. This would allow for reduction of false identifications of targets and as well as improvement of detection of secondary targets that sometimes occurs with invariant thresholds such as a CFAR threshold.
(17) In accordance with an exemplary embodiment,
(18)
(19)
(20) Returning to
(21) In block 116 of
(22) In block 118 of
(23) In block 124, controller 54 may output first target angle 74 as an angle of arrival of a reflected radar signal from a first detection target and output second target angle 82 as an angle of arrival of a reflected radar signal from a second detection target.
(24) It will be understood that method 102 shown in
(25) In block 1220, it is determined whether there is a portion of received beam response spectrum greater than the cumulative threshold envelope. If there is a portion greater than the cumulative threshold envelope (“yes” in block 1220), then the method proceeds to block 1222. A “yes” result in block 1220 corresponds to a situation where there are still additional targets to detect in the beam response spectrum. If there is no portion greater than the cumulative threshold envelope (“no” in block 1224), then the method proceeds to block 1224. A “no” result in block 1220 corresponds to a situation where there are no further targets to detect in the beam response spectrum.
(26) In block 1222, a new maximum value and angle corresponding to the maximum value are detected from the portion of received beam response spectrum greater than the cumulative threshold envelope. The method then returns to block 1216 in which the new maximum value and angle corresponding to the maximum value are used to obtain a new interim threshold envelope. In block 1224, the method ends and all of the detected maximum values and their corresponding angles are reported as angles of arrival of detection targets.
(27)
(28)
(29) Returning to
(30) Similar to the method illustrated in
(31) In block 204, controller 54 (see
(32)
(33) In block 304, a confusion matrix is calculated according to the equation:
ConfusionMat=CalibrationMat.Math.CalibrationMat.sup.H; (1)
where ConfusionMat is the confusion matrix, CalibrationMat is a calibration matrix of the radar system 40 (see
(34)
B.sub.θ=|A.sub.θ|.sup.2; (2)
where A.sub.θ=ConfusionMat[θ,:]. The notation [θ,:] indicates all values of the confusion matrix corresponding to angle of arrival θ. In block 404, a monotonous transformation is performed on arbitrary beam response B.sub.θ. Calculation of the monotonous transformation is described in detail herein. In block 406, envelope widening is performed on the monotonous transformation calculated in block 404. Calculation of the envelope widening is described in detail herein. In block 408, a bias is added to the calculated envelope to account for noise. It will be understood that by adding the bias in block 408, the calculated envelope will have a value greater than a value of the arbitrary beam response at all angles.
(35) The following Table 1 shows an exemplary embodiment of an algorithm for performing a monotonous transformation on arbitrary beam response Bo (Algorithm 1) and an exemplary algorithm for performing envelope widening on the result of the monotonous transformation (Algorithm 2).
(36) TABLE-US-00001 TABLE 1 Algorithm 1 Monotonous Transformation Require: B,I = len(B) E = B M = max(B) for i = 1 : M − 1 do Left Side if E[i + 1] < E[i] then E[i + 1] = E[i] end if end for for i = I : −1 : M + 1 do
Right Side if E[i + 1] > E[i] then E[i] = E[i + 1] end if end for return E Algorithm 2 Envelope Widening Require: E,I = len(E),N for i = N : M − 1 do
Left Side if E[i + 1] == E[i] then E[i − N + 1 : i] = E[i + 1] end if end for for i = I − N : −1 : M + 1 do
Right Side if E[i + 1] == E[i] then E[i + N : −1 : i] = E[i + 1] end if end for E(M − 2N :M+ 2N) = E(M)
Main Lobe return E
(37) In Algorithm 1 in Table 1, left side monotonous transformation is first performed by setting index i equal to 1. Next, it is determined whether E.sub.i+1 is less than E.sub.i. The notation E.sub.x indicates a value of array E at index x. If E.sub.i+1 is less than then Algorithm 1 sets E.sub.i+1 equal to E. If E.sub.i−1 is not less than E.sub.i, then E.sub.i+1 is left unchanged. Next, it is determined whether index i is equal to M−1. If i is equal to M−1, then the left side monotonous transformation ends and Algorithm 1 proceeds to the right side monotonous transformation described herein. If i is not equal to M−1, Algorithm 1 increments index i by 1 and loops back.
(38) In the right side monotonous transformation, Algorithm 1 sets index i is equal to 1, which is a length of array E. Next, it is determined whether E.sub.i+1 is greater than E.sub.i. If E.sub.i+1 is greater than E.sub.i then Algorithm 1 sets E.sub.i equal to E.sub.i+1. If E.sub.i+1 is not greater than then E.sub.i, remains unchanged. Next, it is determined whether index i is equal to M+1. If i is equal to M+1, then the right side monotonous transformation ends. If i is not equal to M+1, then Algorithm 1 decrements index i by 1 and loops back.
(39) In Algorithm 2 of Table 1, a value N is determined, which is a distance between nulls in the arbitrary beam pattern Bo. Algorithm 2 performs left side envelope widening by first setting index i equal to N. Next, it is determined whether E.sub.i+1 is equal to E.sub.i. If E.sub.i+1 is equal to then Algorithm 2 sets all elements of array E from E.sub.i−N+1 to E.sub.i equal to E.sub.i+1. If E.sub.i+1 is not equal to then no change is made. Next, it is determined whether index i is equal to M−1. If i is equal to M−1, left side envelope widening ends and Algorithm 2 proceeds to right side envelope widening. If i is not equal to M−1, Algorithm 2 increments index i by 1 and loops back.
(40) In the right side envelope widening, Algorithm 2 first sets index i to be equal to I−N. Next, it is determined whether E.sub.i+1 is equal to E.sub.i. If E.sub.i+1 is equal to E.sub.i, Algorithm 2 sets all elements of array E from E.sub.i+N to E.sub.i to be equal to E.sub.i+1. If E.sub.i+1 is not equal to E.sub.i then no change is made. Next, it is determined whether index i is equal to M+1. If i is equal to M+1 then the right side envelope widening ends. If i is not equal to M+1, then Algorithm 2 decrements index i by 1 and the method loops back. For the main lobe of array E, Algorithm 2 performs envelope widening by setting values for E.sub.M−2N through E.sub.M+2N equal to E.sub.M.
(41)
(42)
(43)
(44) Determining the simplified threshold envelope 170 requires less processing resources than determining threshold envelope 76 or threshold envelope 84. For example if the plurality of stored threshold envelopes 76.sub.θ1, 76.sub.θ2, . . . 76.sub.θn (see
(45) The exemplary embodiments described above result in significant advantages. For example, the exemplary embodiments use an angle dependent threshold to make it possible to detect multiple targets within a beam response that may otherwise be missed if an invariant CFAR threshold was applied to the beam response. Additionally, the exemplary embodiments help to prevent false identification of side lobes as detection targets. This allows a radar system to provide more accurate information regarding angles of arrival of detection targets, thereby resulting in a more accurate radar system.
(46) While the above disclosure has been described with reference to exemplary embodiments, it will be understood by those skilled in the art that various changes may be made and equivalents may be substituted for elements thereof without departing from its scope. In addition, many modifications may be made to adapt a particular situation or material to the teachings of the disclosure without departing from the essential scope thereof. Therefore, it is intended that the present disclosure not be limited to the particular embodiments disclosed, but will include all embodiments falling within the scope thereof.