Modeling of Monopulse Radar Signals Reflected from Ground Clutter in a Time Domain Considering Doppler Effects

Article information

J. Electromagn. Eng. Sci. 2020;20(3):190-198
Publication date (electronic) : 2020 July 31
doi : https://doi.org/10.26866/jees.2020.20.3.190
1Deparment of Electronic Engineering, Inha University, Incheon, Korea
2Agency for Defense Development, Daejeon, Korea
3Hanwha Systems, Yongin, Korea
*Corresponding Author: Il-Suek Koh (e-mail: ikoh@inha.ac.kr)
Received 2019 May 31; Revised 2019 October 01; Accepted 2020 May 12.

Abstract

To evaluate the performance of a monopulse radar system, it is necessary to accurately model radar return signals from the ground surface in a time domain. In this paper, we propose a numerical method to model these return signals including radar radio frequency specifications, such as the pulse repetition interval, frequency, and polarization, as well as the antenna geometry, the ground clutter backscattering characteristics, and so on. The Doppler effect is also incorporated into the signal generation scheme because of the dynamics of the platform/clutter and the antenna orientation. The Doppler frequency shift caused by the ground clutter is modeled by employing the time correlation of the received signals. In some scenarios, the monopulse signals are generated and numerically examined. For real radar application, the effect of the platform’s roll stabilization on the monopulse signals is investigated based on the proposed signal generation scheme.

I. Introduction

The radar system mounted on an airborne platform operates in an air-to-ground mode as well as in an air-to-air mode [1]. For example, the AN/APG-77 is the F-22’s primary sensor that provides several functions, including terrain mapping, ground-moving target searching/tracking, and air-to-ground ranging (AGR) [2]. The AN/APG-77 radar adopts a monopulse system with an active electronically scanned array (AESA) antenna. The amplitude-or phase-comparison monopulse radar transmits one pulse and then receives three monopulse channels as the sum, elevation, and azimuth difference channels. The channels are calculated as three weighted sums of partial signals received by each antenna element in a receiving mode [3, 4]. Also, in the air-to-ground mode, the clutter backscattering may become much stronger than that of the air-to-air mode, and such strong backscattering can seriously degrade the radar’s performance [5].

The radar echo signals are determined by basic radar radio frequency (RF) specification, such as frequency and pulse repetition interval (PRI), etc., and the radar’s operation environment, such as platform dynamics, ground clutter type, antenna foot-print on the ground plane (mainly determined by the platform orientation), and the ASEA antenna’s beam illumination. Due to the dynamics of the platform and ground clutter, the received signal is slightly distorted, which can be modeled by considering the Doppler frequency shift. Since the range-Doppler (RD) map clearly shows the Doppler frequencies, RD-processing is crucial in many radar applications, including a synthetic aperture radar (SAR) [6] and a moving target indicator [7]. Hence, the modeling of the radar echo signal in a time-domain is required when considering the effect of all the aforementioned factor’s needed to generate an accurate RD map.

Much attention has been paid to the backscattering determined by the ground clutter, which can be used to determine the mean power received by the radar. Since the clutter signal is a typical random signal, several probability density functions (PDF) have been used, such as Gaussian, Weibull, K-distribution, and so on [8]. Except in the case of sea clutter, the Gaussian distribution is a simple and accurate way to develop a time-domain signal model [9]. The correlated Gaussian random signal can be generated in a relatively easy and accurate manner [10, 11], and is a conventional method for incorporating the Doppler effect of the clutter into the modeling of pulse-Doppler radar signals [12].

In Section II, we address the proposed modeling approach of the radar echo signal for a phase-comparison monopulse radar [3], considering the Doppler effects due to the dynamics of the platform and the ground clutter simultaneously. Section III presents the numerical results for operating scenarios that consider the platform dynamics and antenna beam steering direction, etc. Finally, an AGR simulation is examined for a real-world application.

II. Modeling Method

1. Scenario and Coordinate Systems

