I. Introduction
Target detection by scanning a given search volume is a basic surveillance function of radar. The radar scans a search volume by steering beam to angular positions based on a pre-defined beam lattice (or grid) in pre-defined sequences. In mechanically scanning radar, target tracking is generally performed by associating detections in consecutive search frames. The advent of electronically scanning radar made it possible to allocate an independent beam to the target being tracked. Modern radars are composed of an active electronically scanning array (AESA) antenna that offers extreme beam agility [1] and a processor that allocates beam to perform multi-functions, such as search, acquisition, and track for multiple targets. In the case of airborne AESA radar, interleaving of air-to-air and air-to-surface operations is also possible [2]. Thus, to maximize the performance of each function, efficient resource management in AESA radars is critical [3].
Previous researches have focused on optimizing parameters for search task to maximize the radar performance with minimum radar resources [4–7]. In [4, 5], parameters such as search frame time, instrument range Rs (for a case of low pulse repetition frequency [PRF]), beam spacing, and duty factor were investigated to maximize track initiation range Rt with minimum power. With respect to the search frame time, ΔR⁄Rt against the required power to achieve a given Rt was analyzed (where ΔR is a radial distance moved by the target between successive looks) to find the condition that minimizes the required power. With respect to the instrument range, the combination of ΔR⁄Rt and Rs⁄Rt that also minimizes the required power was investigated. With respect to the beam spacing, the optimum beam spacing which minimizes the required power was analyzed by calculating the ratio between pulse integration gain per dwell and beamshape loss depending on beam spacing. In this analysis, a fixed search volume and search frame time were assumed. With respect to the duty factor, eclipse losses for medium and high PRF waveforms were considered (together with the increase in mean power as duty factor increases) to find the optimum value of the duty factor. In [6], similar parameter optimization to maintain the pre-determined Rt for ensuring minimum load for search task was presented. Specifically, a radar search load was defined as an average time required to conduct surveillance over a given search volume during a search frame time, and an optimization (minimization) process of the search load was proposed to satisfy a given Rt for the cumulative probability of 90% by tuning signal-to-noise ratio (SNR), beam overlap angle, and search frame time. In [7], optimizing beamwidth, search dwell time, and beam spacing for pre-defined single probability of detection with minimum search load was discussed. An analytical optimization by expressing the search load as a function of dwell time, beam width, beam spacing, search frame time, and SNR is proposed and compared with numerical optimization results.
While the previous works are interested in parameter optimization for search task, the work presented in this paper aims to provide a guideline for parameter selection when it is required to sacrifice radar resource for search task by reducing search frame time. In modern airborne AESA radars with multi-functionality, each function or task has its own priority, and it is desired to allocate sufficient radar resources to high-priority tasks (such as tracking high-threat targets, or air-to-surface operation in interleaved mode). Generally, search task has relatively low priority, and therefore, the time allocated to the search task is sacrificed first if other tasks with higher priority occur. In case when only search task exists, reducing the time allocated to the search task is equivalent to reducing search frame time. Also, modern airborne AESA radars provide operators with options for customizing radar functions. For instance, a radar operator can select or change (usually, reduce) search frame time either directly by selecting the desired search frame time or indirectly by changing the range-scale on multi-function display.
Reduction in the search frame time can be achieved either by increasing angular area for a single beam with broader beamwidth, by reducing the time allocated to a single beam with shorter dwell time, or by sparsely illuminating overall search volume with wider beam spacing. Regardless of which parameter is controlled, reducing search frame time from optimally designed search parameters results in degradation of detection range performance (which is defined as a cumulative probability of detection [8–12] in this work). In this paper, the detection range performance is analyzed for reduced search frame time by tuning either of the three parameters (beamwidth, dwell time, and beam spacing), and a guideline for selecting the parameter that minimizes detection performance degradation among the three parameters is presented.
II. Search Frame Time Variation
To scan a given search volume, radar illuminates search beam toward beam positions on pre-defined beam lattice in predefined sequence. For radar antenna beam pattern with one-way 3 dB beamwidth θbw, the beam spacing Δθ between the adjacent beam positions on the beam lattice is conventionally defined with a beam overlap ratio ρ as follows:
In the case of radar with an electronically scanned array antenna, θbw increases as beam steering angle with respect to the antenna boresight increases due to the reduced effective antenna aperture [13]. Considering such beam broadening of the antenna, the beam lattice in this work is designed to have uniform beam spacing in sine space, resulting in non-uniform beam spacing but uniform beam overlap ratio in antenna or radar coordinates. Also, instead of rectangular beam lattice, triangular beam lattice [14] is considered in this work.
For the beam lattice with uniform beam spacing in sine space, the number of beams Nb for a single search frame is approximately calculated as follows. Assuming a search volume centered at (az0, el0) with angular span of (azw, elw) in azimuth and elevation, the number of elevation bar Nel is given as
where Δv = ρel sin(θbw,el ), ρel and θbw,el are beam overlap ratio and one-way 3dB beamwidth in elevation, respectively. The upper and lower elevation limits of the search volume in sine-space are denoted by
v l i m + and
v l i m - , respectively, where are
v l i m ± = sin ( e l 0 ± 0.5 e l w ) . Eq. (2) assumes that the elevation offset of the first beam position from the boundary of the search volume is ΔelΔv⁄ρel. Using Nel, Nb is given as follows:
where Δu=ρaz sin(θbw,az), ρaz and θbw,az are beam overlap ratio and one-way 3dB beamwidth in azimuth, respectively. The left and right azimuth limits of the search volume in sine-space are denoted by
u l i m , i + and
u l i m , i - , respectively, where are
u l i m , i ± = sin ( a z 0 ± 0.5 a z w ) cos ( e l i ) . The elevation angle of i-th bar is denoted by eli which is given as follows:
Note that Eq. (3) assumes that the azimuth offset of the first beam position from the boundary of the search volume is ΔazΔu cos(eli )⁄ρaz. In Fig. 1, an example of a beam lattice is shown when (az0, el0) = (10°, −25°), (azw, elw) = (120°, 12°), (ρaz, ρel) = (1, 1), and (θbw,az, θbw,el) = (3°, 4°) for both uv coordinate (sine space) and radar coordinate, respectively. The search volume boundary is shown with the red box, the search volume center is marked with the red-cross symbol, the beam steering angle of each search beam is marked with the black dot symbol, and the beam coverage corresponding to one-way beamwidth is shown with the gray ellipse. Nb calculated from Eq. (3) is 146, while the true total number of search beam Nb,true is 147. To prove the validity of Eq. (3), the error ratio ɛb from Eq. (5) is shown in Fig. 2 for azw and elw up to 120° with a fixed search volume center (az0, el0) = (0°, 0°), respectively, when (θbw,az, θbw,el ) and (ρaz, ρel) are the same as those of the beam lattice in Fig. 1.
It is observed that ɛb is only less than 6.8% for whole values of (azw, elw).
Finally, with a dwell time of a single search beam td, the search frame time Tfr without interrupt of any other task is given as
1. Beamwidth
Assuming θbw,az = θbw,el = θbw for simplicity, it is observed that Tfr is inversely proportional to sin2(θbw) from Eqs. (3) and (6). Fig. 3(a) shows that, as θbw is increased from a reference beamwidth θbw,0 to
2 θ b w , 0 , Tfr is reduced from Tfr,0 to 0.5Tfr,0, where Tfr,0 is a reference search frame time when θbw = θbw,0. In AESA radar, the beamwidth can be increased by tapering the amplitude distribution of the array elements, or simply, by not using a number of elements. The cost of reducing search frame time by broadening beamwidth is loss of SNR. SNR is linearly proportional to the antenna gain G2, which is inversely proportional to
θ b w 2 [15]. The loss of SNR as search frame time is reduced by broadening beamwidth is also shown in Fig. 3(a). As Tfr is reduced from Tfr,0 to 0.5Tfr,0, the loss is increased by 6.02 dB.
Besides the SNR loss, increasing the beamwidth causes degradation of the angular measurement accuracy in the case of monopulse processing, since the accuracy is linearly proportional to the beamwidth [16], i.e., the wider the beamwidth, the poorer the accuracy.
2. Dwell Time
From Eq. (6), it is evident that Tfr is linearly proportional to td. Since SNR is also linearly proportional to td, for a reference dwell time td,0, reduction in search frame time by shortening dwell time also results in an increase in SNR loss, as shown in Fig. 3(b). As td is reduced from a reference dwell time td,0 to 0.5td,0, Tfr is reduced from Tfr,0 to 0.5Tfr,0, and the SNR loss is increased by 3.01 dB.
Dwell time can be shortened either by PRF or by reducing the number of pulses. Since PRF is selected by taking into account various factors, such as unambiguous range or velocity for low or high PRF, visibility and decodability in range-Doppler domain in medium PRF [17, 18], it is recommended not to change PRF for reducing search frame time. However, shortening the dwell time degrades velocity resolution of Doppler processing since the resolution is inversely proportional to the dwell time (more precisely, coherent processing interval) [19].
3. Beam Spacing
Broadening the beamwidth and shortening the dwell time to reduce search frame time both result in SNR loss. In addition, broadening the beamwidth is accompanied by the drawback of angular measurement accuracy degradation in monopulse processing, while shortening the dwell time results in poorer velocity resolution.
Instead of the previous two methods, the search frame time can be reduced by increasing beam spacing without changing the beamwidth or the dwell time. Assuming ρaz = ρel = ρ for simplicity, it is observed that Tfr is inversely proportional to ρ2 from Eqs. (3) and (6). Fig. 3(c) shows that, as ρ is increased from a reference beam overlap ratio ρ0 to
2 ρ 0 , Tfr is reduced from Tfr,0 to 0.5Tfr,0. Although SNR loss is negligible when a target is within the beamwidth, increasing the beam spacing increases the probability of target not being illuminated by the beam. The reduction in the probability can be quantified as SNR loss by calculating beamshape loss [20]. The beamshape loss calculation in [20] is complicated because various parameters, such as the probabilities of detection and false alarms, should be considered. In this work, instead of that, the average beamshape loss &Lmacr;bs is simply modeled as an average loss of two-way antenna gain for a given ρ as follows:
To provide an example of &Lmacr;bs, G is modeled with a far-field pattern of cosine-illuminated rectangular aperture [14, 20], as shown in Fig. 3(d).
where η = 1.1889. In Fig. 3(c), the SNR losses due to increase in &Lmacr;bs calculated from Eqs. (7) and (8) are shown when ρ0 = 0.5 (blue dashed line) and 1.0 (red dashed line), respectively. When Tfr = 0.5Tfr,0, the SNR losses are 0.97 dB and 4.07 dB when ρ0 = 0.5 and 1.0, respectively.
Note that, by applying multi-step interlaced scan (MIS) [21], the beamshape loss can be drastically reduced. The basic concept of MIS for one-dimensional (1D) scanning (i.e., Nel = 1), which was introduced in [21], is expanded to two-dimensional (2D) scanning (i.e., Nel > 1) in this work. The procedure for beam allocation for a given search volume using 2D-MIS is presented as a flow chart in Fig. 4. The steps for interlaced scanning in azimuth and elevation are denoted by Ns,az and Ns,el, respectively. The optimal values of the step Ns depending on ρ are provided in Table 1. The process of optimal step selection was provided in [21], and not repeated in this paper for brevity. The range of ρ is expanded up to 9.9 and beyond in order to take into account for up to 10% reduction in radar resource allocated to search task (i.e., 90% of the radar resource are allocated to other tasks, such as tracking targets or air-toground operation in interleaved mode) for 1D MIS (For 2D-MIS, the range of ρ up to
10 is sufficient to account for up to 10% reduction in radar resource allocated to search task). In the flow chart illustrated in Fig. 4, (u0, v0) and (u, v) refer to the initial and current beam positions of each elevation bar in u and v of sine-space, respectively. The procedure can be summarized as follows:
• Step#A: Search beam allocation begins from (u0, v0), and the beam position shifts by Δu until it reaches
u l i m , i + .
• Step#B. After shifting (u0, v0) to the initial beam position of the next elevation bar, Step#A is repeated.
• Step#C. Once Step#B is repeated for the last elevation bar (i.e., after one search frame is fully scanned), then the values of u0 and u of the previous search frame are shifted by Δu⁄Ns,az, and then Step#A and Step#B are repeated.
• Step#D: Step#A to Step#C are repeated until the accumulated extent of the shift becomes (Ns,az – 1)Δu/Ns,az (i.e., after MIS in azimuth is completed), then the values of v0 and v of the previous search frame are shifted by Δv⁄Ns,el, and then Step#A to Step#C are repeated.
• Step#E. Step#A to Step#D are repeated until the accumulated extent of the shift becomes (Ns,el – 1)Δu/Ns,el (i.e., after MIS in elevation is completed), then whole procedure is repeated from the beginning.
By applying MIS, the SNR losses are reduced to 0.12 dB and 0.86 dB when ρ0 = 0.5 (blue dash-dotted line) and 1.0 (red dash-dotted line), respectively, as shown in Fig. 3(c).
In Fig. 5, the SNR losses for the previously described search frame time reduction methods are compared. Compared to broadening beamwidth, shortening dwell time has an advantage in that the SNR loss is reduced by 50%. Thus, it is recommended to shorten dwell time to reduce search frame time, instead of broadening beamwidth, when sacrificed velocity resolution is not critical. The SNR loss for increasing beam spacing can be either higher or lower than that for shortening dwell time, depending on ρ0. Regardless of ρ0, however, the SNR loss can be drastically reduced by applying MIS when increasing beam spacing. Thus, in the view point of detection range performance which is directly related to SNR, it is recommended to increase beam spacing with MIS for search frame time reduction.
III. Detection Range Performance
In this section, the detection range performance is analyzed for reduced search frame time by tuning either of the three parameters described in the previous section: beamwidth, dwell time, and beam spacing. The detection range performance is defined as a cumulative probability of detection Pc [8–12, 21]. At the n-th search frame scan, Pc (n) is given as
where Pds(m) is a single probability of detection at the m-th scan. Assuming that there are k search beams per single search frame, Pds(m) is given as
where Pd(i, m) is a single probability of detection at the i-th search beam of the m-th scan. Note that, in this work, the detections of the search beams in the same search frame are assumed to be uncorrelated for the calculation of Pds. For a target of Swerling case 1, Pd is given at range R as
where Pfa is a probability of false alarm. In Eq. (11), SNR is calculated as follows:
where Pt is the transmit power of each element in the AESA, Narray is the number of elements, δ is duty factor, σ is the radar cross section (RCS) of the target, λ is wavelength corresponding to RF carrier frequency, k is Boltzmann’s constant, T0 is the standard temperature of 290 K, F is the noise figure, Lscan is scan loss, and Lsys is the overall system loss without Lbs and Lscan. Lbs is instantaneous beamshape loss (or antenna pattern loss [22]) which is calculated at each search beam as G2(θt ), where θt is the angular offset of the target from the beam center.
For detection range performance analysis, the following situation is assumed. The radar is scanning a search volume of (azw, elw) = (120°, 20°) and (az0, el0) = (0°, 0°) with a search beam of (θbw,az, θbw,el ) = (4°, 4°) and td = 50 ms on a beam lattice with (ρaz, ρel) = (1, 1). The number of search beams are 248, and therefore, Tfr = 12.4 seocnds. The rest of the parameters for the calculation of Pc are listed in Table 2.
From this case (Case#1), suppose that it is desired to reduce the search frame time to 50% of the original value, and 50% of the radar resource is allocated to other tasks. This can be accomplished either by broadening the beamwidth by
2 (Case#2), shortening the dwell time by 0.5 (Case#3), or increasing the beam overlap ratio by
2 (Case#4). Specifically, the search frame time of 12.4 seconds in Case#1 is reduced to 50% (thus, Tfr = 6.2 seconds) by broadening the beamwidth to (θbw,az, θbw,el) = (5.64°, 5.64°) in Case#2, by shortening the dwell time to td = 25 ms in Case#3, and by increasing beam overlap ratio to
( ρ a z , ρ e l ) = ( 2 , 2 ) in Case#4. In Case#4, it is possible to apply 2D-MIS with (Ns,az, Ns,el) = (4, 4) as an option. In this work, Case#4 without MIS is denoted by Case#4(a), and Case#4 with MIS is denoted by Case#4(b).
To analyze the detection range performance, it is assumed that the radar platform (for instance, a fighter) is heading north at a constant velocity Vr = 200 m/s, while the target is moving at velocity Vt = 200 m/s. For each of the four cases (Case#1 to #4), Pc is analyzed using Eq. (9) for the target azimuth angle azt spans from 0° to 60°. It is assumed that the target aspect angle aspt with respect to the target’s heading is set as
so that both azt and aspt are kept constant as the target approaches to the radar platform [17]. The target altitude is the same as that of the radar platform.
In Fig. 6, the range R90 corresponding to Pc = 90% is shown against azt. In all of the four cases, ripples are observed as azt changes. Such ripples are originated from variation of Lbs. With a fixed beam lattice, R90 is maximized when a target is at the center of each beam since Lbs = 0, and minimized when it is at the mid-point between the centers of adjacent beams since Lbs is maximized. Besides the ripples, the gradual decrease in R90 as azt increases is due to the increase in Lscan, which is a common characteristic of electronically scanned arrays.
As expected, R90 of Case#1 is superior to those of the other cases, since whole radar resources are allocated to search task in Case#1. Among the four cases, Case#2 shows the worst detection range performance due to the largest SNR loss, as observed in the previous section. R90 of Case#3 and Case#4(a) are longer than R90 of Case#2 for the entire range of azt. However, it is unclear which is superior between Case#3 and Case#4(a). When the target is near the mid-point between two adjacent search beams, R90 of Case#3 is longer than that of Case#4(a). While Lbs of Case#3 remains the same as that of Case#1, since there is no change of either beamwidth or beam spacing, Lbs is drastically increased in Case#4(a) when the target is near the mid-point between two adjacent search beams due to wider beam spacing. On the other hand, R90 of Case#4(a) is superior to that of Case#3 when the target is close to the center of search beams. When the target is at the exact center of search beams, the SNR of Case#3 is half of the SNR of Case#1, while the SNR of Case#4(a) is the same as that of Case#1. Consequently, detection range performance of Case#4(a) highly fluctuates depending on azt, as observed in Fig. 6. Such fluctuation is not desirable since the detection range performance varies drastically depending on the angular position of a target [18], and can be smoothed by applying MIS.
When MIS is applied, beam positions of beam lattice gradually shift as search frame is repeated, and therefore, Lbs at a fixed azt within Ns,azNs,el (= 16) search frames changes continuously. In Fig. 7(a) and 7(b), the values of Lbs at azt = 0.20° (where Lbs is minimized in Case#1) and 2.15° (where Lbs is maximized in Case#1) within 2Ns,azNs,el search frames are shown for Case#1, Case#4(a), and Case#4(b), respectively. It is observed that Lbs is constant in both Case#1 and Case#4(a) where the beam lattices are fixed. Specifically, Lbs = 0.89 dB at azt = 0.20° in both Case#1 and Case#4(a). At azt = 2.15°, Lbs = 2.26 dB and 5.39 dB in Case#1 and Case#4(a), respectively. Unlike Case#1 or Case#4(a), Lbs changes continuously with a period of Ns,azNs,el in Case#4(b). Note that the average values of Lbs in Case#4(b) are higher or lower than those in Case#4(a) at azt = 0.20° or 2.15°. As a result, the fluctuation observed in Case#4(a) is mitigated in Case#4(b). Note that R90 of Case#4(b) is superior to that of Case#3 for the entire range of azt, which shows that increasing beam spacing with MIS minimizes detection range performance degradation.
Although the results in Fig. 6 provides insights on why tuning beam spacing with MIS has advantages over controlling the other two parameters, the scenario of a target approaching to the radar platform at a constant azimuth and aspect angles are not a practical scenario. Thus, in the rest of this section, further analysis on detection range performance is presented for various realistic scenarios. Specifically, the following four scenarios with different target trajectories are considered for variety of the analysis, as shown in Fig. 8.
In the 1st scenario (Sc#1), the target starts from is an azimuth angle of +5° and moves southward (−x direction in Fig. 8(a)) with a constant heading of −15° referenced to north, as shown in Fig. 8(a). Thus, the target azimuth angle decreases as the target gets closer to the radar, as shown in Fig. 8(b). The target altitude is the same as that of the radar platform.
In the 2nd scenario (Sc#2), starting from an azimuth angle of +45°, the target moves with a constant heading of −80° referenced to north. The target altitude is 4,000 ft lower than that of the radar platform.
In the 3rd scenario (Sc#3), starting from an azimuth angle of +5°, the target moves with a constant heading of +15° referenced to north. The target altitude is 4,000 ft higher than that of the radar platform.
In the 4th scenario (Sc#4), the target moves southward while maneuvering with S-turn. The target altitude is the same as that of the radar platform.
The cumulative probabilities of detection from the four scenarios are shown in Fig. 9. The values of R90 from each scenario are listed in Table 3. For all of the scenarios, it is observed that the detection range performance degradation is most severe for Case#2 due to the highest SNR loss, as analyzed in the previous section (Note that Pc does not reaches 90% in Sc#2). The second-worst performance is observed for Case#3, and the degradation is minimized in Case#4. In Sc#1, #2, and #4, Case#4(b) outperforms Case#4(a), while Case#4(a) shows slightly better performance in Sc#3, thus highlighting the advantages of applying MIS.
IV. Conclusion
Among the three parameters (beamwidth, dwell time, and beam spacing), it has been shown that increasing beam spacing has the advantage over controlling the other two parameters for reducing search frame time in terms of detection range performance. However, increasing beam spacing may suffer large beamshape loss when a target is near the mid-point between two adjacent search beams. Applying MIS has an effect of averaging beamshape loss, and therefore increase in beamshape loss corresponding to increased beam spacing is mitigated. On the other hand, as shown with the SNR loss analysis in Section II, broadening beamwidth and shortening dwell time suffer larger SNR loss which degrades the detection range performance. With various scenarios, the detection range performance, which is defined as the cumulative probability of detection, has been analyzed for the four cases (Case#1, #2, #3, and #4), respectively, in Section III. The results verify that increasing beam spacing provides the superior detection range performance compared to controlling the other two parameters, and that applying MIS has an advantage of mitigating the large fluctuation of detection range performance depending on target angle.








