The latest developments in Computer Aided Engineering, by means of the remarkable progress that several softwares have been made in terms of computational speed and of enlarged ranges of phaenomena that can be simulated and accurately analyzed, represent an immense span to inventiveness. For this particular instance, it is also a great opportunity for who wants to reach or already has an expertise in the Antenna Design field. In this article, among a vast amount of possibilities, some of the Matlab
toolboxes dedicated to antenna designing have been adopted in such a way as to create the structure of a coaxial-fed folded patch antenna combined with a low performance dielectric substrate (i.e. FR4). The main purpose was to obtain an antenna resonating at 2.1 GHz and with an input impedance matched at 50 Ohms. Subsequently, this designed radiating object has been introduced in a broader structure, that is a linear array of five antennas, as its representative single element. The Tchebyshev array synthesis has been used in order to reach particular design and radiating behaviour. Specifically, the beam-steering in the broadside and at 45ยฐ off the boresight direction have been analyzed by means of the overall gain and the electric and magnetic field behaviour close and far from the array. All the mathematical models which have been used together with the software feedback represent the structure of an overall process of development, optimization and analysis which have been very carefully performed and shown in the following paragraphs.
Authors:
- Blasi Luca
- Mastrofini Alessandro
- Mucenica Stefan Leonard
Wireless Electromagnetic Technologies Professor: Marrocco Gaetano University of Rome Tor Vergata 2022

Download full pdf report

Introduction
All the parameters necessary to generate the array structure are shown in this paper together with some considerations on the implications of running certain simulations and of the deriving results. The array of folded patch antennas has been designed entirely in Matlab
. Some flaws in the tools that have been used emerged in terms of lacking of full compatibility between them and some limitations in creating and analyzing the exact structure intended for this design. For example, the single folded patch antenna couldnโt be used as an array element in the Antenna Toolbox (which would have allowed studying the overall array behaviour in terms of gain and electromagnetic fields), while the Phased Sensor Array Analyzer
made possible creating the array structure but it couldnโt give results describing the electromagnetic field. This puzzling situation made it necessary to reach a compromise by using a similar structure to that of the folded patch antenna: a PCB stack shorting the patch with a series of small cylindrical pins instead of using a rectangular wall of metal as in the original design recommendation. Some results could be optained both with the original structure and with its approximation: in those cases several comparations shown their radiating behaviour is very similar, thus the approximation represents an acceptable compromise.
Starting from the theoretical modeling based on mathematical formulations [1], a set of parameters related to the single antenna (the reflection coefficient, the size of the patch and the ground, the substrate thickness, the feed location) have been optimized based upon a mesh density refinement process.
Generally, the study of an antenna array starts with the design of the array factor. The array factor is the complex-valued far field radiation pattern obtained by considering the elements to be punctiform sources. The total field is both influenced by the angular filtering effect of the array factor (which disregards the particular structure of the array element, but taking into account informations as spatial arrangement and currents circulating on the array alements) and the single element radiation. This gives an idea of how the behavior of the single antenna is affected by the other antennas of the array. All the elements are participating to those mutual interactions, bringing out a different behaviour in all.
Subsequently, several tests were conducted in order to verify that the the overall behavior of the array complies with the project requirements. A flowchart of all the steps of the project is available at the end.
Tchebyshev array factor design
A Tchebyshev array factor will be designed in this part, so that the specific structure of the single element of the array is disregarded for the moment. The array is represented by a linear distribution of five elements (๐๐๐ = 5), which are supposed to be uniformly spaced and with a non-uniform amplitude for the feed arrangement, that anyway is going to be symmetrical (starting from the element in the middle, located in the origin of the geometrical axes). Also the “mean lobe to side lobet ratio” (๐ ) will be one of the input design variables. Thus, the variables of this design part are the element inter-spacing (๐๐๐๐ก, which needs to be optimal such that the beamwidth is minimized), the feed amplitude of the single antenna (each one represented by a ๐ถ๐, with ๐ โ {โ2, โ1, 0, 1, 2}, the ๐๐กโ current coefficient – these unknowns are actually just three and not five, due to the symmetrical amplitude distribution hyphotesis), the operating resonant frequency required for the single antenna (๐) and some other quantities related to them, such as the tapering efficiency (๐๐ ) and the beamwidth (๐ต๐fn). Moreover, the difference between the beamwidth of the Tchebyshev array (Uniform Spacing Non Uniform Amplitude case, shorten ๐๐๐ด) and that of a uniform array having the same number of elements and inter-spacing value but a uniform amplitude feed distribution (Uniform Spacing Uniform Amplitude, shorten ๐๐ด) will be discussed. The input design variables related to the Tchebyshev array factor (i.e. ๐๐๐, ๐ and ๐)are listed in table 1.


