# Effective Antenna Modellings for NF-FF Transformations

Authors: Francesco D'Agostino; Flaminio Ferrara; Claudio Gennarelli; Rocco Guerriero; Massimo Migliozzi
Journal: International Journal of Antennas and Propagation
 (2011)

---

## Abstract

Two efficient probe-compensated near-field-far-field transformations with spherical scanning for antennas having two dimensions very different from the third one are here developed. They rely on the nonredundant sampling representations of the electromagnetic fields and on the optimal sampling interpolation expansions, and use effective antenna modellings. In particular, an antenna with a predominant dimension is no longer considered as enclosed in a sphere but in a cylinder ended in two half spheres, whereas a surface formed by two circular “bowls” with the same aperture diameter but different lateral bends is adopted to shape an antenna with two predominant dimensions. These modellings are able to fit very well a lot of antennas by properly setting their geometric parameters. It is so possible to remarkably lower the number of data to be acquired, thus significantly reducing the measurement time. Numerical tests assessing the accuracy and the robustness of the techniques are reported.


---

## Body

## 1. Introduction

As is well known, near-field-far-field (NF-FF) transformation techniques play a significant role in modern antenna measurements [1–4]. As a matter of fact, the pattern evaluation from NF measurements allows one to overcome those drawbacks which, for electrically large radiating systems, make it unpractical to measure the pattern in a conventional FF range. Two main techniques for the FF pattern reconstruction from NF measurements can be found in literature. The former relies on the wave-mode expansion method and uses canonical scanning surfaces (sphere, cylinder, plane) [1–4]. The latter is based on the equivalence theorem [5, 6] and can be applied to arbitrary scanning surfaces although at a remarkably increased computational cost.Among the NF-FF transformation techniques, the employing of the spherical scanning has attracted considerable attention [7–15], since it allows one to reconstruct, from a single set of NF measurements, the complete radiation pattern of the antenna under test (AUT). In the classical NF-FF transformation with spherical scanning [7–13], the choice of the highest spherical wave and the data spacings on the measurement sphere are determined according to a rule of thumb, related to the smallest sphere able to contain the AUT.In [14], such an NF-FF transformation has been properly modified by taking into account the spatial band-limitation properties of the electromagnetic (EM) fields [16]. In particular, the choice of the highest spherical wave has been rigorously fixed by the band-limitation properties and the number of data to be acquired on the rings decreases towards the poles. Moreover, the nonredundant sampling representations of the EM field [17] have been applied in the same paper [14] to remarkably reduce the number of required NF data when antennas with one or two predominant dimensions are considered. These results have been obtained by considering the AUT as enclosed in a prolate or oblate ellipsoid, respectively. Then, an optimal sampling interpolation (OSI) formula has been developed to recover the data required by the aforementioned NF-FF transformation from the knowledge of a nonredundant number of NF samples. However, since the case of an ideal probe has been considered in [14], the effects of the measurement probe have not been taken into account in the evaluation of the antenna FF pattern. The ideal probe assumption has been removed in [15], wherein an efficient probe-compensated NF-FF transformation technique with spherical scanning tailored for elongated or quasiplanar antennas has been developed, by properly extending the nonredundant sampling representation to the voltage measured by a nondirective probe. Then, the nonredundant voltage samples are interpolated via an efficient OSI algorithm for reconstructing the NF data required by the probe-compensated NF-FF transformation technique with spherical scanning [7], again properly modified to account for the spatial bandlimitation properties.It must be stressed that the reduction in the number of the NF data to be acquired reflects in a decrease of the measurement time, and this is a very important issue for the antenna measurement community, since such a time is very much greater than that needed to carry out the NF-FF transformation. A measurement time reduction can also be obtained by performing a fast electronic scanning via the modulated scattering technique [18], which employs arrays of scattering probes. However, antenna testing facilities based on such a technique are less flexible than those using mechanical scans. Accordingly, it appears more convenient the reduction of the needed NF data attainable by using the nonredundant sampling representations, which does not require to change the positioning systems, but only the software for controlling them.The aim of this paper is to develop even more effective NF-FF transformations with spherical scanning tailored for antennas having one or two predominant dimensions. To this end, very flexible antenna modellings will be adopted. In particular, an electrically long antenna will be considered as enclosed in a rounded cylinder, namely, a cylinder which ends in two half spheres, whereas a surface formed by two circular “bowls” with the same aperture diameter but different lateral bends (double bowl) will be adopted to shape a quasiplanar antenna. The so-obtained NF-FF transformation techniques result to be more effective from the data reduction viewpoint than those employing the prolate or oblate ellipsoidal modelling. In fact, these flexible modellings allow one to fit the shape of a lot of actual antennas better by properly setting their geometric parameters. Moreover, they remain quite general and contain the spherical modelling as a particular case.


## 2. Nonredundant Sampling Representation of the Voltage on a Sphere

Let us consider an AUT enclosed in a convex domain bounded by a rotational surfaceΣ and a nondirective probe scanning a spherical surface of radius d in the antenna NF region. The spherical coordinate system (r,ϑ,φ) will be adopted in the following to denote an observation point P in the NF region. As shown in [19], the voltage measured by such a kind of probe has the same effective spatial bandwidth of the AUT field, and, accordingly, the nonredundant sampling representations of EM fields [17] can be applied to it. Moreover, since the scanning sphere can be represented by meridians and parallels, in the following we deal with the voltage representation on an observation curve C described by an analytical parameterization r̲=r̲(ξ). According to [17], let us introduce the “reduced voltage”(1)Ṽ(ξ)=V(ξ)ejγ(ξ),
where γ(ξ) is a proper phase function to be determined. The error, occurring when Ṽ is approximated by a spatially bandlimited function, becomes negligible as the bandwidth exceeds a critical value Wξ [17]. In fact, it exhibits a steplike behaviour whose transition occurs at Wξ [16]. Accordingly, such an error can be effectively controlled [20] by choosing a bandwidth equal to χ′Wξ, where χ′ is an excess bandwidth factor, slightly greater than unity for electrically large antennas.WhenC is a meridian, by adopting(2)Wξ=βl′2π
(β being the wavenumber and ℓ′ the length of the curve C′ obtained as intersection between the meridian plane passing through the observation point P and Σ), we get (3)γ=β2[R1+R2+s1′-s2′],ξ=πl′[R1-R2+s1′+s2′],
where R1,2 are the distances from P to the tangency points P1,2 on C′, and s1,2′ are their arclength coordinates.When the observation curve is a parallel, due to the involved symmetry, the phase function is constant and it is convenient to use the azimuthal angleφ as parameter. The corresponding bandwidth is given [17] by(4)Wφ=β2max⁡z′⁡(R+-R-)=β2max⁡z′((z-z′)2+(ρ+ρ′(z′))2-(z-z′)2+(ρ-ρ′(z′))2),
wherein ρ=dsinϑ,  ρ′(z′) is the equation of Σ in cylindrical coordinates, and the maximum is attained [17] on that zone of the surface Σ lying on the same side of the observation circumference with respect to its maximum transverse circle.In the light of these results, the reduced voltage atP(ξ,φ) on the meridian fixed by φ can be evaluated via the OSI expansion(5)Ṽ(ξ,φ)=∑m=m0-q+1m0+qṼ(ξm,φ)ΩM(ξ-ξm)DM′′(ξ-ξm),
where m0=Int[ξ/Δξ] is the index of sample nearest (on the left) to P,2q is the number of the retained intermediate samples Ṽ(ξm,φ), that is, the reduced voltage values at the intersection points between the sampling parallels and the considered meridian, and(6)ξm=mΔξ=2πm(2M′′+1),M′′=Int(χM′)+1,M′=Int(χ′Wξ)+1,χ being a factor that controls the truncation error [17] and Int(x) denoting the integer part of x. Moreover,(7)DM′′(ξ)=sin[(2M′′+1)ξ/2](2M′′+1)sin(ξ/2),ΩM(ξ)=TM[2cos⁡2(ξ/2)/cos⁡2(ξ0/2)-1]TM[2/cos⁡2(ξ0/2)-1]
are the Dirichlet and Tschebyscheff Sampling functions, respectively, TM(·) being the Tschebyscheff polynomial of degree M=M′′-M′ and ξ0=qΔξ.The intermediate samplesṼ(ξm,φ) can be evaluated via the OSI formula(8)Ṽ(ξm,φ)=∑n=n0-p+1n0+pṼ(ξm,φn,m)ΩNm(φ-φn,m)×DNm′′(φ-φn,m),
where Ṽ(ξm,φn,m) are the reduced samples on the parallel fixed by ξm,2p is the retained samples number, n0=Int(φ/Δφm), and(9)φn,m=nΔφm=2πn2Nm′′+1,Nm′′=Int(χNm′)+1,Nm′=Int[χ*Wφ(ξm)]+1,Nm=Nm′′-Nm′,χ*=1+(χ′-1)[sinϑ(ξm)]-2/3.The variation ofχ* with ξ is required to ensure a band-limitation error constant with respect to ξ.By properly matching (5) and (8), the following two-dimensional OSI expansion results: (10)Ṽ(ξ,φ)=∑m=m0-q+1m0+q{∑n=n0-p+1n0+pṼ(ξm,φn,m)ΩNm(φ-φn,m)DNm"(φ-φn,m)ΩM(ξ-ξm)DM′′(ξ-ξm)∑n=n0-p+1n0+p×Ṽ(ξm,φn,m)ΩNm(φ-φn,m)DNm′′(φ-φn,m)ΩM(ξ-ξm)DM′′(ξ-ξm)⋅∑n=n0-p+1n0+p},
which makes it possible to recover the NF data needed to carry out the spherical NF-FF transformation [7] as modified in [15].


### 2.1. Double Bowl Modelling

An effective and flexible modelling for an AUT having a quasiplanar geometry is obtained by choosingΣ coincident with the smallest surface formed by two circular “bowls” with the same aperture diameter 2a (see Figures 1 and 2). Note that their lateral surfaces have not the same bend because they are generally determined by rotating two different circular arcs, each equal to a quarter of circumference (with radius c and c′). This surface allows one to fit very well a lot of actual antennas by properly choosing the values of the parameters c, c′, and a. For instance, Σ coincides with a sphere if c=c′=a, and it becomes a half sphere if c=0 and c′=a, and it reduces to a circular dish for c=c′=0. It can be easily verified that, for such a modelling, the length of the curve C′ is ℓ′=2[b+b′+(c+c′)π/2], wherein b=a-c and b′=a-c′. The phase function and the optimal parameter relevant to a meridian can be determined by substituting in relations (3) the appropriate values of R1,2 and s1,2′, whose expressions change depending on the location of the tangency points P1,2. For such a modelling, five cases must be considered for ϑ ranging in [0,π] (see Figure 2).Figure 1
Spherical scanning for a quasiplanar antenna.Figure 2
Relevant to the double bowl modelling.When0≤ϑ≤ϑA=sin-1(a/d), it results that (11)R1=d2+b2+2bdsinϑ-c2,s1′=-(b+cα1),α1=tan-1(R1c)-tan-1[(b+dsinϑ)dcos⁡ϑ],R2=d2+b2-2bdsinϑ-c2,s2′=b+cα2,α2=tan-1(R2c)-tan-1[(b-dsinϑ)dcos⁡ϑ].WhenϑA<ϑ≤ϑB=cos-1(c/d), R1, s1′ and α1 are again given by relations (11), whereas it results that(12)R2=d2+b′2-2b′dsinϑ-c′2,s2′=b+c(π2)+c′α2,α2=tan-1(R2c′)-tan-1[dcos⁡ϑ(dsinϑ-b′)].WhenϑB<ϑ≤ϑc=π-cos-1(c′/d), R2, s2′, and α2 are again given by relations (12), whereas it results that(13)R1=d2+b2-2bdsinϑ-c2,s1′=b+c(α1+π2),α1=-tan-1(R1c)-tan-1[dcos⁡ϑ(dsinϑ-b)].WhenϑC<ϑ≤ϑD=π-sin-1(a/d), R1, s1′, and α1 are again given by relations (13), whereas it results that(14)R2=d2+b′2+2b′dsinϑ-c′2,s2′=b+2b′+(c+c′)(π2)+c′α2,α2=tan-1(R2c′)-tan-1[(dsinϑ+b′)|dcos⁡ϑ|].WhenϑD<ϑ≤π, R2, s2′, and α2 are again given by relations (14), whereas it results that(15)R1=d2+b′2-2b′dsinϑ-c′2,s1′=b+c(π2)+c′(π2-α1),α1=tan-1(R1c′)-tan-1[(b′-dsinϑ)|dcos⁡ϑ|].For what concerns the evaluation of the maximum in (4) which allows the determination of the azimuthal bandwidth Wφ, it is convenient (when ϑ≤π/2) to introduce the angular coordinate η (Figure 2) such that z′=ccosη and ρ′=b+csinη. It can be easily verified through simple geometrical considerations or analytical manipulations that the sought maximum corresponds to the η value that zeroes the derivative of R+-R- with respect to η. A quite analogous result holds when π/2<ϑ≤π.


### 2.2. Rounded Cylinder Modelling

An effective modelling for an elongated AUT is got by considering it as enclosed in a cylinder of heighth′ ended in two half spheres of radius a′ (see Figures 3 and 4). In such a case, it results that ℓ′=2(h′+πa′). The phase function and the optimal parameter relevant to a meridian can be again obtained by substituting the appropriate values of R1,2 and s1,2′ in relations (3). It can be easily verified that their expressions change depending on the location of the tangency points P1,2. In particular, three cases must be considered for ϑ ranging in [0,π] (see Figure 4).Figure 3
Spherical scanning for an elongated antenna.Figure 4
Relevant to the rounded cylinder modelling.When0≤ϑ≤sin-1(a′/d), it results that(16)R1=(dsinϑ)2+(dcos⁡ϑ-h′2)2-a′2,s1′=a′sin-1(a′dsinϑ+R1(h′/2-dcos⁡ϑ)R12+a′2),R2=R1,s2′=a′sin-1(a′dsinϑ-R2(h′/2-dcos⁡ϑ)R22+a′2).Whensin-1(a′/d)<ϑ≤π-sin-1(a′/d), the expressions of R1 and s1′ are again given by relations (16), whereas it results that(17)R2=(dsinϑ)2+(dcos⁡ϑ+h′2)2-a′2,s2′=h′+a′[π-sin-1(a′dsinϑ+R2(h′/2+dcos⁡ϑ)R22+a′2)].Whenπ-sin-1(a′/d)<ϑ≤π, the expressions of R2 and s2′ are again given by relations (17), whereas it results that(18)R1=(dsinϑ)2+(dcos⁡ϑ+h′2)2-a′2,s1′=h′+a′[π2-sin-1(R1dsinϑ+a′(h′/2+dcos⁡ϑ)R12+a′2)].As regards the azimuthal bandwidthWφ, it can be shown that the maximum in (4) is attained at(19)z′={z,|z|≤h′2,[h′2+(|z|-h′/2)a′2(dsinϑ)2+(|z|-h′/2)2]sgn⁡(z),|z|>h′2,
where sgn(·) is the sign function.


## 2.1. Double Bowl Modelling

An effective and flexible modelling for an AUT having a quasiplanar geometry is obtained by choosingΣ coincident with the smallest surface formed by two circular “bowls” with the same aperture diameter 2a (see Figures 1 and 2). Note that their lateral surfaces have not the same bend because they are generally determined by rotating two different circular arcs, each equal to a quarter of circumference (with radius c and c′). This surface allows one to fit very well a lot of actual antennas by properly choosing the values of the parameters c, c′, and a. For instance, Σ coincides with a sphere if c=c′=a, and it becomes a half sphere if c=0 and c′=a, and it reduces to a circular dish for c=c′=0. It can be easily verified that, for such a modelling, the length of the curve C′ is ℓ′=2[b+b′+(c+c′)π/2], wherein b=a-c and b′=a-c′. The phase function and the optimal parameter relevant to a meridian can be determined by substituting in relations (3) the appropriate values of R1,2 and s1,2′, whose expressions change depending on the location of the tangency points P1,2. For such a modelling, five cases must be considered for ϑ ranging in [0,π] (see Figure 2).Figure 1
Spherical scanning for a quasiplanar antenna.Figure 2
Relevant to the double bowl modelling.When0≤ϑ≤ϑA=sin-1(a/d), it results that (11)R1=d2+b2+2bdsinϑ-c2,s1′=-(b+cα1),α1=tan-1(R1c)-tan-1[(b+dsinϑ)dcos⁡ϑ],R2=d2+b2-2bdsinϑ-c2,s2′=b+cα2,α2=tan-1(R2c)-tan-1[(b-dsinϑ)dcos⁡ϑ].WhenϑA<ϑ≤ϑB=cos-1(c/d), R1, s1′ and α1 are again given by relations (11), whereas it results that(12)R2=d2+b′2-2b′dsinϑ-c′2,s2′=b+c(π2)+c′α2,α2=tan-1(R2c′)-tan-1[dcos⁡ϑ(dsinϑ-b′)].WhenϑB<ϑ≤ϑc=π-cos-1(c′/d), R2, s2′, and α2 are again given by relations (12), whereas it results that(13)R1=d2+b2-2bdsinϑ-c2,s1′=b+c(α1+π2),α1=-tan-1(R1c)-tan-1[dcos⁡ϑ(dsinϑ-b)].WhenϑC<ϑ≤ϑD=π-sin-1(a/d), R1, s1′, and α1 are again given by relations (13), whereas it results that(14)R2=d2+b′2+2b′dsinϑ-c′2,s2′=b+2b′+(c+c′)(π2)+c′α2,α2=tan-1(R2c′)-tan-1[(dsinϑ+b′)|dcos⁡ϑ|].WhenϑD<ϑ≤π, R2, s2′, and α2 are again given by relations (14), whereas it results that(15)R1=d2+b′2-2b′dsinϑ-c′2,s1′=b+c(π2)+c′(π2-α1),α1=tan-1(R1c′)-tan-1[(b′-dsinϑ)|dcos⁡ϑ|].For what concerns the evaluation of the maximum in (4) which allows the determination of the azimuthal bandwidth Wφ, it is convenient (when ϑ≤π/2) to introduce the angular coordinate η (Figure 2) such that z′=ccosη and ρ′=b+csinη. It can be easily verified through simple geometrical considerations or analytical manipulations that the sought maximum corresponds to the η value that zeroes the derivative of R+-R- with respect to η. A quite analogous result holds when π/2<ϑ≤π.


## 2.2. Rounded Cylinder Modelling

An effective modelling for an elongated AUT is got by considering it as enclosed in a cylinder of heighth′ ended in two half spheres of radius a′ (see Figures 3 and 4). In such a case, it results that ℓ′=2(h′+πa′). The phase function and the optimal parameter relevant to a meridian can be again obtained by substituting the appropriate values of R1,2 and s1,2′ in relations (3). It can be easily verified that their expressions change depending on the location of the tangency points P1,2. In particular, three cases must be considered for ϑ ranging in [0,π] (see Figure 4).Figure 3
Spherical scanning for an elongated antenna.Figure 4
Relevant to the rounded cylinder modelling.When0≤ϑ≤sin-1(a′/d), it results that(16)R1=(dsinϑ)2+(dcos⁡ϑ-h′2)2-a′2,s1′=a′sin-1(a′dsinϑ+R1(h′/2-dcos⁡ϑ)R12+a′2),R2=R1,s2′=a′sin-1(a′dsinϑ-R2(h′/2-dcos⁡ϑ)R22+a′2).Whensin-1(a′/d)<ϑ≤π-sin-1(a′/d), the expressions of R1 and s1′ are again given by relations (16), whereas it results that(17)R2=(dsinϑ)2+(dcos⁡ϑ+h′2)2-a′2,s2′=h′+a′[π-sin-1(a′dsinϑ+R2(h′/2+dcos⁡ϑ)R22+a′2)].Whenπ-sin-1(a′/d)<ϑ≤π, the expressions of R2 and s2′ are again given by relations (17), whereas it results that(18)R1=(dsinϑ)2+(dcos⁡ϑ+h′2)2-a′2,s1′=h′+a′[π2-sin-1(R1dsinϑ+a′(h′/2+dcos⁡ϑ)R12+a′2)].As regards the azimuthal bandwidthWφ, it can be shown that the maximum in (4) is attained at(19)z′={z,|z|≤h′2,[h′2+(|z|-h′/2)a′2(dsinϑ)2+(|z|-h′/2)2]sgn⁡(z),|z|>h′2,
where sgn(·) is the sign function.


## 3. NF-FF Transformation with  Spherical Scanning

For reader’s convenience, the key steps of the probe-compensated NF-FF transformation with spherical scanning [7] as modified in [15] by properly taking into account the spatial band-limitation properties of the EM field [16] are reported in this section.According to [7], the electric field in the FF region can be expressed in terms of the following truncated spherical wave expansion:(20)E̲(R,Θ,Φ)=e-jβRR∑n=1Nmax⁡×∑m=-nn[jn+1bnm1f̲̃1nm(Θ)+jnbnm2f̲̃2nm(Θ)]×ejmΦ,
where the spherical coordinate system (R,Θ,Φ) has been adopted to specify the FF observation point.The choice of the highest spherical wave to be considered is determined in the classical approach according to the following rule of thumb:(21)Nmax⁡=Int(βa)+10,
where a is the radius of the smallest sphere enclosing the AUT. Whereas in the approach [15], the highest spherical harmonic is rigorously fixed by the aforementioned band-limitation properties of the radiated EM fields [16]. Accordingly,(22)Nmax⁡=Int(χ′βa)+1.The vectorial functionsf̲̃1,2nm(Θ) are given by(23)f̲̃1nm(Θ)=(-m|m|)m12πn(n+1)×[jmsinΘP¯n|m|(cos⁡Θ)Θ̂-ddΘP¯n|m|(cos⁡Θ)Φ̂],f̲̃2nm(Θ)=(-m|m|)m12πn(n+1)×[ddΘP¯n|m|(cos⁡Θ)Θ̂+jmsinΘP¯n|m|(cos⁡Θ)Φ̂],P¯n|m|(cosΘ) being the normalized associated Legendre functions as defined in [21].As shown in [7], when using a probe with a first-order azimuthal dependence pattern (e.g., an open ended circular waveguide), the modal expansion coefficients in (20) can be determined from the knowledge of the voltages V and V′ measured by the probe and rotated probe, respectively, and are given by(24)bnm1=2n+116π∫0π∫02π[I1(ϑ)V-jI2(ϑ)V′]e-jmφsinϑdφdϑ∑υ=1υmax⁡[b′υ11Aυ1n(βd)-b′υ12Bυ1n(βd)],bnm2=2n+116π∫0π∫02π[I2(ϑ)V-jI1(ϑ)V′]e-jmφsinϑdφdϑ∑υ=1υmax⁡[b′υ11Bυ1n(βd)-b′υ12Aυ1n(βd)],
where υmax is the highest spherical wave of the probe,(25)I1(ϑ)=[d1mn(ϑ)-d-1mn(ϑ)],I2(ϑ)=[d1mn(ϑ)+d-1mn(ϑ)],dμmn(ϑ)=(n+μ)!(n-μ)!(n+m)!(n-m)!∑σ(n+mn-μ-σ)×(n-mσ)(-1)n-μ-σ×(cos⁡ϑ2)2σ+μ+m(sinϑ2)2n-2σ-μ-m
are the rotation coefficients [7, 22], the summation over σ involving all terms in which the binomial coefficients do not lead to negative arguments for the factorials. Moreover,(26)Aυμn(βd)=(2n+1)(2υ+1)n(n+1)υ(υ+1)(υ+μ)!(n-μ)!(υ-μ)!(n+μ)!×(-1)μ12jn-υ∑p=|n-υ|n+υ×{τ(μ,n,-μ,υ,p)hp(2)(βd)j-p[n(n+1)+υ(υ+1)-p(p+1)]×τ(μ,n,-μ,υ,p)hp(2)(βd)},Bυμn(βd)=(2n+1)(2υ+1)n(n+1)υ(υ+1)(υ+μ)!(n-μ)!(υ-μ)!(n+μ)!×(-1)μ12jn-υ∑p=|n-υ|n+υ×[j-p(2jμβd)τ(μ,n,-μ,υ,p)hp(2)(βd)]
are the translation coefficients [7, 11, 23], wherein hp(2)(x) is the spherical Hankel function of second kind, and τ(μ,n,-μ,υ,p) are the linearization coefficients defined by the expansion of the product of two unnormalized associated Legendre functions(27)Pnμ(x)Pυ-μ(x)=∑p=|n-υ|n+υτ(μ,n,-μ,υ,p)Pp(x).The probe’s expansion coefficientsb′υ11,2 in relations (24) can be evaluated from the knowledge of its tangential electric field on a scanning sphere when it is used as transmitting antenna(28)b′υ11,2=1βg1,2υ(βd)∫0π∫-ππE̲t′(d,ϑ,φ)⋅f̲̃1,2υ1*(ϑ)×e-jφsinϑdφdϑ,(·) and (*) denoting the inner product and the complex conjugate, respectively, and(29)g1,υ(x)=hυ(2)(x),g2,υ(x)=1xddx[xhυ(2)(x)].According to (28), the probe’s expansion coefficients can be determined from the simulated measurements of the field radiated by an open-ended circular waveguide by using the uncompensated NF-FF transformation [14].The integration overφ in relations (24) can be carried out [14] by expanding the voltage measured by the probe and rotated probe on each parallel in Fourier series with respect to φ, and this can be efficiently performed by employing the Fast Fourier Transform (FFT) algorithm. In a quite similar way, the remaining integration over ϑ can be efficiently performed through the following steps: (i) the aforementioned Fourier series coefficients are extended [7] from the range [0,π] to [0,2π] and then expanded in Fourier series via FFT, (ii) the rotation coefficients are expanded in Fourier series too [7, 22]. The double integrals are so transformed in double summations, of indexes ℓ and i, involving the Fourier series coefficients (in ϑ) and the integrals(30)∫0πej(l-i)ϑsinϑdϑ.It must be stressed that, to take advantage of the numerical efficiency of the FFT, the number of NF parallels to be considered in the NF-FF transformation and the number of samples on them must be the first power of two greater or equal toNmax and 2Mφ, respectively, where Mφ=Int(χ*βasinϑ)+1.It is worth noting that, by inverting the summations order, the spherical wave expansion (20) can be rewritten as(31)E̲(R,Θ,Φ)=e-jβRR∑m=-MΦMΦ∑n=|m|(n≠0)Nmax⁡×[jn+1bnm1f̲̃1nm(Θ)+jnbnm2f̲̃2nm(Θ)]ejmΦ=e-jβRR∑m=-MΦMΦc̲m(Θ)ejmΦ,
where MΦ=Int(χ*βasinΘ)+1. Relation (31) allows an efficient evaluation of the antenna far field at the considered elevation angle Θ by carrying out the summation again via the FFT algorithm.From a numerical viewpoint, it is convenient to apply the described FF reconstruction process to evaluate only the FF samples required by the OSI expansion in [24], properly modified to deal with an even number of samples on the FF parallels:
(32)F̲(Θ,Φ)=∑m=m0-q+1m0+q{∑n=n0-p+1n0+pF̲(Θm,Φn,m)ΩLΦ(Φ-Φn,m)ΩN(Θ-Θm)DN′′(Θ-Θm)×2LΦ′′-12LΦ′′∑n=n0-p+1n0+pF̲(Θm,Φn,m)×ΩLΦ(Φ-Φn,m)DLΦ′′-1(Φ-Φn,m)∑n=n0-p+1n0+pF̲(Θm,Φn,m)ΩLΦ(Φ-Φn,m)⋅2LΦ′′-12LΦ′′}
wherein (33)Θm=mΔΘ=2πm(2N′′+1),N′′=Int(χNmax⁡)+1,N=N′′-Nmax⁡,Φn,m=nΔΦm=nπLΦ′′,LΦ′′=2i≥Int(χMΦ)+1,LΦ=LΦ′′-MΦ,F̲(Θ,Φ)=RejβRE̲(R,Θ,Φ), and the other symbols have the same or analogous meanings as in (10).It is useful to note that the expansion (32) is based on a spherical modelling of the antenna. OSI expansions requiring a lower number of FF samples can be obviously derived when using more effective modellings tailored to the actual antenna geometry.


## 4. Numerical Tests

Many numerical simulations have been performed in order to assess the effectiveness and robustness of the proposed NF-FF transformation techniques with spherical scanning for antennas having two dimensions very different from the third one. Two sets of simulations are reported in the following. The former (from Figure5 to 11) refers to a quasiplanar antenna, whereas the latter (from Figure 12 to 16) is relevant to the case of an elongated one. In both the cases, an open-ended circular waveguide having radius 0.338λ, λ being the wavelength is chosen as measurement probe.Figure 5
Amplitude of the output voltageV′ on the meridian at φ=90°. Solid line: exact. Crosses: interpolated.In the first set, the numerical tests refer to a scanning sphere having radiusd=25λ and to three uniform planar circular arrays placed at z=-2λ,1λ,4λ, having radius equal to 18λ,19.5λ, and 16λ, respectively. Their elements, symmetrically placed with respect to the plane y=0, are elementary Huygens sources linearly polarized along y and spaced of 0.5λ along the radial and azimuthal lines. Such an AUT has been fitted by the double bowl modelling with a=20λ,c=4.5λ, and c′=2.5λ. Figures 5 and 6 show representative reconstruction examples of the amplitude and phase of the rotated probe voltage V′ on the meridian at φ=90°. As can be seen, the exact and recovered voltages are practically indistinguishable. To assess in a more quantitative way the performances of the algorithm, the maximum and mean-square reconstruction errors (normalized to the voltage maximum value over the sphere) have been evaluated for χ′=1.20, χ=1.10,1.15,1.20,1.25, and p=q ranging from 3 to 10 (see Figures. 7 and 8). They have been computed by comparing the reconstructed and exact voltage V′ on a spherical grid with Δϑ=Δφ=1°. As can be seen, the errors decrease up to very low values on increasing the oversampling factor and/or the number of retained samples, thus assessing the effectiveness of the interpolation procedure. The algorithm robustness has been verified by adding random errors to the exact samples. These errors simulate a background noise (bounded to Δa in amplitude and with arbitrary phase) and an uncertainty on the voltage samples of ±Δar in amplitude and ±Δα in phase. As shown in Figures 9 and 10, the algorithm works well also when dealing with error-affected data. At last, the developed interpolation algorithm has been applied to recover the NF data required to carry out the spherical NF-FF transformation [7], as modified in [15]. The reconstruction of the antenna FF pattern in the principal plane E is shown in Figure 11. As can be seen, the FF reconstruction, in spite of the significantly reduced number of used NF samples (see Table 1), is everywhere very accurate thus validating the effectiveness of the proposed approach.Table 1NF-FF transformation techniqueNdataStandard approach [7] (same number of data on each parallel)130 562Standard approach as modified in [15] (number of data on rings decreasing towards the poles)80 610Here proposed approach26 467Figure 6
Phase of the output voltageV′ on the meridian at φ=90°. Solid line: exact. Crosses: interpolated.Figure 7
Normalized maximum errors in the reconstruction of the voltageV′.Figure 8
Normalized mean-square errors in the reconstruction of the voltageV′.Figure 9
Amplitude of the output voltageV′ on the meridian at φ=90°. Solid line: exact. Crosses: interpolated from error-affected data.Figure 10
Phase of the output voltageV′ on the meridian at φ=90°. Solid line: exact. Crosses: interpolated from error-affected data.Figure 11
E-plane pattern. Solid line: exact. Crosses: reconstructed from NF data.Figure 12
Amplitude of the output voltageV′ on the meridian at φ=90°. Solid line: exact. Crosses: interpolated.The second set of figures refers to a scanning sphere having radiusd=35λ and to a uniform planar array of elementary Huygens sources polarized along the z axis, spaced by 0.5λ. These sources cover a zone in the plane y=0 formed by a 16λ×44λ rectangle ended in two half-circles. Accordingly, such an antenna has been modelled by a rounded cylinder with h′=44λ and a′=8λ. An example of voltage reconstruction relevant to the meridian at φ=90° is reported in Figure 12. Again, a very accurate reconstruction results. Figure 13 confirms the algorithm accuracy by showing the normalized maximum and mean-square errors for p=q ranging from 3 to 10 and χ′=χ=1.20. Also in such a case, the algorithm is resulted to be very stable with respect to random errors affecting the exact samples as shown in Figure 14. At last, the reconstructions of the antenna FF pattern in the E and H planes are shown in Figures 15 and 16, respectively. As can be seen, the exact and recovered fields are practically indistinguishable, thus assessing the effectiveness of the developed NF-FF transformation. Also in such a case, it can be interesting to compare the number of the employed NF data with those required by the other techniques (see Table 2).Table 2NF-FF transformation techniqueNdataStandard approach [7] (same number of data on each parallel)130 562Standard approach as modified in [15] (number of data on rings decreasing towards the poles)102 786Here proposed approach27 231Figure 13
Normalized errors in the reconstruction of the voltageV′.Figure 14
Amplitude of the output voltageV′ on the meridian at φ=90°. Solid line: exact. Crosses: interpolated from error-affected data.Figure 15
E-plane pattern. Solid line: exact. Crosses: reconstructed from NF data.Figure 16
H-plane pattern. Solid line: exact. Crosses: reconstructed from NF data.


## 5. Conclusions

Fast and accurate probe-compensated NF-FF transformations with spherical scanning for nonspherical antennas have been here developed by employing flexible AUT modellings, which allow one to fit very well the shape of a lot of actual antennas. In particular, an electrically long antenna has been considered as enclosed in a rounded cylinder, whereas a double bowl has been adopted to shape a quasiplanar antenna. The use of these modellings makes possible to remarkably lower the number of data to be acquired, thus significantly reducing the measurement time. The developed techniques work very well as widely assessed by the numerical simulations.



---

*Source: 936781-2011-03-21.xml*