Preprint

Article

Analysis, Design, and Experimental Validation of a High-Isolation, LowCross-Polarization Antenna Array Demonstrator for Software-Defined-Radar Applications

Altmetrics

Downloads

Views

Comments

A peer-reviewed article of this preprint also exists.

Nicholas Ricciardella,

Walter Fuscaldo

,Tito Mattei,Anna Maria Fiorello,Leopoldo Infante,

Alessandro Galli^*

Nicholas Ricciardella,

Walter Fuscaldo

,Tito Mattei,Anna Maria Fiorello,Leopoldo Infante,

Alessandro Galli^*

This version is not peer-reviewed

Submitted:

24 May 2024

Posted:

28 May 2024

You are already at the latest version

Alerts

Abstract

In a software-defined radar (SDR) system, most of the signal processing usually implemented in hardware is implemented by software, thus allowing for higher flexibility and modularity compared to conventional radar systems. However, the majority of SDR demonstrators and proofs-of-concept reported in the open literature so far have been based on simple antenna systems. As a result, the full potentialities of an SDR approach have not been completely exploited yet. In this work, we propose a flexilbe antenna module to be integrated into an active electronically scanning array (AESA) with controlled sidelobe level over a wide angular range, exhibiting polarization reconfigurability with low cross-polarization level and high isolation. For this purpose, analytical and numerically efficient techniques for the synthesis of the aperture distribution and the correct evaluation of the radiating features (e.g., beamwidth, pointing angle, sidelobe levels, etc.) are presented in order to grant a real-time control of the digital beamforming network. A sub-array module demonstrator is fabricated and measured to corroborate the concept.

Keywords:

Subject: Engineering - Electrical and Electronic Engineering

1. Introduction

In the last few years, the attention on software-defined radar (SDR) has grown exponentially. As a matter of fact, this approach offers a series of advantages that goes from digital beamforming to hardware flexibility [1]. The possibility to drive an active electronically scanning array (AESA) with an SDR platform would pave the way for a fully-digital radiofrequency (RF) front-end with unprecedented capabilities in terms of reconfigurability, re-use of hardware, and modularity. In the past, the growth in SDR platforms has been limited mostly by the low-speed interface and the low sampling frequency of the analog-to-digital converters (ADC) and digital-to-analog converters (DAC) [2,3]. However, the recent developments in the state of the art for DAC–ADC components allowed for developing fully digital L-band products, and nowadays fully digital X-band architectures represent the next step on the SDR framework [4,5].

Recently, SDR platforms have been configured and tested to drive passive bistatic radar (PBR) systems based on linear arrays to perform detection and tracking functions. However, these works [6,7,8,9] were mainly focused on the software implementation of the calibration network and of the Digital Signal Processing (DSP) architecture, whereas the software implementation of a fully digital beamformer was not analyzed in detail. Moreover, in order to fully exploit the potentialities of the SDR approach, the design flexibility has to be ensured at the hardware level. Antenna elements that allow for dynamic polarization and beam steering are thus required, especially when used in large planar AESA to reach the high-gain and high-coverage requirements of radar tracking systems.

To this aim, we discuss here the design of a

40 \times 40

planar AESA capable of scanning a very directive beam over a wide angular range with controlled sidelobe in C band. Efficient analytical and numerical methods for both the realization of a Tseng–Cheng beamforming network [10] and the evaluation of the main radiating properties are proposed to easily interface the AESA with an SDR platform. In order to respond to the needs of design flexibility and high performance, inspired by the recent works on multifunction phased array radar (MPAR) [11,12,13], we propose an attractive unit cell fed with two ports that can be independently controlled in amplitude and phase, and capable of providing low cross-polarization and high isolation. In order to corroborate the antenna performance, a

4 \times 6

array demonstrator is fabricated and measured, showing good agreement with the expected theoretical results.

The paper is organized as follows. In Section 2 we describe the targeted antenna performance of the AESA and present analytical methods for both the design of the digital beamformer and the evaluation of the AESA radiating features. In Section 3 the antenna element is described and designed to ensure polarization flexibility, low cross-polarization, and high isolation. In Section 4 the fabrication and experimental validation of a

4 \times 6

demonstrator is shown. Finally, conclusions are drawn in Section V.

2. Array Design

We aim at designing a planar AESA capable of producing a pencil beam (about

3^{\circ}

of half-power beamwidth (HPBW), at broadside) over an angular range of

60^{\circ}

on the elevation plane and a full coverage on the azimuthal plane with a sidelobe level (SLL) of

- 30

dB to be maintained for the whole scanning range and over the design frequency range, also exhibiting polarization reconfigurability with a low cross-polarization level and high isolation. C band is considered for the design, with a bandwidth of

10 %

. The main frequencies in the design band will be called

f_{min}

for the lower frequency,

f_{max}

for the higher frequency and

f_{0}

for the central frequency. In particular,

f_{min} = 0.95 f_{0}

and

f_{max} = 1.05 f_{0}

. The latter properties can be demanded to the element factor and will be discussed in Section III. In the next subsections, we focus on the design of the Array Factor (AF) with analytical techniques to meet the constraints of having a pencil beam with controlled SLL over the entire scanning range.