3
Starting from the optimal inter-element spacing, it requires some secondary variables in orderer to be evaluated. Firstly, a coefficient (๐พ, depending on ๐ and indirectly on the number of array elements) neeeds to be calculated, among with the wavelength in the free space (๐, as the ratio between the light speed in the free space – ๐ – and the operating frequency ๐). Having those quantities, ๐๐๐๐ก minimizing the ๐ต๐fn can be computed by using the eq. (1) group.
\begin{gathered} d_{opt}\,{\leadsto}\,\min\{BW_{\operatorname{fn}}\}\\ \\ d_{opt}\,{=}\,\lambda\,\left[1\,-\,\frac{\arccos\left(\frac{1}{\gamma}\right)}{\pi}\right] \quad (1)\\ \\ \gamma\,{=}\,\cosh\left[\frac{1}{2N}\,\ln\left(R+\sqrt{R^2-1}\right)\right] \end{gathered}
The array factor evaluation through the Tchebyshev polynomial comes as the next step. The Tchebyshev polynomial approximation to the second order (๐2(๐ฅ)) is sufficient so the current coefficients can be obtained, remembering that a symmetrical amplitude distribution hypothesis has been made. In this case the current coefficients ๐ถ๐ , ๐ โ {โ2, โ1, 0, 1, 2} follow the rule ๐ถ๐ = ๐ถโ๐, thus their symbolic representation can be simplified as follows: ๐ถ๐ , ๐ โ {0, 1, 2} (or just ๐ = 0, 2). Since the Riblet variation to the Dolph- Tchebyshev synthesis model has been used, the variable of ๐2(๐ฅ) becomes ๐ฅ = ๐ + ๐ cos(๐ข) and the formulation for ๐๐๐๐ก โ (๐โ2, ๐] differs from the standard model. The coefficients ๐ and ๐ are related to the maximum value selected in the sub-domain of ๐2(๐ฅ) (the window of visible radiation lobes). This maximum, called for example ๐ฅ1, corresponds to the main lobe in which it can be convertend by evaluating the array factor |๐2(๐ฅ1)| (some additional references to the specific formulas that need to be used in order to compute can be found in [1] and [6]). That said, the current coefficients can be extracted from ๐2(๐ฅ) (so by using eq. (2)). Resulting uniquely out of a ๐ถ๐ dependence, the tapering efficiency is obtained by using eq. (3).
Next, both non uniform amplitude (Tchebyshev array, ๐๐๐ด) and uniform amplitude (๐๐ด) cases are compared. The comparison shows how the ๐ต๐fn in the ๐๐ด case (i.e ๐ต๐[๐๐ด] fn ) is narrower than that of the ๐๐๐ด case (i.e. ๐ต๐[๐๐๐ด] fn ). This result (see eq. (4) and table 2) is generally effective and the comparison has been made to show that this particular design case confirms the general condition.
All the numerical results related to the Tchebyshev array design are gathered together in table 2. After these quantities are calculated, some important design considerations can be made about the array efficiency. A Non-Uniform Amplitude Array Factor has been designed by using the Riblet variation of the Dolph-Tchebyshev array synthesis model. Considering the maximum to minimum feed ratio:
r_{\max/\min}={C_{\max}\over C_{\min}}
the less ๐ max โ min is, the more efficient distribution of current is reached. In this particular design, the requirement was to to find the ๐๐๐๐ก which minimizes the beamwidth, starting from the input design variables. Thus, the ๐max โ min is a straight consequence of the current coefficients and its optimal value has not been the seek of this project. Anyway, for this design, ๐max โ min โ 4.39 meaning that if a damage of the element with the ๐ถmax (i.e. ๐ถ0) level of feed occurs, most part of the efficiency will be lost. In any case, the tapering efficiency shows how it will not be possible to take advantage of 21% of the array in an ideal situation, remembering that this design model can be discerned by the real circumstance in terms of the Tchebyshev error (see [1]).
\begin{split} \begin{aligned} T_2\left[x=a+b\cos\,u\right]=\dots\quad\\ =\,C_0\,+\,2\,C_1\,\cos\,u\,+\,C_2\cos\,2u\quad (2)\\ =\,(2a^2+b^2-1)\,+\,4ab\,\cos\,u+b^2\,\cos\,2u\quad \end{aligned} \end{split}\\
\eta_T\,=\,\frac{1}{2N+1}\,\frac{||C_0+2C_1+2_2||^2}{C_0^2+2C_1^2+2C_2^2} \quad (3)
\begin{gathered} BW_{\operatorname{fn}}^{[UA]}\,{<}\,BW_{\operatorname{fn}}^{[NUA]}\\ \\ BW_{\operatorname{fn}}^{[NUA]}\,{=}\,2\frac{180}{\pi}\left[\frac{\pi}{2}-\arccos\left(\frac{\arccos\left(\frac{\cos\left(\frac{\pi}{2N}-a\right)}{b}\right)}{k_0d}\right)\right]\quad (4)\\ \\ BW_{\operatorname{fn}}^{[UA]}\,{=}\,\frac{2\lambda}{N\,d}\,\frac{180}{\pi} \end{gathered}
In the end, the 2D array patterns resulting by the use of the calculated parameters are shown in fig. 1. Two polar patterns and thei corresponding rectangular ones have been plotted (in the azimuth cut plane and in the elevation cut one).
Rectangular folded patch antenna design
The main components of a rectangular folded patch are: the patch, the substrate (generally accessory, but used in this project), the ground, the rectangular shorting wall between the patch and the ground, and the feed. More details about them will be presented in a short while. Before that, some other remarks are necessary: this antenna will be the element of the array, which will be designed starting actually from a PIFA (Planar Inverted F Antenna). A general PIFA realized with a dielectric substrate is shown in fig. 2. A particular condition is imposed to the PIFA, by which the width of the rectangular shorting wall (๐ค๐ ๐ ) equals the patch width (๐patch). so that the PIFA and the folded patch antenna will be two equivalent structures (see eq. (5)). This remark on the PIFA is necessary because generally its structure is not equivalent to that of the folded patch antenna because of the possible variability of the shorting width (๐ค๐ ๐ ), which doesnโt always satisfy the above imposed condition. This condition on ๐ค๐ ๐ represents one of the input variables of the PIFA design, such as the matched input resistance (๐ ๐๐), the resonant frequency (๐) and the substrate characteristics (see table 3). A preliminary evaluation of the patch parameters have been realized by leaning on a theoretical set of formulas (see [1]).

