High-Performance Anti-Retransmission Deception Jamming Utilizing Range Direction Multiple Input and Multiple Output (MIMO) Synthetic Aperture Radar (SAR)

Retransmission deception jamming seriously degrades the Synthetic Aperture Radar (SAR) detection efficiency and can mislead SAR image interpretation by forming false targets. In order to suppress retransmission deception jamming, this paper proposes a novel multiple input and multiple output (MIMO) SAR structure range direction MIMO SAR, whose multiple channel antennas are vertical to the azimuth. First, based on the multiple channels of range direction MIMO SAR, the orthogonal frequency division multiplexing (OFDM) linear frequency modulation (LFM) signal was adopted as the transmission signal of each channel, which is defined as a sub-band signal. This sub-band signal corresponds to the transmission channel. Then, all of the sub-band signals are modulated with random initial phases and concurrently transmitted. The signal form is more complex and difficult to intercept. Next, the echoes of the sub-band signal are utilized to synthesize a wide band signal after preprocessing. The proposed method will increase the signal to interference ratio and peak amplitude ratio of the signal to resist retransmission deception jamming. Finally, well-focused SAR imagery is obtained using a conventional imaging method where the retransmission deception jamming strength is degraded and defocused. Simulations demonstrated the effectiveness of the proposed method.


Introduction
Synthetic aperture radar (SAR) has been widely used in civil exploration and military surveillance due to its all-weather, long-time and large-range detection capabilities [1]. However, rapidly developed electronic countermeasures (ECMs) can degrade the SAR imaging quality and system efficiency significantly [2][3][4]. The ECMs for SAR can be classified into barrage jamming and deception jamming, according to the signal characteristics and purpose of the ECM. Barrage jamming uses radio frequency noise with enough power to cover and corrupt the target signal. Deception jamming modulates fake target information on an intercepted signal retransmitted to the SAR via the direct-path or multipath routes [5]. Due to the high power consumption and easy identification and suppression, barrage jamming is less utilized compared to deception jamming. The false target formation by deception jamming will mislead the SAR image interpretation, which is difficult to recognize and mitigate, making deception jamming more flexible and efficient than barrage jamming.
As a result, anti-jamming and jamming suppression methods have become an appealing research topic in recent years [6][7][8]. Until now, most of the research and literature has been focused on the problem of narrow band interference (NBI) suppression, which is a classical unintentional interference caused by communication sources and broadcast and TV signals [9][10][11][12][13]. The relevant research on deception jamming suppression, wide band interference, is insufficient and not comprehensive. Miller proposed a parametric method based on building a mathematical model of the interference. Nevertheless, the effectiveness of jamming suppression is determined by the model's degree of matching with the jamming. In [14], an adaptive two-dimensional filter, relying on estimating the statistical characteristics of the SAR image, was designed to suppress the NBI. Meanwhile, as for a classical non-parametric method, it is arduous to design an appropriate filter for the echo signal, which is not only blurred by NBI but also interfered with by wide band jamming. Methods such as Independent Component Analysis (ICA) and Independent Subspace Analysis (ISA) [12,13], were adopted for NBI suppression. These can extract the statistically independent jamming signal relying on the significant variance between jamming and target echoes in the data domain. ISA improves signal projection from the one-dimensional time domain to the two-dimensional time-frequency domain and reduces the constraints of ICA. However, when the jamming signal is not independent of the target echoes, the jamming suppression effect will be degraded. The aforementioned NBI suppression methods are not suitable for wideband interference, such as deception jamming. Tao et al. utilized the short-time Fourier transform (STFT) to convert the wideband interference (WBI) into a series of instantaneous NBI spectrum mitigation problems [15]. The premise of Tao's method is that WBI has distinct energy variations in the time-localized spectra that can be captured by the STFT. However, deception jamming can intercept the SAR signal and retransmit it to the SAR without energy variations in localized time.
In essence, the main approach to deception jamming suppression is to use a complex signal during transmission and increase the signal to interference ratio (SIR) at the receiving terminal. The MIMO SAR [16,17] realizes this idea based on its multiple channel receiving and transmitting ability. In [18], Gebert proposed the structure of an azimuth MIMO SAR in order to improve the resolution. In [17], the OFDM signal, a typical independent waveform for MIMO radar, is used in MIMO SAR to reduce the range ambiguity. MIMO SAR can not only achieve high-resolution and wide-swath imaging but also has great potential in anti-jamming and jamming suppression. In [19], Rosenberg made use of the multiple channel signal receiving ability, a vital advantage of MIMO SAR, to cancel hot-clutter jamming, one type of multipath jamming. However, direct-path wide band jamming, such as retransmission deception jamming, was not mentioned.
To overcome the limitations of the aforementioned methods, we propose herein a novel MIMO structure, called range direction MIMO SAR, to suppress retransmission deception jamming. The multi-channel antenna of the MIMO SAR is arranged along the range direction. First, at the transmitting terminal, the orthogonal frequency division multiplexing (OFDM) LFM signal is adopted as the sub-band signal [20], which has a one-to-one correspondence with the transmitting channel of the antenna array. Then, random phases were modulated in each sub-band signal. After that, the modulated sub-band signal will be concurrently transmitted by the corresponding antenna. The concurrently transmitted sub-band signals will arrive at the jammer at the same time, which makes them hard to intercept and be recognized by the jammer. The jammer can only retransmit the signals with frequencies located in the jammer's working frequency band. Hence, the jammer's retransmission time delay will increase more than a Pulse Repetition Time (PRT), which means the phase of the retransmitted deceptive jamming signal will not match the target echo in the azimuth. Second, at the receiving terminal, the sub-band signals will be pre-processed, including the phase compensation, time shift and frequency shift. Then, the sub-band signals can be utilized to synthesize a wide band signal in the range direction. Meanwhile, part of the frequency band of the synthesized wide band signal will be blurred by the retransmission deceiving jamming. Finally, the traditional methods could be applied to get the SAR image of the synthesized wide band signal. After using the proposed methods, the jammer signal could not get the signal processing gain and will defocus in the SAR image because of the mismatch with the target echo. Meanwhile, the sub-band signal synthesis will increase the target signal processing gain and the signal to interference ratio will increase. This paper is organized as follows: Section 2 presents the principle of range direction MIMO SAR. Section 3 proposes the pre-processing method of the concurrently transmitted sub-band signals of range direction MIMO SAR. Section 4 analyzes the efficiency of the proposed retransmission deception jamming suppression. Section 5 shows the results of our simulations. Section 6 concludes the paper.

