Next Article in Journal
Medical Devices with Embedded Sensor Systems: Design and Development Methodology for Start-Ups
Previous Article in Journal
An Online Hashing Algorithm for Image Retrieval Based on Optical-Sensor Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bistatic Radar Cooperative Imaging Based on Complementary Random Waveform

1
School of Information Science and Engineering, Zhejiang Sci-Tech University, Hangzhou 310018, China
2
School of Computer Science and Technology, Zhejiang Sci-Tech University, Hangzhou 310018, China
3
Institute of Land Aviation, Beijing 101121, China
4
College of Electronic Science and Engineering, National University of Defense Technology, Changsha 410073, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(5), 2577; https://doi.org/10.3390/s23052577
Submission received: 27 December 2022 / Revised: 28 January 2023 / Accepted: 23 February 2023 / Published: 25 February 2023
(This article belongs to the Section Communications)

Abstract

:
The cooperative imaging of the bistatic radar is an important research topic for missile-borne radar detection. The existing missile-borne radar detection system is mainly based on the target plot information separately extracted by each radar for data level fusion, without considering the gain brought by the cooperative processing of the radar target echo signal. In this paper, a random frequency-hopping waveform is designed for the bistatic radar to achieve efficient motion compensation. A coherent processing algorithm for bistatic echo signals is designed to achieve band fusion and improve the signal quality and range resolution of the radar. Simulation and high-frequency electromagnetic calculation data were used to verify the effectiveness of the proposed method.

1. Introduction

Radars with all-weather, all-day, long-range detection ability, are one of the important sensors in the onboard platforms for target imaging [1]. However, with the progress of target stealth technology and the rapid technical development of jamming countermeasures, such as target RCS reduction, which requires a radar to adopt a large antenna aperture and signal system with a large bandwidth, which brings great challenges to single radar detection [2]. Compared to single radar detection, multi-radar collaborative detection can not only obtain more observation periods but also obtain the characteristics of different dimensions, such as frequency character [3], polarization character [4], space character [5], etc., which is an important development direction for stealth target detection in a complex electromagnetic environment. According to the varied detection angles and radar working frequency bands, multi-radar cooperative detection mode can be divided into spatial cooperation and frequency band cooperation, to obtain the benefits of spatial diversity [6,7] and frequency diversity [8,9]. The former realizes target detection and recognition by jointly processing information for multiple detection angles to reduce the probability of misjudgment in single-view detection. The latter processes varied radar frequencies coherently to obtain higher-range resolution and gain by pulse compression.
For the application scenario of the missile-borne radar, multiple missile-borne radars accompany the flight, and the difference in detection perspective is small. Frequency band coordination is usually adopted to improve detection performance. Moreover, the missile-borne radar usually works in complex electromagnetic environments with a large number of jammers. In order to improve the detection performance, the missile-borne radar usually adopts a random frequency-hopping signal with excellent anti-jamming performance [10,11,12,13]. Therefore, based on the bistatic radar detection framework and the random frequency-hopping signal system, this paper studies the bistatic radar cooperative imaging method. Due to the movement of different radars with high speed and the distribution of echo data in discontinuous multiple frequency bands, it is necessary to eliminate the phase error introduced by the platform motion and compensate for the phase difference between the varied radars. Furthermore, the data fusion of different frequency bands needs to be realized. The process involves motion compensation, coherent registration, and data fusion [14,15,16]. In terms of motion compensation, reference [17] used the cross-correlation of adjacent range directions for velocity estimation. The compensation accuracy is determined by the bandwidth, which cannot be applied to radar waveforms sensitive to motion, such as random frequency-hopping. In [18], the extremum of the searching cost function determined the target velocity. However, the performance of the algorithm is sensitive to the step size and range of the searching step. In [19], a velocity estimation method based on the complementary code modulation was proposed to achieve precise velocity estimation by transmitting adjacent pulse trains with complementary relations. However, the radar data rate is reduced as two pulse trains are needed. In terms of coherent registration, [20] divided the compensated phase between subbands into a linear phase and a fixed phase and used the all-pole model to solve the problem. However, the pole order is difficult to determine. In [21], a cross-correlation method was used to solve the linear phase term, and FFT transformation was used to solve the fixed phase term, which avoided the pole order problem. In terms of frequency band data fusion, in [22,23], the defective frequency band was filled in a wideband signal for imaging. However, the band estimation was based on the existing frequency band without additional information. In [24], a probabilistic model was established for sparse-frequency band echo signals. After Doppler focusing on the target, a Bayesian learning algorithm was used for range fusion imaging to avoid the errors introduced by the frequency band filling algorithm. However, Doppler focusing requires accurate target motion compensation, and the compensation accuracy of the sub-frequency band needs further analysis.
In this paper, we studied the bistatic radar cooperative imaging based on a complementary random waveform. Firstly, we constructed a dual-platform cooperative detection geometric model of bi-transmitting and bi-receiving bistatic radars. Secondly, a random frequency-hopping signal waveform pair with complementary characteristics is proposed and analyzed. Then, the motion compensation of the target is carried out to realize the phase-coherent registration of the echo signal. Through the complementary characteristics of frequency-hopping coefficients, the distance and velocity terms of the target in the phase of the echo signal are separated, and then the target velocity is accurately estimated. Finally, a band fusion algorithm is designed to achieve two-dimensional imaging of the target. The radar imaging method studied in this paper has the advantages of the bistatic radar system and a random frequency-hopping signal [3,4,5,10,11,12,13]. It improves the detection performance and realizes the efficient estimation of target velocity, which lays a foundation for the future study of missile-borne radar collaborative detection imaging.
This paper is organized as follows. In Section 2, the waveform pairs of the bistatic random frequency-hopping signal with complementary characteristics are designed. A model of a bistatic radar echo is derived. In Section 3, the compensation phase between echoes is accurately calculated using the echoes of the homologous transmitted signals. The coherent registrations of echoes are effectively realized. Furthermore, an efficient velocity estimation method is designed based on the complementary characteristics of echoes to fulfill the relative motion compensation. Section 4 provides that the band fusion algorithm is based on tight constraint rearrangement and zero complement for bistatic radar cooperative imaging. The point-target simulation data and the high-frequency electromagnetic calculation data are used to verify the algorithms in Section 5. Finally, the conclusions are summarized in Section 6.

