US 12,235,344 B2
Radar target detection method based on estimation before detection
Benzhou Jin, Jiangsu (CN); Yutong Shen, Jiangsu (CN); Jianfeng Li, Jiangsu (CN); Xiaofei Zhang, Jiangsu (CN); and Qihui Wu, Jiangsu (CN)
Assigned to NANJING UNIVERSITY OF AERONAUTICS AND ASTRONAUTICS, Jiangsu (CN)
Filed by NANJING UNIVERSITY OF AERONAUTICS AND ASTRONAUTICS, Jiangsu (CN)
Filed on Apr. 29, 2022, as Appl. No. 17/732,531.
Application 17/732,531 is a continuation of application No. PCT/CN2021/104677, filed on Jul. 6, 2021.
Claims priority of application No. 202010736841.4 (CN), filed on Jul. 28, 2020.
Prior Publication US 2023/0059515 A1, Feb. 23, 2023
Int. Cl. G01S 13/50 (2006.01); G01S 13/58 (2006.01)
CPC G01S 13/505 (2013.01) [G01S 13/583 (2013.01)] 4 Claims
OG exemplary drawing
 
1. A radar target detection method based on estimation before detection (EBD), comprising:
1) obtaining baseband signal by a receiver, and performing pulse compression and coherent integration on the baseband signal to obtain a range-Doppler map, performing pre-detection based on the range-Doppler map to obtain interested pre-detect targets (PDTs), wherein the corresponding ranges and Doppler frequencies of cells, wherein the PDTs are pre-detected, are represented by rζ and fζ, respectively;
2) estimating ranges and the Doppler frequencies of the PDTs, wherein the estimates are represented by custom character and custom character;
3) establishing a dimension-reduction observation model of a received signal based on custom character and custom character;
4) reconstructing a target vector in the dimension-reduction observation model based on a sparse recovery algorithm; and
5) adopting a generalized likelihood ratio detector for target detection based on a reconstruction result and outputting target detection results and their parameters, and determining whether the PDTs are true targets based on the target vector,
wherein the step 2) comprises:
adopting custom character and custom character to represent true ranges and Doppler frequencies of the PDTs, respectively, and letting custom character=[custom character; custom character] and θζ=[rζ;fζ], wherein the symbol; in the square brackets represents connecting two vectors, the received signal is represented as:
y=custom character(custom character)β+n  (7)
y in the equation (7) represents a received signal of one coherent processing interval, wherein β is a dimension-reduction target vector, an ith element βi represents a true complex amplitude of an ith PDT, custom character is an observation matrix, and n is an additive white Gaussian noise vector;
based on a maximum likelihood criterion, estimates of custom character and β are given as:

OG Complex Work Unit Math
wherein θ=[r;f], a minimum over β is attained for:
β=(A(θ)HA(θ))−1A(θ)Hy  (9)
then, minimizing a cost function in (8) is equivalent to minimizing the function:
g(θ)=∥y−A(θ)(A(θ)HA(θ))−1A(θ)Hy∥22,  (10)
obviously, a minimum value of the equation (10) is obtained at θ=custom character, the derivative of g(θ) at custom character is evaluated by the first-order Taylor series:
θg(θ)≈∇θg(custom character)+∇θ2g(custom character)(θ−custom character),  (11)
obviously, ∇θg(custom character)=0, then:
custom character≈θ−(∇θ2g(custom character))−1(∇θg(θ))  (12)
θ is replaced with θζ to obtain an estimate of custom character:
custom character≈θζ−(∇θ2gζ))−1(∇θgζ))  (13).