diagrams
Thatโs the reason why the characteristics of the PIFA, shown in table 3, are distinguished between pre-optimal and post-optimal (same thing applies to the ground component). Thus, an optimization process of all these and other parameters will be performed in some following steps. Just before that, now the theoretical model will be shortly pointed out: it furnishes a provisional set of results (that will be adjusted afterward) . Firstly, the length and the width of the patch starting from the ๐ค๐ ๐ = ๐patch imposed condition, the substrate wavelength and thickness (i.e. ๐FR4 and โFR4) are calculated (eq. (6)). Then, other parameters, such as the effective permittivity ๐๐๐๐ and length (๐ฟ๐๐๐ ) in eq. (7), the radiation resistance (๐ ๐) in eq. (8), the half-power beamwidth in the E-cut (ฮ๐ธ) and the H-cut (ฮ๐ป) in eq. (9) and eventually the feed location, defined by ๐feed in the patch length direction with respect to the free edge (eq. (10)), have been calculated.
Refinement with Matlab Method of Moments

Substrate thickness selection
Three thickness levels were available for the ๐น๐ 4 substrate required in this project (see table 3). In this part it will be explained why choosing a thinner substrate (if the ๐น๐ 4 is used) is more convienient. This choice will be motivated both in relation to the physical behaviour (radiation efficiency) and to the limits of the Matlab tools that have been used (in terms of the mesh density level selection in a range that gives more reliable results).
W_{\operatorname{patch}}\,=\,w_{sc} \quad(5)
\begin{aligned}L_{\operatorname{patch}}\,+\,W_{\operatorname{patch}}\,-\,w_{sc}\,&{=}\,\frac{\lambda_{\operatorname{FR4}}}{4}\,+\,h_{\operatorname{FR4}}\quad (6)\\ W_{\operatorname{patch}}\,&{=}\,\frac{\lambda}{2}\,\sqrt{\frac{2}{\varepsilon_{\operatorname{FR4}}\,+\,1}}\\ \end{aligned}
\begin{aligned} \varepsilon_{eff}\,{=}\,\frac{\varepsilon_{\operatorname{FR4}}+1}{2}\,+\,\frac{\varepsilon_{\operatorname{FR4}}-1}{2}\,\left(1+12\frac{h_{\operatorname{FR4}}}{W_{\operatorname{patch}}}\right)^{-\frac{1}{2}}\\ L_{eff}\,{=}\,\frac{\lambda_{\operatorname{FR4}}}{4}\\ \Delta L{=}0.412\,h\,\left[\frac{(\varepsilon_{eff}+0.3)\,\left(\frac{W_{\operatorname{patch}}}{h_{\operatorname{FR4}}}+0.268\right)}{(\varepsilon_{eff}-0.258)\,\left(\frac{W_{\operatorname{patch}}}{h_{\operatorname{FR4}}}+0.8\right)}\right]\quad (7)\\ L\,=\,L_{eff}-2\Delta L \end{aligned}
R_r\,{=}\,\frac{120\,\lambda}{W_{\operatorname{patch}}}\,\left[1-\frac{1}{24}\,\left(2\pi\,\frac{h_{\operatorname{FR4}}}{\lambda}\right)^{2}\right]^{-1}\quad (8)
\begin{aligned} \Theta_E\,&{=}\,2\,\arccos\,\sqrt{\frac{7.03\,\lambda^2}{4\,(3\,L_e^2+h_{\operatorname{FR4}}^2)\,\pi^2}} \quad(9)\\ \Theta_H\,&{=}\,2\,\arccos\sqrt{\frac{1}{2\,+\,2\pi\,\frac{W_{\operatorname{patch}}}{\lambda}}}\\ \end{aligned}
\ell_{\operatorname{feed}}\,{=}\,\frac{L_{\operatorname{patch}}}{\pi}\,\arccos\,\sqrt{\frac{R_{in}}{R_r}} \quad(10)
The Antenna Toolbox
gives specific information about the mesh density level that should be adopted for the design of the patch antenna components. The only issue is that these details are given only for particular ranges of the ratio indicator called relative thickness or electrical thickness โ๐ (see [2] and [4]). The electrical thickness depends on the ratio between the substrate thickness (โFR4) and the wavelength related to the substrate medium (๐FR4). When a mesh is configured in the Antenna Toolbox environment, a specific parameter needs to be adjusted: the maximum edge length of the generic triangle covering the geometry of the antenna (๐max). In the case of a relative length โ๐ comparable to 1โ10, itโs recommended to select an ๐max โ
๐โ10. A substrate thickness respecting this relationship is called a thick substrate. None of the available substrates verifies this condition. Among them, only the thinnest substrate and the second to last one (thus โFR4 = 0.8๐๐ and โFR4 = 1.0๐๐) are part of a range which the Antenna Toolbox provides instructions of. Itโs the thin substrate range: the automatic mesh mode should be adopted for a thin substrate, namely having a relative thickness less or equal than one fiftieth (โ๐ โค 1โ50). This leads to a further explanation about the actual substrate thickness chosen for this project. The thinner substrate choice rationale starts with the consideration of its quality factor (๐) depending on the loss tangent (๐ก๐๐๐ฟ) and being generally low in the FR4 substrate case (being the two quantities inversely proportional, i.e. ๐ โ 1โ๐ก๐๐๐ฟ, and being (tan ๐ฟ) very high compared to that of other more efficient substrates). This means the FR4 is a big power dispersor. Since increasing โFR4 will provoke just more losses in terms of a radiation efficiency drop and since the only thickness values of 0.8๐๐ and 1.0๐๐ would give reliable/accurate results in the Antenna Toolbox
simulations, the 0.8๐๐ thickness level will be adopted.

Mesh density refinement
Although a mesh density choice has already been made by selecting the best maximum edge length ๐max, the accuracy achievable by using the mesh automatic mode in the case of substrates belonging to the thin substrate range will be proved hereafter. An initial study of the mesh density level influence on the reflection coefficient (ฮ in ๐๐ต) evaluated at the resonant frequency (๐ = 2.1๐บ๐ป๐ง) has been realized, thus a ฮ2.1๐บ๐ป๐ง = ๐น(๐max) function has been plotted with an initial step of ฮ๐max = 2.5.10โ4 ๐ between every two mesh densities related to their specific ๐max. This first simulation considered a broader range of ๐max variation: [2.5.10โ4 ๐, 6.0.10โ4 ๐].
![Figure 3: Minimum of the reflection coefficient ฮ [๐๐ต] in the frequency range 2.0 รท 2.2๐บ๐ป๐ง depending on the varying mesh
density level](https://alessandromastrofini.it/wp-content/uploads/2022/01/image-5-1024x549.png)
Since the resulting plot (fig. 3) has shown big uncertainty of the reflection coefficient value at the resonant frequency (ฮ2.1๐บ๐ป๐ง) at almost every mesh ๐max level (primarily due to the big step selected between one density level and another), some more detailed tests have been run, by considering a slightly narrower range ([2.5.10โ4 ๐, 5.0โ4 ๐]) and a tchicker evaluation of the maximum edge values (so that the mesh variation step has been remarkably reduced to 1.0.10โ4 ๐). Specifically, the step between two mesh density levels in terms of the maximum edge length of each one has been reduced from a ฮ๐๐ = 2.5.10โ4 to ฮ๐๐ = 1.0.10โ4. In all those simulations, an important fact needs to be noted. Even very small variations on the maximum edge length value involved considerable inconsistencies in almost every part of the mesh range in terms of considerable variations of the frequency at wich ฮ reaches its minimum (min(ฮ), therefore the resonant frequency of the antenna changes very easily). Thus, considering the frequency ๐โ at which min(ฮ) is actually obtained, instead of evaluationg it at the theoretical resonating frequency value at every mesh level, not only the standard test comparing ๐max and ฮ has been run, but also some mesh refinement plots representing the relationship between ๐โ and ๐max (fig. 8c), ฮ๐โ and ๐max (fig. 8b) and also min(ฮ) and ๐max (fig. 8a) have been taken into account (where ฮ๐โ is the difference between ๐โ and the resonant frequency ๐ = 2.1๐บ๐ป๐ง). More parameter relationships have been collected and this led to the setting of the mesh density choice in terms of ๐max that has been selected inside the most stable region (i.e. showing the smallest deviation of the reflection coefficient minimum from the resonant frequency). In the end itโs been specifically taken the โautomaticโ ๐max (= 3.5โ4 ๐) suggested by the Antenna Toolbox, since this value belongs to the stable region and seems to give the most accurate results. The ๐max values belonging to the stable region ([3.1 . 10โ4 ๐, 3.7.10โ4 ๐]) exibit slight deviations from the resonant frequency (ฮ๐โ โ [0.01๐บ๐ป๐ง, 0.03๐บ๐ป๐ง]) and the minimum of the reflection coefficient varies in the range [โ24 ๐๐ต , โ33 ๐๐ต].

Patch parameters refinement
After the selection of the maximum edge length of the mesh (consequently, of its density level), a more refined computation of the reflection coefficient will be made, depending on the patch size (i.e. on its length and width), but also on the feed location. Firstly, only a parametric variation of the feed position across the patch length direction has been considered, depending on variations of ๐ฟpatch and ๐patch. This means that the first refinement of the feed position has been evaluated starting from its theoretical equation (depending indirectly by ๐patch). The change of the feed location has been taken into account in the computation of every step of the simulation, thus in every evaluation of the reflection coefficient, keeping ๐max as a constant (the previously selected length). The patch size variations provoked wide modifications of the reflection coefficient, which values depending on that have been represented by an initial contour plot (with variations of the patch size in a broader range and with a larger step between one value and another, see fig. 5a). After that, another simulation (fig. 5b) in a narrower range of the patch size variation has been run in order to choose from there a set of ฮ values (related to a set of coupled values (๐ฟpatch,๐patch)) that should put the patch antenna in the best resonant condition (which means that a ฮ as close as possible to zero in the linear scale and as negative as possible in the ๐๐ต, logarithmic scale is sought). The logarithmic (๐๐ต scale) was used for the contour plots because this way the resulting ฮ values are easier to distinguish from one another even grafically and mathematically.


A set of 20 values of ฮ and respective coupled values (๐ฟpatch,๐patch). has been selected from the second simulation range so that a more specific simulation could be run. In this third case (fig. 15), the reflection coefficient has been plotted in a range around the resonant frequency ([2.0๐บ๐ป๐ง, 2.2๐บ๐ป๐ง]) in order to find which is the best combination for the patch size that makes actually resonate the antenna at the project frequency. Another determining and discriminating factor was the input impedance (๐๐๐ = ๐ ๐๐ + ๐ฅ๐๐๐, where the real part of ๐๐๐ is the input resistance, while the imaginary one is the input reactance), because a impedance matching (at 50ฮฉ) needed to be achieved for the project. In the ideal case, of course, a reflection coefficient ฮ(๐๐) = 0.00 โ โโ๐๐ต would be required in order to reach the perfect impedance matching (perfect matching with input resistance at 50.00ฮฉ and null reactance). As an actual result, before seeking a better matching, the ฮ value related to all the couple candidates ((๐ฟpatch,๐patch)) ranged from โ24 ๐๐ต to โ30 ๐๐ต. The two best candidates of the ฮ value related to the coupled size (๐ฟpatch,๐patch) are plotted in fig. 15. Two refined plots (with narrower frequency steps in the same range) of ฮ and ๐ for one of the candidates are represented in fig. 4. An additional design strategy contributed to the final patch size choice: the feed location varying across the patch width direction (๐คfeed). Thus, a parametric impedance matching study depending on that has been run on the best couple candidates and a few others from the set of 20. The resulting values are listed in table 4 and the plots are shown in fig. 6.


antenna

From the PIFA to a PCB stack approximation
Since the Antenna Array Designer
(AAD) is not able to generate an array of PIFAs , the Sensor Array Analyzer
(SAA) has been used to create that type of structure, starting from the Tchebyshev array factor design made in the first part and from the previous optimization process of a single PIFA. In the last part the overall array of PIFAs needs to be analyzed by means of beamsteering for particular angles, total gain, electric and magnetic field patterns. Unfortunately, not all of this information can be extracted by using uniquely the the SAA tool. This fact led to the use of another tool, the PCB Antenna Designer
, so that all the parameters required for the project could be analyzed. With this last tool, a similar structure to that of the PIFA has been realized and introduced as the antenna element of the array distribution. Just before moving to the array analysis, a comparison between the single antenna being a PIFA structure and the antenna created by using a PCB stack was performed.

A single limitation in the use of the PCB Antenna Designer has been encountered, but it has been overcome with an approximation strategy. Thus, the PIFA and the PCB stack designed are not perfectly identical, but their single and overall performance inside the array had shown a slight quantifiable difference. The limitation consisted in the absence of a specific option or combination of commands that would allow to realize the rectangular shorting wall between the patch and the ground. The only thing the designer can do is to replace this kind of shorting pin with a series of small diameter (e.g. 0.4๐๐) cylindrical shorting pins (see fig. 16) and fig. 9 close to one edge of the patch, across the patch width direction. With this design choice, a very similar behaviour to that of the PIFA structure was possible to simulate. To support this results, a comparison of the 2D gain patterns (elevation and azymuth cut) of the PIFA to those of the PCB stack antenna has been realized in terms of the mean squared error values, shorten ๐๐๐ธ, related to the two specific two-dimensional patterns considered (๐๐๐ธ๐๐ = 0.55๐๐ต in the elevation cut directivity pattern and ๐๐๐ธ๐๐ง = 0.16๐๐ต in the azimuth cut directivity pattern). Furthermore, since the SAA allows the plotting of some patterns such as the directivity (but not of the gain), a comparison between the 2D patterns of the antenna array of PIFAs and the array of PCBs has been presented. In that case the error is presented as a comparison between the main lobe levels in the PIFA array and the PCB array cases, but also in terms of first first side lobe levels of the two arrays (see table 5). This discussion was necessary, because moving to the last part of the project required a deeper study of the array of antennas, that will be displayed by using the PCB antenna as the element of the array. A further discussion about the differences between the tools can be consulted in the Methods.
Overall array performance
In this last part of the project, the array factor and the the element antenna (PIFA or PCB stack) designs will be combined so that their total effect will be examined. As itโs been already mentioned in the previous paragraph, some of the information describing the overall array performance can be asked from both the SAA (where the array of PIFAs was designed) and the PCB/AAD (where the single patch antenna made starting from the PCB stack and also the whole array can be constructed). The performance of the overall array will be evaluated in two cases: in the broadside case (90โฆ) but also at 45โฆ off the boresight direction. First of all, it is very simple to identify the phase shift coefficients in the broadside case because they all equal 0โฆ (thereโs no actual phase shift between antennas). In the second case, some manual calculation can be made in order to insert the phase shifts between the single antennas (this is the case of the PCB stack array), and it can be also automatically computed by using the array of PIFAs. The generic procedure is shown in eq. (12). Just before moving to the last part of the analysis, where only the characteristics of the array of PCBs will be considered, a last comparison between the array of PCBs and the array of PIFAs will be made in terms of 2D directivity patterns (both in broadside and in 45โฆ off the boreside direction, see fig. 19). An error analysis based on the difference between the PCBs array main lobe and that of the PIFAs array, but also on the difference between the PCBs array main side lobe and that of the PIFAs array will be shown (see table 5). The patterns and the numerical results according to the numerical error values illustrate how the array of PCBs represent a good enough approximation of the array of PIFAs (see fig. 10 and fig. 19).


Total array gain
In order to obtain and describe the last plots of the project, it is not possible to recur to the SAA because some of the commands, such as and EHfields()
and patternMultiply()
, are not available in this tool. Thus, in this last part, only the array of PCBs will be used, because itโs compatible with these commands. The total array gain will be computed and displayed in two different cases: by means of a fullwave model (with no approximation) and by using the pattern multiplication principle. This second case is idealistic because it makes an approximation of the real patterns by not considering the mutual induced interactions between the single antennas of the array; however, a “mathematical” interaction is considered in this model because it comes out of a formulation that mixes the single element effect (which gives an element gain factor ๐บ0) on the gain pattern with the angular filtering influence of the array factor, weighted by a combination of the current coefficients (which gives an array gain factor ๐บ๐น ). The overall idealistic case is based on the pattern multiplication principle, which (in a case such as the Tchebyshev array synthesis model, thus with non uniform amplitude feed), will take into account also the effect of the tapering efficiency (๐๐ ) in its formulation (see eq. (13)). Thus, the pattern multiplication principle can be applied angle by angle and this will provide a set of values that can be computed and inserted in matrices by using MatLab. Starting from the matrices, the 2D patterns can be represented. On the other hand, an easier way is to use the patternMultiply()
command and get the ideal patterns.

Electric and magnetic near fields
The near field characteristics (distinguishing between its inductive and radiative regions) have been analyzed in both the broadside and by beamsteering to 45โฆ off the boresight conditions. Two situations have been analyzed: the general near field behaviour (in its electric and magnetic field components) but also the specific orientation of the fields in the center position of each of the five elements composing the antenna (so, considering the ๐ฅ and ๐ฆ coordinates of each patch corresponding to the length and width directions of its planar development and the ๐ง coordinate representing the vertical component of the distance between the observation point and the array). The vertical distance ๐ง has been considered in three different cases. First, the Fresnel region (inductive near field region) has been analyzed at ๐ง๐๐๐ = 0.62 โ ๐ฟ3 maxโ๐ (where ๐ฟmax is the maximum size of the array, thus ๐ฟmax = 4๐๐๐๐ก + ๐ฟpatch), while the radiative near field region was measured, starting from the Fraunhofer limit expression (๐ง๐๐๐๐ข๐ = 2๐ท2โ๐), at ๐ง๐๐๐๐ข๐โ2. Moreover, itโs been considered an observation point at a generic distance from the array in the far field region (4๐ง๐๐๐๐ข๐). The three cases can were compared both in the broadside and in the 45โฆ off the boresight direction examples. In the broadside case, the inductive and radiative near electric and magnetic fields can be ordered by following the increasing of their intensities as a criterion. It can be noted that the order of these intensities (calling them ๐ธ๐, ๐ป๐ , with ๐, ๐ = โ2, 2) is the same of that of the feed currents (the relationship is shown or the ๐ธ field intensities – see eq. (11) – but it works analogously for the ๐ป field). In the far field, a totally different phaenomenon is observable (at least in the broadside case and with good approximation in the 45โฆ direction too): the ๐ธ field components are all equal, as those relative to the magnetic field ๐ป. This can be interpreted as follows: in the near field, the elements of the array still donโt behave completely as one radiating object so that their individual radiative effect can be observed, prevailing on the mutual interactions that will make them work as a whole at a further distance (i.e. in the far field). This fact can be explicitly analyzed by looking at the electric and magnetic fields registered at the specific distances: everything is shown in table 6. The facts described above can be also observed in fig. 11.

\begin{aligned} &\text{If}\quad &C_0\,>\,C_{1}\,>\,C_2\\ &\text{and} \quad &C_0>C_{-1}>C_{-2}\\ &\text{then} \quad &E_0\,>\,E_{1}\,>\,E_2\\ &\text{then} \quad &E_0>E_{-1}>E_{-2}\\ \end{aligned} \quad (11)
Methods
In this conclusive paragraph, a list of the MatLab
commands that are not available for all the tools that have been used will be pointed out. Next, the computational costs related to the different parts of the project will be discussed. In the end, a comparison between the array of PIFAs realized in Matlab
and the CST Studio Suite (Student edition)
will be shown. The theoretical formulations extracted from [1] and adopted for this project regarding the pre-optimized PIFA model and the Tchebyshev array factor design have been also implemented in Wolfram Mathematica
, which enables the user to write and compute formulas in a very fast and clear fashion (for this purpose the software acts as a Computer Algebra System
). Explaining what are the limitations of every tool and our previous suggestions about how to combine them (this made it possible to get conclusive results) may help the reader who wants to engage himself into doing some similar studies about array of antennas. Four different MatLab
tools belonging to two different toolboxes have been used for this project. By using the Antenna Designer
tool, the PIFA structure has been generated and analyzed, while the PCB stack has been created and studied through the PCB Antenna Designer
and we acquired the array of PCBs from the AAD. All these tools belong to the Antenna Toolbox
([3]), while the SAA (by which the study of the array of PIFAs is possible) is part of the Phased SAA ([5]). The differences in terms of the informations than can be directly extracted by using commands from the two main toolboxes used for this projects are shown in table 7. We dedicate also a short paragraph to discuss about the computational costs. In terms of random access memory (RAM) occupied by the different elements during simulations, the single PIFA weights 960 MB, while the single PCB 1.2 GB; the array of PIFAs occupies 4 GB while the array of PCBs engages 5.6 GB of memory. The RAM occupation is not enough in order to describe the computational costs, but it represents just a starting point (it gives information about the “stone” that needs to be carried, but not about how long). Thus, the informations about the the time spent by the workstation adopted for the project to complete the calculations of different simulations are collected in table 8. The running time for one the radiation pattern simulation (in both single PIFA and PCB antennas and arrays) and the sparameters() function (which calculates a list of parameters which ฮ belongs to) for both the single PIFA and the PCB is shown in fig. 12. We realized a PIFA model also in CST Studio Suite and plotted its 2D radiation pattern diagrams. For the seek of a further comparison, the data related to the gain pattern in the azimuth cut have been exported as a .txt file which have been converted into an array in Matlab so that the gain pattern in this cut could be overlapped to that of the PIFA realized with the Antenna Designer (see).




\begin{aligned} u_0&=\alpha d_{opt}\,=\,2\pi\frac{d_{opt}}{\lambda}\,\cos(\theta_0)\\ \alpha_n&=n\,\alpha\,d_{opt}=nu_0\qquad \left(n=\overline{1,5}\right) \end{aligned} \quad (12)
G_{total}=\,\eta_T\,G_0\,G_T\quad (13)
Conclusion
The different analysis conducted in this paper have fulfilled the project requirements (none of them has been omitted). The toolbox I toolbox which have been used made possible to reach the final results (related to electric and magnetic fields) using the PCB structure as the antenna element of the array since the array of PIFAs isnโt compatible for this type of study. Using the array of PCBs introduced an admissible set of errors. The use of the array of PCBs led to a notable increase of the computational costs and this may represent an obstacle to users that would like to replicate the results but having machines who donโt have enough memory resources.
Data availability
All the Matlab
code is available in the online repository (file final_wet_project.m
) at the following link: https://github. com/mastroalex/antenna-design. There are also available more plot and information about antenna and array properties. There are all the raw data for plot and calculation. In the appendix there are different patterns (for PCB and PIFA comparison), 3D gain patterns, Smith diagram, and different antenna / array diagrams. There is also a project workflow and some code snippets.
References
- [1] Constantine A. Balanis. Antenna Theory: Analysis and Design. John Wiley & Sons Inc, 2016.
- [2] Sergey N. Makarov. Antenna and EM modeling with MATLAB. A John Wiley & Sons, Inc., publication, 202.
- [3] MathWorksยฎ. Antenna Toolbox. 2022. url: https://it.mathworks.com/products/antenna.html (visited on 11/01/2022).
- [4] MathWorksยฎ. Meshing. 2022. url: https://it.mathworks.com/help/antenna/ug/meshing.html (visited on 11/01/2022).
- [5] MathWorksยฎ. Phased Array System Toolbox. 2022. url: https : / / it . mathworks . com / products / phased – array.html (visited on 11/01/2022).
- [6] Sophocles J. Orfandis. Electromagnetic Waves and Antennas. The MathWorks, Inc., 1999-2016.
Appendix








for i=1:length(v)
mesh(p,'MaxEdgeLength',v(i));
s = sparameters(p,freqRange);
slog=20*log10(abs(s.Parameters));
ssi=slog(:,:);
freq(i)=s.Frequencies(ssi==min(ssi));
smin(i)=min(ssintax);
disp('Step: ');
disp(length(v)+1-i);
close all
end
u0 = -k0*dopt*cos(pi/4);
an = zeros(5,1);
for i=0:4
an(i+1,1) = i*u0;
end
kfa_ster=kfa;
kfa_ster.PhaseShift = an*(180/pi);