The dynamic scenario for the air-to-ground radar operation is shown in Fig. 1. An airborne platform moves at a constant velocity, v⃗ = vxx̂ + vyŷ + vz, while repeatedly transmitting a signal to the ground surface at every PRI. Given that an antenna beam is illuminated toward the line of sight (LOS), the return signals may be reflected primarily from the main-lobe beam projected on the ground that is truncated by the half-power beam width (HPBW) for a numerical simulation. The return signals of multiple pulses are aligned along the fast time (t) and the slow time (η), generating 2D raw signals [13]. The mean reflection power from the ground clutter inside the projected main-lobe beam can be calculated using the backscattering coefficient [14]. Here, we assume the ground clutter to be a rough surface [15] for which a theoretical backscattering model, such as an Integral Equation Model (IEM) may be accurate [16] for two co-polarized scatterings. The return signal can be assumed to be a Gaussian random signal whose backscattering coefficient is used as a variance of the random signal. The movement of the platform and ground clutter can be modeled using the Doppler frequency shift in the received signals.

Fig. 1

Dynamics scenario for air-to-ground radar.

Note that the shape of the projected main-lobe beam in Fig. 1 is dependent on the orientations of the platform and the antenna. Hence, the coordinate systems of the platform and the antenna are important. The relationship between the global and local coordinate systems of the platform is shown in Fig. 2. The parameters related to the coordinate systems in Fig. 2 are summarized as follows:

Fig. 2

Local and global coordinate systems of the platform and antenna.

As the local coordinate of the antenna coincides with that of the platform, ŝl is defined as

s^l=sin θcos ϕx^p+sin θsin ϕy^p+cos θz^p.

Considering the dynamics of the platform, the LOS vector of the antenna in the global coordinates is rotated depending on θr θy and θp as seen in Fig. 2. Based on θr θy and θp: the antenna’s LOS vector in the global coordinate, ŝg, is represented as

s^g=Rx(θr)Ry(θy)Rz(θp)s^l

where Rx, Ry, and Rz are rotation matrices along the global x, y, and z axes, respectively [17]. In Fig. 3, the signs of the antenna weights for three monopulse channels are defined. The V and W axes of the monopulse radar coincide with the yg and zg axes in the global coordinates, respectively [18]. In Figs. 4 and 5, the projected main-lobe beam and the V and W axes are found on the ground plane using a ray-tracing technique: first, the normal vector of the antenna’s VW plane is aligned in the direction of ŝg. Then, rays are launched from the antenna’s center to the rim of the HPBW of the antenna pattern or the V and W lines in the antenna’s VW plane. Finally, the intercept points between the rays and the ground plane are calculated.

Fig. 3

Sign of antenna weight for phase-comparison monopulse sum and difference patterns: (a) sum (w), (b) elevation difference (wΔEL ), and (c) azimuth difference ( wΔAZ ).

Fig. 4

Simulation results according to the azimuth angle ϕ. (a) Projected antenna footprint and V and W axes, and (b) analytic and simulation results of the V and W axis widths.

Fig. 5

Simulation results according to the rolling angle θr. (a) Projected antenna footprint and the V and W axes, and (b) analytic and simulation results of γ.

The simulation parameters are presented in Table 1. In Fig. 4, the platform orientation and antenna beam steering are given as θr = 0º, θy = 0º, θp = 0º, and θ = 105º, ϕ = 0º–60º, respectively. As observed in Fig. 4(a), it is well-known that as the azimuth angle (ϕ) increases, the HPBW of the V axis broadens while the HPBW of the W axis remains constant [19]. The simulation and theoretical estimations for the HPBW of the V and W axes on the ground are compared in Fig. 4(b), where an excellent agreement can be observed. As shown in Fig. 5(a), if the platform is rolled, the monopulse axes rotate by a small angle, γ, with respect to the major axes parallel to the LOS vector [14]. The estimations of γ by simulation and by the known theoretical formula [14] are in excellent agreement, as seen in Fig. 5(b). For this calculation, the platform orientation angles are assumed as θr = 0º–60º, θy = 0º, and θp = 0º. The beam illumination angles of the antenna are fixed at θ = 130º and ϕ = 0º. Therefore, the projected antenna footprint on the ground plane can be accurately found for any platform and antenna orientation.