2.1. Analytical Techniques

2.1.1. Sidelobe Level

We first examine the abovementioned specifications to find the most suitable technique for the realization of an efficient digital beamformer. The requirement of an invariant SLL considerably restricts the available beamforming techniques. As a matter of fact, apart from optimization algorithms (see, e.g., the excellent reviews in [14,15] and refs. therein), only the Tseng–Cheng technique [10] guarantees a constant SLL at arbitrary azimuthal cuts and for any scan angle. Moreover, the Tseng–Cheng technique has the following advantages compared to other techniques: i) it is fully analytical, ii) it is computed once, without requiring to refresh the amplitude coefficients as the scan angle is changing, iii) being based on the Dolph–Chebyshev technique, it provides the minimum beamwidth for a given SLL [10]. All these aspects make the Tseng–Cheng technique particularly attractive for its implementation in an SDR-driven beamforming network.

In order to implement the Tseng–Cheng technique we assume an even number of elements along the two main directions, i.e.,

N_{x} = N_{y} = M = 2 N

and quadrantal symmetry of the amplitude coefficients, i.e.,

a_{m n} = a_{- m, n} = a_{m, - n} = a_{- m, - n}

(the extension to rectangular lattices has been proposed in [16]). Under this hypothesis the AF can conveniently be recast as [10,17]:

AF (θ, ϕ) = 4 \sum_{m = 1}^{N} \sum_{n = 1}^{N} a_{m n} cos [(2 m - 1) u] cos [(2 n - 1) v]

(1)

where

u = 0.5 k d_{x} (sin θ cos ϕ - sin θ_{0} cos ϕ_{0})

and

v = 0.5 k d_{y} (sin θ sin ϕ - sin θ_{0} sin ϕ_{0})

d_{x} (d_{y})

are the spacing along the

x (y)

-axis,

θ

and

ϕ

are the elevation and azimuthal angles of the spherical coordinate system,

k = 2 π / λ

is the free-space wavenumber,

λ

being the free-space wavelength. As shown in [10], an invariant SLL over the entire angular range can be achieved if the amplitude coefficients are set equal to

\begin{matrix} a_{m n} & = \frac{1}{N^{2}} \sum_{p = 1}^{N} \sum_{q = 1}^{N} T_{M - 1} [w_{0} cos \frac{π (p - \frac{1}{2})}{2 N} cos \frac{π (q - \frac{1}{2})}{2 N}] \\ cos (\frac{π (m - \frac{1}{2}) (p - \frac{1}{2})}{N}) cos (\frac{π (n - \frac{1}{2}) (q - \frac{1}{2})}{N}) \end{matrix}

(2)

where

T_{M - 1} [\cdot]

is the Chebyshev polynomial of the first kind and order

M - 1

w_{0}

is related to the SLL through the expression

w_{0} = cosh [1 / (M - 1) arccosh (PSR)]

, where PSR is the peak-to-sidelobe ratio, i.e.,

PSR = \sqrt{1 / SLL}

(assuming all quantities in linear scale). In order to avoid grating lobes the usual condition for untapered distributions modifies as follows:

\frac{d_{x (y)}}{λ} < \frac{1 - [arccos (1 / w_{0})] / π}{1 + sin θ_{\max}}

(3)

where

θ_{\max}

is the maximum scan angle. As is clear from the expression (2), its straight implementation would require four nested loops to evaluate all the coefficients (viz., two for the m and n elements indices, and two for the p and q summation indices), thus resulting time-consuming computations. To take advantage of the MATLAB potentialities in handling matrix calculations, it is therefore important to recast (2) to avoid as much as possible the use of nested loops. Indeed, the two summations in (2) can conveniently be expressed as a quadratic form:

a_{m n} = N^{- 2} {\underset{̲}{t}}_{p} (m) {\underset{̲}{\underset{̲}{T}}}_{p q} {\underset{̲}{t}}_{q} (n)

(4)

where

{\underset{̲}{t}}_{p} (m) = cos [π (m - 0.5) (\underset{̲}{p} - 0.5)] / N

{\underset{̲}{t}}_{q} (n) = cos [π (n - 0.5) (\underset{̲}{q} - 0.5)] / N

, and

{\underset{̲}{\underset{̲}{T}}}_{p q} = T_{2 N - 1} {w_{0} cos [π (p - 0.5) / 2 N] cos [π (q - 0.5) / 2 N]}

with

\underset{̲}{p} = 1, 2 \dots N

and

\underset{̲}{q} = 1, 2 \dots N

row and column vectors of N elements, respectively.

2.1.2. Beamwidth

As is known, the HPBW (

Θ_{x (y) 0}

) at broadside along the principal planes of a uniform planar array is well approximated by the formula:

Θ_{x (y) 0} = π - 2 arccos (\frac{0.886 λ}{2 N_{x (y)} d_{x (y)}}) ≃ \frac{0.886 λ}{N_{x (y)} d_{x (y)}}

(5)

