# Shaping Single Offset Reflector Antennas Using Local Axis-Displaced Confocal Quadrics

Authors: Rafael A. Penchel; José R. Bergmann; Fernando J. S. Moreira
Journal: International Journal of Antennas and Propagation
 (2016)

---

## Abstract

This work investigates a novel numerical procedure for the solution of an exact formulation for the Geometrical Optics synthesis of a single reflector antenna by simultaneously imposing Snell’s Law and Conservation of Energy in a tube of rays, yielding a second-order nonlinear partial differential equation of Monge-Ampère type, which can be solved as a boundary value problem. The investigation explores the interpolating properties of confocal quadrics to locally represent the shaped reflector surface. It allows the partial derivatives involved in the formulation to be analytically expressed. To illustrate the method, two examples of offset single reflectors shaped to radiate a Gaussian power density within a superelliptical contoured beam are presented. The results are validated by Physical Optics analysis with equivalent edge currents.


---

## Body

## 1. Introduction

Shaped reflector antennas can be employed to provide an optimized coverage with high gains and low cochannel interference. Several techniques based on Geometrical Optics (GO) have been used for the shaping of single and dual reflector antenna configurations. Among these techniques, [1–4] presented exact formulations for the GO synthesis of single and dual offset shaped reflector configurations by simultaneously imposing Snell’s Law and Conservation of Energy in a tube of rays. They lead to second-order nonlinear partial differential equations of Monge-Ampère (MA) type, which can be numerically solved as boundary value problems [1–4]. For the solution, the authors employed a finite-difference scheme associated to regular grids to approximate the derivatives. Then, the linearized form of the MA equation was iteratively solved [4–7]. Alternatively, [8] has proposed a GO synthesis technique for multiple reflectors where the surfaces are locally approximated by a biparabolic expansion and their parameters are obtained by solving linear equations. It has been effectively applied in designs of multiple reflectors to produce a uniform phase in front of an antenna aperture [8]. For circular symmetric geometries, [9, 10] presented a GO synthesis scheme based on the combination of conic sections to locally represent the reflector generatrices. Similarly, the conic parameters were also obtained by solving linear equations.In this work, motivated by the effectiveness of the local representation for shaped reflectors, we explore the use of axis-displaced confocal quadrics into the numerical solution of the MA equation. Instead of using finite-difference schemes to calculate the derivatives at grid points, the surface is locally represented by the confocal quadrics, where their parameters are obtained by interpolation, allowing the analytical representation of the derivatives, together with the simplification of the MA equation. It also allows some flexibility in the definition of the grid points that are used to improve the representation of the reflector surface. To illustrate the numerical scheme, it has been applied to single reflector synthesis, where the reflector surface yields a prescribed far-field radiation pattern when illuminated by a point source (feed horn), as shown in Figure1.Figure 1
Geometry of the offset reflector antenna.This work is organized as follows: Section2 presents the several steps involved in the synthesis technique. The shaped reflector surface is defined by a set of grid points and confocal quadrics that are employed to locally interpolate the surface. The analytical surface derivatives resulting from the interpolation process are replaced into the MA equation and the boundary condition at the reflector’s edge leads to a system of nonlinear equations. To solve the system, an iterative numerical technique is employed. Section 3 describes the application of the numerical scheme in the synthesis of two examples of single offset reflector antennas designed to generate a far-field beam with a superelliptical contour. A Physical Optics (PO) analysis plus edge currents technique is employed to verify the synthesized reflector antenna [11, 12]. The conclusions are presented in Section 4.


## 2. GO Shaping Technique