Simulation parameters

2. Radar Return Signal in a Time-Domain

The configuration of the signal transmission/reception at the radar and the corresponding projected main-lobe beam needed to calculate the radar return signals in a time-domain are shown in Fig. 6. First, the projected main-lobe beam is divided into uniform rectangular grids (cells) based on a preset range resolution. The range resolution is determined using the sampling time interval of the fast time (t) [13]. The projected area can be approximated as an ellipse [20] whose major axis is parallel to the LOS direction from the radar to the center of the ellipse. Since the minor axis is perpendicular to the major axis, the ellipse can be uniformly discretized using a rectangle along the major and minor axes, as seen in Fig. 6. The ellipse is located in the far-field region from the platform, and so the return signals from the cells along a line parallel to the minor axis are coherently added into one time-domain response in each sampling time (range bin). The received monopulse signals from each q-th ground cell is represented as

Fig. 6

Signal transmission/reception of antenna to/from surface-projected main-lobe beam.

(1) ΨqΣ(tn,η)=AFΣAFΣG2σq(tn,η)A(q)Rc2(tn,η)δ[t-2Rc(tn,η)c],ΨqΔEL(tn,η)=AFΣAFΔELG2σq(tn,η)A(q)Rc2(tn,η)δ[t-2Rc(tn,η)c],ΨqΔAZ(tn,η)=AFΣAFΔAZG2σq(tn,η)A(q)Rc2(tn,η)δ[t-2Rc(tn,η)c],

where σq and A(q) are the backscattering coefficient and area of the q-th cell, respectively, G is the gain of the antenna element, and δ (·) is the Dirac delta function. Note that the q-th cell belongs to the sampling time (tn), as shown in Fig. 6. The sampling time can be computed simply as tn = 2Rc (tn, η)/c, where c is the speed of light. Rc (tn, η) is the distance from the platform to the q-th cells corresponding to the major axis, as shown in Fig. 6. AF, AFΔEL, and AFΔAZ are the array factors of the antenna for the sum, elevation difference, azimuth difference channels, respectively, on the q-th cell, which is given as

(2) AFΣ=p=1NpwpΣexp (-jk0[R^q(tn,η)-s^g]·rp),AFΔEL=p=1NpwpΔELexp (-jk0[R^q(tn,η)-s^g]·rp),AFΔAZ=p=1NpwpΔAZexp (-jk0[R^q(tn,η)-s^g]·rp),

where Np is the number of the antenna elements and k0 is the free-space propagation constant. q (tn, η) is the direction vector from the phase center of the antenna to the q-th cell at t = tn · p is the displacement vector from the phase center of the antenna to the p-th antenna element. wpΣ,wpΔEL, and wpΔAZ are the weight of the p-th antenna element for three monopulse channels, as seen in Fig. 3. In monopulse systems, a sum channel is used at the transmitting (Tx) mode and three channels, such as the sum, elevation difference, and azimuth difference channels, are employed in the receiving (Rx) mode [4]. The total return signals sampled at t=tn are a summation of the return signal from specific ground cells. The calculated cells at t=tn are mathematically denoted as Cell (q) ∈ tn. Adding all responses from the ground cells belonging to Cell(q) ∈ tn, three monopulse signals can be formulated without the consideration of the Doppler effects as

(3) ΨΣ,ΔEL,ΔAZ(tn,η)=Cell(q)tnTerrainΨqΣ,ΔEL,ΔAZ(tn,η)

where Ψ, ΨΔEL, and ΨΔAZ are the total radar return signals for the monopulse sum, elevation difference, and aizumuth difference channels, respectively.