Geometry Model
The antenna is staggered in a direction that is vertical with the azimuth direction, which is called the range direction. There are M transmitting and receiving sub antennas comprising the whole linear MIMO antenna array. The OFDM-LFM [18,20] is used in the Range MIMO SAR, and each sub-band has a one-to-one correspondence with the sub-antenna element, which means that the number of OFDM sub-band signals is determined by the sub-antenna element count.
The geometry model is illustrated in Figure 1. We suppose that there are M sub-antenna elements, that the distance between them is equal, and that each antenna element can receive the transmitted signals from the others. As shown in Figure 1, we assume that the MIMO antenna linear array is along the y-axis. Its length is L, and the interval between each element is d. The MIMO SAR travels along with the x-axis with a speed of v. The height of the MIMO SAR is H. The projection of the SAR platform in the xoy plane is the coordinate origin. The range direction is along the y-axis. The middle point of the antenna E mid is located at (0, y m , 0) when the azimuth time is zero. Likewise, when the azimuth time is t a , the k-th antenna's coordinate E k can be indicated as (vt a , y k , H). Hence, an arbitrary element E k 's y-axis coordinate is as follows: After using the proposed methods, the jammer signal could not get the signal processing gain and will defocus in the SAR image because of the mismatch with the target echo. Meanwhile, the sub-band signal synthesis will increase the target signal processing gain and the signal to interference ratio will increase. This paper is organized as follows: Section 2 presents the principle of range direction MIMO SAR. Section 3 proposes the pre-processing method of the concurrently transmitted sub-band signals of range direction MIMO SAR. Section 4 analyzes the efficiency of the proposed retransmission deception jamming suppression. Section 5 shows the results of our simulations. Section 6 concludes the paper.

Geometry Model
The antenna is staggered in a direction that is vertical with the azimuth direction, which is called the range direction. There are M transmitting and receiving sub antennas comprising the whole linear MIMO antenna array. The OFDM-LFM [18,20] is used in the Range MIMO SAR, and each sub-band has a one-to-one correspondence with the sub-antenna element, which means that the number of OFDM sub-band signals is determined by the sub-antenna element count.
The geometry model is illustrated in Figure 1. We suppose that there are M sub-antenna elements, that the distance between them is equal, and that each antenna element can receive the transmitted signals from the others. As shown in Figure 1, we assume that the MIMO antenna linear array is along the y-axis. Its length is L, and the interval between each element is d. The MIMO SAR travels along with the x-axis with a speed of v. The height of the MIMO SAR is H. The projection of the SAR platform in the xoy plane is the coordinate origin. The range direction is along the y-axis. The middle point of the antenna Emid is located at (0, ym, 0) when the azimuth time is zero. Likewise, when the azimuth time is ta, the k-th antenna's coordinate Ek can be indicated as (vta, yk, H). Hence, an arbitrary element Ek's y-axis coordinate is as follows: There is a scatter target and a jammer in the observation scene. Meanwhile, if the point scatter target P is located at (x0, y0, z0) and the jammer is located at (xJ, yJ, zJ), the double trip instantaneous slant range of target Rkl, which presents the slant range from the k-th transmitting antenna element to the l-th receiving antenna element, can be indicated as: There is a scatter target and a jammer in the observation scene. Meanwhile, if the point scatter target P is located at (x 0 , y 0 , z 0 ) and the jammer is located at (x J , y J , z J ), the double trip instantaneous slant range of target R kl , which presents the slant range from the k-th transmitting antenna element to the l-th receiving antenna element, can be indicated as: where, t a is the azimuth time, and R kP represents the single trip slant range from the k-th transmitting antenna to the target P. Similarly, R Pl is the single trip slant range from target P to the l-th receiving antenna. Likewise, the double trip instantaneous slant range of the jammer is:

Signal Model
Usually, the sub-band signal can be formed by dividing the wide band signal in the frequency domain. This is illustrated in Figure 2. The bandwidth of the original wide band signal is B w , and it is divided by the bandwidth B. After a time shift, the sub-band signal will be concurrently transmitted in the corresponding transmission channel.
where, ta is the azimuth time, and RkP represents the single trip slant range from the k-th transmitting antenna to the target P. Similarly, RPl is the single trip slant range from target P to the l-th receiving antenna. Likewise, the double trip instantaneous slant range of the jammer is:

Signal Model
Usually, the sub-band signal can be formed by dividing the wide band signal in the frequency domain. This is illustrated in Figure 2. The bandwidth of the original wide band signal is Bw, and it is divided by the bandwidth B. After a time shift, the sub-band signal will be concurrently transmitted in the corresponding transmission channel. Supposed that the pulse width is Tw, the bandwidth is Bw, the center frequency is fc, and the center time Tc is 0. Therefore, the modulated frequency ratio Kr is equal to Bw/Tw. At is the signal amplitude. Then, the wideband signal can be indicated as: If the wideband signal Sw is divided into M sub-bands, the bandwidth of the sub-band signal B is Bw/M and the pulse width of the sub-band signal T is Tw/M. Therefore, the center frequency of the k-th sub-band signal is as follows: Then, the center time of the k-th sub-band signal is indicated as follows: So that the k-th sub-band transmitting signal with random phase modulation can be indicated as: Without consideration of the time delay, the concurrent transmission of all sub-band signals is as follows: Supposed that the pulse width is T w , the bandwidth is B w , the center frequency is f c , and the center time T c is 0. Therefore, the modulated frequency ratio K r is equal to B w /T w . A t is the signal amplitude. Then, the wideband signal can be indicated as: If the wideband signal S w is divided into M sub-bands, the bandwidth of the sub-band signal B is B w /M and the pulse width of the sub-band signal T is T w /M. Therefore, the center frequency of the k-th sub-band signal is as follows: Then, the center time of the k-th sub-band signal is indicated as follows: So that the k-th sub-band transmitting signal with random phase modulation can be indicated as: Without consideration of the time delay, the concurrent transmission of all sub-band signals is as follows: We suppose that there are n c targets in the scene. Then, the target echoes, which are transmitted by the k-th antenna and received by the l-th antenna, can be indicated as: where t r and t a are the time of range direction and azimuth, respectively. T is the range direction time width of the sub-band signal. T a is the azimuth synthesis time. f k is the k-th working center frequency of the k th transmitting antenna. δ n c kl is the n c -th target radar cross-section (RCS). A r is the receiving signal amplitude which is attenuated by the radar equation. Theoretically, the antenna elements corresponding with the sub-band will receive other antenna-transmitted sub-band signals because the frequency band of the antenna elements is equal to the wide band signal that is synthesized by the sub-band signal. Hence, the whole echoes of target signal are as follows: Usually, we assume that the jammer can only intercept the signal located in its working frequency channel and accurately estimate all the parameters. According to the theoretical model of retransmission deception jamming [3,5], the major procedures are time delay and false information modulation. If the jammer intercepts the I-th sub-band signal, the deception jamming signal model is as follows: where the ∆R JP n (t a ) is the instantaneous single-trip slant range from point J to the n-th false target P, and τ s is the system inherent time delay of the jammer. A J is the jamming amplitude. The operator * represents the convolution, which means the time delay of the intercepted signal. The false information is included in the additional phase modulation. Generally, the system time delay is beyond the PRT. Therefore, τ s is utilized to modify the system time delay as follows: so that the echoes, which include the target signal and the retransmitted deception jamming signal, can be represented as:

Range MIMO Signal Processing
Due to the one-to-one correspondence between the sub-band signal and the antenna element, the range MIMO SAR antenna with M channels is equivalent to having the capability of using M sub-band signals [20]. Moreover, the geometry model of the range MIMO SAR antenna will determine the position of the equivalent phase center [21], which is established to solve the MIMO signal processing.
The M antenna elements ranging along the y-axis will produce M × M equivalent phase centers, where (2 × M − 1)'s positions are different from each other. Hence, as Figure 3 shows, the 3-channel MIMO SAR has nine phase centers in total, and only five positions are independent. The equivalent phase center of the k-th transmitting antenna and l-th receiving antenna is located at y kl , and indicated by: y kl = (y k + y l )/2 (14) phase center position of three transmitting elements is overlapped and situated in the 2nd antenna element position. In the overlapped phase center, the wide band signal will be synthesized without any additional conditions because the OFDM signal is utilized as the sub-band signal. The bandwidth of the synthesized signal in the 2nd antenna elements is triple that of the sub-band signal when the bandwidth of the sub-band signal is equal. The adjacent phase center will synthesize twice the original band width of the sub-band signal. Simultaneously, the range direction resolution of each phase center is diverse, which is positively related to the bandwidth of the synthesized wide band signal. Meanwhile, for the range direction MIMO SAR with M transmitting and receiving antenna channels, the bandwidth range of the synthesized signal is from 1 to M times the sub-band signal. To get the highest resolution and largest bandwidth, the middle equivalent phase center is utilized as an effective jamming suppression terminal. Range direction MIMO SAR has advantages compared to azimuth MIMO [20,22]. The merit of range MIMO is that the azimuth processing is equivalent to a traditional single channel SAR without consideration of non-uniform sampling and any Doppler ambiguity in azimuth direction [23,24].
The movement of every equivalent phase center is translational after an azimuth interval. Hence, the azimuth sample will not be staggered. According to the operation principle of equivalent phase center, the sum of the single-trip slant range RkP and RPl, which is the actual range Rkl, is approximate to the double-trip slant range kl R from the equivalent phase center Ekl with coordinate (vta, ykl, H) to the target P: The slant range Rkl from the k-th transmitting antenna element Ek to the l-th receiving antenna element El is approximately equal to the double-trip slant range kl R from the equivalent phase center Ekl to the target P. For any antenna elements E with coordinate (x, y, H), the Taylor series expansion of the single trip slant range R(ta, y) around arbitrary point yT is: Usually, the slant range R is expanded around Emid, i.e., the center point of the linear antenna array. Hence, Equation (2) can be written as follows: From the equation, the location of y kl is the same with y lk . As the figure shows, the equivalent phase center position of three transmitting elements is overlapped and situated in the 2nd antenna element position. In the overlapped phase center, the wide band signal will be synthesized without any additional conditions because the OFDM signal is utilized as the sub-band signal.
The bandwidth of the synthesized signal in the 2nd antenna elements is triple that of the sub-band signal when the bandwidth of the sub-band signal is equal. The adjacent phase center will synthesize twice the original band width of the sub-band signal. Simultaneously, the range direction resolution of each phase center is diverse, which is positively related to the bandwidth of the synthesized wide band signal.
Meanwhile, for the range direction MIMO SAR with M transmitting and receiving antenna channels, the bandwidth range of the synthesized signal is from 1 to M times the sub-band signal. To get the highest resolution and largest bandwidth, the middle equivalent phase center is utilized as an effective jamming suppression terminal. Range direction MIMO SAR has advantages compared to azimuth MIMO [20,22]. The merit of range MIMO is that the azimuth processing is equivalent to a traditional single channel SAR without consideration of non-uniform sampling and any Doppler ambiguity in azimuth direction [23,24].
The movement of every equivalent phase center is translational after an azimuth interval. Hence, the azimuth sample will not be staggered. According to the operation principle of equivalent phase center, the sum of the single-trip slant range R kP and R Pl , which is the actual range R kl , is approximate to the double-trip slant range R kl from the equivalent phase center E kl with coordinate (vt a , y kl , H) to the target P: The slant range R kl from the k-th transmitting antenna element E k to the l-th receiving antenna element E l is approximately equal to the double-trip slant range R kl from the equivalent phase center E kl to the target P. For any antenna elements E with coordinate (x, y, H), the Taylor series expansion of the single trip slant range R(t a , y) around arbitrary point y T is: Usually, the slant range R is expanded around E mid , i.e., the center point of the linear antenna array. Hence, Equation (2) can be written as follows: The position of the equivalent phase center is the center between the transmitting and receiving antenna elements. The Taylor series expansion of the equivalent phase center can be indicated as follows: Meanwhile, the slant range difference is as follows:

Phase Error Correction
The phase difference, which is caused by the range difference and has some negative effects on the SAR imaging, needs to be corrected in any imaging process. As Equation (16) shows, B can be indicated as: where R c (y m ) presents the instantaneous slant range between the center point of the antenna array and the target point. R c⊥ (y m ) is the projection of R c (y m ) in the xoz plane. Hence, the phase difference is as follows: where d kl = y k − y l . It presents the absolute interval between element E k and element E l . On the one hand, the phase error is related to the synthesis aperture time for the same target. When the azimuth changed, the phase error for the same target P with coordinate (x 0 , y 0 , z 0 ) can be indicated as: We suppose that the velocity of the SAR platform is 300 m/s, d ij is 5 m, λ is 0.056 m, and the center point P is located at (0, 16010, 0). When t a1 is zero and t a2 is from −2.5/s to 2.5/s, the phase error change is as shown in Figure 4.
The position of the equivalent phase center is the center between the transmitting and receiving antenna elements. The Taylor series expansion of the equivalent phase center can be indicated as follows: Meanwhile, the slant range difference is as follows:

Phase Error Correction
The phase difference, which is caused by the range difference and has some negative effects on the SAR imaging, needs to be corrected in any imaging process. As Equation (16) shows, B can be indicated as: where Rc(ym) presents the instantaneous slant range between the center point of the antenna array and the target point.
is the projection of Rc(ym) in the xoz plane. Hence, the phase difference is as follows: where dkl = yk − yl. It presents the absolute interval between element Ek and element El.
On the one hand, the phase error is related to the synthesis aperture time for the same target. When the azimuth changed, the phase error for the same target P with coordinate (x0, y0, z0) can be indicated as: We suppose that the velocity of the SAR platform is 300 m/s, dij is 5 m,  is 0.056 m, and the center point P is located at (0, 16010, 0). When ta1 is zero and ta2 is from −2.5/s to 2.5/s, the phase error change is as shown in Figure 4.  The phase error will increase when t a2 is far from t a1 . Thus, the phase error needs to be corrected in each azimuth sampling time. On the other hand, the phase error is related to the positions of the different targets in the main-lobe beam coverage. Usually, we use the slant range of the center point in the beam coverage to replace the slant ranges of the other targets in the scene. Meanwhile, the phase error will be imported. If there are two targets P 1 and P 2 with coordinates (x 1 , y 1 , z 1 ) and (x 2 , y 2 , z 2 ), the phase error between P 1 and P 2 is: When the wave length and antenna interval are fixed, the phase error is only affected by the target position. We suppose that the position of target P 2 is moving along the x-axis and y-axis. Then, P 1 is located at (0, 0). Accordingly, for different intervals of the antenna channel d, the phase error is as shown by Figure 5. The phase error will increase when ta2 is far from ta1. Thus, the phase error needs to be corrected in each azimuth sampling time.
On the other hand, the phase error is related to the positions of the different targets in the main-lobe beam coverage. Usually, we use the slant range of the center point in the beam coverage to replace the slant ranges of the other targets in the scene. Meanwhile, the phase error will be imported. If there are two targets P1 and P2 with coordinates (x1, y1, z1) and (x2, y2, z2), the phase error between P1 and P2 is: When the wave length and antenna interval are fixed, the phase error is only affected by the target position. We suppose that the position of target P2 is moving along the x-axis and y-axis. Then, P1 is located at (0, 0). Accordingly, for different intervals of the antenna channel d, the phase error is as shown by Figure 5. If target P2 is at any point in the scene and d is 1.5 m, the phase error of the beam coverage will be obtained as in Figure 6.  If target P 2 is at any point in the scene and d is 1.5 m, the phase error of the beam coverage will be obtained as in Figure 6. The phase error will increase when ta2 is far from ta1. Thus, the phase error needs to be corrected in each azimuth sampling time.
On the other hand, the phase error is related to the positions of the different targets in the main-lobe beam coverage. Usually, we use the slant range of the center point in the beam coverage to replace the slant ranges of the other targets in the scene. Meanwhile, the phase error will be imported. If there are two targets P1 and P2 with coordinates (x1, y1, z1) and (x2, y2, z2), the phase error between P1 and P2 is: When the wave length and antenna interval are fixed, the phase error is only affected by the target position. We suppose that the position of target P2 is moving along the x-axis and y-axis. Then, P1 is located at (0, 0). Accordingly, for different intervals of the antenna channel d, the phase error is as shown by Figure 5. If target P2 is at any point in the scene and d is 1.5 m, the phase error of the beam coverage will be obtained as in Figure 6.   To obtain a high resolution image, the phase difference needs to be modified by using the center point of the beam coverage. The simulation proved that when the antenna interval d is smaller, the phase error is acceptable for image processing.

Signal Processing Flowchart
The signal processing includes two primary parts: preprocessing and imaging processing. After preprocessing, the traditional imaging algorithm can be utilized. The procedure of preprocessing is as follows: the first step in multiple sub-band signal processing is using the corresponding center frequency to demodulate the sub-band signals. Then, the sub-band signals are filtered to reduce the out-band noise and wide band jamming. The second step of the filtered sub-band signal is modifying the random phase of the sub-band signal. After this step, the target signal will be accurately matched. However, the system time delay of the deception jamming permits the jamming signal mismatch. The third step of the echo is the phase difference compensation for the multiple band signals. The fourth step of the compensated signal is the synthesis of the wide band signal using the sub-band signals. We supposed that there are M receiving and transmitting channels in the range MIMO SAR. Hence, the receiving signal can be indicated as: According to the position of the signal equivalent phase center, the sub-band signal is settled in the matrix. Every column presents a phase center. Meanwhile, the Q-th equivalent phase center signal is: The band width is determined by the position of the equivalent phase center. B Q , the band width of S E Q , is as follows: Thus, the widest band width B W = MB. The number and bandwidth of the sub-band signal are two decisive components of B W . To analyze the preprocessing easily, we ignore the influence of the target number and RCS. Let τ a = R n c ij (t a )/c, the equation can be indicated as: The S OB , echo of the transmitted wide band signal S w , is as follows: Simultaneously, it can be represented by the sub-band signal S E Q .
where ∆T k is equal to T k − T c . The sub-band signal will be separated according to different carrying frequencies. Therefore, the sub-band signal is demodulated by the center frequency and compensated random phase: Then, every sub-band signal in the same equivalent phase center is time shifted. Let t r = t r − ∆T k , so the process is as follows: Comparing this equation with Equation (30), the time shifted sub-band signal lacks a linear phase and a constant phase. The frequency shift and phase compensation will recover the wide band by using the time-shifting sub-band signals. In conclusion, the signal preprocessing flow chart of range MIMO SAR which adopts the OFDM signal as the sub-band signal and random phase modulation, is as seen in Figure 7.
Then, every sub-band signal in the same equivalent phase center is time shifted. Let r r k t t T   , so the process is as follows: Comparing this equation with Equation (30), the time shifted sub-band signal lacks a linear phase and a constant phase. The frequency shift and phase compensation will recover the wide band by using the time-shifting sub-band signals. In conclusion, the signal preprocessing flow chart of range MIMO SAR which adopts the OFDM signal as the sub-band signal and random phase modulation, is as seen in Figure 7. After preprocessing, the sub-band signals have to be synthesized into a wide band signal that can be utilized to obtain SAR images by a traditional imaging algorithm. Moreover, the preprocessing will suppress the jamming signals.

Range Direction MIMO SAR Anti-Deception Jamming Efficiency Analysis
Considering the retransmision deception jamming in the echoes as Figure 8, the anti-jamming efficiency of the range MIMO SAR is analyzed after signal pre-processing and imaging. Because the range direction MIMO SAR using the concurrently transmitted sub-band signal is modulated by a random phase, the jammer processing time will increase and even cannot intercept and recognize concurrently transmitted signals.  After preprocessing, the sub-band signals have to be synthesized into a wide band signal that can be utilized to obtain SAR images by a traditional imaging algorithm. Moreover, the preprocessing will suppress the jamming signals.

Range Direction MIMO SAR Anti-Deception Jamming Efficiency Analysis
Considering the retransmision deception jamming in the echoes as Figure 8, the anti-jamming efficiency of the range MIMO SAR is analyzed after signal pre-processing and imaging. After preprocessing, the sub-band signals have to be synthesized into a wide band signal that can be utilized to obtain SAR images by a traditional imaging algorithm. Moreover, the preprocessing will suppress the jamming signals.

Range Direction MIMO SAR Anti-Deception Jamming Efficiency Analysis
Considering the retransmision deception jamming in the echoes as Figure 8, the anti-jamming efficiency of the range MIMO SAR is analyzed after signal pre-processing and imaging. Because the range direction MIMO SAR using the concurrently transmitted sub-band signal is modulated by a random phase, the jammer processing time will increase and even cannot intercept and recognize concurrently transmitted signals.  Because the range direction MIMO SAR using the concurrently transmitted sub-band signal is modulated by a random phase, the jammer processing time will increase and even cannot intercept and recognize concurrently transmitted signals. The inherent system time delay τ s can be indicated as follows: where represents selecting the lesser round number, and mod means selecting the remainder. The integer n indicates that the jamming signal will fall behind n PRTs in the SAR receiving terminal. As in Equation (11), we suppose that the jammer intercepts the I-th sub-band signal: After carrying out frequency demodulation and random phase demodulation, the deception jammer signal is as follows: where A m means the jammer signal gain [4,25]. τ is the time delay and is given as follows: After range compression, the jammer signal is as follows: where ∆ϕ(t a ) means the phase difference caused by the jammer system time delay and is equal to the minus phase: Next, we suppose that the slant range of the deception target generated by the retransmit jamming is approximate with the real target: Hence, after the Range Cell Migration Correction, the false target P will be corrected into the shortest slant range R P Il . The jamming signal is as follows: If the azimuth compression response function is h a (t a ), the azimuth compression process is as follows: Based on the integral mean value theorem, the azimuth gain of the jamming signal is in the range given by the expression: where J a is the azimuth signal of the entire jamming signal J. Equation (41) shows that the jamming signal defocuses in the azimuth direction due to the random phase mismatch. The greater the phase difference randomness, the lower the jamming signal pulse compression processing gain. The shorter the integrating range, the lower the jamming signal compression processing gain. We suppose that the jamming signal azimuth processing gain is K (T a ), which is less than the real target azimuth processing gain T a . As Equation (40) shows, the jamming signal will achieve sub-band range compression gain. Hence, according to whole signal processing flowchart, the total gain of the jamming signal is proportional to the following equation: As for the single sub-band target signal, after the range and azimuth compression processing, the total gain is as follows: Then, the sub-band signal will be synthesized to a wide band signal, which will increase the range processing gain. The total gain will be as follows: Hence, the peak aptitude ratio of the target to jamming without sub-band synthesizing is indicated as: The peak aptitude ratio of target to jamming after sub-band synthesizing is indicated as:

Simulation
Retransmission deception jamming is simulated in this part. We supposed that the range direction MIMO SAR has 5-channel antenna elements. The geometry model is as shown in Figure 1. The target is located in a 3 × 3 matrix that has a middle point P in the center of the scene with coordinates (0, 16010, 0). Then, the interval in the range direction of the point matrix is 10 m, and the interval in the azimuth is 30 m. The jammer J is located in the adjacent point in the same row with coordinates (0, 16000, 0) in the point matrix. The other simulation conditions are given the Table 1. First, the random phase modulation is simulated. We suppose that the modulation random phase obeys Gaussian noise N(0, 2π) and that the jammer system time delay is 2 PRTs. When the azimuth synthesis time and frequency modulation rate are fixed such that the azimuth time is 1 s and the chirp rate in the azimuth direction is 500 Hz, the result of the jammer signal and target signal azimuth compression processing is as follows: Figure 9 shows that the deception jamming signal will defocus after the azimuth compression because of the phase residue in random phase compensation. The phase residue is associated with the time delay of retransmitting jamming signal. Next, we analyze the effect of the azimuth time to the peak amplitude ratio (PAR) of the target signal to the jamming signal without sub-band synthesis. The Monte-Carlo experiment is simulated when the azimuth chirp rate is 50 kHz/s and the sampling rate is 60 kHz/s. The azimuth time is from 0.01 s to 1 s, and the random phase still obeys N(0, 2π). The peak amplitude ratio simulation is as seen in Figure 10. the azimuth chirp rate is 50 kHz/s and the sampling rate is 60 kHz/s. The azimuth time is from 0.01 s to 1 s, and the random phase still obeys (0, 2 ) N  . The peak amplitude ratio simulation is as seen in Figure 10.  It is noticeable that the PAR will increase when the azimuth time is increasing, which is consistent with the theoretical analysis that rA1 is proportional to the azimuth time. The effect of the random phase modulation relies on the azimuth time. The longer the azimuth time, the better the jamming suppression result. Second, we suppose that the inherent jammer delay time is 2 PRT and that the jammer intercepts the third band signal and forms the retransmitted jamming signal. When the jammer J forms two false point targets PF1 and PF2 located at (15, 16010, 0) and (−15, 16010, 0), respectively, the imaging result is as seen in Figure 11.
When the interference to signal ratio is from 0 dB to 40 dB and the step is 20 dB, which means that the interference power is increasing, the sub-band synthesis processing of range MIMO SAR the azimuth chirp rate is 50 kHz/s and the sampling rate is 60 kHz/s. The azimuth time is from 0.01 s to 1 s, and the random phase still obeys (0, 2 ) N  . The peak amplitude ratio simulation is as seen in Figure 10.  It is noticeable that the PAR will increase when the azimuth time is increasing, which is consistent with the theoretical analysis that rA1 is proportional to the azimuth time. The effect of the random phase modulation relies on the azimuth time. The longer the azimuth time, the better the jamming suppression result. Second, we suppose that the inherent jammer delay time is 2 PRT and that the jammer intercepts the third band signal and forms the retransmitted jamming signal. When the jammer J forms two false point targets PF1 and PF2 located at (15, 16010, 0) and (−15, 16010, 0), respectively, the imaging result is as seen in Figure 11.
When the interference to signal ratio is from 0 dB to 40 dB and the step is 20 dB, which means that the interference power is increasing, the sub-band synthesis processing of range MIMO SAR It is noticeable that the PAR will increase when the azimuth time is increasing, which is consistent with the theoretical analysis that r A1 is proportional to the azimuth time. The effect of the random phase modulation relies on the azimuth time. The longer the azimuth time, the better the jamming suppression result. Second, we suppose that the inherent jammer delay time is 2 PRT and that the jammer intercepts the third band signal and forms the retransmitted jamming signal. When the jammer J forms two false point targets P F1 and P F2 located at (15,16010,0) and (−15, 16010, 0), respectively, the imaging result is as seen in Figure 11. As Figures 12-14 show, the deception jamming signal is obviously defocusing in the azimuth in both sub-band synthesis and single band processing so that there is a jamming stripe forming along the azimuth. Owing to the random phase modulation, the deception jamming only gets the range direction compression gain. Meanwhile, the jamming signal is located in the right range bin with the range coordinate 16010.
As the jamming power increases from 0 dB to 40 dB, the amplitude of the jamming stripe also increases. With the same jamming power, the jamming suppression effect of sub-band synthesis processing is better than that of the single sub-band processing. As Figure 14b shows, the amplitude of jamming is greater than that of the real target in single sub-band processing, which has negative effects on the image interpretation. Nevertheless, as Figure 14a shows, the amplitude of the jamming signal is still less than that of the real target in sub-band synthesis processing, which almost has no negative effects on the image interpretation. The quantitative results for the SIR and PAR of the signal to interference is indicated in Table 2.  When the interference to signal ratio is from 0 dB to 40 dB and the step is 20 dB, which means that the interference power is increasing, the sub-band synthesis processing of range MIMO SAR and single band processing are simulated, respectively.
As Figures 12-14 show, the deception jamming signal is obviously defocusing in the azimuth in both sub-band synthesis and single band processing so that there is a jamming stripe forming along the azimuth. Owing to the random phase modulation, the deception jamming only gets the range direction compression gain. Meanwhile, the jamming signal is located in the right range bin with the range coordinate 16010. As Figures 12-14 show, the deception jamming signal is obviously defocusing in the azimuth in both sub-band synthesis and single band processing so that there is a jamming stripe forming along the azimuth. Owing to the random phase modulation, the deception jamming only gets the range direction compression gain. Meanwhile, the jamming signal is located in the right range bin with the range coordinate 16010.
As the jamming power increases from 0 dB to 40 dB, the amplitude of the jamming stripe also increases. With the same jamming power, the jamming suppression effect of sub-band synthesis processing is better than that of the single sub-band processing. As Figure 14b shows, the amplitude of jamming is greater than that of the real target in single sub-band processing, which has negative effects on the image interpretation. Nevertheless, as Figure 14a shows, the amplitude of the jamming signal is still less than that of the real target in sub-band synthesis processing, which almost has no negative effects on the image interpretation. The quantitative results for the SIR and PAR of the signal to interference is indicated in Table 2.    As Table 2 shows, after image processing, the SIR and PAR, including the sub-band synthesis process and the single sub-band process, will decrease when the ISR in the echoes is increasing. However, the SIR and PAR of the sub-band synthesis process are higher than those of the single band process, which has also been illustrated in Figures 12-14. Even when the deception jamming signal power is 40 dB, the PAR of the sub-band signal synthesis is 2.9. Nevertheless, the PAR of the single sub-band process is 0.59, which means the amplitude of the jamming signal is larger than what the target shows in Figure 14b. After sub-band synthesis is used, the SIR and PAR will have  As Table 2 shows, after image processing, the SIR and PAR, including the sub-band synthesis process and the single sub-band process, will decrease when the ISR in the echoes is increasing. However, the SIR and PAR of the sub-band synthesis process are higher than those of the single band process, which has also been illustrated in Figures 12-14. Even when the deception jamming signal power is 40 dB, the PAR of the sub-band signal synthesis is 2.9. Nevertheless, the PAR of the single sub-band process is 0.59, which means the amplitude of the jamming signal is larger than what the target shows in Figure 14b. After sub-band synthesis is used, the SIR and PAR will have As the jamming power increases from 0 dB to 40 dB, the amplitude of the jamming stripe also increases. With the same jamming power, the jamming suppression effect of sub-band synthesis processing is better than that of the single sub-band processing. As Figure 14b shows, the amplitude of jamming is greater than that of the real target in single sub-band processing, which has negative effects on the image interpretation. Nevertheless, as Figure 14a shows, the amplitude of the jamming signal is still less than that of the real target in sub-band synthesis processing, which almost has no negative effects on the image interpretation. The quantitative results for the SIR and PAR of the signal to interference is indicated in Table 2.
As Table 2 shows, after image processing, the SIR and PAR, including the sub-band synthesis process and the single sub-band process, will decrease when the ISR in the echoes is increasing. However, the SIR and PAR of the sub-band synthesis process are higher than those of the single band process, which has also been illustrated in Figures 12-14. Even when the deception jamming signal power is 40 dB, the PAR of the sub-band signal synthesis is 2.9. Nevertheless, the PAR of the single sub-band process is 0.59, which means the amplitude of the jamming signal is larger than what the target shows in Figure 14b. After sub-band synthesis is used, the SIR and PAR will have increases in the sub-band number of times in the theoretical analysis. The simulation result is as follows: As Figure 15 shows, the PAR promotion r A2 /r A1 and SIR promotion are both approximately 5, which is the sub-band number in the simulation. Hence, by comparing the results of sub-band synthesis processing and single sub-band processing, one obvious feature is that the range direction MIMO SAR can effectively suppress deception jamming. increases in the sub-band number of times in the theoretical analysis. The simulation result is as follows: As Figure 15 shows, the PAR promotion rA2/rA1 and SIR promotion are both approximately 5, which is the sub-band number in the simulation. Hence, by comparing the results of sub-band synthesis processing and single sub-band processing, one obvious feature is that the range direction MIMO SAR can effectively suppress deception jamming.