2. Echo Model of Bistatic Radar Cooperative Imaging

The bistatic radar cooperative imaging is shown in Figure 1. Suppose that both radars 1 and 2 adopt random frequency-hopping waveforms, and both coherent processing pulse trains contain N sub-pulses with the same frequency-hopping rules. The pulse repetition time (PRT) is T r . The pulse width is T . The random modulation frequency distribution is within a given bandwidth B . The minimum frequency-hopping (FH) step is Δ f = B / N . The FH coefficient is a random integer sequence of length N in [ 0 , N 1 ] . In order to make full use of the frequency band, the frequency-hopping coefficient usually traverses the interval and is not repeated. Suppose that the FH coefficients of radars 1 and 2 are c 1 ( n ) and c 2 ( n ) , respectively. The transmitted signal models of the two radars are expressed as follows:
s 1 ( t ) = n = 0 N 1 rect t n T r T exp j 2 π f 1 + c 1 ( n ) Δ f t n T r
s 2 ( t ) = n = 0 N 1 rect t n T r T exp j 2 π f 2 + c 2 ( n ) Δ f t n T r
where f 1 and f 2 are the carrier frequencies of radar 1 and 2, respectively.
In order to separate the varied radar signal components by the corresponding band pass filter, the carrier frequency difference of f 1 and f 2 for both radars should be greater than the radar signal modulation bandwidth. In Figure 1, the separated set of echoes is denoted as s 11 ( t ) , s 12 ( t ) , s 21 ( t ) , s 22 ( t ) , where s p q represents the pulse transmitted by radar p   ( p = 1 , 2 ) and the echo in radar q   ( q = 1 , 2 ) receiving mode. Suppose that the target contains K scattering centers. At the initial time, the distance between the k t h scattering center and radar 1 is r 1 k . The radial velocity of the target is v . When the n t h sub-pulse is emitted, the distance between the scattering center and radar 1 can be expressed as r 1 k ( n ) = r 1 k + v n T r . When radar 1 transmits and radar 1 receives, the n t h sub-pulse echo can be expressed as
s 11 n = k = 0 K 1 σ k exp j 4 π c ( f 1 + c 1 ( n ) Δ f ) r 1 k n = k = 0 K 1 σ k exp j 4 π c ( f 1 + c 1 ( n ) Δ f ) r 1 k + v n T r
where the σ k is the intensity of the k t h scattering center. For the scene of collaborative imaging of multiple missile-borne radars, the lines of sight from varied radars are almost the same. In other words, the intensity and velocity of a scattering center can be regarded as constant for the two radars.
Suppose that, at the initial time, the distance between the k t h scattering center and radar 2 is r 2 k . When radar 2 transmits and radar 2 receives, the n t h sub-pulse echo can be expressed as
s 22 n = k = 0 K 1 σ k exp j 4 π c ( f 2 + c 2 ( n ) Δ f ) r 2 k n = k = 0 K 1 σ k exp j 4 π c ( f 2 + c 2 ( n ) Δ f ) r 1 k + v n T r + Δ r
where Δ r = r 2 k r 1 k represents the electromagnetic wave transmission time corresponding to the distance difference between the two radars.
Similarly, the echo s 21 n received by radar 1 transmitted by radar 2 and echo s 12 n received by radar 2 transmitted by radar 1 are, respectively, expressed as
s 21 n = k = 0 K 1 σ k exp j 4 π c ( f 2 + c 2 ( n ) Δ f ) r 2 k + r 1 k 2 + v n T r = k = 0 K 1 σ k exp j 4 π c ( f 2 + c 2 ( n ) Δ f ) r 1 k + v n T r + Δ r 2
s 12 n = k = 0 K 1 σ k exp j 4 π c ( f 1 + c 1 ( n ) Δ f ) r 1 k + v n T r + Δ r 2