In Fig. 7, the return signal is calculated from a flat ground whose backscattering is assumed to be a constant at a slow time, η = 0. In Fig. 7(a) and 7(b), the platform orientation and antenna beam steering are assumed as θr = 60º, θy = 0º, θp = 0º, and θ = 130º, ϕ = 0º, respectively. The “/w γ” and “/wo γ” indicate “with γ-correction” and “without γ-correction”, respectively. In Fig. 7(a), the roll angle is not zero, so it can be observed that the null of ΨΔEL locates off the boresight and the monopulse response of ΨΔAZ is not zero. But, as shown in Fig. 7(b), when the return signal is compensated with the known γ-correction scheme, the null point of ΨΔEL is shifted to the boresight and the level of ΨΔAZ is decreased to almost zero [4]. The exact and simulated monopulse ratios, Im(ΨΔEL )/Re(Ψ), are compared in Fig. 7(c) for θr = 0º [14]. The proposed signal-generation scheme can provide an accurate monopulse ratio. The Re(Ψ) is approximately 2 × 10−19, which is essentially zero. Here, Re(z) and Im(z) are the real and imaginary parts of z, respectively. Therefore, the generated time-domain signal can provide accurate monopulse signal behaviors.

Fig. 7

Monopulse response vs angle-off boresight: (a) monopulse signal before γ-correction, (b) monopulse signal after γ-correction, and (c) monopulse ratio of Im(ΨΔEL ) to Re(Ψ).

3. Doppler Effects

There are two dominant Doppler effects in the dynamic situation considered here because of to the movements of the platform and ground clutter. The Doppler shift of the clutter can be included simply by generating the correlated Gaussian random variable of σq (tn, η) in the slow time in (1) for a proper correlation function, such as the Gaussian and exponential functions. Here, the Cholesky decomposition method is used [11]. The correlation time (tc) is inversely proportional to the Doppler frequency of the ground clutter [21, 22]. After computing (1) with the correlated σq (tn, η)along the slow time, the Doppler shift due to the platform dynamic is considered. The time-dependent distance between the platform and the q-th cell is expressed at the slow time (ηm) as

(4) Rq(tn,ηm)=[xf+vx(tn+ηm)-xq]2+[yf+vy(tn+ηm)-yq]2+[zf+vz(tn+ηm)-zq]2

where (xf, yf, zf ) and (xq, yq, zq ) are the coordinates of the platform and the center of the q-th cell at (t = 0,η = 0) of the global coordinates, respectively. By using the Taylor series expansion, the Doppler frequencies along the fast-time, fD, tn, and slow-time, fD, ηm axes, can be formulated as

(5) fD,tn=-Rq(0,ηm)·vRq(0,ηm)λ0,fD,ηm=-Rq(0,0)·vRq(0,0)λ0

where λ0 is the free-space wavelength. Finally, the total return signal, including the Doppler effect due to platform movement, is given as

(6) ΨΣ,ΔEL,ΔAZ(tn,ηm)=Cell(q)tnTerrain[ΨqΣ,ΔEL,ΔAZ(tn,η)×exp(-j2πfD,tntn)×exp(-j2πfD,ηmηm)].

If two Doppler effects are considered at the same time, it is physically expected that the correlation time of the generated signal may be similar to the shorter one between the correlation times due to the platform and clutter dynamics.

III. Simulation Results

The simulation parameters are summarized in Table 1. The parameters of correlation length, root mean square (RMS), and dielectric constant of the rough surface are found in [5]. The pulse repetition frequency (PRF) is the reciprocal of PRI. A complete signal generation is considered, including two Doppler effects. For the rest of the simulations, the platform orientation and antenna beam steering are set as θr = 45º, θy = 0º, θp = 0º, and θ = 130º, ϕ = 0º, respectively. In Fig. 8, the RD maps are shown for the sum and elevation difference signals. The RD map is calculated by taking the fast Fourier transform (FFT) along the slow-time axis, which clearly shows the signal Doppler frequencies. The correlation time of the clutter is assumed as 0.03 seconds for an X-band ground clutter [25]. The Doppler frequency can be exactly computed as 9,192.5 Hz at the center of the projected ellipse. The signal Doppler frequencies on the projected antenna footprint distributes over the wrapped Doppler frequency range (−800 < fD < −700) ∪ (344.6 < fD < 800) [19]. Hence, the correct Doppler frequency generation is observed in Fig. 8. In Fig. 8(b), the null of the elevation signal can be observed. In the Tx mode, only the sum pattern is transmitted, and so the Doppler frequency distribution of the difference pattern is similar to that of the sum pattern. Over the projected antenna footprint on the ground, the Doppler frequencies exit over a 500-Hz band, which indicates that the 500-Hz frequency band is generated due to the platform dynamics. Hence, the correlation time of the return signal due to the platform movement is around 1/500≈0.002 seconds [22, 24]. Thus, two correaltion times of the ground clutter is chosen: a longer correlation time, 0.03 seconds, and a shorter one, 0.001 seconds. As explained earlier, the correlation time of the generated radar signal should be similar to the shorter one between the correlation times due to the platform and clutter dynamics. In Fig. 9, the correlation of the generated signal is computed along the slow time to produce two correlation times of the clutter. It shows that the correlation time becomes the shorter one between the two values of tc of the gound clutter and the platform.