Conclusions
In this paper, a range direction MIMO SAR system is proposed as a novel retransmission deception jamming suppression SAR. Based on concurrent transmission of sub-band signals with random initial phase modulation and sub-band synthesis in the equivalent phase center of the MIMO SAR, the jamming suppression capability for retransmission deception jamming is enhanced.
First, the sub-band signal synthesis processing model is built. Then, the efficiency of deception jamming suppression is analyzed. Because of the combination of the sub-band synthesis and random phase modulation, the retransmission deception jamming will defocus in the azimuth, and the interference to signal ratio will decrease. Finally, the simulation of deception jamming proves the validity of the jamming suppression method. Hence, the MIMO SAR will not only improve the spatial resolution and wider swath coverage, but will also enhance the jamming suppression capability. In the future, we will investigate non-coherent integration of different equivalent phase centers in the image domain to further improve the anti-jamming capability of the proposed method.

Acknowledgments:
The research is supported in part by the National Natural Science Foundation of China (Nos. 61132006, 61301187 and 61371133).
Author Contributions: Ruijia Wang designed the method and performed the experiment and simulation and wrote the manuscript; Jie Chen supervised the study and modified the manuscript; Xing Wang and Bing Sun supervised the study and wrote part of the manuscript.

Conclusions
In this paper, a range direction MIMO SAR system is proposed as a novel retransmission deception jamming suppression SAR. Based on concurrent transmission of sub-band signals with random initial phase modulation and sub-band synthesis in the equivalent phase center of the MIMO SAR, the jamming suppression capability for retransmission deception jamming is enhanced.
First, the sub-band signal synthesis processing model is built. Then, the efficiency of deception jamming suppression is analyzed. Because of the combination of the sub-band synthesis and random phase modulation, the retransmission deception jamming will defocus in the azimuth, and the interference to signal ratio will decrease. Finally, the simulation of deception jamming proves the validity of the jamming suppression method. Hence, the MIMO SAR will not only improve the spatial resolution and wider swath coverage, but will also enhance the jamming suppression capability. In the future, we will investigate non-coherent integration of different equivalent phase centers in the image domain to further improve the anti-jamming capability of the proposed method.