3. Motion Compensation Method Based on Complementary Random Waveform

3.1. Distance Difference Compensation

Comparing Equation (3) with Equation (6), the echo signals transmitted by one radar and received by another radar have phase modulation terms introduced by different distance delay. s 12 n can be represented by s 11 n as
s 12 n = exp j 2 π Δ r c ( f 1 + c 1 ( n ) Δ f ) k = 0 K 1 σ k exp j 4 π c ( f 1 + c 1 ( n ) Δ f ) r 1 k + v n T r = exp j 2 π Δ r c ( f 1 + c 1 ( n ) Δ f ) s 11 n = exp j ϕ ( n ) s 11 n
where ϕ ( n ) = 2 π Δ r ( f 1 + c 1 ( n ) Δ f ) / c and ϕ ( n ) can be calculated from s 12 n and s 11 n , and it is used to compensate the range difference of s 12 ( n ) . ϕ ( n ) can be expressed as
ϕ ( n ) = a n g l e s 12 ( n ) · c o n j s 11 ( n )
The compensation process of s 22 ( n ) and s 21 ( n ) is similar. The echo after range compensation can be expressed as
s 11 n = s 12 n = k = 0 K 1 σ k exp j 4 π c ( f 1 + c 1 ( n ) Δ f ) r 1 k + v n T r
s 21 n = s 22 n = k = 0 K 1 σ k exp j 4 π c ( f 2 + c 2 ( n ) Δ f ) r 1 k + v n T r
After the radar range difference is compensated, the echoes emitted by varied radars have the same law. Therefore, the collaborative imaging method is studied based on echoes s 11 n and s 21 n .

3.2. Echo Motion Compensation

It can be seen from the above echo model that there is a coupling term between the frequency variation and range variation of the random frequency-hopping waveforms, which will lead to the introduction of a higher order phase term in the echo. The high-resolution range profile obtained by direct FFT has a serious defocus. It is necessary to compensate for the target motion. By comparing Equations (9) and (10), it can be seen that s 11 n and s 21 n have the same range delay and different signal frequencies. If the random frequency-hopping mode of the two radars is designed reasonably, the influence of frequency modulation can be eliminated to achieve the rapid estimation of target motion parameters.
The random frequency-hopping modulation of the two radars is shown in Figure 2, where f 0 = f 1 + f 2 / 2 , c ( n ) = c 1 ( n ) + f 2 f 1 / 2 and c 2 ( n ) = c 1 ( n ) . Both random frequency-hopping signals from radars comply with the frequency rules. The frequency-hopping coefficient c ( n ) of the pulse transmitted by radar 1 equals to the one of radar 2, and the direction is opposite. It can be obtained by substituting the above equations into Equations (9) and (10).
s 11 n = exp j 4 π v n T r c ( f 0 + c ( n ) Δ f ) k = 0 K 1 σ k exp j 4 π c ( f 0 + c ( n ) Δ f ) r 1 k
s 21 n = exp j 4 π v n T r c ( f 0 c ( n ) Δ f ) k = 0 K 1 σ k exp j 4 π c ( f 0 c ( n ) Δ f ) r 1 k
In order to estimate the velocity of the target, the echo components s 11 and s 21 are multiplied and the following results were obtained.
s n = s 11 n · s 21 n = exp j 8 π c v n T r f 0 k = 0 K 1 l = 0 K 1 σ k σ l exp j 4 π c f 0 ( r 1 k + r 1 l ) + c ( n ) Δ f ( r 1 k r 1 l ) = exp j 8 π c v n T r f 0 k = 0 K 1 σ k 2 exp j 8 π f 0 c r 1 k + k = 0 K 1 l = 0 , l k K 1 σ k σ l exp j 4 π c f 0 ( r 1 k + r 1 l ) + c ( n ) Δ f ( r 1 k r 1 l ) = exp j 8 π c v n T r f 0 s 0 + s n
where k = 0 K 1 σ k 2 exp j 8 π f 0 c r 1 k is a constant value and s n varies with the number of sub-pulses n . Therefore, s n comprises a single-frequency signal component s 0 and a variable component s n . After the FFT transformation, the energy of the single-frequency signal component is coherently accumulated to form a peak, while s n is still in the defocus state. As the result, the target velocity can be estimated by the peak position of the sequence after FFT transformation.
v = k 0 c 4 N T r f 0
where k 0 is the peak position, N = 128 , T r = 10 4 s , f 0 = 15 G H Z . By Equation (14), the target velocity estimation precision is Δ v = c 8 N T r f 0 , with the above parameters, the target velocity estimation precision can be calculated as Δ v = 0.1953 m / s . To ensure the focus of the range imaging, the phase introduced by the velocity error during the imaging time should be no more than π / 4 , namely,
2 π Δ v N T r f 0 c π 4 Δ v c 8 f 0 N T r
The target velocity estimation accuracy meets the imaging requirements of the random frequency-hopping signal. After the target velocity is estimated, the target translational velocity can be compensated for Equations (11) and (12). The compensated echo can be expressed as
s 11 n = s 11 n exp [ j 4 π v n T r c ( f 0 + c ( n ) Δ f ) ] = exp j 4 π v n T r c ( f 0 + c ( n ) Δ f ) k = 0 K 1 σ k exp j 4 π c ( f 0 + c ( n ) Δ f ) r 1 k exp [ j 4 π v n T r c ( f 0 + c ( n ) Δ f ) ] = k = 0 K 1 σ k exp j 4 π c ( f 0 + c ( n ) Δ f ) r 1 k
s 21 n = s 21 n exp j 4 π v n T r c ( f 0 c ( n ) Δ f ) = exp j 4 π v n T r c ( f 0 c ( n ) Δ f ) k = 0 K 1 σ k exp j 4 π c ( f 0 c ( n ) Δ f ) r 1 k exp j 4 π v n T r c ( f 0 c ( n ) Δ f ) = k = 0 K 1 σ k exp j 4 π c ( f 0 c ( n ) Δ f ) r 1 k
Equations (16) and (17) describe the echo of target random frequency-hopping signals in different frequency bands. In the following, the coherent fusion high-resolution imaging algorithm of random frequency-hopping echoes in different frequency bands will be studied.

