Research Article | Open Access
Property Analysis of MIMO-Based Missile-Borne Forward-Looking SAR
As a special multiple-input multiple-output (MIMO) radar networking mode, missile-borne forward-looking synthetic aperture radar (MFL-SAR) has many potential applications. This paper describes and analyzes properties of this new configuration. Range history and Doppler history are analyzed and derived using the designed geometric configuration. Then the expressions of range and Doppler resolution are determined based on the validity of two-dimensional (2D) resolution imaging capability. To help to design the proper system and motion parameters of this configuration, key parameters affecting the imaging ability are found out. Due to high velocities and accelerations of both transmitter and receiver, high-order terms in the slant range equation should be kept to reduce the approximation error. The range resolution and Doppler resolution of MFL-SAR are both space-variant and time-variant owing to the complexity of this configuration. The tiny changes of 2D resolution during the synthetic aperture time should be considered when designing the imaging algorithm of MFL-SAR.
Multiple-input multiple-output (MIMO) [1–4] radar networking system is a new radar mode in which many transmitters synchronize to transmit signals and many other receivers get target echo signals simultaneously. Proposed at the beginning of 21th century, it has been attracting much attention because of its good performance in antistealth, low probability intercept and antijamming capability. Besides, application of MIMO radar networking to radar imaging  could improve radar imaging resolution and thus can be further applied to enhance properties of targets detection, recognition, and tracking [6, 7]. As one special MIMO radar networking mode, missile-borne multistatic synthetic aperture radar (SAR) has potential applications on ground detection, battlefield reconnaissance, active homing guidance, and so forth. Electromagnetic wave transmitted by many transmitters impinges receivers after reflecting from target area and gets centralized processing [8, 9]. Such technology can achieve high-resolution imaging of high speed platform for those targets in the straight-ahead position and realize tracking and attacking with high precision, so it can be used in terminal guidance of missile [10, 11] especially.
In order to confirm the validity of application of such MIMO radar networking system to missile-borne platform, the mode of “single-transmitting and single-receiving” is assumed, that is, missile-borne bistatic forward-looking SAR (MBFL-SAR). MBFL-SAR is a new bistatic SAR imaging mode , which has advantage of 2D imaging ability in the forward-looking mode over monostatic SAR . Owing to the “far-transmitting and near-receiving” mode and forward-looking mode of such configuration, transmitter can be placed in safe area far away from target area, while attack missile receives target echo signals in electromagnet silence state, acting as receiver. MBFL-SAR can work to image and seek targets actively during the whole terminal attacking stage because of its low probability intercept and good stealth performance. As one mode of bistatic SARs, MBFL-SAR is also subject to the problems and special requirements of synchronization of the transmitter and receiver. Some investigations were conducted and appealing approaches were suggested on synchronization of bistatic SAR [14–16]. This paper concentrates on description and analysis of the properties of this new configuration and ignores the synchronization problem.
Configuration and resolution analysis are the basis of MBFL-SAR imaging study, and some researches have been published. Three different bistatic forward-looking SAR configurations have been analyzed, and the optimal geometry configurations have been proposed in . Qiu and Hu designed bistatic forward-looking geometry with stationary transmitter and forward-looking receiver and analyzed its 2D resolution in . Experiments have already been carried out to test the bistatic forward-looking imaging ability [19, 20]. Some imaging algorithms of bistatic forward-looking configuration have been suggested [21–23]. But all investigations above are focused on airborne or space-borne bistatic forward-looking configurations; a few studies about missile-borne platform have been achieved.
In this paper, we describe and analyze this special configuration. Geometric configuration and signal model of MBFL-SAR are introduced in Section 2. High-order terms appearing in the range equation cannot be ignored, because of the fact that both transmitter and receiver have high velocities and accelerations. Doppler parameters are found out and analyzed based on the Taylor series expansion of the range equation. In Section 3, we give the derivation of 2D resolution of MBFL-SAR and find out key parameters affecting the imaging performance. Several examples are provided to illustrate the resolution ability in the MBFL-SAR in Section 4. Finally, conclusions are drawn in Section 5.
2. Geometric Configuration and Signal Model
2.1. Geometric Configuration
MBFL-SAR configuration is shown in Figure 1, in which both transmitter and receiver travel in curvilinear descending motion. (here, the height of the scene is assumed to be zero) is a target point at any position of the dashed area, that is, the imaging area. Transmitter is descending along the curve in the plane which forms an included angle with plane , in which receiver is travelling along . is the origin of coordinates, and is the projection of transmitter when the slow time . For convenience of description, all motion parameters of transmitter are described in coordinates , and the ones of receiver are described in coordinates . Assume that, at , transmitter and receiver are at and , respectively. Velocity vectors of transmitter and receiver at that moment are and , respectively; acceleration vectors are and , respectively. The components in directions and are and , respectively; the components in directions and are and , respectively.
The locations of transmitter and receiver at are and ; then we have
The position of transmitter in coordinates , according to the rotation relationship between coordinates and , can be expressed as
Then we can obtain the expression of instantaneous bistatic slant range as where is the position of in coordinates ; it can be obtained using the method similar to (2) and with
2.2. Signal Model
Suppose that the transmitted waveform is the linear frequency modulation (LFM), and scattering from can be written as where is the range time, is the wavelength, is the speed of light, and is the chirp rate. and are the range and azimuth envelopes, respectively.
To design imaging algorithms for MBFL-SAR more conveniently, it is important to take efficient approximation of ; some simulations are given using the parameters listed in Table 1 and the results are shown in Figure 2.
(a) Up to quadratic term
(b) Up to the third term
It can be seen that the maximum approximation error when keeping the terms up to the quadratic term is about 0.02 m, nearly the value of wavelength, so it is necessary to keep the higher term. When keeping the terms up to the third term, the maximum approximation error is only 2.44 × 10−6 m, which is much less than the value of wavelength. Therefore, expanding (3) as a Taylor series of needs to keep the terms up to the third-order with
2.2.1. Doppler Frequency and Doppler Centroid
Doppler frequency can be obtained through derivation of slant range and is expressed as
Doppler centroid represents echo Doppler frequency at the moment beam center cross target, given by
Doppler centroid represents correlativity only to motion parameters and irrelevance to location of targets in the traditional case where monostatic SAR moves with invariant velocity along a straight line or the case where transmitter and receiver move with equal velocities along parallel trajectories in traditional airborne bistatic SAR system. However, it can be seen from (9) that Doppler centroid in MBFL-SAR conducts different performance that varies with such parameters correlative to location of targets as , , , and , so is space-variant in MBFL-SAR configuration.
2.2.2. Doppler Frequency Rate
Doppler frequency rate is the change rate of Doppler frequency and can be expressed as
We can find from (2.1) and (10) that Doppler frequency rate not only depends on the coordinates of targets but also relates strongly to original velocities and accelerations of platforms in this special configuration. Given such consumption that transmitter and receiver have equivalent velocities and invariant moving heights and accelerations are ignored and that the included angle is consumed as 0°, we can get the traditional bistatic parallel and velocity-equivalent mode. Therefore, the traditional bistatic parallel and velocity-equivalent mode is the special case of MBFL-SAR to some extent.
2.2.3. Range Walk Ratio (RWR)
Range walk ratio represents the range walk increment of echo signal in per time unit; it can be obtained as
Similar to Doppler centroid, range walk ratio in MBFL-SAR configuration is also space-variant. Such space-variance leads to the result that MBFL-SAR does not have azimuth-invariant property, which is important to conventional algorithms in monostatic SAR and airborne bistatic SAR in parallel and velocity-equivalent mode. So space-variance must be considered when the range cell migration correction is achieved.
3. Analysis of 2D Resolution
In this section, 2D imagingresolution of MBFL-SAR is analyzed. The expressions of range resolution and Doppler resolution are determined using gradient method . Then key parameters affecting the 2D imaging ability are found out, which is helpful to design system parameters and motion parameters.
3.1. Range Resolution
Assume that and represent slant range vectors from transmitter and receiver to target at , respectively, and can be written as where , , and are unit vectors in the directions axis, axis, and axis.
According to the geometrical relationship in Figure 1, we can determine the coefficients in (12) such that where and are squint angle and look angle of transmitter, respectively, and are angle deviating from the beam center and look angle of receiver, respectively, and the angles are as follows:
According to gradient method in , we can get the range resolution of MBFL-SAR: where is bandwidth of transmitted signal and is the range resolution which reflects the capability of distinguishing two targets in horizontal plane of imaging area. Given a certain bandwidth , is related to such angles between radars and target as , , , and and included angle (collectively called “space angles”). Although heights of transmitter and receiver do not appear in (16) explicitly, they can also affect through such “space angles.” In addition, the values of are different because of different space angles caused by different coordinates of targets and different slow time . Thus, is not only space-variant but also time-variant.
3.2. Doppler Resolution
As shown in the coordinate system in Figure 1, velocity vectors of transmitter and receiver platforms are expressed as where and represent speed components of transmitter and receiver in horizontal plane, respectively, and and denote speed components of transmitter and receiver in the vertical plane. Because the directions of motion of transmitter and receiver in the vertical plane are downward, speed components in vertical direction are negative.
Again, using the gradient method, the Doppler resolution can be obtained as with
Similarly, is also dependent on “space angles.” Besides, synthetic aperture time and signal wavelength can also have influence on . Similar to , is space-variant and time-variant too.
4. Numerical Results
In this section, several examples are provided to illustrate the resolution ability in the MBFL-SAR. Below are some simulation results of the resolution capability and 2D resolution. Parameters used in the simulations are listed in Table 1. Due to the features of space-variance and time-variance, we give the results at a certain time (here, ). Figure 3 depicts MBFL-SAR range contours and Doppler frequency contours, and MBFL-SAR range resolution and Doppler resolution are given in Figure 4. To test the time-variant property, some simulations of the 2D resolution of the point target located at (0, 4500, 0) at different slow time are also implemented, and the results are shown in Figure 5.
(a) Range resolution
(b) Doppler resolution
(a) Range resolution
(b) Doppler resolution
To further analyze the signal properties of MBFL-SAR, several simulations about airborne bistatic forward-looking SAR (ABFL-SAR) are also realized, where the velocities of the transmitter and receiver are 150 m/s and 100 m/s with included angle of 20 degrees, and the heights are 5 km and 2 km, respectively. Doppler frequency and RWR comparisons of five-point targets 150 m apart from each other in the imaging area are conducted between ABFL-SAR and MBFL-SAR during the same synthetic aperture time. The results are given in Figures 6 and 7, respectively. Besides, experiments on range cell curvature (RCC) of one-point target are finished and results are in Figure 8.
MBFL-SAR range contours and Doppler frequency contours when are shown in Figure 3, where range contours and Doppler frequency contours are intersectant but all of them are not perpendicular. It can be seen that MBFL-SAR configuration has the capability of 2D imaging resolution. In addition, range contours and Doppler frequency contours are nearly perpendicular in some places of the imaging area but are nearly parallel in some other places, which describe the space-variant feature of 2D resolution as discussed previously.
Figure 4 represents the results of MBFL-SAR range resolution and Doppler resolution when . Range resolution is shown in Figure 4(a), where we can find that the values of range resolution vary with the location of the targets in the imaging area. And because of the particularity of MBFL-SAR configuration, there exist isolines among different targets located in different location and location; for example, two targets located in (1350, 1510) and (2220, 3810) have equal value of 6.5759 m. Unlike range resolution, the value of Doppler resolution is very low, which can be seen from Figure 4(b). The space-variant property of Doppler resolution can also be found from the same value of 0.9113 m of two targets located in (810, 3130) and (3640, 2720), respectively. In addition, targets located in smaller location conduct better performance than ones located in greater location in MBFL-SAR configuration.
Some simulations are also achieved to analyze range resolution and Doppler resolution of MBFL-SAR at different , and the results are like Figure 5. Range resolution and Doppler resolution vary with but the changes are tiny during the synthetic aperture time, 0.044 m of range resolution and 0.016 m of Doppler resolution, respectively. This characteristic can be considered when designing the imaging algorithms of this configuration, and some approximations could be taken.
Because of the presence of high velocities and accelerations in MBFL-SAR configuration, Doppler frequencies of targets in the imaging area have greater changes than ABFL-SAR. And range cell migrations (RCMs) are also more severe during the same synthetic aperture time; that is to say, RWRs are greater. As shown in Figure 6, Doppler frequency curves in both configurations are nearly parallel linear lines, but the maximum changes of Doppler frequency in ABFL-SAR are only 131 Hz, while the one in MBFL-SAR is much greater, about 9600 Hz. Results given in Figure 7 show that RWRs in both configurations are space-variant, but the one in MBFL-SAR is more grievous. The maximum RCM in ABFL-SAR is about 16 m, while the one during the same synthetic aperture time in MBFL-SAR is 45 m, which results from the fact that the range in MBFL-SAR is much further and changes more rapidly than ABFL-SAR. Figure 8 depicts RCC comparison between ABFL-SAR and MBFL-SAR, where we can find that RCC in MBFL-SAR is more serious than ABFL-SAR, about 20 m and 2 m, respectively. The greater RCC could not be neglected in MBFL-SAR when designing imaging algorithms, while the one in ABFL-SAR can be ignored if the influence on the resolution is small enough.
Applications of MIMO radar networking technology to radar imaging can improve the performance of radar imaging resolution, targets detection, recognition, and tracking. As a special MIMO radar networking system, MBFL-SAR is a special mode of bistatic forward-looking imaging by combining geometry configuration of bistatic motion with forward-looking mode and has potential of being applied to missile precision terminal guidance. Properties of MBFL-SAR are analyzed in this paper. The high velocities and the two square-root terms lead to the presence of high-order terms in slant range history and Doppler history which cannot be ignored. Simulations have been done to illustrate the resolution ability of this configuration, and both range resolution and Doppler resolution are space-variant and time-variant. The tiny changes of the values of 2D resolution in this configuration during short synthetic aperture time and comparisons with ABFL-SAR can be used to design imaging algorithms in this special configuration. Research on the properties of MBFL-SAR configuration validates the efficiency of the application of MIMO radar networking to radar imaging, which also provides theoretical groundwork for the whole system parameter design and imaging algorithms of MBFL-SAR.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This work was supported by the National Natural Science Foundation of China under Grant nos. 61001211 and 61303035, the Fundamental Research Funds for the Central Universities (K5051202016), and the Science Foundation for Navigation (20110181004).
- E. Fishler, A. Haimovich, R. Blum, D. Chizhik, L. Cimini, and R. Valenzuela, “MIMO radar: an idea whose time has come,” in Proceedings of the IEEE Radar Conference, pp. 71–78, April 2004.
- X. Li, Z. Zhang, W. Mao, X. Wang, J. Lu, and W. Wang, “A study of frequency diversity MIMO radar beamforming,” in Proceedings of the IEEE 10th International Conference on Signal Processing (ICSP '10), pp. 352–356, October 2010.
- R. Sharma, “Analysis of MIMO radar ambiguity functions and implications on clear region,” in Proceedings of the IEEE International Radar Conference (RADAR '10), pp. 544–548, June 2010.
- J. Li and P. Stoica, MIMO Radar Signal Processing, John Wiley & Sons, Hoboken, NJ, USA, 2009.
- I. Cumming and F. H. Wong, Digital Processing of Synthetic Aperture Radar Data: Algorithms and Implementation, Artech House, Norwood, Mass, USA, 2004.
- J. Li, G. Liao, and H. Griffiths, “Bistatic MIMO radar space-time adaptive processing,” in Proceedings of the IEEE Radar Conference: In the Eye of the Storm (RadarCon '11), pp. 498–501, May 2011.
- X. H. Wu, A. A. Kishk, and A. W. Glisson, “MIMO-OFDM radar for direction estimation,” IET Radar, Sonar & Navigation, vol. 4, no. 1, pp. 28–36, 2010.
- J. Li and P. Stoica, “MIMO radar with colocated antennas,” IEEE Signal Processing Magazine, vol. 24, no. 5, pp. 106–114, 2007.
- A. Haimovich, R. Blum, and L. Cimini, “MIMO radar with widely separated antennas,” IEEE Signal Processing Magazine, vol. 25, no. 1, pp. 116–129, 2008.
- Z. Shi, H. Wang, P. Zhang, X. Tang, B. Wu, and C. Wang, “Study on the probability of successful handoff of missile trajectory from midcourse guidance to terminal guidance,” in Proceedings of the International Conference on Computational and Information Sciences (ICCIS '10), pp. 1005–1008, Chengdu, China, December 2010.
- Q. Zhao, X. Zhang, W. Song et al., “Rotational analysis of certain terminal-guidance projectile,” in Proceedings of the International Conference on Electronic Measurement and Instruments (ICEMI '07), pp. 415–418, July 2007.
- F. Comblet, A. Khenchaf, A. Baussard et al., “Bistatic synthetic aperture radar imaging: theory, simulations, and validations,” IEEE Transactions on Antennas and Propagation, vol. 54, no. 11, pp. 3529–3540, 2006.
- I. Walterscheid, T. Espeter, J. Klare, A. R. Brenner, and J. H. G. Ender, “Potential and limitations of forward-looking bistatic SAR,” in Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS '10), pp. 216–219, Honolulu, Hawaii, USA, July 2010.
- W. Wang, “Approach of adaptive synchronization for bistatic SAR real-time imaging,” IEEE Transactions on Geoscience and Remote Sensing, vol. 45, no. 9, pp. 2695–2700, 2007.
- W. Wang, C. Ding, and X. Liang, “Time and phase synchronisation via direct-path signal for bistatic synthetic aperture radar systems,” IET Radar, Sonar and Navigation, vol. 2, no. 1, pp. 1–11, 2008.
- M. Weib, “Synchronisation of bistatic radar systems,” in Proceedings of IEEE International Geoscience and Remote Sensing Symposium (IGARSS '04), pp. 1750–1753, September 2004.
- J. Wu, J. Yang, H. Yang, and Y. Huang, “Optimal geometry configuration of bistatic forward-looking SAR,” in Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '09), pp. 1117–1120, 2009.
- X. Qiu and D. Hu, “Some reflections on bistatic SAR of forward-looking configuration,” IEEE Geoscience and Remote Sensing Letters, vol. 5, no. 4, pp. 735–739, 2008.
- T. Espeter, I. Walterscheid, J. Klare et al., “Bistatic forward-looking SAR: results of a spaceborne-airborne experiment,” IEEE Geoscience and Remote Sensing Letters, vol. 8, no. 4, pp. 765–768, 2011.
- J. Balke, “Field test of bistatic forward-looking synthetic aperture radar,” in Proceedings of IEEE International Radar Conference, pp. 424–429, May 2005.
- Z. Li, D. Yao, and T. Long, “SPECAN algorithm for forward-looking bistatic SAR,” in Proceedings of 9th International Conference on Signal Processing ( ICSP '08), pp. 2517–2520, October 2008.
- R. Wang, O. Loffeld, H. Nies, and V. Peters, “Image formation algorithm for bistatic forward-looking SAR,” in Proceedings of the 30th IEEE International Geoscience and Remote Sensing Symposium (IGARSS '10), pp. 4091–4094, Honolulu, Hawaii, USA, July 2010.
- J. Wu, Z. Li, and Y. Huang, “Focusing bistatic forward-looking SAR with stationary transmitter based on keystone transform and nonlinear chirp scaling,” IEEE Geoscience and Remote Sensing Letters, vol. 11, no. 1, pp. 148–152, 2014.
- G. P. Cardillo, “On the use of gradient to determine bistatic SAR resolution,” in Proceedings of Antennas and Propagation Society International Symposium, vol. 2, pp. 1032–1035, May 1990.
Copyright © 2014 Yachao Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.