In the problem illustrated in Figure1, the cone of rays within the contour Ω′ and with radiation pattern Iη is reflected within a beam with contour Ω′′ where the power density is represented by Gζ. It is simpler to use the complex coordinates (stereographic projection) adopted in [2–4] to describe directions. η and ζ are such coordinates and here they represent the incident and reflected ray directions, respectively, defined as(1)η=cot⁡θ2eiϕ,ζ=cot⁡α2eiβ,where θ, ϕ and α, β are the spherical angular directions of the incident and reflected rays, respectively, as shown in Figure 2. Two auxiliary coordinate systems are used to describe ray directions: primed (′) and double primed (′′), related to feeder and far-field, respectively. The function r(η) describes the distance of the reflector point at η to the origin and Lη is a real function that represents the reflector’s surface [2–4]:(2)rη=eLηη2+1.The transformation τ:η→ζ is obtained by imposing Snell’s Law on the reflector and is expressed by(3)ζ=η+1Lη,where Lη is the derivative of L with respect to η [4].Figure 2
Feeder and far-field auxiliary coordinate systems.From Conservation of Energy, the relation between the feed radiation and the specified far-field power pattern is given by a nonlinear MA equation of elliptic (− sign) or hyperbolic (+) type in terms ofLη [4] as follows: (4)Lηη-Lη22-Lηη-2=±IηGζ1+ζ21+η22,where Lηη is the second partial derivative of L with respect to η and Lηη- is the second partial derivative of L with respect to η and its complex conjugate η-.


### 2.1. Reflector Surface Representation

Several iterative numerical techniques have been employed to solve the elliptical type of the MA equation, leading to a discrete representation of the surface at a grid of points. In [5, 7] the authors employed a 9-point cell finite-difference numerical scheme to approximate the derivatives of Lη in (3) and (4).Here, alternatively, the numerical solution employs axis-displaced confocal quadrics to locally describe the reflector’s surface, which allows an analytical representation of the derivatives and the simplification of (3). The solution yields a finite set of points on the reflector surface associated to a grid on ηη1,1,…,ηj,k,…,ηJ,K, where j and k describe the node position at the grid. Although the technique could be applied to any type of grid, in this work we adopt a polar grid for comparative purposes, where J and K are the numbers of rings and radials, respectively. The polar grid is defined in terms of θ′ and ϕ′ referred to the z′-axis (feed axis), as shown in Figure 2. The external ring of the polar grid defines the feed cone edge with semiangle θc. The reflector point associated to the center of the grid is assumed to be known for the numerical solution and can be used in a postsynthesis adjustment of the reflector dimensions, as the GO formulation involves only derivatives of the real function L employed to represent the reflector surface.For each ray directionηj,k corresponding to node j,k the surface point is defined by Lj,k with the help of (2). The quadric j,k that describes the surface around this point, illustrated in Figure 3, is given by [4] as follows:(5)aj,ke-L=η-+ηbj,k+iη--ηcj,k+η2-1dj,k-η2-1,where(6a)aj,k=-ej,kvj,k,(6b)bj,k=ej,ksin⁡γj,kcos⁡φj,k,(6c)cj,k=ej,ksin⁡γj,ksin⁡φj,k,(6d)dj,k=ej,kcos⁡γj,k;aj,k is thesemilatus rectum, ej,k is the eccentricity, vj,k is the distance between O and the quadric vertex Vj,k, and φj,k and γj,k are the azimuth and elevation angles of the quadric axis (see Figure 3). The four parameters aj,k, bj,k, cj,k, and dj,k are calculated by imposing that the quadric matches the reflector at the points associated with the 4-point cell around the grid node j,k, as depicted in Figure 4. This leads to analytical expressions for the parameters involving the values of L at the 4-point cell (Lj,k, Lj-1,k, Lj+1,k-1, and Lj+1,k+1). The polar grid in Figure 4 is generated with the help of (1) and incrementing θ1,k=0 to θJ,k=θc, with θj,k-θj-1,k=Δθ, and ϕj,1=0 to ϕj,K=2π-Δϕ, with ϕj,k-ϕj,k-1=Δϕ=2π/K.Figure 3
Local axis-displaced confocal quadric surface.Figure 4
Node scheme for the inner 4-point cells.By differentiating (5) with respect to η and its conjugate η¯ and replacing the derivatives of L into (4), one obtains(7)Iηj,kGζj,kζj,k2+1ηj,k2+12=bj,k2+cj,k2+dj,k2-1dj,k-1η-+bj,k-icj,k22,where Snell’s Law on the reflector imposes the relation between ηj,k and ζj,k, expressed by the following equation [4]:(8)ζj,k=dj,k+1-bj,k+icj,kη-j,kbj,k-icj,k+dj,k-1η-j,k.It should be noted that, for a quadric, Lηη-Lη2=0, yielding the simplification of (4).From (7) it is possible to define the operator Γj,k for all inner nodes of the grid as follows:(9)Γj,k=-Iηj,kGζj,kζj,k2+1ηj,k2+12+bj,k2+cj,k2+dj,k2-1dj,k-1η-j,k+bj,k-icj,k22which represents a nonlinear equation in terms of the values of L in the 4-point cell (Lj,k, Lj-1,k, Lj+1,k-1, and Lj+1,k+1) centered at node j,k. For j=1, Lj-1,k corresponds to the central point of the surface, that is defined prior to the numerical solution and is kept constant along the iterations.


### 2.2. Boundary Conditions

To apply the synthesis technique, we consider that the feed source radiates a conical tube of rays with semivertex angleθc and axis tilted by θo from the z-axis, as illustrated in Figure 2. After reflection, the rays are confined within a solid angle with perimeter Ω′′ in the plane ζ′′ given by [4] as follows:(10)ζ′′=ζoζ-1ζo-ζ.In the examples investigated in Section3, we have considered an approximately square perimeter Ω′′ described by a superellipse as follows: (11)12ζuζoζ+1ζo-ζ+ζoζ-+1ζo-ζ-2σ+12iζvζoζ+1ζo-ζ-ζoζ-+1ζo-ζ-2σ=1.Asσ→1, Ω′′ tends to be an ellipse with axes 2ζu=cot⁡αu/2 and 2ζv=cot⁡αv/2, and center at ζo=cot⁡αo/2, where 2αu and 2αv are the far-field beamwidths in the vertical and horizontal planes; respectively, and αo is the beam offset angle with respect to the z-axis, as shown in Figure 2. As σ→∞,Ω′′ tends to be a rectangle with sides 2ζv and 2ζu.The boundary condition is imposed by enforcing the rays at the feed cone perimeterΩ′ to map into the far-field superelliptical perimeter Ω′′ after reflection. For the numerical solution, at each ray direction ηJ,k corresponding to a boundary node J,k (see Figure 5), the real function L is approximated by a quadric given by (5). The parameters aJ,k, bJ,k, cJ,k, and dJ,k are calculated by imposing that the quadric passes at the reflector points associated with the 4-point cell around node J,k, as depicted in Figure 5. This leads to analytical expressions for the four parameters in terms of the values of L at the cell nodes (LJ,k, LJ-1,k, LJ,k-1, and LJ,k+1).Figure 5
Node scheme for the boundary 4-point cells.Once the expressions for the parameters are found, the mappingτ:ηJ,k→ζJ,k can be expressed with the help of (8) and replaced into (11). It allows one to define an operator for the grid nodes at the perimeter Ω′ as(12)ΓJ,k=-1+12ζuζoζJ,k+1ζo-ζJ,k+ζoζ-J,k+1ζo-ζ-J,k2σ+12iζvζoζJ,k+1ζo-ζJ,k-ζoζ-J,k+1ζo-ζ-J,k2σ.Equation (12) represents a nonlinear equation in terms of L at the 4-point cell around node J,k.


### 2.3. Numerical Iterative Scheme

The application of (9) or (12) at each grid node leads to a system of J×K nonlinear equations in terms of LL1,1,…,Lj,k,…,LJ,K. For the numerical solution of the system of nonlinear equations, Newton’s method has been employed [4–7]. It is based on an iterative scheme where convergence and stability depend on the choice of an adequate initial solution for the reflector’s surface. In this work, a single quadric was adopted as initial solution at all grid nodes, chosen to approximate the far-field beamwidth in the vertical plane. The iterative scheme continues until Γj,k<ε at all nodes, where ε is a specified maximum value for the operator residues.


## 2.1. Reflector Surface Representation

Several iterative numerical techniques have been employed to solve the elliptical type of the MA equation, leading to a discrete representation of the surface at a grid of points. In [5, 7] the authors employed a 9-point cell finite-difference numerical scheme to approximate the derivatives of Lη in (3) and (4).Here, alternatively, the numerical solution employs axis-displaced confocal quadrics to locally describe the reflector’s surface, which allows an analytical representation of the derivatives and the simplification of (3). The solution yields a finite set of points on the reflector surface associated to a grid on ηη1,1,…,ηj,k,…,ηJ,K, where j and k describe the node position at the grid. Although the technique could be applied to any type of grid, in this work we adopt a polar grid for comparative purposes, where J and K are the numbers of rings and radials, respectively. The polar grid is defined in terms of θ′ and ϕ′ referred to the z′-axis (feed axis), as shown in Figure 2. The external ring of the polar grid defines the feed cone edge with semiangle θc. The reflector point associated to the center of the grid is assumed to be known for the numerical solution and can be used in a postsynthesis adjustment of the reflector dimensions, as the GO formulation involves only derivatives of the real function L employed to represent the reflector surface.For each ray directionηj,k corresponding to node j,k the surface point is defined by Lj,k with the help of (2). The quadric j,k that describes the surface around this point, illustrated in Figure 3, is given by [4] as follows:(5)aj,ke-L=η-+ηbj,k+iη--ηcj,k+η2-1dj,k-η2-1,where(6a)aj,k=-ej,kvj,k,(6b)bj,k=ej,ksin⁡γj,kcos⁡φj,k,(6c)cj,k=ej,ksin⁡γj,ksin⁡φj,k,(6d)dj,k=ej,kcos⁡γj,k;aj,k is thesemilatus rectum, ej,k is the eccentricity, vj,k is the distance between O and the quadric vertex Vj,k, and φj,k and γj,k are the azimuth and elevation angles of the quadric axis (see Figure 3). The four parameters aj,k, bj,k, cj,k, and dj,k are calculated by imposing that the quadric matches the reflector at the points associated with the 4-point cell around the grid node j,k, as depicted in Figure 4. This leads to analytical expressions for the parameters involving the values of L at the 4-point cell (Lj,k, Lj-1,k, Lj+1,k-1, and Lj+1,k+1). The polar grid in Figure 4 is generated with the help of (1) and incrementing θ1,k=0 to θJ,k=θc, with θj,k-θj-1,k=Δθ, and ϕj,1=0 to ϕj,K=2π-Δϕ, with ϕj,k-ϕj,k-1=Δϕ=2π/K.Figure 3
Local axis-displaced confocal quadric surface.Figure 4
Node scheme for the inner 4-point cells.By differentiating (5) with respect to η and its conjugate η¯ and replacing the derivatives of L into (4), one obtains(7)Iηj,kGζj,kζj,k2+1ηj,k2+12=bj,k2+cj,k2+dj,k2-1dj,k-1η-+bj,k-icj,k22,where Snell’s Law on the reflector imposes the relation between ηj,k and ζj,k, expressed by the following equation [4]:(8)ζj,k=dj,k+1-bj,k+icj,kη-j,kbj,k-icj,k+dj,k-1η-j,k.It should be noted that, for a quadric, Lηη-Lη2=0, yielding the simplification of (4).From (7) it is possible to define the operator Γj,k for all inner nodes of the grid as follows:(9)Γj,k=-Iηj,kGζj,kζj,k2+1ηj,k2+12+bj,k2+cj,k2+dj,k2-1dj,k-1η-j,k+bj,k-icj,k22which represents a nonlinear equation in terms of the values of L in the 4-point cell (Lj,k, Lj-1,k, Lj+1,k-1, and Lj+1,k+1) centered at node j,k. For j=1, Lj-1,k corresponds to the central point of the surface, that is defined prior to the numerical solution and is kept constant along the iterations.


## 2.2. Boundary Conditions

To apply the synthesis technique, we consider that the feed source radiates a conical tube of rays with semivertex angleθc and axis tilted by θo from the z-axis, as illustrated in Figure 2. After reflection, the rays are confined within a solid angle with perimeter Ω′′ in the plane ζ′′ given by [4] as follows:(10)ζ′′=ζoζ-1ζo-ζ.In the examples investigated in Section3, we have considered an approximately square perimeter Ω′′ described by a superellipse as follows: (11)12ζuζoζ+1ζo-ζ+ζoζ-+1ζo-ζ-2σ+12iζvζoζ+1ζo-ζ-ζoζ-+1ζo-ζ-2σ=1.Asσ→1, Ω′′ tends to be an ellipse with axes 2ζu=cot⁡αu/2 and 2ζv=cot⁡αv/2, and center at ζo=cot⁡αo/2, where 2αu and 2αv are the far-field beamwidths in the vertical and horizontal planes; respectively, and αo is the beam offset angle with respect to the z-axis, as shown in Figure 2. As σ→∞,Ω′′ tends to be a rectangle with sides 2ζv and 2ζu.The boundary condition is imposed by enforcing the rays at the feed cone perimeterΩ′ to map into the far-field superelliptical perimeter Ω′′ after reflection. For the numerical solution, at each ray direction ηJ,k corresponding to a boundary node J,k (see Figure 5), the real function L is approximated by a quadric given by (5). The parameters aJ,k, bJ,k, cJ,k, and dJ,k are calculated by imposing that the quadric passes at the reflector points associated with the 4-point cell around node J,k, as depicted in Figure 5. This leads to analytical expressions for the four parameters in terms of the values of L at the cell nodes (LJ,k, LJ-1,k, LJ,k-1, and LJ,k+1).Figure 5
Node scheme for the boundary 4-point cells.Once the expressions for the parameters are found, the mappingτ:ηJ,k→ζJ,k can be expressed with the help of (8) and replaced into (11). It allows one to define an operator for the grid nodes at the perimeter Ω′ as(12)ΓJ,k=-1+12ζuζoζJ,k+1ζo-ζJ,k+ζoζ-J,k+1ζo-ζ-J,k2σ+12iζvζoζJ,k+1ζo-ζJ,k-ζoζ-J,k+1ζo-ζ-J,k2σ.Equation (12) represents a nonlinear equation in terms of L at the 4-point cell around node J,k.


## 2.3. Numerical Iterative Scheme

The application of (9) or (12) at each grid node leads to a system of J×K nonlinear equations in terms of LL1,1,…,Lj,k,…,LJ,K. For the numerical solution of the system of nonlinear equations, Newton’s method has been employed [4–7]. It is based on an iterative scheme where convergence and stability depend on the choice of an adequate initial solution for the reflector’s surface. In this work, a single quadric was adopted as initial solution at all grid nodes, chosen to approximate the far-field beamwidth in the vertical plane. The iterative scheme continues until Γj,k<ε at all nodes, where ε is a specified maximum value for the operator residues.


## 3. Design Examples

