Aiming to solve the problem that the joint range and azimuth super-resolution algorithm of vehicle millimeter wave radar is too complex to be implemented quickly, a low complexity joint super-resolution method based on direct selection of frequency domain data is presented. The algorithm first transforms the space-time range-domain joint data into frequency domain by fast Fourier transformation, and stores and processes the two-dimensional frequency domain data of the area of interest. Based on the equivalence between Fourier transformation and beam space transformation based on DFT transformation, the range-azimuth joint MUSIC super-resolution in frequency domain data is achieved, and the fast joint estimation of target information is completed. The orthogonality of frequency domain subspace and the theory of frequency domain beam dimension reduction super-resolution algorithm are deduced. The relationship between the resolution and estimation performance of distance and azimuth of the algorithm and signal-to-noise ratio is simulated. The simulation results show that the accuracy and resolution of the algorithm are much higher than traditional FFT, and the computational complexity of the algorithm is greatly reduced compared with traditional MUSIC.
accuracy and resolution of the traditional Fast Fourier Transform (FFT) algorithm is insufficient. And super-resolution algorithms such as Multiple Signal Classification (MUSIC) algorithm and Estimating Signal Parameter via Rotational Invariance Technologies (ESPRIT) algorithm have high accuracy and resolution but huge computation.
algorithm is an urgent problem to be solved at present, which can not be avoided in engineering. Bienvenu G[1] et al. proposed a high-resolution target bearing estimation method to improve accuracy and resolution and even the statistical stability[3]-[4].
range azimuth joint super-resolution method based on frequency domain beam dimensionality reduction. Based on the equivalence of range angle FFT and multi-dimensional Discrete Fourier Transform (DFT) beam space transformation of space-time data, the corresponding frequency domain data area of the target area is selected according to the prior information of the target, MUSIC joint super-resolution based on beam dimension reduction. The dimension of beam space data is greatly reduced, which makes it possible to realize joint super resolution engineering and realize fast joint estimation of multiple parameters of vehicle borne radar. Compared with the traditional FFT algorithm, the resolution of target parameter estimation is significantly improved, which is a fast joint super-resolution algorithm feasible in engineering. 2. TDM-MIMO LFMCW Radar Signal Model
（TDM-MIMO）radar as LTX , and the number of receiving antenna elements as LRX . The simplified array model is shown in Figure 1(a). Set the spacing of receiving antenna elements as dr =λ 2 , the spacing of transmitting antenna elements as dt , and meet the requirements of dt =LRX × dr . Assume that the spacing between any receiving antenna lRX and the first receiving antenna in the receiving antenna array is drlRX = (lRX −1) dr .
Receiving array RxLRx Tx1
dt Tx1 ...
TxLTx
Rx1 dr ...
LRx ...
TxLTx ... ...
LRx ... f B f fIF
Tm τ
t
Virtual array：L=LTx×LRx t (a)TDM-MIMO radar array model. (b) Sawtooth time-frequency diagram.
Under the above model, obtain TDM-MIMO LFMCW radar sawtooth beat signal model[5] x(t) as x (t ) = Ampe j2π 2vcfc + 2μcR + 2μvmTm t+ 2Rcfc + 2 fcvmTm +ldr sinθ
c c λ + G (t )
In the formula, Amp represent the signal amplitudes, fc represents the signal carrier frequency, μ = B Tm represents the FM slope, where B represents the signal bandwidth, Tm represents the signal repetition period, m = 0,1...., M −1 represents the sequence number of repetition period, l = 0,1...., L −1 represents the array element sequence number of virtual array receiving antenna, dr = λ 2 is the virtual antenna spacing. G(t) is additive white Gaussian noise.
aR ( R) = 1, e− j2π 2cμfsR ,, e− j2π 2cμfsR(N −1)
T
aθ (θ ) = 1, e − j2π d sinθ λ ,, e − j2π d sinθ (L−1) T λ 3. Beam space transformation in frequency domain based on time-frequency equivalence
The beam space dimension reduction super-resolution algorithm based on DFT transform is to
obtain the beam space data by multiplying the original sampling data and the beam space conversion
matrix. At the same time, the steering vector also reduces the dimension according to the beam
selection. In this method, the multi-dimensional data are first FFT transformed, and only the
frequency domain data corresponding to the parameter region of interest are stored; Only when the
(
According to reference [6] and space-time equivalence, the beam conversion matrix in the beam
space array flow pattern of each dimension is defined. The beam transformation matrices defining the
angle and range dimensions are WθH ∈ Lb×L ,WRH ∈ Nb×N respectively, and the number of effective
beam selections in the angle and range dimensions are Lb , Nb respectively, that is, the reduced
dimension data length. Then the l and n elements are
wθH = 1 e− jl 2Lπb e− j(Lb −1)l 2Lπb wRH = 1 e− jn 2Nπb ,, e− j(Nb −1)n 2Nπb (
In the formula, l, 0 ≤ l ≤ ( Lb −1) , Lb ≤ L , when Lb = L , that is, WθH is the angle dimension full beam transformation matrix. 0 £ n £ (Nb - 1), Nb £ N , when Nb = N , WRH is the range dimension full beam transformation matrix.
Assuming that the matrices of one-dimensional Fourier transform are FN and FL respectively, where N、L respectively represent the length of the vector to be operated in the distance dimension and angle dimension, and data X NL is the distance angle two-dimensional raw data of N × L , then the frequency domain data matrix YNL obtained by two-dimensional Fourier transform is expressed as:
YNL = FN X NL FLT
equivalent to the two-dimensional Fourier transform. Eq. (
The above proves the equivalence of beam space dimensionality reduction and direct selection of
frequency domain data. For the beam space super-resolution algorithm, literature [6] has completed
the one-dimensional MUSIC proof based on the beam space. The beam space based MUSIC
algorithm can be expanded from one-dimensional to two-dimensional through equation (
same to do the left first and the right first. Two dimensional Fourier transforms are two onedimensional Fourier transforms, and they are independent of order. The same is true for higher dimensional Fourier transforms. If x1,, xL is the column vector of X NL , X NL can be written as
L X NL = xieiT transformation, where ei is the ith column of the identity matrix IL , so i
L L L vec (YNL ) = vec ( FN X NL FLT ) = vec FN xieiT FLT = vec ( FN xieiT FLT ) = vec{( FN xi )( FLei )T } i=1 i=1 i=1
L L = FLei ⊗ FN xi = ( FL ⊗ FN ) (ei ⊗ xi ) = ( FL ⊗ FN ) vec ( X NL )
i=1 i=1
vec ( Z ) = (WRH ⊗WθH ) vec ( X NL )
(
b (θ, R ) = (WRH aR ) ⊗ (WθH aθ ) = bR ⊗ bθ
Rbeam−2MUSIC =
1 Nb Lb
ZZ H
(
After eigenvalue decomposition, the noise subspace Un−beam−2MUSIC can be obtained, and the spectral peak search function is obtained as follows:
Pbeam−2MUSIC = bH (θ, R )Un−beam−2MUSICUnH−beam−2MUSICb (θ, R) 1
The spectral peak is obtained by searching the spectral peak of the above equation, and the corresponding information of the peak is the two-dimensional information of the range azimuth angle of the target. Compared with the traditional two-dimensional MUSIC algorithm, the MUSIC algorithm after the beam space reduces the data dimension in terms of covariance estimation, eigenvalue decomposition and spectral peak search, greatly reducing the calculation time.Based on the joint super-resolution in the frequency domain, the data storage pressure and computational complexity have been significantly reduced. According to existing research, the resolution and parameter estimation performance of the beam space dimension reduction algorithm under the condition of reasonable beam selection.
platform transmits FMCW signals with a bandwidth of 150MHz. The period of a Chirp is 16 us, the number of snapshots is 300, the number of transmitting and receiving antennas is 2 and 15 respectively. The search points of azimuth dimension and distance dimension are ls = 89, ns = 236 . The number of wave beams in angle dimension and distance dimension is Lb = 8 and Nb = 13 respectively. Set the parameters as Target1: (110 m, 12 °); Target2: the (110.5 m ,15 °). 4.1 Effective Estimation Diagram of Target Information
When plotting SNR = 0 dB , the range velocity dimension effectiveness estimation diagram of the range azimuth joint frequency domain beam reduction MUSIC algorithm is compared with the effectiveness estimation diagram of the traditional 2DFFT algorithm, as shown in Figure 2 and Figure
) B (edud lti p m
A (a) The proposed algorithm spectral peak estimation (b) Effective estimation of the proposed algorithm (a) 2DFFT algorithm joint information estimation spectral peak (b) 2DFFT algorithm joint information effective estimation
effectively realize the effective estimation of target information. The comparison between Figure 3 and Figure 2 shows that under this condition, 2DFFT cannot distinguish two targets. The algorithm proposed in this paper can achieve effective resolution of two targets and more accurate parameter estimation. 4.2 Target Information Distribution Map
Through 50 Monte Carlo experiments, parameters distribution diagram of the proposed algorithm in this paper and 2DFFT algorithm are compared as in Figure 4. 14.5 14 13.5 13 12.5 120 Azimuth true value 1 Azimuth true value 2 Proposed algorithm azimuth estimate1 Proposed algorithm azimuth estimate2 2DFFT azimuth estimate
Azimuth true value 1 Azimuth true value 2 Proposed algorithm azimuth estimate1 Proposed algorithm azimuth estimate2 2DFFT azimuth estimate (a) Distance dimension target information distribution.
(b) Azimuth dimension target information distribution.
4.3 Performance Analysis - RMSE Statistics
Set the echo signal to noise ratio variation range of the signal as −30 dB:10 dB:20 dB , and the number of Monte Carlo is 50. The traditional 2DFFT algorithm and the proposed algorithm proposed in this paper are used to analyze the estimation error of Target1. The experimental results are as follows. 1.2 1 0 -30 -25 -20 -15 -10 SNR-5(dB) 0 5 10 15 20 0 -30 -25 -20 -15 -10 SNR-5(dB) 0 5 10 15 20 (a) Distance dimension RMSE.
(b) Azimuth dimension RMSE.
proposed algorithm proposed in this paper has higher estimation accuracy than the traditional 2DFFT algorithm in terms of target range and angle dimensions, and its estimation performance improves with the increase of signal-to-noise ratio, and the range dimension estimation accuracy is slightly higher than the angle dimension estimation accuracy.
amount is (( NL)) after data acquisition. The algorithm proposed in this paper can reduce the data storage amount to (( Nb Lb )) , which can reduce the data storage amount by (30 × 200) (8×13) = 57.7 . ae 3 ö Based on MUSIC joint super-resolution algorithm, the computational complexity includes ççç(NL) ÷÷ è ÷ø of eigenvalue decomposition processing part and ((NL) nl ) of spectral peak search part. The total ae 3 ö complexity of the algorithm is ççç(NL) +(NL) nl÷÷÷ . The complexity of the proposed algorithm, è ø ae 3 ö eigenvalue processing and spectral peak search is reduced to ççç(NbLb ) ÷÷÷ø and ((NbLb ) nl ) è respectively, and t Total computational complexity reduced to 6.5×104 . When the selected target area is small, the number of beams can also be smaller to further reduce the complexity.
In this paper, a range azimuth joint super-resolution algorithm based on frequency domain beam dimensionality reduction for vehicle borne radar is proposed. Compared with the traditional 2DFFT algorithm, the algorithm can effectively improve the resolution and estimation accuracy. Compared with the general MUSIC algorithm, the dimensionality reduction of the frequency domain space-time beam greatly reduces the data storage pressure and the scale of data processing, realizes the multilevel reduction of data storage and computational complexity, facilitates the realization of engineering applications, and provides effective technical support and reasonable solutions for the solution of the vehicle borne radar target accurate detection problem.