Fig. 8

RD maps of monopulse signals for tc = 0.03 seconds: (a) sum signal (Ψ) and (b) elevation signal (ΨΔEL ).

Fig. 9

Correlation of received signals: (a) tc = 0.03 seconds and (b) tc = 0.001 seconds.

Next, an AGR simulation is considered for a real radar application. Here, we use the coherent integration technique to combine the 2D time-domain signal [26]. The coherent technique truncates the RD map of the sum signal around the Doppler center. Also, the RD map of the elevation difference signal is chosen from the same range as that for the sum pattern. Then, the monopulse ratio is calculated and averaged to one ratio [14, 26]. In Fig. 10, the ratio is shown with and without a γ-correction scheme. The γ-correction scheme can retrieve the correct monopulse ratio, which enhances the accuracy of the range estimation. The null point of the monopulse ratio is around the range estimation [26]. Without the γ-correction scheme, the range estimation error is about 2.06% but with the γ-correction scheme, the error is drastically reduced to about 0.07% which is within the known error range(< ≈ 1%) of the phase-comparison monopulse AGR scheme [14].

Fig. 10

Monopulse ratio for AGR simulation with and without γ-correction.

IV. Conclusion

A numerical algorithm was proposed to generate the ground return signal in a time domain for a monopulse radar system, which can accurately consider the radar operation scenario and the Doppler effects due to the movements of the platform and ground clutter. To truncate the return area, the main beam up to the HPBW is projected on the ground surface. Since the projected antenna footprint is estimated using a ray-tracing technique, the accuracy of the footprint can be very high for any radar operation scenario. To simulate the monopulse system, the sum signal is transmitted and then, using the appropriate array weights, the sum and two difference signals are synthesized. This method guarantees the correlation relationship between the sum and difference signals from an identical ground cell. The backscattering coefficient of the ground clutter is assumed to be a Gaussian random process. Thus, the Doppler effect of the ground clutter can be efficiently modeled using a correlated random number generation. For an application example, an AGR simulation is considered. The proposed signal-generation method is applied to a monopulse AGR simulation. It is numerically shown that with the known γ-correction scheme, the signal distortion caused by rolling of the platform can be effectively removed. The range estimation error is within the known error boundary for the phase-comparison monopulse AGR scheme.

Acknowledgments

This work was supported by a grant-in-aid of Hanwha Systems (No. U-18-015).

References