4. Band Fusion Imaging Based on Tight Constrained Rearrangement and Zero Complement

4.1. Coherent Processing

The existing frequency band fusion algorithms are mainly based on two ideas. One is to use band extrapolation to fill the defect band after imaging, which is based on the existing frequency band information with no additional information. Another idea is to use a modern spectrum, sparse representation, or the Bayesian method to extract scattering center parameters for super-resolution range imaging. Because the original information of the radar data is lost in the process of parameter extraction, it is difficult to carry out two-dimensional imaging.
In this paper, a tight constrained rearrangement and zero complement frequency band fusion algorithm is proposed. The main idea is to rearrange the random frequency-hopping in different frequency bands. Zeros are made up for the echo data at the defect frequency points. Moreover, tight constraints are set so that the echo signals of the two radars can be used for coherent fusion imaging by FFT transformation. The process of zero rearrangement and complement is shown in Figure 3. The echo sequences of Equations (16) and (17) are sorted from small to large according to the minimum FH step. The obtained echo sequence is expressed as
s a l l n = s a l l 0 , s a l l 1 , s a l l 2 , , s a l l L 1
where L = M + 2 N is the length of echo sequence after zero padding. The first N components of s a l l n are the results of the s 21 n rearrangement. The last N components are the results of the s 11 n rearrangement, which are all valid observations. The defect frequency band is f 2 , f 1 . In order to ensure the phase-coherence of the two radar FH signals, tight constraints are set: the defect bandwidth f 1 f 2 is an multiple of the FH step Δ f , and the corresponding number of zero-padding frequency points is M = f 1 f 2 / Δ f 1 .
The zero-padding position of s a l l n locates in the middle of defect frequency band. The range resolution of s a l l n is consistent with the performance of the echo sequence with continuous bandwidth 2B. The high-resolution range image obtained by FFT of s a l l n can be expressed as
h r r p k = F F T S a l l n
Thus, the range imaging results with high-resolution under bistatic cooperative imaging can be obtained. The two-dimensional range image of the target with high-resolution can be obtained by further combining the envelope alignment and the initial phase correction algorithm.

4.2. Imaging Process

The process of bistatic radar cooperative imaging based on complementary random waveform is shown in Figure 4. First, the echo signals of two radars that meet the complementary characteristics and tight constraints are designed. Then varied band pass filters are set according to the difference of carrier frequencies to obtain the separated target echo sequences s 11 ( t ) , s 12 ( t ) , s 21 ( t ) , s 22 ( t ) .
After obtaining the echo components of the four paths, the modulation phase term introduced by the difference in distance delay is calculated under different paths according to Equation (7). It is used to compensate for the additional phase introduced by the distance difference under different paths. The distance difference correction of different paths is completed. Four-path echo signals with the same distance delay are obtained (As shown in Equations (9) and (10)).
With four-path echo signals, the same range delays are obtained. After multiplication, combined with FFT, the target speed can be effectively estimated. The target echo is compensated based on the estimated velocity, and the random frequency-hopping echo sequences of different frequency bands are obtained. Then, the range imaging of random frequency-hopping echoes in different frequency bands with high-resolution is realized by rearranging and zero padding. Finally, combined with the envelope alignment and the initial phase correction methods in the ISAR imaging process, the two-dimensional imaging of the bistatic cooperative imaging with high-resolution is completed.
The distance difference correction process requires 2 M multiplication and N addition. The target motion compensation process requires 2 M multiplication, N addition and one F F T . The frequency band fusion imaging requires M multiplication and 2 F F T . Thus, the complexity of the proposed imaging algorithm is o ( 3 M N log ( M N ) + 5 M + 2 N ) . Where, N is the number of pulses and M is the number of sampling points.