To illustrate the numerical scheme described in Section2, two examples of offset shaped reflector antennas are presented. They are shaped to provide, according to GO principles, a Gaussian far-field radiation pattern Gζ in a solid angle defined by a superelliptical perimeter Ω′′ in the ζ′′ plane and described by (13)Gζ=G0e-gρ′′,where the parameter g is chosen to specify the attenuation at the far-field boundary, requiring equal coverage at the superelliptical contour defined by ρ′′ as follows: (14)ρ′′=12ζuζoζ+1ζo-ζ+ζoζ-+1ζo-ζ-2σ+12iζvζoζ+1ζo-ζ-ζoζ-+1ζo-ζ-2σ.For both examples, the parameter g=1.382 is set to provide 6 dB attenuation at the contour Ω′′. A feed with a raised-cosine model Iη=I0cos2n⁡θ′ illuminates the reflector antenna, with n=9.605 to provide 12 dB attenuation at the reflector edge. The feed cone has an offset angle θo=130∘and half-angle θc=30∘ (see Figure 2) and it is vertically polarized.For the first example, the far-field coverage region is specified byαo=-18∘, αu=8∘, αv=12∘, and σ=1.6 (see Figure 2). The initial solution L0 for the iterative numerical scheme is an ellipsoid with a0=-40.893cm, b0=-0.077795, c0=0, and d0=0.62708, which provides a beamwidth equal to the specification in the symmetry plane, αu=8∘. To ensure convergence, the far-field perimeter Ω′′ is gradually deformed from a circle to the specified superelliptical perimeter. For the synthesis, we first employed a regular polar grid with J=15 rings, K=70 radials (i.e., 980 interior nodes and 70 boundary nodes), and center at the feed axis, as illustrated in Figure 6(a). Consequently, 1,050 local quadrics were used to attain the shaped reflector. A total of m=16 iterations were necessary to obtain Γj,k<ε=10-5 at all grid nodes, where Γmean=7.683×10-8. During the iterative process, we have observed that the largest residues occur close to the reflector’s edge, probably due to the superelliptical perimeter imposed at Ω′′. The shaped reflector has approximately an elliptical projection in xy plane. By adding a proper constant to the final solution Lη, the projected diameters in the symmetry plane ϕ=0 and orthogonal plane are 25 cm and 29 cm, respectively. Consequently, the distance from O to the final reflector along the ray at center of the polar grid becomes r0=27.958cm. Figures 6(a) and 6(b) illustrate the grids resulting from the synthesis in the complex planes η′ and ζ′′, respectively, for the initial and final reflectors.Figure 6
Uniform grid representation of the ray directions at the (a) feed (complex planeη′) and (b) far-field (complex plane ζ′′).
(a) (b)As observed in Figure6(b), the superelliptical perimeter is achieved at Ω′′. To compensate for the feed attenuation (12 dB) at the reflector edge, the synthesis moves far-field rings closer to each other toward the coverage perimeter. The increase in the concentration of points indicates the reflector area where the shaping is more intense when compared with the initial solution and some distortion may be observed in areas close to the symmetry plane Im⁡[1/ζ′′]=0. On the other hand, the regular polar grid intrinsically generates a denser grid of points at the most inner ring along the asymmetry plane which may bring some numerical instabilities along the iterative procedure. As mentioned before, the numerical scheme described in Section 2 presents flexibility to handle different spacings in the polar grid. To compensate these effects and provide a better description of the reflector surface, we alternatively employed a nonuniform polar grid in the η′ plane where the spacing between the rings follows approximately a Gaussian rule with larger Δθ at the center and shorter one at the outer rings, thus reducing the density of points at the center and increasing it close to the boundary [8]. Figures 7(a) and 7(b) present the new distributions of grid points in the complex planes η′ and ζ′′, respectively, and, as one can observe, the mapped points in the plane ζ′′ show a more coherent distribution, especially close to the boundary and the symmetry plane.Figure 7
Nonuniform polar grid representation of the ray directions at the (a) feed (complex planeη′) and (b) far-field (complex plane ζ′′).
(a) (b)Figure8 shows a cut in the symmetry plane of the final reflector and the corresponding rays obtained from the nonuniform grid arrangement of Figure 7. To validate the results, the shaped reflector was analysed by Physical Optics (PO) with equivalent edge currents [11]. To obtain a smooth description of the reflector, the reflector points obtained from the GO shaping were interpolated by a polynomial plus Fourier series [12]. Figure 9 compares the GO and PO radiation patterns in the symmetry plane (ϕ=0) at 28 and 38 GHz. A good agreement was achieved and the differences observed in Figure 9 can be attributed to diffractive effects not taken into account by the GO synthesis. Figures 10(a) and 10(b) show the contour plots in the uv-plane of the PO radiation patterns at 28 GHz and 38 GHz, respectively, which are compared with the desired GO superelliptical beam contour described by the dashed contour.Figure 8
Shaped reflector surface and rays at the symmetry plane (ϕ=0) of the first example.Figure 9
Radiation patterns in the symmetry plane (ϕ=0) of the first example.Figure 10
Radiation pattern (dBi) in theuv-plane at (a) 28 GHz and (b) 38 GHz of the first example.
(a) (b)In the second design example, the far-field coverage region has a superelliptical contour, with the longer axis along the symmetry plane, that is,αu=12∘ and αv=8∘. Different from the first case, now the stretching of the far-field grid is more intense in the asymmetry plane close to the boundary. The initial solution L0 for the iterative numerical scheme is an ellipsoid with a0=-37.469cm, b0=-0.007115, c0=0, and d0=0.510939, which provides a beamwidth coverage equal to the specified in the symmetry plane with αu=12∘. A total of m=15 iterations were necessary to obtain Γj,k<ε=10-5 at all grid nodes, where Γmean=3.706×10-8. Figure 11 illustrates the shaped reflector and the rays from the feed point source at the origin to the far-field in the symmetry plane. The shaped reflector has approximately an elliptical projection in plane xy and, by adding a constant to the synthesized solution Lη, the diameter in the symmetry plane is again adjusted to 25 cm, yielding the diameter of 29.75 cm in the asymmetry plane and the distance r0=28.32cm from O to the reflector along the ray at center of the polar grid. For this example, the same nonuniform grid in the complex plane η′ was employed (as illustrated in Figure 7(a)) and the resulting mapping in the complex plane ζ′′ is illustrated in Figure 12.Figure 11
Shaped reflector surface and rays at the symmetry plane (ϕ=0) of the second example.Figure 12
Representation of the ray directions at the far-field (complex planeζ′′) of the second example.Once again, the shaped reflector was interpolated with a polynomial plus Fourier series and analysed by PO [11, 12]. Figure 13 shows the GO and PO radiation patterns in the symmetry plane (ϕ=0) at 28 GHz and 38 GHz. A good agreement between the desired GO and PO radiation pattern results is observed. Figures 14(a) and 14(b) show the PO radiation pattern at 28 GHz and 38 GHz, respectively, in the uv-plane, where one observes the desired superelliptical beam contour achieved for the main lobe (the dashed contour represents the desired GO coverage).Figure 13
Radiation patterns in the symmetry plane (ϕ=0) of the second example.Figure 14
Radiation pattern (dBi) in theuv-plane at (a) 28 GHz and (b) 38 GHz of the second example.
(a) (b)


## 4. Conclusions

This work presented an alternative numerical scheme for the solution of the MA equation resulting from the exact formulation of the GO synthesis of offset reflector antennas. It explores the local interpolating properties of confocal quadrics to represent shaped reflector surfaces, which allows the partial derivatives involved in the formulation to be analytically expressed and provides a significant simplification of the MA equation. The numerical scheme was used to shape two reflectors radiating Gaussian patterns within superelliptical cones. The shaped reflectors were analysed by the PO approximation with equivalent edge currents and the results showed good agreement with the GO radiation pattern specifications. The method proved to be very fast and stable, achieving convergence with relatively few iterations. This behaviour is directly linked with the analytical determination of the derivatives∂Γj,k(m)/∂Lp,q. Different from the finite-difference approach, the present method can be extended to other types of grids and contours Ω′ with minor efforts.



---

*Source: 4715681-2016-12-29.xml*