1. Skolnik MI. Radar Handbook New York, NY: McGraw-Hill; 1970.
2. Malas JA. F-22 radar development. In : Proceedings of the IEEE 1997 National Aerospace and Electronics Conference (NAECON). Dayton, OH; 1997; p. 831–839.
3. Willis NJ, Griffiths HD. Advances in Bistatic Radar Raleigh, NC: SciTech Publishing Inc; 2007.
4. Sherman SM, Barton DK. Monopulse Principles and Techniques 2nd edth ed. Boston, MA: Artech House; 2011.
5. Oh Y, Sarabandi K, Ulaby FT. An empirical model and an inversion technique for radar scattering from bare soil surfaces. IEEE Transactions on Geoscience and Remote Sensing 30(2):370–381. 1992;
6. Curlander JC, McDonough RN. Synthetic Aperture Radar: Systems and Signal Processing New York, NY: John Wiley & Sons; 1991.
7. Li XM, Feng DZ, Liu HW, Xing MD, Luo D. Two-dimensional pulse-to-pulse canceller of ground clutter in airborne radar. IET Radar, Sonar & Navigation 3(2):133–143. 2009;
8. Zhu G, Chen Y, Yin H. Analysis of typical ground clutter statistical characteristics. In : Proceedings of 2017 International Applied Computational Electromagnetics Society Symposium (ACES). Suzhou, China; 2017; p. 1–2.
9. Gurbuz AC. Determination of background distribution for ground-penetrating radar data. IEEE Geoscience and Remote Sensing Letters 9(4):544–548. 2011;
10. Zhang Y, Hu H, Luo J. Distributed Antenna Systems: Open Architecture for Future Wireless Communications Boca Raton, FL: Auerbach Publications; 2007.
11. Sorooshyari S, Daut DG. Generation of correlated Rayleigh fading envelopes for accurate performance analysis of diversity systems. In : Proceedings of the 14th IEEE Proceedings on Personal, Indoor and Mobile Radio Communications (PIMRC). Beijing, China; 2003; p. 1800–1804.
12. Szabo AP. Clutter simulation for airborne pulse-Doppler radar. In : Proceedings of the International Conference on Radar (IEEE Cat No 03EX695). Adelaide, Australia; 2003; p. 608–613.
13. Mahafza BR. Introduction to Radar Analysis 2nd edth ed. Boca Raton, FL: CRC Press; 2017.
14. Zwagerman PK. Air-to-ground ranging using electronic roll stabilization of monopulse data. In : Proceedings of the 1988 IEEE National Radar Conference. Ann Arbor, MI; 1988; p. 151–156.
15. Fung AK. Microwave Scattering and Emission Models and Their Applications Boston, MA: Artech House; 1994.
16. Laur H, Bally P, Meadows P, Sanchez J, Schaettler B, Lopinto E, Esteban D. Derivation of the backscattering coefficient σ0 in ESA ERS SAR PRI products. European Space Agency. Noordwijk, The Netherlands: ESA Document ES-TN-RE-PM-HL09 2002.
17. Anton H, Busby RC. Contemporary Linear Algebra Hoboken, NJ: John Wiley & Sons; 2003.
18. Hiroshige JJ, Kennedy TA. Error analysis of electronic roll stabilization for electronically scanned antennas. In : Proceedings of the 1991 IEEE National Radar Conference. Los Angeles, LA; 1991; p. 71–75.
19. Mahafza BR. Radar Systems Analysis and Design Using MATLAB 3rd edth ed. Boca Raton, FL: CRC Press; 2013.
20. Zheng Y, Tseng SM, Yu KB. Closed-form four-channel monopulse two-target resolution. IEEE Transactions on Aerospace and Electronic Systems 39(3):1083–1089. 2003;
21. Curry GR. Radar System Performance Modeling 2nd edth ed. Boston, MA: Artech House; 2004.
22. Tse D, Viswanath P. Fundamentals of Wireless Communication Cambridge, MA: Cambridge University Press; 2005.
23. Dietrich CR, Newsam GN. Fast and exact simulation of stationary Gaussian processes through circulant embedding of the covariance matrix. SIAM Journal on Scientific Computing 18(4):1088–1107. 1997;
24. Morris GV, Harkness L. Airborne Pulsed Doppler Radar 2nd edth ed. Boston, MA: Artech House; 1996.
25. Xu J, Peng Y, Wan Q, Wang X, Xia X. Doppler distributed clutter model of airborne radar and its parameters estimation. Science in China Series F: Information Sciences 47(5):577–586. 2004;
26. Nam JH, Lee HS, Koh IS, Song JH, Bae CS. Performance assessment of mono-pulse radar air-to-ground ranging algorithm for various clutter environments. The Journal of Korean Institute of Electromagnetic Engineering and Science 31(7):631–639. 2020;

Biography

Jeong-Hun Nam