5. Experimental Results

5.1. Validation of Simulation Data

In order to verify the effectiveness of the dual base cooperative imaging method, a simulation is adopted. The point scattering model of a target ship is shown in Figure 5, which is composed of 367 scattering points. The ship size is l = 120 , b = 30 . The carrier frequency of radar 1 is f 1 = 15 GHz . The carrier frequency of radar 2 is f 2 = 15.255 GHz . The signal bandwidth is B = 128 MHz and Δ f = 1 MHz . The modulated signal c 1 ( n ) from radar 1 is a sequence of non-repeating random integers of length N in interval 0 , N 1 . The modulated rule from radar 2 is c 2 ( n ) = c 1 ( n ) . The geometric scene of the bistatic cooperative imaging is shown in Figure 6, with the ship located at the origin of the coordinate system O X Y Z . The bistatic imaging platform flies forward along the Y-axis at altitude H = 1.73 k m and has the flight velocity v r = 225 m / s . The distance between radar 1 and the center of the ship at the initial moment is 10 km. The target is in a maneuvering state, and the projected component of the moving speed on the radar line of sight is v t L O S = 3 m / s . At the initial moment, the distance between radar 2 and the center of the ship is 10.017 km. The azimuth angle θ and pitch angle β in the target coordinate system are θ = 80 o and β = 10 o , respectively.
In this paper, the complementary random frequency-hopping waveform designed is used to detect the target. The target echo sequence after separation using the band pass filter is expressed as s 11 ( t ) , s 12 ( t ) , s 21 ( t ) , s 22 ( t ) . Equation (8) is used to calculate the phase introduced by radar range difference, as shown by the black asterisk in Figure 7a, and the theoretical calculation results of this phase are also given in the figure by red circles. In Figure 7a, the proposed method can effectively eliminate the phase term introduced by the platform distance difference. Figure 7b shows the FFT transformation result of s n . By Equation (14), the target velocity estimated according to the peak position is 3.099 m/s. It meets the requirements of Equation (15) Δ v 0.1953 m / s . Therefore, target range imaging has good focusing performance. Figure 7c and d show the range imaging results of the two radars, respectively. Due to the limited bandwidth, the scattering centers are distributed in a small number of range units. Few target details can be obtained. The range image with high resolution fused by the proposed method is shown in Figure 7e. Meanwhile, the high-resolution range image with ideal compensation of target velocity is given in Figure 7e. In Figure 7c and e compared to, for example, the fused range image provides richer target details than the single radar range image, and the estimated speed accuracy can meet the focusing requirements of high-resolution range imaging.
The two-dimensional radar imaging results obtained by the envelope alignment and initial phase correction using the envelope cross-correlation method and PGA method are shown in Figure 8. Among them, Figure 8a and b show the independent imaging results of radars 1 and 2, respectively. Due to the limitation of the range resolution, the targets are distributed in about 15 range units in the range direction. The scattering center distribution is fuzzy, and the target contour is unclear. The fusion imaging results adopting the model are shown in Figure 8c. As the signal band grows, the target scattering center distribution area becomes larger, about 50 cells, providing more target details. ISAR imaging under ideal compensation of speed is presented to verify the focusing performance of the image by contrast [25]. Among them, the contrast of the actual fused ISAR image is 16.38, while the one with the ideal compensated ISAR image is 17.26. Figure 9a and b show the comparison of imaging results in a noisy environment (SNR = 0 d B ). The contrast of the actual fused image in a noisy environment is 7.88. Moreover, the contrast of the ideal compensated ISAR image in a noisy environment is 8.45. Using AMD Ryzen 7 5800H with Radeon Graphics 3.20 GHz CPU and MATLAB 2020a, the average computation time for the proposed imaging algorithm is 9.18 s. The calculation time of the single base radar imaging algorithm is 8.82 s. In practical applications, the computing speed can be greatly improved by improving hardware configuration and the GPU parallel acceleration. Thus, the experimental results show that the proposed method can effectively realize the motion compensation and obtain the fusion images of the target.

5.2. Data Validation of High-Frequency Electromagnetic Calculations

