The problems of the “Three Crossings” motions and compensations for the moving target of the high speed moving MIMO radar are studied. Firstly, the space model is established to describe the problems by educing the equations of delay and transmitting pattern, which are changing with time. The echo characteristics are analyzed and the formulas of the “Three Crossings” are given. Secondly, a compensation method, which uses the fractional Fourier transform (FrFT) to compensate the crossing Doppler-cell and the preprocess to compensate the crossing range-cell and crossing-beam, is proposed. This method could compensate the “Three Crossings” simultaneously with little complexity of calculation because the preprocess deals only with the transmitting signals. Lastly, the phantom antenna array is employed in the simulations, the working model of “sparse transmitting dense receiving” is applied to the simulation of crossing-beam, and the results of simulations demonstrate the validity and the practicability of the method of combining the preprocess and the FrFT to compensate “Three Crossings.”
1. Introduction
The aircrafts are moving faster and faster as the development of the aviation technology, some of them could attain a velocity more than mach 10, and the maneuverability is becoming stronger too, which make the moving targets much more difficult to be searched, detected, and tracked by radar. Long time accumulation could increase the SNR output and the performance of detection, but adding accumulation time may cause “Three Crossings”, which means crossing Doppler-cell, crossing range-cell, and crossing-beam. The “Three Crossings” ruin the coherence of the echo. In order to accumulate coherently to get the best detection ability, the “Three Crossings” must be compensated before the accumulation.
Some valuable studies have been proposed for the crossing Doppler-cell and crossing range-cell [1–8], but most of them could compensate only one of these crossings. The narrow beam of the phase array (PA) radar limited the study about the crossing-beam, for PA radar cannot compensate the crossing-beam; so far there is no reference for this study. The MIMO radar transmits orthogonal signals to form a low amplitude but wide beam to cover the whole surveillance area [9–11]. Even though the “Three Crossings” happen simultaneously, the MIMO radar can receive all the echoes without any energy loss. So, the MIMO radar produces the qualification for the study of the “Three Crossings,” especially the crossing-beam.
In this paper, the “Three Crossings” of the target detected by a high-speed moving MIMO radar are studied. A space model is established firstly to educe the equations about how the delay and the transmitting pattern change with the time. Secondly, an analysis of the echo is proposed to show the “Three Crossings” of the target. Thirdly, the preprocessing method [12] is employed to compensate the crossing range-cell and the crossing-beam, while the fractional Fourier transform (FrFT) [13] is employed to compensate the crossing Doppler-cell. A performance brief analysis about the compensation mismatching is proposed after all this, and the simulation results are produced at last.
2. Space Model
At time t, the radar lies at point P with the height of H and moves along the Y-axis with a high speed vr, while the target lies at point Q and moves with a low speed vt, assuming that vr is much larger than vt. Let the distance of the radar and the target PQ=R0; the angle between the direction of the radar and PQ is θ. After t1, the transmitting signal propagates to the point where the target lies; assume that the point is Q′, the radar moves to point P′, and the distance becomes R(t1). Because the target moves very slowly, the range QQ′ is very small, PQ and PQ′ are nearly the same, and angle ∠P′PQ and angle ∠P′PQ′ are almost equal to each other, then we can assume that PQ′=PQ=R0 and ∠P′PQ′=∠P′PQ=θ. The space geometry relationship of the MIMO radar and the detected target is shown in Figure 1.
The space geometry relationship of the MIMO radar and the detected target.
According to the cosine theorem, R(t1) can be calculated; after Taylor expandedness, the delay of the echo should be
(1)τ(t)=2R0c+-2vc-v(t-R0c)+ac2(c-v)3(t-R0c)2,
where v is the relative velocity, c is the velocity of light, and a is the acceleration. Educe a further step; (2) can be gotten as follows:
(2)t-τ(t)=t-2R0c+2vc-v(t-R0c)-ac2(c-v)3(t-R0c)2=c+vc-vt-c+vc-v2R0c+v-ac2(c-v)3t2=kt-kt0-pt2,
where
(3)k=c+vc-v,t0=2R0c+v,p=ac2(c-v)3.
Extracting the triangle PP′Q′, in the whole surveillance period, the changing angle from the target to the radar is equal to α, as is shown in Figure 2.
The changing angle.
Assuming that α is constant in one pulse repeat interval (PRI) and variable between each PRI, the MIMO radar has an M elements transmitting uniform linear array (ULA) and an N elements receiving ULA; L pulses are transmitted in the accumulation period. So, the phase matrixes of the steering vectors of the transmitting and receiving arrays are
(4)ΨT=[ψ(t1),ψ(t2),…,ψ(tL)]T=[ψ1(t1)ψ1(t2)⋯ψ1(tL)ψ2(t1)ψ2(t2)⋯ψ2(tL)⋮⋮⋱⋮ψM(t1)ψM(t2)⋯ψM(tL)]T,(5)ΨR=[Φ(t1),Φ(t2),…,Φ(tL)]T=[ϕ1(t1)ϕ1(t2)⋯ϕ1(tL)ϕ2(t1)ϕ2(t2)⋯ϕ2(tL)⋮⋮⋱⋮ϕN(t1)ϕN(t2)⋯ϕN(tL)]T,
where the phase of the mth transmitting channel is
(6)ψm(tl)=2π(m-1)dtλsin[α1(tl)],m=1~M,l=1~L
and the phase of the nth receiving channel is
(7)ϕn(tl)=2π(n-1)drλsin[α2(tl)],n=1~N,l=1~L.
The angles α1 and α2 are the α in Figure 2. Assume that, after b PRI, the angle changes, so tl should be
(8)tl=⌈lb⌉T,l=1~L,
where T could be any moment in the radar period Tr. According to the cosine theorem, the sinα can be calculated in the triangle PP′Q′; put sinα into (6); then we can get
(9)ψm(tl)=(1+vrR0⌈lb⌉Tcosθ)vrR0⌈lb⌉T·2π(m-1)dsinθλ,m=1~M,l=1~L.
If the transmitting array and the receiving array are the same, sin[α2(tl)] would have the same formula with sin[α1(tl)], so (9) can be used for ϕn(tl) too.
3. Echo Analysis
The transmitting signal can be expressed as
(10)S=[s1(t),…,sL(t)]=[s11(t),…,s1M(t)s21(t),…,s2M(t)⋮sL1(t),…,sLM(t)]T,
where the signal transmitted by the mth element at the lth period is
(11)slm(t)=rect(t-lTrTp)sbm(t-lTr)ej2πfct,
where sbm(t) denotes the transmitted baseband signal. The echo received by the nth channel after mixing should be
(12)sr_ln(t)=∑m=1Mrect[t-τ(t)-lTrTp]sbm[t-τ(t)-lTr]×e-j2πfcτ(t)ejψm(tl)·ejϕn(tl)+vn(t).
In this paper, the clutter is ignored to focus on the “Three Crossings” and their compensations; the clutter compression in this case of the airborne MIMO radar, including the clutter diffusing and agglomerating, will be discussed in another paper, so vn(t) in (12) only stands for the noise.
Name the rectangle function in (12) as f(t), so
(13)f(t)=rect[t-τ(t)-lTrTp]=rect(kt-kt0-pt2-lTrTp).
Let the numerator in (13) be 0; the front edge in the fast-slow time domain of the echo can be calculated as
(14)t^(l)=pTr2k3l2+(1k-1-2pTrt0k2)lTr+(2k-1)t0+pt02k.
The echo envelope will move in the fast-slow time domain as the number of the pulse increases, according to (14), so the range migration would happen in this case.
Name the baseband signal and the first exponential function in (12) as g(t), so
(15)g(t)=∑m=1Msbm[t-τ(t)-lTr]e-j2πfcτ(t).
Here, the delay τ(t) is a quadric polynomial and the Doppler-frequency of the echo fd(t) would be a function of time, which would make the coherence of the echo be destroyed; that is to say, the Doppler migration would happen in this case.
Name the last two exponential functions in (12) as h(t), so
(16)h(t)=∑m=1Mejψm(tl)·ejϕn(tl).
According to (4)~(9), we have
(17)ψm(tl1)≠ψm(tl1+b),ϕn(tl1)≠ϕn(tl1+b).
That is to say, after b PRI, the target moves into the next beam and the crossing-beam would happen in this case.
4. Compensation Methods
In this paper, the preprocess method and the FrFT are employed to compensate the “Three Crossings.” The preprocess method compensates the range migration by combining the direct digital synthesis (DDS) technology to control the transmit signal elaborately, that is, to regulate the pulse width of the transmitting signal to make the echo envelopes be aligned automatically in the fast-slow time domain. The beam-crossing is compensated by preprocessing the transmitting signal too. Based on the algorithm of the preprocess, the conjugate of the transmitting pattern is added in the transmitting signals to make the transmitting angles be compensated at the target.
From p=ac2/(c-v)3≈a/c, p provides little effect on the range migration, so f(t) can be approximated as
(18)f(t)=rect(kt-kt0-lTrTP).
So, the new transmitting signal can be written as
(19)slm(t)=rect(t/k-lTrTp)sbm(tk-lTr)e-jψm(tl)ej2πfct,m=1~M.
Then, the added signal at the target should be
(20)sT_lm(t)=∑m=1Mrect[(t-t1)/k-lTrTp]sbm×[(t-t1)k-lTr]ej2πfc(t-t1),
where t1 is the transmitting delay from radar to the target. So, the echo after mixing should be
(21)sr_ln(t)=∑m=1Mrect(t^-t0TP)sbm[t-τ(t)k-lTr]×e-j2πfcτ(t)ejϕn(tl)+vn(t).
The front edge of the echo envelopes in the fast-slow time domain is a constant t0, which means the echo envelopes are aligned, so the range-cell crossing is compensated.
After the receiving beam forming, we have
(22)sr_l(t)=∑n=1Nsr_ln(t)e-jϕn(tl)=∑n=1N∑m=1Mrect(t^-t0TP)sbm[t-τ(t)k-lTr]e-j2πfcτ(t)+∑n=1Nvn(t).
In (22), the transmitting pattern is canceled by the receiving pattern, so the beam-crossing is compensated. Next, the FrFT is employed to compensate the Doppler-cell crossing, which would make the phases of the echo contain only constant and one-order polynomial; the result should be
(23)sr_l(t)=∑n=1N∑m=1Mrect(t^-t0TP)sbm(λ1t+λ0)×e-j2πfc(ξ1t+ξ0)+∑n=1Nvn(t),
where λ1 and ξ1 are the coefficients of the one-order polynomial of t and λ0 and ξ0 are constants.
In (23), the echo envelopes are aligned in the fast-slow time domain, the orders of the echo phases are less than 2, and the pattern is constant, so after the processes of the preprocess and the FrFT, the “Three Crossings” are corrected, echo pulses are coherent, and the highest accumulation output can be gotten after the matched filter and FFT.
5. Simulations
The LFM signal is employed in the simulations. The conditions are as follows: the PRI is 1 ms, the pulse width is 0.1 ms, the original range between the radar and the target is 100 km, the velocity of the radar platform is 1 km/s, θ=π/6, the number of the transmitting antennas is 2, the bandwidth is 0.2 MHz, the frequency interval of the channels is 0.2 MHz, the number of the pulses is 500, the sample frequency is 1 MHz, the carrier frequency is 2 GHz, and the zero intermediate frequency is 0.2 MHz. After a little simple calculation, the original relative velocity is vr≈866 m/s, the acceleration is 5 m/s2, during the accumulation period, and the range between the radar and the target shifts ΔR≈433.6 m.
The relative range shifts about 450 m in Figure 3(a); it is quite close to the number calculated by equation. The MTD result expands in the velocity dimension because of the Doppler shift.
The matched filter (MF) and MTD results without any compensation.
After the range compensation, the peaks of the MF results are superposition, but the MTD result is still expanding. From Figure 4, the range-cell crossing has been compensated, but the Doppler-cell crossing still exists.
The matched filter and MTD results after range compensation.
After the range and Doppler compensations are done, the MF results are totally superposition and the MTD result is an ideal thumb pin; then we can gain the highest accumulation output (see Figure 5).
The matched filter and MTD results after range and Doppler compensations.
When the beam-crossing happens, three conditions are satisfied generally; those are as follows: (1) the number of the transmitting elements is large, (2) the relative range between the radar and the target is not too far away, and (3) the vertical velocity is high. So, the relative range is reset as 20 km, the velocity of the radar is reset as 3000 m/s, the number of the transmitting elements is 10, and the PRI changes into 2 ms. The phantom element [14] is employed in the transmitting array and the interval distance is equal to the product of the number of the transmitting elements and the half of the wavelength, which could expand the transmitting caliber to 100 transmitting elements to decrease the width of the transmitting main lobe. It can be calculated that the width of the main lobe is about 0.731°, the shift of the angle in one PRI is 0.0744°, after 100 PRI, and the target moves across 2.0362 beams.
The accumulations before and after compensating the beam-crossing, on condition that the range-cell crossing and the Doppler-cell crossing have been compensated, are shown in Figure 6, and the accumulations are normalized by the accumulation result before compensating the beam-crossing. It is obvious to see that the peak after compensation is twice higher than that before, which is consistent with the case that the target moves across about 2 beams.
Accumulation results before and after the beam-crossing compensations.
To prove the validity of the precompensation method, simulations are taken under the main conditions that the numbers of the transmitting and receiving channels are M=N=4, the number of the transmitting pulses is L=200, the radar PRI is Tr=1 ms, the pulse width is Tp=50μs, the original distance of the radar and the target is R0=100 km, and the real relative velocity is v0=4 km/s. Consider that it is difficult to know the real relative velocity in practice, so in order to check the performance of the precompensation when v≠v0, four different possible relative velocities are selected as the compensation velocities and v=4,5,6,7 km/s; for the symmetry, the performances of these velocities should be the same, respectively, to v=4,3,2,1 km/s.
Firstly, the accumulation results of different compensation velocities are shown in Figure 7. It is easy to see that when the compensation velocity is more different from the real relative velocity, the accumulation result is wider and the peak is lower, which means the compensation performance is worse. When the compensation velocity is 7 km/s, the peak of the accumulation is nearly 3 dB less than the accumulation of v= 4 km/s. So, when the difference between the real relative velocity and the compensation velocity is less than 3 km/s, the accumulation performance is acceptable. At the same time, one could see that there is a good velocity difference tolerance.
Accumulations of four different compensation velocities.
Secondly, the detection probabilities of the four different compensation velocities are shown in Figure 8. Corresponding to the results in Figure 7, at the same SNR value, when the compensation velocity is more different from the real relative velocity, the detection probability is lower. If there is no compensation, the scenario is the same to the v=8 km/s or v=0 km/s and the detection performance would be even worse than the one when v=7 km/s.
Detection probabilities of four different compensation velocities.
From these two groups of simulation results, one can figure out that the precompensation can compensate the range migration valid in a tolerance of the difference between the real relative velocity and the compensation velocity. The detection probability at 90% after compensation when v=v0 is about 7 dB or 8 dB better than the one when these two velocities are 3 km/s different.
6. Conclusion
As the developments of the velocity and the maneuverability of aircrafts, the “Three Crossings” are becoming bigger problems for the MTD. Traditional radars cannot compensate the beam-crossing, so no article discussed this problem. MIMO radar, which transmits a wide beam, provides the probability to achieve the compensation of the beam-crossing. In this paper, the “Three Crossings” are discussed simultaneously based on the airborne MIMO radar system. To detect the moving target, the characteristic of the echo is analyzed, the “Three Crossings” model is proposed, and then the preprocessing method and the FrFT are employed to compensate the “Three Crossings.” The simulation results show that these two methods can compensate the “Three Crossings” validly; the accumulation peak is significantly higher after the compensations. Even when the compensation velocity is not exactly the same with the real one, the compensation performance is also acceptable within a large velocity tolerance.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors would like to thank all the editors and referees for their suggestions and advices. This work is supported by the National Natural Science Foundation of China (no. 11076006, no. 61032010, and no. 61101173).
PerryR. P.DipietroR. C.FanteR. L.SAR imaging of moving targets1999321188200LiuY.MengH.LiG.WangX.Velocity estimation and range shift compensation for high range resolution profiling in stepped-frequency radar2010747917952-s2.0-7795798442510.1109/LGRS.2010.2047492BoersY.DriessenJ. N.Multi-target particle filter track before detect application20041516351357SuJ.XingM.WangG.BaoZ.High-speed multi-target detection with narrowband radar2010445956032-s2.0-7795650563110.1049/iet-rsn.2008.0160LiuY.MengH.LiG.WangX.Velocity estimation and range shift compensation for high range resolution profiling in stepped-frequency radar2010747917952-s2.0-7795798442510.1109/LGRS.2010.2047492BoersY.DriessenJ. N.Multi-target particle filter track before detect application20041516351357GeroleoF. G.Brandt-PearceM.Detection and estimation of LFMCW radar signals20124814054182-s2.0-8485618998210.1109/TAES.2012.6129644GuanJ.ChenX. -L.HuangY.HeY.Adaptive fractional Fourier transform-based detection algorithm for moving target in heavy sea clutter201265389401RabideauD. J.ParkerP.Ubiquitous MIMO multifunction digital array radar1Conference Record of the 37th Asilomar Conference on Signals, Systems and ComputersNovember 2003105710642-s2.0-4143151023MooP. W.Multiple-input multiple-output radar search strategies for high-velocity targets2011532562652-s2.0-7925162275510.1049/iet-rsn.2010.0044HassanienA.VorobyovS. A.GershmanA. B.Moving target parameters estimation in noncoherent MIMO radar systems2012605235423612-s2.0-8485996182010.1109/TSP.2012.2187290LuoC.LiJ.LiuH. M.HeZ. S.Compensation method for envelop migration of MIMO radar high speed moving target based on transmit signal pre-processInternational Workshop on Microwave and Millimeter Wave Circuits and System Technology (MMWCST '12)2012Cheng Du, China14SunH.-B.LiuG.-S.GuH.SuW.-M.Application of the fractional Fourier transform to moving target detection in airborne SAR2002384141614242-s2.0-003682149510.1109/TAES.2002.1145767DonnetB. J.LongstaffI. D.MIMO radar, techniques and opportunitiesProceedings of the 3rd European Radar ConferenceSeptember 2006Manchester, UK1121152-s2.0-4174912546910.1109/EURAD.2006.280286