received a B.S. degree in electronic engineering from Gangneung-Wonju National University, Gangneung, Korea, in 2017. He is currently working toward his M.S. degree at the department of electronic engineering, Inha University, Incheon, Korea. His research interests include radar system modeling and radar clutter signal modeling and analysis.

Jae-Won Rim

received B.S. and M.S. degrees in electronic engineering from Inha University, Incheon, Korea, in 2014 and 2016, respectively. He is currently working toward a Ph.D. at the department of electronic engineering, Inha University, Incheon, Korea. His research interests include modeling and simulation of radar systems, electronic warfare, and numerical methods in electromagnetic fields.

Hyunsoo Lee

received his M.S. and Ph.D. degrees in electronics engineering from Inha University, Incheon, Korea, in 2014 and 2019, respectively. He is currently working at the agency for defense development, Daejeon, Korea, as a research engineer. His research interests include electromagnetic numerical analysis and diffraction by discontinuous wedge structure.

Il-Suek Koh

received B.S. and M.S. degrees in electronics engineering from Yonsei University, Seoul, Korea, in 1992 and 1994, respectively, and a Ph.D. degree from the University of Michigan, Ann Arbor, MI, USA, in 2002. In 1994, he joined LG electronics Ltd., Seoul, as a research engineer. Currently, he works at Inha Univerity, Incheon, Korea, as a professor. His research interests include wireless communication channel modeling and numerical and analytical methods for electromagnetic fields.

Jong-Hwa Song

received B.S. and M.S. degrees in electronics engineering from Ajou University, Suwon, Korea, in 2005 and 2007, respectively. He received his Ph.D. from Konkuk University, Korea, in 2016. Currently, he works as a senior engineer on a Avionics Radar Team at Hanwha Systems. His main research interests are radar system design, radar signal processing, and navigation systems for airborne radar.

Article information Continued

(p, ŷp, p): Orientation vector of the platform in its local coordinates,
(g, ŷg, g): Orientation vector of the platform in the global coordinates,
ŝl and ŝg: LOS vectors of the antenna in the platform’s local coordinates and in the global coordinates, respectively,
θ and ϕ: Elevation and azimuth angles in the antenna’s local coordinates,
θr θy and θp: Roll, yaw, and pitch angles in the platform’s local coordinates.

Fig. 1

Dynamics scenario for air-to-ground radar.

Fig. 2

Local and global coordinate systems of the platform and antenna.

Fig. 3

Sign of antenna weight for phase-comparison monopulse sum and difference patterns: (a) sum (w), (b) elevation difference (wΔEL ), and (c) azimuth difference ( wΔAZ ).

Fig. 4

Simulation results according to the azimuth angle ϕ. (a) Projected antenna footprint and V and W axes, and (b) analytic and simulation results of the V and W axis widths.

Fig. 5

Simulation results according to the rolling angle θr. (a) Projected antenna footprint and the V and W axes, and (b) analytic and simulation results of γ.

Fig. 6

Signal transmission/reception of antenna to/from surface-projected main-lobe beam.

Fig. 7

Monopulse response vs angle-off boresight: (a) monopulse signal before γ-correction, (b) monopulse signal after γ-correction, and (c) monopulse ratio of Im(ΨΔEL ) to Re(Ψ).

Fig. 8

RD maps of monopulse signals for tc = 0.03 seconds: (a) sum signal (Ψ) and (b) elevation signal (ΨΔEL ).

Fig. 9

Correlation of received signals: (a) tc = 0.03 seconds and (b) tc = 0.001 seconds.

Fig. 10

Monopulse ratio for AGR simulation with and without γ-correction.

Table 1

Simulation parameters

Parameter Value
Polarization VV
Frequency X band
Number of pulses 50
Platform height (feet) 30,000
Platform speed (m/s) 180
PRF (Hz) 1,670
Rectangle size (m) 30 × 30
Number of antenna array
V axis 20
W axis 20
Antenna array spacing
V axis λ / 4
W axis λ / 4
Correlation function Gaussian [23]
Correlation length (cm) 8.4 [5]
rms height (cm) 0.4 [5]
Dielectric constant 15.4 – j2.15 [5]