The experimental data are the ship model shown in Figure 9. The length and width of the ship are l = 214.6 m and b = 28.04 m, respectively. The carrier frequencies of radars 1 and 2 are f 1 = 10.1 GHz and f 2 = 10.2 GHz , respectively. The signal bandwidths are B = 100 MHz , Δ f = 0.3125 MHz and the number of frequency-hopping points is N = 320 . The modulated signal from radar 1 c 1 ( n ) is a non-repeated random integer sequence of length N in the interval 0 , N 1 . The modulated rule of radar 2 can be expressed as c 2 ( n ) = c 1 ( n ) . The azimuth angle θ and the pitch angle β of the line of sight of the radar in the ship body coordinates system are 30° and 10°, respectively. The initial ranges of radars 1 and 2 are 12 km and 12.02 km , respectively. Suppose that the speed of the radar platform has been compensated by GPS-assisted data. The target velocity is v t L O S = 4 m / s . The accumulated number of pulses of two-dimensional imaging is 960. The corresponding azimuth synthetic aperture is 1.72.
Figure 10 is the model of ship. Figure 11 is comparison of imaging results. Figure 11a,b show the estimated results of the phase introduced by the range difference and the target velocity, respectively. The target velocity estimated from the peak position is −3.925 m/s. It meets the requirements of Equation (15). Figure 11c and d show the independent imaging results of radars 1 and 2, respectively. Due to the limitation of range resolution, the target is distributed around 120 range units in the range direction. The sidelobe of the scattering center has a great influence, resulting in a long trailing. The results obtained by the fusion imaging mode are shown in Figure 11e. The target is distributed around 360 range units in the range direction. With the increase in the number of effective coherent frequency points, the sidelobe of the target scattering center decreases significantly, and the contour becomes clearer. Figure 12a shows the local magnification of the radar 2 ISAR imaging, and Figure 12b shows the local magnification of the fusion ISAR imaging after ideal compensation. By comparing the local magnification area of the two images, the image with ideal compensation has a better focusing effect of scattering center than that of the single base radar imaging. Compared with the ISAR imaging results under ideal velocity compensation (Figure 11f), the estimated contrast of the fused ISAR image after velocity compensation is 104.8591, while the contrast of the ideal compensation ISAR image is 106.4266. The proposed method can realize bistatic cooperative imaging provides high-quality target information for target recognition.

6. Conclusions

Bistatic radar cooperation can effectively improve the detection performance of target characteristics. In this paper, a signal-level cooperative imaging method based on random complementary waveforms was studied based on the two-transmitter and two-receiver cooperative detection system of the bistatic radar. Firstly, the echo signals of the two radars satisfy the complementary characteristics and tight constraints are designed. The echo signal model is derived. The compensation phase between echoes is accurately calculated by using the echoes of the homologous transmitted signals. The coherent compensation of the echoes is completed. Furthermore, an efficient velocity estimation method is designed based on the complementary characteristics of the frequency modulation. The estimated result of the target velocity is within the range of velocity accuracy. Finally, a band fusion algorithm based on tight constraint rearrangement and zero complement is used to realize two-dimensional imaging of the bistatic radar with high-resolution. The point-target simulation data and high-frequency electromagnetic calculation data are used to verify the effectiveness of the proposed method. By comparing with the imaging results of the single base radar, the proposed imaging algorithm can provide more target details. The proposed method is expected to provide beneficial exploration and basic algorithm support for the cooperative imaging of missile-borne, dual-based radars.

Author Contributions

Software, Z.S.; Formal analysis, B.Z. and Y.Z.; Writing—original draft, X.L.; Writing—review & editing, J.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Key Laboratory of National Defense Science and Technology Fund, grant number WDZC20205500208 and HTKJ2022KL504016.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data is unavailable due to privacy restrictions.

Acknowledgments

We are very grateful for the support and help from Long Wu of Zhejiang Sci-Tech University. Thanks for Ziyi Liu of ZSTU help.

Conflicts of Interest

The authors declare no conflict of interest.

Correction Statement

This article has been republished with a minor correction to the Funding statement. This change does not affect the scientific content of the article.