where the last approximation is obtained in the limit of large arrays (we recall that

arccos (x) ≃ π / 2 - x

for

x \to 0

As the beam is scanned at a generic angle (

θ_{0}

ϕ_{0}

), the beamwidth over the two principal planes changes according to the following equations [18]

\begin{matrix} Θ_{h} = {(\sqrt{{cos}^{2} θ_{0} [Θ_{x 0}^{- 2} {cos}^{2} ϕ_{0} + Θ_{y 0}^{- 2} {sin}^{2} θ_{0}, ϕ_{0}]})}^{- 1} \end{matrix}

(6)

\begin{matrix} Ψ_{h} = {(\sqrt{Θ_{x 0}^{- 2} {sin}^{2} ϕ_{0} + Θ_{y 0}^{- 2} {cos}^{2} ϕ_{0}})}^{- 1} \end{matrix}

(7)

where

Θ_{h}

and

Ψ_{h}

represent the HPBW on the elevation plane, and its orthogonal plane, respectively. When the array has the same number of elements and the same spacing over the x and y axes (as in the case shown discussed here),

Θ_{x 0} = Θ_{y 0} = Θ_{0 h}

, thus

Ψ_{h} = Θ_{0 h}

is constant, whereas

Θ_{h} = Θ_{0 h} sec θ_{0}

, giving rise to an unavoidable pattern distortion in the scanning towards endfire (viz., for

θ_{0} = 60^{\circ}

sec θ_{0} = 2

, hence

Θ_{h}

is doubled). While the numerical evaluation of the beamwidth along the plane

ϕ = ϕ_{0}

is quite straightforward, a few words should be spent on the evaluation along the orthogonal plane. Indeed, the AF is represented as a function of

θ

and

ϕ

, hence the HPBW over the

ϕ_{0}

-plane is simply given by the distance between the

- 3

dB points along

θ

, while the HPBW over the orthogonal plane needs an explicit definition of the angle

ψ

. This definition can be recovered from simple geometrical considerations and leads to the following expression

ψ = arcsin [sin \tilde{θ} sin (\tilde{ϕ} - ϕ_{0})]

(8)

where

\tilde{θ}

and

\tilde{ϕ}

are the zeros of the equation

tan θ cot θ_{0} = sec (ϕ - ϕ_{0})

, that is but the equation of the plane orthogonal to the unit vector

\hat{θ}

in (

θ_{0}

ϕ_{0}

), i.e., the plane defined by

\hat{r} \cdot \hat{θ} |_{θ_{0}, ϕ_{0}} = 0

, where

\hat{r}

is the position unit vector. We should note that

\tilde{θ}

and

\tilde{ϕ}

should be sought by means of a root-finding algorithm. Once all the

\tilde{θ}, \tilde{ϕ}

pairs are found, the AF calculated at (

\tilde{θ}, \tilde{ϕ}

) must be rearranged according to the ordering given by

- π / 2 \leq ψ (\tilde{θ}, \tilde{ϕ}) \leq π / 2

. Then, the

- 3

dB points along

ψ

can finally be found.

However, when tapering techniques for SLL reduction are taken, the HPBWs are notably affected and accounted for by the spreading factor

s_{f}

. In particular, if one applies a Dolph–Chebyshev technique (as that of Tseng–Cheng) to achieve an SLL of

- 30

dB, the

s_{f}

is given by

s_{f} = 1 + 0.636 {\{\frac{2}{PSR} cosh [\sqrt{{(arccosh PSR)}^{2} - π^{2}}]\}}^{2} .

(9)

Therefore,

Θ_{x 0}

and

Θ_{y 0}

as they appear in (5) have to be multiplied by

s_{f}

as given by (9).

2.1.3. Bandwidth

As is well known from array theory [17,19], the array bandwidth is severely limited by the phase shifters unless one resorts to true time-delay devices. The latter usually require expensive and complex architectures that can only be justified when the array needs to guarantee an excellent performance over a considerably large bandwidth. Conversely, for relatively small fractional bandwiths as those considered here, phase shifters are usually preferred. As a result, the beam is exactly steered at a point in space

(θ_{0}, ϕ_{0})

, only at the center frequency

f_{0}

, i.e., when the u and v terms appearing in (1) are evaluated at

k_{0} = k (f_{0})

. As a result, for

f \neq f_{0}

, the beam is squinted by an angular deviation

Δ θ

from the theoretical pointing angle

θ_{0}

given by the formula

Δ θ = | arcsin (sin θ_{0} f_{0} / f) - θ_{0} | .

(10)

The main consequence of the beam squint is that the

- 3

dB fractional bandwidth

Δ f / f_{0}

is affected. In particular, for a linear array of length L it is determined by the expression [19]

\frac{Δ f}{f_{0}} = 0.886 s_{f} \frac{λ_{0}}{L sin θ_{0}} ≃ s_{f} Θ_{0 h} csc θ_{0}

(11)

λ_{0}

being the center wavelength. Again, the last expression is obtained in the limit of large arrays using the right-hand side of (5).

2.2. Numerical Results

The analytical techniques presented in the previous paragraphs are exploited here in a proof-of-concept to test their effectiveness and accuracy in the array synthesis and its performance evaluation.

The proof-of-concept is chosen to fit the constraints mentioned at the beginning of Section 2, namely

SLL < - 30

dB,

θ_{0, max} = 60^{\circ}

Θ_{0 h} = 3^{\circ}

over the

f_{min}

–

f_{max}

range. In order to fulfill the SLL requirement, we have to accept a spreading factor of about

1.3

, but only a slight modification to the condition for avoiding grating lobes (

w_{0} ≃ 1.0057

). We thus set

d_{x} = d_{y} = λ_{0} / 2

which mostly satisfies (3). With these numbers at hand, it is found from (5) that

N_{x} = N_{y} = 40

provides for

Θ_{0 h}

around

2 . 5^{\circ}

that, after taking into account the spreading factor (viz.,

s_{f} ≃ 1.3

), yields

3 . 3^{\circ}

. Under this condition the

- 3

dB fractional bandwidth at the maximum scan angle would be around

6.6 %

. Using (11), it is found that the targeted bandwidth performance would be reached by either limiting the maximum scan angle to

40^{\circ}

, or broadening the broadside beamwidth to about

4 . 5^{\circ}

. In the following we consider the original setting, thus accepting a reduction in the operating bandwidth.

Therefore, we consider a

40 \times 40

planar array with

d / λ = 0.5

spacing among the elements and apply the Tseng–Cheng technique to obtain an invariant SLL of

- 30

dB. The distribution of the amplitude coefficients is reported in Figure 1. As shown, the coefficients at the four corners are almost zero, thus the corresponding elements may not be fed, with minor impact on the radiating performance. In particular, in Figure 2(a)–(d) a comparison is shown between the normalized power radiation patterns produced with an ideal (i.e., keeping all the coefficients) Tseng–Cheng distribution (see Figure 2(a)), and those produced by setting to 0 the coefficients below

0.01

(see Figure 2(b)),

0.05

(see Figure 2(c)), and

0.1

(see Figure 2(d)), when the array is scanned at

(θ_{0} = 60^{\circ}, ϕ_{0} = 45^{\circ})

. For a threshold of

0.01

no appreciable difference is seen with the ideal case, but in the former case one would achieve a reduced filling factor of

83 %

. As the threshold raises, the filling factor reduces further (viz.,

62 %

and

43 %

for thresholds of

0.05

and

0.1

, respectively), but the SLL raises as well (viz., peaks of

- 26

dB and

- 21

dB for thresholds of

0.05

and

0.01

, respectively). As a result, further reductions of the number of elements should not be considered as this would require optimization techniques for sparse arrays, an aspect that goes beyond the scope of this work.

The effect of the elements reduction is even more appreciable looking at the 1-D normalized power patterns along the two principal planes (solid black and gray dashed lines in Figure 2(e)–(h)). It is noted that the patterns along the elevation plane are more affected than those along the orthogonal plane in all cases. Interestingly, the beamwidths are not affected by the elements reduction, and their numerical evaluation agree well (a percent error less than

10 %

) with the formulas provided in (5)–(7).

With respect to the bandwidth performance, the beam squint

Δ θ

and the gain loss are evaluated as a function of the frequency in Figure 3(a), (b). It is seen that the numerical evaluation of the beam squint (solid black line in Figure 3(a)) compares well with the theoretical prediction of (10) (dashed grey line in Figure 3(a)), and reveals a beam squint of about

3^{\circ}

within the

- 3

dB bandwidth shown in Figure 3(b). The beamwidths along the principal planes and as a function of frequency are also shown in Figure 4 where their numerical evaluation is compared against the analytical expressions provided in Eqs. (6)–(7) accounting for the spreading factor. As commented before, the agreement is good, showing a percent error lower than

10 %

which remains rather constant with respect to the frequency. Interestingly, it is noted that

Ψ_{h}

is considerably less affected by the frequency with respect to

Θ_{h}

3. Element Design

3.1. Antenna Structure

In the previous Section 2, the array has been designed to obtain a radiation pattern featuring a pencil beam with a HPBW of about

3^{\circ}

and a SLL lower than

- 30

dB to be scanned over an angular range of

60^{\circ}

on the elevation plane and

360^{\circ}

on the azimuthal plane. At this stage, it is necessary to design the antenna element to meet the requirements of low cross-polarization and isolation over the

f_{min}

–

f_{max}

band, and ensure polarization reconfigurability. In this regard, the low cross-polarization requirement calls for elements that are symmetric along the principal planes, such as square patch antennas [18]. The latter, however, rarely maintain such performance over fractional bandwidths as large as

10 %

, unless stacked configurations (typically with the top substrate thicker than the bottom one) are considered [20,21]. As per the polarization reconfigurability, two ports are considered and fed with aperture-coupled techniques to guarantee the required isolation [11,22]. These aspects have led to consider the antenna architecture proposed in [11] as a reference structure to be suitably modified and optimized to meet our requirements in the

f_{min}

–

f_{max}

band.

3.2. Full-Wave Results

The aforementioned antenna, that is shown in Figure 5 (a), consists of a dual-pol aperture-coupled stacked patch made of three dielectric substrates. On the bottom of the lower substrate (Rogers Kappa438) two printed microstrip lines are used to excite two slots on the top of the same substrate. The plane where the slots are located is also used as a ground plane. On the lower part of the upper substrate there is a parasite patch slightly smaller than the main patch printed on the top of the same layer. Both upper and middle dielectric layers are realized with a Taconic TLX-8 laminate.

The radiating element has been simulated in CST [23] in a 2-D periodic environment, in order to take into account the coupling effects that occur in a 2-D array. The S-parameters results obtained from the simulation of the optimized antenna are shown in Figure 5 (b). The isolation between ports is rather satisfactory, with values below

- 30

dB, as well as the reflection coefficient that is well below

- 10

dB in the frequency range of interest.

In Figure 6 the radiation patterns are shown on the principal planes, and as expected from the reflection coefficients a low cross polarization level and a symmetric co-polar pattern is obtained for both the H- and V-port excitations. To complete the picture, full-wave 3-D radiation patterns for co-polar and cross-polar components of both H and V port are shown in Figure 7. A detailed view of the resulting antenna is also shown in Figure 8, whereas in Table 1 a list of the various design parameters is reported.

4. Sub-Array Demonstrator

4.1. Full-Wave Results

The radiating element analyzed in the previous Section is here used on a

2 \times 8

array to realize a demonstrator. This analysis is particularly useful to understand if the unit-cell simulation results (that emulate an infinite array environment) are consistent with a finite array simulation setup.

In Figure 9, the F-parameters obtained from the

2 \times 8

array simulation are shown. F-parameters are used to calculate antenna coupling coefficients sources in the case of simultaneous excitation, when is not possible to apply the general S-Parameter definition anymore, and can be therefore considered as active scattering parameters. In particular,

F_{hh (vv)}

parameters refers to the reflection coefficient of a port V(H) when all the other ports H(V) are excited, while the

F_{hv (vh)}

parameters refers to the isolation between the ports H(V) when all the others V(H) are excited. In particular, in the

f_{min}

–

f_{max}

band

F_{hh}

and

F_{vv}

parameters are well below

- 10

dB, and have a similar behavior to the reflection coefficients obtained from the unit cell simulation for both H and V ports. The same holds for the

F_{hv}

and

F_{vh}

parameters, that are similar to the transmission coefficients values obtained from the unit-cell simulation. The good agreement between the F parameters of the array simulation with the S parameters of the unit-cell simulation demonstrate that the latter takes into account very well the mutual coupling effects.

Moreover, in Figure 10 and Figure 11 are shown the 1-D and 3-D far-field radiation patterns, respectively. It can be noticed that the radiation pattern of the array is squeezed along the array axis due to the array factor. On the transverse axis, the total pattern remains almost unaffected with respect to the element pattern provided by the unit-cell. In all cases, the cross-polar values remain mostly below

- 20

dB, as expected from the F parameters.

4.2. Experimental Results

Two prototypes of this

2 \times 8

array have been fabricated and have been used to implement a total

4 \times 8

array for the measurement setup. The prototypes are placed on a metallic plane properly shaped, and two plexiglass bars are used to secure the two prototypes on the metallic support. For that reason, two columns of radiating elements are unavailable, and the array under test is a

4 \times 6

. The measurement setup is shown in Figure 12 (a), with all the measurements performed in an anechoic chamber.

The measured

S_{i j}

parameters are shown in Figure 12 (b) for both V and H port, respectively signed as port 1 and port 2. The measurements show a frequency shift of the resonance, that is centered to

{\tilde{f}}_{0} = f_{0} - 0.09 f_{0}

instead of

f_{0}

. However, a fractional bandwidth of

10 %

is achieved between

{\tilde{f}}_{min} = f_{min}

0.09 f_{0}

and

{\tilde{f}}_{max} = f_{max}

0.09 f_{0}

, and the isolation between ports is lower than

- 30

dB. The frequency shift can be attributed to different aspects:

An important PCB curvature has been observed. This is due to the manufacturing process, and at the design frequency this issue can affect the antenna gain and the $S_{i j}$ parameters;
Possible variations of the dielectric constant respect to the one shown in the data-sheet of the materials;
Possible interactions between the metallic support and the radiating elements. These metallic parts are very near to the radiating board, and their effect could be not negligible.

With this setup, it is possible to evaluate the total pattern, from the active reflection coefficient. As is known, the total pattern can be obtained by pattern multiplication between the array factor and the active element pattern

G_{e, act}

, which is in turn related to the active reflection coefficient

Γ_{act}

for every scan angle

θ_{0}, ϕ_{0}

through the relation [24]

G_{e, act} (u_{0}, v_{0}) = \frac{4 π cos θ_{0} p_{x} p_{y}}{λ^{2}} (1 - | Γ_{act} (u_{0}, v_{0}) |^{2})

(12)

where

p_{x}

and

p_{y}

represent the periodicity of the array along the x- and y-axis, respectively, with

u_{0} = sin θ_{0} cos ϕ_{0}

and

v_{0} = sin θ_{0} sin ϕ_{0}

The results obtained with (12) from the measurement of the scattering parameters are shown in Figure 13, whereas in Figure 14 are shown the principal planes of the active element pattern for both polarizations, and the gain plot versus frequency at broadside for both polarizations. From these results it can be seen that the active element gain has a quite flat and regular shape in the

u v

plane and a peak value near to the theoretical expected value in the

{\tilde{f}}_{min}

–

{\tilde{f}}_{max}

band. The main reason for this behavior is due to the frequency shift reported in Figure 12 (b) and in Figure 14, where the peak gain and the minimum of the reflection coefficient are located around

{\tilde{f}}_{0}

5. Conclusions

A planar AESA capable of producing a pencil beam with invariant sidelobe level over a wide angular range and at any frequency in C band has been designed with fully analytical techniques. An efficient MATLAB implementation of the Tseng–Cheng technique is presented that looks particularly attractive for realizing a Software-Defined-Radar (SDR) interface for the real-time control of the beamforming network. The radiating performance of a case study is analyzed in detail showing the effects of a reduction of the filling factor. A smaller

4 \times 6

demonstrator has been fabricated and measured to test the main radiating performance of the proposed antenna architecture. The flexibility of the polarization control, the low level of cross-polarization, and the good isolation between ports over a

10 %

fractional bandwidth make this prototype attractive for future SDR applications.

Acknowledgments

This work has been developed under the contract

COLB - CTR - 2017 - 009 - A

between Leonardo S.p.A. and Sapienza University of Rome.

References

Fulton, C.; Yeary, M.; Thompson, D.; Lake, J.; Mitchell, A. Digital Phased Arrays: Challenges and Opportunities. Proc. IEEE 2016, 104, 487–503. [Google Scholar] [CrossRef]
Kwag, Y.K.; Jung, J.S.; Woo, I.S.; Park, M.S. Modern software defined radar (SDR) technology and its trends. J. Electromagn. Eng. Sci. 2014, 14, 321–328. [Google Scholar] [CrossRef]
Debatty, T. Software defined radar a state of the art. 2nd International Workshop on Cognitive Information Processing (CIP); , 2010; pp. 253–257.
Ciociola, A.; Infante, L.; Ricciardella, N.; Solimene, R.; Felaco, M.; Pellegrini, G. Digitally Synthesized Antenna Test Bench for Next Generation Phased Array Systems. IEEE Int. Conf. Phased Array Syst. Technol., 2022, pp. 1–6. [CrossRef]
Castillo-Rubio, C.F.; Pascual, J.M. Current Full Digital Phased-Array Radar developments for Naval applications. IEEE Int. Conf. Phased Array Syst. Technol., 2019, pp. 1–6. [CrossRef]
Capria, A.; Petri, D.; Moscardini, C.; Conti, M.; Forti, A.C.; Massini, R.; Cerretelli, M.; Ledda, S.; Tesi, V.; Dalle Mese, E.; others. Software-defined Multiband Array Passive Radar (SMARP) demonstrator: A test and evaluation perspective. OCEANS 2015; IEEE, , 2015; pp. 1–6.
Capria, A.; Giusti, E.; Moscardini, C.; Conti, M.; Petri, D.; Martorella, M.; Berizzi, F. Multifunction imaging passive radar for harbour protection and navigation safety. IEEE Aerosp. Electron. Syst. Mag. 2017, 32, 30–38. [Google Scholar] [CrossRef]
Jamil, K.; Alam, M.; Hadi, M.A.; Alhekail, Z.O. A multi-band multi-beam software-defined passive radar. Part I: System design. Radar 2012; , 2012; pp. 64–67.
Jamil, K.; Alam, M.; Hadi, M.A.; Alhekail, Z.O. A multi-band multi-beam software-defined passive radar. Part II: Signal processing. Radar 2012; , 2012; pp. 72–75.
Tseng, F.I.; Cheng, D.K. Optimum scannable planar arrays with an invariant sidelobe level. Proc. IEEE 1968, 56, 1771–1778. [Google Scholar] [CrossRef]
Saeidi-Manesh, H.; Karimkashi, S.; Zhang, G.; Doviak, R.J. High-isolation low cross-polarization phased-array antenna for MPAR application. Radio Sci. 2017, 52, 1544–1557. [Google Scholar] [CrossRef]
Saeidi-Manesh, H.; Zhang, G. Low cross-polarization, high-isolation microstrip patch antenna array for multi-mission applications. IEEE Access 2018, 7, 5026–5033. [Google Scholar] [CrossRef]
Saeidi-Manesh, H.; Zhang, G. High-isolation, low cross-polarization, dual-polarization, hybrid feed microstrip patch array antenna for MPAR application. IEEE Trans. Antennas Propag. 2018, 66, 2326–2332. [Google Scholar] [CrossRef]
Rocca, P.; Oliveri, G.; Mailloux, R.J.; Massa, A. Unconventional phased array architectures and design methodologies – A review. Proc. IEEE 2016, 104, 544–560. [Google Scholar] [CrossRef]
Cheng, D.K. Optimization techniques for antenna arrays. Proc. IEEE 1971, 59, 1664–1674. [Google Scholar] [CrossRef]
Kim, Y.; Elliott, R. Extensions of the Tseng–Cheng pattern synthesis technique. J. Electromagn. Waves Appl. 1988, 2, 255–268. [Google Scholar] [CrossRef]
Elliot, R.S. Antenna Theoryand Design; John Wiley & Sons, 2006.
Balanis, C.A. Antenna Theory–Analysis and Design; John Wiley & Sons, Inc, 2005.
Mailloux, R.J. Phased Array Antenna Handbook; Artech House, 2005.
Waterhouse, R. Microstrip Patch Antennas: A Designer’s Guide; Springer Science & Business Media, 2013.
Mishra, P.K.; Jahagirdar, D.R.; Kumar, G. A Review of Broadband Dual Linearly Polarized Microstrip Antenna Designs with High Isolation [Education Column]. IEEE Antennas Propag Mag 2014, 56, 238–251. [Google Scholar] [CrossRef]
Ghorbani, K.; Waterhouse, R. Dual polarized wide-band aperture stacked patch antennas. IEEE Trans. Antennas Propag. 2004, 52, 2171–2175. [Google Scholar] [CrossRef]
SIMULIA CST Studio Suite, Dassault Systemes. (2023).
Pozar, D. The active element pattern. IEEE Trans. Antennas Propag. 1994, 42, 1176–1178. [Google Scholar] [CrossRef]

Figure 1. Map of the amplitude coefficients of a

40 \times 40

array for the design specifications on the beam given in the text through the Tseng–Cheng technique.

Figure 1. Map of the amplitude coefficients of a

40 \times 40

array for the design specifications on the beam given in the text through the Tseng–Cheng technique.

Figure 2. Color maps of

{| AF |}^{2} / {| AF |}_{\max}^{2}

vs.

θ

and

ϕ

, for (a) an ideal Tseng–Cheng technique, (b)–(d) a Tseng–Cheng technique neglecting coefficients with amplitudes below

0.01

0.05

, and

0.1

, respectively, when

θ_{0} = 60^{\circ}, ϕ_{0} = 45^{\circ}

. The dotted-dashed lines represent the

- 30

dB iso-lines of the contour plot. Parameters of the Tseng–Cheng technique and of the array are reported in the text. (e)–(h) The radiated power density normalized to the absolute maximum (i.e.,

{| AF |}^{2} / {| AF |}_{\max}^{2}

) along the principal planes (in black solid lines and dashed gray lines, the plane

ϕ = ϕ_{0}

, and the orthogonal one, respectively), for the same cases analyzed in (a)–(d), respectively.

Figure 2. Color maps of

{| AF |}^{2} / {| AF |}_{\max}^{2}

vs.

θ

and

ϕ

, for (a) an ideal Tseng–Cheng technique, (b)–(d) a Tseng–Cheng technique neglecting coefficients with amplitudes below

0.01

0.05

, and

0.1

, respectively, when

θ_{0} = 60^{\circ}, ϕ_{0} = 45^{\circ}

. The dotted-dashed lines represent the

- 30

dB iso-lines of the contour plot. Parameters of the Tseng–Cheng technique and of the array are reported in the text. (e)–(h) The radiated power density normalized to the absolute maximum (i.e.,

{| AF |}^{2} / {| AF |}_{\max}^{2}

) along the principal planes (in black solid lines and dashed gray lines, the plane

ϕ = ϕ_{0}

, and the orthogonal one, respectively), for the same cases analyzed in (a)–(d), respectively.

Figure 3. (a) Beam squint

Δ θ

and (b) gain loss vs. f in the range

f_{min} = 0.95 f_{0}

f_{max} = 1.05 f_{0}

. In (a) numerical (theoretical) results are reported with solid black (dashed grey) lines.

Figure 3. (a) Beam squint

Δ θ

and (b) gain loss vs. f in the range

f_{min} = 0.95 f_{0}

f_{max} = 1.05 f_{0}

. In (a) numerical (theoretical) results are reported with solid black (dashed grey) lines.

Figure 4. The beamwidths along the principal planes

Θ_{h}

(in black) and

Ψ_{h}

(in grey) vs. f in the range

f_{min}

–

f_{m a x}

. Numerical (theoretical) results are reported with solid (dashed) lines.

Figure 4. The beamwidths along the principal planes

Θ_{h}

(in black) and

Ψ_{h}

(in grey) vs. f in the range

f_{min}

–

f_{m a x}

. Numerical (theoretical) results are reported with solid (dashed) lines.

Figure 5. (a) Proposed antenna, (b) S parameters of the optimized radiating element.

Figure 6. (a)–(d) Main cuts of the radiation pattern for the optimized radiating element. (a)–(b) refer to H-port excitation; (c)–(d) refer to V-port excitation. (a) and (c) show the

ϕ = 0^{\circ}

cuts, whereas (b) and (d) show the

ϕ = 90^{\circ}

cuts.

Figure 6. (a)–(d) Main cuts of the radiation pattern for the optimized radiating element. (a)–(b) refer to H-port excitation; (c)–(d) refer to V-port excitation. (a) and (c) show the

ϕ = 0^{\circ}

cuts, whereas (b) and (d) show the

ϕ = 90^{\circ}

cuts.

Figure 7. (a)-(d) Radiation patterns for (a)–(b) horizontal polarization, (c)–(d) vertical polarization at

f_{0}

. In (a)-(c) co-polar component, in (b)–(d) cross-polar component.

Figure 7. (a)-(d) Radiation patterns for (a)–(b) horizontal polarization, (c)–(d) vertical polarization at

f_{0}

. In (a)-(c) co-polar component, in (b)–(d) cross-polar component.

Figure 8. (a)–(c) Detailed antenna structure. (a) Front view and cross sections at different heights. (b) Back view from the ground plane and feeding structure. (c) Back view from SMA connectors.

Figure 9. 2x8 array F parameters of the optimized radiating element.

Figure 10. (a)-(d) Main cut of the radiation patterns for the 2x8 array. (a)-(b) show results while the H port is excited, (c)-(d) while the V port is excited. In (a)-(c) the

ϕ = 0^{\circ}

cut, in (b)-(d) the

ϕ = 90^{\circ}

cut is plotted.

Figure 10. (a)-(d) Main cut of the radiation patterns for the 2x8 array. (a)-(b) show results while the H port is excited, (c)-(d) while the V port is excited. In (a)-(c) the

ϕ = 0^{\circ}

cut, in (b)-(d) the

ϕ = 90^{\circ}

cut is plotted.

Figure 11. (a)-(d) Radiation patterns for (a)-(b) horizontal polarization, (c)-(d) vertical polarization at

f_{0}

. In (a)-(c) co-polar component, in (b)-(d) cross-polar component.

Figure 11. (a)-(d) Radiation patterns for (a)-(b) horizontal polarization, (c)-(d) vertical polarization at

f_{0}

. In (a)-(c) co-polar component, in (b)-(d) cross-polar component.

Figure 12. (a) Array demonstrator, (b) Measured S Parameters.

Figure 13. (a)-(c) Gain measurements on the

u v

plane for H polarization at (a)

{\tilde{f}}_{min}

, (b)

{\tilde{f}}_{0}

, (c)

{\tilde{f}}_{max}

. (d)–(f) As in (a)–(c) but for V polarization.

Figure 13. (a)-(c) Gain measurements on the

u v

plane for H polarization at (a)

{\tilde{f}}_{min}

, (b)

{\tilde{f}}_{0}

, (c)

{\tilde{f}}_{max}

. (d)–(f) As in (a)–(c) but for V polarization.

Figure 14. Gain obtained from the

Γ_{a c t}

measurements for the three main frequencies (

{\tilde{f}}_{min}

{\tilde{f}}_{0}

{\tilde{f}}_{max}

) on the plane (a)

ϕ = 0^{\circ}

and (b)

ϕ = 90^{\circ}

for H polarization. In (c) and (d), as in (a) and (b) but for V polarization. In (e) and (f) it is shown the gain in dB obtained from

Γ_{a c t}

measurements for

θ = 0^{\circ}

versus the expected gain, respectively for H and V polarization

Figure 14. Gain obtained from the

Γ_{a c t}

measurements for the three main frequencies (

{\tilde{f}}_{min}

{\tilde{f}}_{0}

{\tilde{f}}_{max}

) on the plane (a)

ϕ = 0^{\circ}

and (b)

ϕ = 90^{\circ}

for H polarization. In (c) and (d), as in (a) and (b) but for V polarization. In (e) and (f) it is shown the gain in dB obtained from

Γ_{a c t}

measurements for

θ = 0^{\circ}

versus the expected gain, respectively for H and V polarization

Table 1. Design parameters for the radiating element

Substrates/Patch/Connectors		Slots		Strips
Parameter	Value [mm]	Parameter	Value [mm]	Parameter	Value [mm]
p	26.98	$d_{x Hslot}$	12.42	$d_{x Hstrip}$	4.80
$h_{1}$	0.51	$d_{y Hslot}$	10.00	$d_{y Hstrip}$	13.05
$h_{2}$	0.79	$d_{x Vslot}$	6.02	$d_{x Vstrip}$	8.55
$h_{3}$	2.36	$d_{y Vslot}$	12.56	$d_{y Vstrip}$	4.80
$w_{p 1}$	12.9	$w_{Hslot}$	1.84	$w_{Hstrip}$	0.9
$h_{p 1}$	12.9	$h_{Hslot}$	1.84	$l_{Hstrip}$	8.87
$w_{p 2}$	14.5	$w_{Hslit}$	0.26	$w_{Tstrip}$	0.55
$h_{p 2}$	14.5	$l_{Hslit}$	3.29	$l_{Tstrip}$	6.00
$d_{x H}$	5.50	$w_{Vslot}$	1.85	$w_{Vstrip}$	0.90
$d_{y H}$	13.49	$h_{Vslot}$	1.85	$l_{Vstrip}$	11.37
$d_{x V}$	9.00	$w_{Vslit}$	0.15	—	—
$d_{y V}$	5.50	$l_{Vslit}$	2.25	—	—

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.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.