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.


  • Blasi Luca
  • Mastrofini Alessandro
  • Mucenica Stefan Leonard

Wireless Electromagnetic Technologies 
Professor: Marrocco Gaetano 
University of Rome Tor Vergata 
Download full pdf report
Development and Analysis of an Array of Folded Patch Antennas Using Matlab


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.

Table 1: Array factor input design variables
Table 1: Array factor input design variables

Table 2: Tchebyshev array design results
Table 2: Tchebyshev array design results

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.

			d_{opt}\,{=}\,\lambda\,\left[1\,-\,\frac{\arccos\left(\frac{1}{\gamma}\right)}{\pi}\right] \quad (1)\\ 

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]).

			=\,C_0\,+\,2\,C_1\,\cos\,u\,+\,C_2\cos\,2u\quad (2)\\ 
	\quad (3)
		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)\\

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]).

Figure 1: Array factor polar, in azimuth (a) and elevation (c) cut, and rectangular in azimuth (b) and elevation (d) cut,
Figure 1: Array factor polar, in azimuth (a) and elevation (c) cut, and rectangular in azimuth (b) and elevation (d) cut,

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

Figure 2: PIFA realized with a dielectric substrate
Figure 2: PIFA realized with a dielectric substrate

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)\\
	\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
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)
	\Theta_E\,&{=}\,2\,\arccos\,\sqrt{\frac{7.03\,\lambda^2}{4\,(3\,L_e^2+h_{\operatorname{FR4}}^2)\,\pi^2}} \quad(9)\\
\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.

Table 3: Project input parameters (frequency, matching and substrate features), theoretical (pre-optimized, shorted "preopt") and optimized ("post-opt") features of the patch and ground PIFA components
Table 3: Project input parameters (frequency, matching and substrate features), theoretical (pre-optimized, shorted “preopt”) and optimized (“post-opt”) features of the patch and ground PIFA components

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
Figure 3: Minimum of the reflection coefficient ฮ“ [๐‘‘๐ต] in the frequency range 2.0 รท 2.2๐บ๐ป๐‘ง depending on the varying mesh density level

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 ๐‘‘๐ต].

Table 4: Final ฮ“ and impedance matching values after last optimization choice: the ๐‘คfeed change
Table 4: Final ฮ“ and impedance matching values after last optimization choice: the ๐‘คfeed change

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.

Figure 4: Reflection coefficient (a) and impedance (b) plots depending on ๐‘“ โˆˆ [2.0, 2.1๐บ๐ป๐‘ง]
Figure 5: First (a) and second (b) contour plots depending on the patch size (๐ฟpatch and ๐‘Špatch variations)

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.

Figure 6: Final ฮ“ (a) and impedance matching (b) plots after further refinement including ๐‘คfeed change.
Figure 7: Gain patterns (a) in the nulle elevation plane, (b) in the null azimuth plane and (c) 3D current plot on the patch
Figure 8: Comparison of the meshing parameter ๐‘’max with minimum of reflection coefficient (a), variation of the resonant frequency with respect to the design one (the dotted line represents the value with respect to automatic meshing) (b), resonant frequency identified by the minimum of ฮ“ (c)

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.

Figure 9: Cropped zoom over PCB via holes. For an overall picture see fig. 16

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).

Table 5: Comparison between the array of PIFAs and the array of PCBs in terms of the directivity evaluated in their main lobes and first side lobes, both in the broadside and 45โ—ฆ off the boresight direction, both in the azimuth cut (๐‘Ž๐‘ง) and elevation cut (๐‘’๐‘™) planes. The comparison has been made in terms of the directivity difference (ฮ”๐ท# = ๐ท๐‘ƒ ๐ผ๐น๐ด๐‘  โˆ’ ๐ท๐‘ƒ๐ถ๐ต๐‘  in ๐‘‘๐ต) in the corresponding positions and also as a relative percentage error (ฮ”๐ท๐‘Ÿ # %, the ratio between ฮ”๐ท# and the directivity in the array of PIFAs case (๐ท๐‘ƒ ๐ผ๐น๐ด๐‘ )
Table 6: Particular electric (๐ธ) and magnetic (๐ป) field values measured in both broadside (90โ—ฆ) and in the 45โ—ฆ off the boresight directions. The ๐ธ-field is measured in [ ๐‘‰ ๐‘š ] while the ๐ป-field in [ ๐ด ๐‘š ] . The fields have been measured with respect to the centers of each element antenna vertically off the array plane at different height levels of ๐‘ง๐‘“๐‘Ÿ๐‘’ (inductive region), ๐‘ง๐‘“๐‘Ÿ๐‘Ž๐‘ข๐‘› (radiative region) and 4.๐‘ง๐‘“๐‘Ÿ๐‘Ž๐‘ข๐‘› (particular coordinate in the far field region).

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.

Figure 10: 2D gain patterns (ideal โ–  and real โ–  ) in both the broadside case ((a) azimuth cut and (b) elevation cut) and 45โ—ฆ off the boresight direction ((c) azimuth cut and (d) elevation cut).

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.

Figure 11: Electric and magnetic fields in Fresnel region in both broadside (a) and 45ยฐ off the boresight direction (b). Same fields calculated in the center of the position of the antennas in both broadside (c) and 45ยฐ off the boresight direction (d). The contribution of the various coefficients is evident.
	&\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}\\
\quad (11)


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).

Figure 12: Analysis of the computational cost of PIFA vs PCB in the calculation of the radiation pattern (single antenna and array) and of the sparameters
Figure 13: 2D gain patterns in the azimuth cut: comparison between the PIFA created in Matlab (โ– ) and that designed in CST Studio Suite (โ– ).
Table 7: This table is listing some of the commands which have been used for this project. Some of them are partly or fully compatible with both the tools used for the array analysus and others have been implemented by Matlab to work only in one of them. The main issue of the simulation part was the array of PIFAs could be realized with SAA but using that tool didnโ€™t allow to study the electromagnetic field. Instead, the AAD permitted the realization of an array of PCBs, which electromagnetic field could be represented and discussed.
Table 8: The optimization process went through different steps, as described in the previous sections. Here, some of the computational duration required to completely execute some of the simulations (related to one of the workstations which have been used) are listed so the reader may have an idea of the computational costs for this project array
		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)


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.


  • [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: (visited on 11/01/2022).
  • [4] MathWorksยฎ. Meshing. 2022. url: (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.


Figure 14: EH Fields in radiative near field region for broadside (a) and 45ยฐ off boresight direction (b) and in far field region for broadside (c) and 45ยฐ off boresight direction (d)
Figure 15: The two best ฮ“ plots depending on the specific patch size expressed coupled variations (๐ฟpatch,๐‘Špatch) and represented in a frequency range of [2.0, 2.1]๐บ๐ป๐‘ง
Figure 16: 3D view and top view of the PIFA approximation realized with a PCB stack and cylindrical shorting pins
Figure 17: PCB array
Figure 18: 2D gain patterns (PIFA (โ– ) and PCB (โ– )) in both the broadside case (azimuth (a) cut and elevation (b) cut) and 45โ—ฆ off the boresight direction (azimuth (c) cut and elevation (d) cut). The curve with the dashed line has not been calculated regularly but mirrored, due to the limitation in the maximum elevation angle range which can be asked for in the SAA patterns.
Figure 19: PIFA (โ– ) and PCB (โ– ) single antenna directivity patterns (dB) in the Azimuth cut (๐œƒ๐‘’๐‘™ = 0โ—ฆ, (a)) and in the Elevation cut (๐œ™๐‘Ž๐‘ง = 0โ—ฆ, (b)). PIFA and PCB arrays directivity patterns (dB) in the Azimuth cut (๐œƒ๐‘’๐‘™ = 0โ—ฆ, (c)) and in the Elevation cut (๐œ™๐‘Ž๐‘ง = 0โ—ฆ, (d))
Figure 20: Gain pattern for PIFA single antenna (a) and PCB array (b). Smith chart (c)
Figure 21: Flowchart
for i=1:length(v)
	s = sparameters(p,freqRange);
	disp('Step: ');
	close all
	u0 = -k0*dopt*cos(pi/4);
	an = zeros(5,1);
	for i=0:4
		an(i+1,1) = i*u0; 
	kfa_ster.PhaseShift = an*(180/pi);