References

  1. Huang, P.; Yang, H.; Zou, Z.; Xia, X.-G.; Liao, G. Multichannel Clutter Modeling, Analysis, and Suppression for Missile-Borne Radar Systems. IEEE Trans. Aerosp. Electron. Syst. 2022, 58, 3236–3260. [Google Scholar] [CrossRef]
  2. Sun, W.; Zhang, Y. Performance Analysis on Coherent Receiving Radar Based on Electromagn-etic Scattering Characteristics of Target. In Proceedings of the 2019 IEEE International Conference on Computational Electromagnetics (ICCEM), Shanghai, China, 20–22 March 2019; pp. 1–3. [Google Scholar] [CrossRef]
  3. Mei, H.; Li, Y.; Xing, M.; Quan, Y.; Wu, C. A Frequency-Domain Imaging Algorithm for Transla-tional Variant Bistatic Forward-Looking SAR. IEEE Trans. Geosci. Remote Sens. 2020, 58, 1502–1515. [Google Scholar] [CrossRef]
  4. Wang, Y.; Chen, W.; Mao, X.; Lei, W. Terrain Classification of Polarimetric SAR Images Based on Optimal Polarization Features. In Proceedings of the 2022 International Conference on Computer Engineering and Artificial Intelligence (ICCEAI), Shijiazhuang, China, 22–24 July 2022; pp. 400–403. [Google Scholar] [CrossRef]
  5. Zhong, H.; Liu, H.; Qi, G. Analysis of Terminal Area Airspace Operation Status Based on Traje-ctory Characteristic Point Clustering. IEEE Access 2021, 9, 16642–16648. [Google Scholar] [CrossRef]
  6. Strelnikov, A.P.; Volkov, A.S.; Bakhtin, A.A.; Gorelik, A.V.; Kobzev, V.A. Alamouti Scheme a-nd Spatial Diversity MIMO Algorithms. In Proceedings of the 2021 IEEE East-West Design & Test Symposium (EWDT-S), Batumi, Georgia, 10–13 September 2021; pp. 1–6. [Google Scholar] [CrossRef]
  7. Haimovich, A.M.; Blum, R.S.; Cimini, L.J. MIMO Radar with widely separated antennas. IEEE Signal Process. Mag. 2008, 25, 116–129. [Google Scholar] [CrossRef]
  8. Liu, H.; Zhou, S.; Su, H.; Yu, Y. Detection performance of spatial-frequency diversity MIMO radar. IEEE Trans. Aerosp. Electron. Syst. 2014, 50, 3137–3155. [Google Scholar]
  9. Zhang, Y.; Wang, T.-J.; Zhao, H.-P.; Zhang, Y.; Zhao, H. Multiple radar subbands fusion algorithm based on support vector regression in complex noise environment. IEEE Trans. Antennas Propag. 2018, 66, 381–392. [Google Scholar] [CrossRef]
  10. Liu, S.; Jia, Y.; Liu, Y.; Zhang, X. Research on Ultra-Wideband NLFM Waveform Synthesis and Grating Lobe Suppression. Sensors 2022, 22, 9829. [Google Scholar] [CrossRef] [PubMed]
  11. Zhang, J.; Hou, C.; Lin, Y.; Zhang, J.; Xu, Y.; Chen, S. Frequency Hopping Signal Modulation Recognition Based on Time-Frequency Analysis. In Proceedings of the 2021 IEEE 18th International Conference on Mobile Ad Hoc and Smart Systems (MASS), Denver, CO, USA, 4–7 October 2021; pp. 46–52. [Google Scholar] [CrossRef]
  12. Zhang, J.; Yuan, B.; Jiang, Z.; Guo, Y.; Wang, D. Microwave Staring Correlated Imaging Method Based on Steady Radiation Fields Sequence. Sensors 2020, 20, 6859. [Google Scholar] [CrossRef] [PubMed]
  13. He, S.; Wu, X. Research of a Radar Imaging Algorithm Based on High Pulse Repetition Random Frequency Hopping Synthetic Wideband Waveform. Sensors 2019, 19, 5424. [Google Scholar] [CrossRef]
  14. Chen, Y.; Chen, R.; Liu, H.; Guo, J.; Wang, Y.; Zhang, J. Research on Configuration Constraints of Airborne Bistatic SARs. Sensors 2022, 22, 6534. [Google Scholar] [CrossRef]
  15. Zhang, J.; Ding, T.; Zhang, L. Longtime Coherent Integration Algorithm for High-Speed Maneuvering Target Detection Using Space-Based Bistatic Radar. IEEE Trans. Geosci. Remote Sens. 2022, 60, 1–16. [Google Scholar] [CrossRef]
  16. Hong, S.; Moon, W.M.; Paik, H.Y.; Choi, G.H. Data fusion of multiple polarimetric SAR images using discrete wavelet transform (DWT). In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Toronto, ON, Canada, 24–28 June 2002; IEEE: New York, NY, USA, 2002; Volume 6, pp. 3323–3325. [Google Scholar]
  17. Liao, Z.-K.; Hu, J.-M.; Lu, D.-W.; Zhang, J. Motion analysis and compensation method for random stepped frequency radar using the Pseudorandom code. IEEE Access 2018, 6, 57643–57654. [Google Scholar] [CrossRef]
  18. Li, X.; Liu, G.-S.; Ni, J. Autofocusing of ISAR images based on entropy minimization. IEEE Trans. Aerosp. Electron. Syst. 1999, 35, 1240–1252. [Google Scholar] [CrossRef]
  19. Liao, Z.-K.; Lu, D.-W.; Hu, J.-M.; Zhang, J. Waveform design for random stepped frequency radar to estimate object velocity. Electron. Lett. 2018, 54, 894–896. [Google Scholar] [CrossRef]
  20. Tian, B.; Chen, Z.-P.; Xu, S.-Y. Sparse sub-band fusion imaging based on parameter estimation of geometrical theory of diffraction model. IET Radar Sonar Navig. 2014, 8, 318–326. [Google Scholar] [CrossRef]
  21. Tian, J.-H.; Sun, J.-P.; Wang, G.-H.; Wang, Y.; Tan, W. Multiband radar signal coherent fusion processing with IAA and apFFT. IEEE Signal Process. Lett. 2013, 20, 463–466. [Google Scholar] [CrossRef]
  22. Bai, X.-R.; Zhou, F.; Wang, Q.; Xing, M.; Bao, Z. Sparse sub-band imaging of space targets in high-speed motion. IEEE Trans. Geosci. Remote Sens. 2013, 51, 4144–4154. [Google Scholar] [CrossRef]
  23. Hu, P.-J.; Xu, S.Y.; Wu, W.-Z.; Chen, Z. Sparse sub-band ISAR imaging based on autoregressive model and smoothed ℓ0 algorithm. IEEE Sens. J. 2018, 18, 9315–9323. [Google Scholar] [CrossRef]
  24. Zhou, F.; Bai, X.-R. High-resolution sparse sub-band imaging based on Bayesian learning with hierarchical priors. IEEE Trans. Geosci. Remote Sens. 2018, 56, 4568–4580. [Google Scholar] [CrossRef]
  25. Yang, T.; Shi, H.; Guo, J.; Liu, M. FMCW ISAR Autofocus Imaging Algorithm for High-Speed Maneuvering Targets Based on Image Contrast-Based Autofocus and Phase Retrieval. IEEE Sens. J. 2020, 20, 1259–1267. [Google Scholar] [CrossRef]
Figure 1. Bistatic radar collaborative imaging.
Figure 1. Bistatic radar collaborative imaging.
Sensors 23 02577 g001
Figure 2. Bistatic radar random frequency-hopping complementary modulation.
Figure 2. Bistatic radar random frequency-hopping complementary modulation.
Sensors 23 02577 g002
Figure 3. Random frequency-hopping echo rearrangement fills zero.
Figure 3. Random frequency-hopping echo rearrangement fills zero.
Sensors 23 02577 g003
Figure 4. Collaborative imaging flow of the bistatic radar based on complementary random waveform.
Figure 4. Collaborative imaging flow of the bistatic radar based on complementary random waveform.
Sensors 23 02577 g004
Figure 5. Target point scattering model.
Figure 5. Target point scattering model.
Sensors 23 02577 g005
Figure 6. Schematic diagram of the geometric scene of dual-base cooperative imaging.
Figure 6. Schematic diagram of the geometric scene of dual-base cooperative imaging.
Sensors 23 02577 g006
Figure 7. Comparison of fusion distance imaging effect.
Figure 7. Comparison of fusion distance imaging effect.
Sensors 23 02577 g007
Figure 8. Comparison of ISAR imaging results.
Figure 8. Comparison of ISAR imaging results.
Sensors 23 02577 g008
Figure 9. Imaging results in noisy environment.
Figure 9. Imaging results in noisy environment.
Sensors 23 02577 g009
Figure 10. Model of ship.
Figure 10. Model of ship.
Sensors 23 02577 g010
Figure 11. Comparison of imaging results.
Figure 11. Comparison of imaging results.
Sensors 23 02577 g011
Figure 12. Locally magnified contrast images.
Figure 12. Locally magnified contrast images.
Sensors 23 02577 g012
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, X.; Hu, J.; Zou, B.; Zhu, Y.; Song, Z. Bistatic Radar Cooperative Imaging Based on Complementary Random Waveform. Sensors 2023, 23, 2577. https://doi.org/10.3390/s23052577

AMA Style

Li X, Hu J, Zou B, Zhu Y, Song Z. Bistatic Radar Cooperative Imaging Based on Complementary Random Waveform. Sensors. 2023; 23(5):2577. https://doi.org/10.3390/s23052577

Chicago/Turabian Style

Li, Xin, Jiemin Hu, Bo Zou, Yongfeng Zhu, and Zhiyong Song. 2023. "Bistatic Radar Cooperative Imaging Based on Complementary Random Waveform" Sensors 23, no. 5: 2577. https://doi.org/10.3390/s23052577

APA Style

Li, X., Hu, J., Zou, B., Zhu, Y., & Song, Z. (2023). Bistatic Radar Cooperative Imaging Based on Complementary Random Waveform. Sensors, 23(5), 2577. https://doi.org/10.3390/s23052577

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop