Planar array antennas are many antennas arranged and connected in a matrix format on a plane in order to reduce the inherent environmental pollution by reducing the grating lobes and levels   . Directivity in the other hand as a fundamental antenna parameter is a measure of how “directional” an antenna’s radiation pattern is  . It can also be seen as a 4π times the ratio of the far field power density to the total radiated power  .
Many array antennas were developed, ranging from linear, circular, phased, conformal and planar array antennas. Planar array of monopole antenna method is chosen for this research  . The total field radiated as well as the directivity by this antenna depends largely on the elemental parameters, such as inter-element spacing.
Many researchers referenced in our research like    and  had researched on different ways and methods of improving the radiation pattern synthesis and attempts to enhance the directivity of the radiated field using different methods and have different aims. But this research, in which no researcher has done before, provides an in-depth study on the effect of varying inter-element spacing on the directivity of the total field radiated by planar array of monopole antenna.
2. Design of Planar Array Using Hadamard Matrix
The “on” and “off” states in the array arrangement correspond to “1” and “−1” of the Hadamard matrix  . The inter-element spacing in both directions of the array was varied. In determining the configuration of the planar array in a matrix form, a special matrix arrangement technique was adopted. This matrix technique or method is known as Hadamard matrix which was introduced by Jacques Hadamard in 1893.
An matrix is an Hadamard matrix of order n if the entries of F are either +1 or −1 and such that , where is the transpose of F and I is the order n identity matrix.
Hadamard matrix is one of the mathematical conjectures in nature. Although many associated ideas have been developed, the very existence of these matrices has extensive consequences in many fields of research, such as optimal design theory, information theory and graph theory. For instance, a Hadamard matrix can be interpreted directly as a weighing design. They can be used in forming optimal fractional factorial designs, orthogonal arrays (It consists of two conventionally summed linear sub-arrays that are oriented at right angle to each other, such that the outputs of the sub-arrays are cross-correlated to form an antenna pattern.) F-square designs  .
From Hadamard states that if is Hadamard matrix of order m and is a Hadamard matrix of n then is a hadamard matrix of order mn.
The symbol denotes direct product or convolution of matrices: if A is a matrix with typical entry , then
A matrix is a Hadamard matrix of order n if the entries of F are either +1 or −1 and such that
That is let
For a 4 × 4 matrix arrangement, that is using Equation (3) above, that product of gives 
Equation (5) implies that F is nonsingular, and has an inverse ; consequently
But in this paper, we considered a uniform planar array antenna with all elements energised. From Hadamard, let
According to Hadamard matrix, is a Hadamard matrix of order N × M if the entries of F are either +1 or −1. In this paper we considered in Equation (7) where all the elements will to be +1, that is,
From the above Equations (8)-(11) it therefore implies that, in a uniform planar array antenna with all the elements in the array matrix energised and represented by “1” we have;
For clarity and easy placement of these elements in the planar array form, the nodal elements are numbered from 1 - 16 according the element’s position in the matrix array as shown in Figure 1.
From Figure 1 below, the element of each node is the element’s position number in the uniform planar array matrix. The inter-element spacings and are spacings along x and y-axis respectively.
10 elements with (+) sign before them, that is, elements 1, 2, 3, 4, 5, 9, 11, 13 are “on” while 6 elements with (−) sign before them, that is, elements 7, 8, 10, 12, 14, and 15 are “off” of the 4 × 4 array. The total elements involved in this planar array evaluation were 16 elements.
3. Derivation of the Formula for Planar Array Factor AFplanar from Linear Array Antenna
Linear antenna array provides a good background in array theory because of the insights that leads into formation of beam and it relates to the array excitation functions and the radiation pattern obtained  .
Consider an isotropic radiator located at the origin of the linear array, the E-field may be written as (assuming θ-polarization)
is the complex excitation of the isotropic radiator, k is the free space wave number and d is distance of the observation point from the origin.
Assume that the N elements of the array are uniformly spaced with a separation
Figure 1. Placement of element in planar array form  .
distance d. The uniformly spaced linear array antenna is presented in Figure 2 below  .
This arrangement is necessary showing that the configuration of a planar array antenna is made up of two linear antennas placed at orthogonally to each other. The nodes, 1 - 4 represent the four elements of the linear array and a, is the inter-element spacing. While d1 - d4 is the far field approximation.
From Figure 2 below, the far field approximations of the array is
The magnitudes of current of the array elements was assumed to be equal and the current on the array element located at the origin was used as the phase reference (zero phase)  see Figure 1.
Therefore, the currents .
The far E-field of the individual elements is:
The total array far field could be found by using superposition:
Figure 2. Uniform linear array antenna  .
Note that the array factor (AF) is independent of the antenna type, assuming all of the elements are identical. The array factor for a uniformly-spaced N-element linear array is  :
We consider a uniform array defined by uniformly-spaced identical elements of equal magnitude with a linearly progressive phase from element to element:
Inserting this linear phase progression into the formula for the general N-element array gives
where is a function which defined as the array phase function and is a function of  the element spacing, phase shift, frequency and elevation angle. If the array factor of Equation (25) is multiplied by the result is
The subtraction of the array factor from Equation (26) above gives
Therefore, according to  and 
The complex exponential term in Equation (28) is the phase shift (array phase factor) of the array phase center relative to the origin. If the position of the array is shifted so that the center of the array is located at the origin, this phase term goes to zero. That is if ,
then Equation (23) becomes;
The array factor is normalized so that the maximum value for any value of N is unity  . The normalized array factor therefore is given as
The antenna is rectangular array antenna as shown in Figure 3. The interest is the arrangement of a linear array which is of N × N matrix with M number of element in x equals N number elements in y  .
The designed principles for planar arrays are similar to those elements placed in two dimensions (Figure 3), the array factor of a planar array can be expressed as the multiplication of the array factors of two linear arrays of Equation (30) above, one along the x-axis and the other one along the y-axis. Therefore, planar array factor can be simply expressed as:
where, and is the array phase function along x-and-y axis while and is the progressive phase shift along x-and-y axis respectively. and , is the inter-element spacing x-and-y axis respectively when and can be obtained as
Figure 3. Geometry of the N × M planar array antenna  .
The nulls of the array function are found by determining the zeros of the numerator term where the denominator is not simultaneously zero  . That is:
The peaks of the array function can be found by determining the zeros of the numerator term where the denominator is simultaneously zero.
4. The Wavelength (Lambda, 𝜆)
It has been established that the speed of electromagnetic wave on air is 3 × 108 m/s, and can be obtained from the formula,
f = the frequency of the electromagnetic wave.
= Lambda, the wavelength.
For the purpose of this investigation, the frequency is 2.5 GHz, therefore, the wavelength, can be obtained as follows;
The (AF) in this research is calculated using the normalized (AF) obtained in Equation (28)
For a beam pattern broadside
Let for a uniform array and . , .
5. Determination of the Directivity Dx and Dy of a Linear Array
The directivity of linear array antenna D, is determined by the formula,
where, is the voltage ratio, L is the length of the linear array, f is the broadening factor of broad array and a is the distance between elements in the linear array  . Let voltage ratio, .
6. Determination of the Directivity D0 of Planar Array
Since the configuration of a planar array antenna is consists of two linear antennas placed orthogonally to each other, the directivity of a planar array antenna, can be obtained as  :
where, = The directivity of a linear array along x-axis.
= The directivity of a linear array along y-axis.
Table 1 is the result of the computation of varying the Inter-Element spacing
Table 1. Result of varying the inter-element spacing.
of the Planar Array of Monopole Antenna in both x and y-directions keeping progressive phase shift, theta/phi = 45, number of elements N = M = 4 constant throughout the variation. The table contain the values of number elements along both x and y-directions, progressive phase-shift and directivity of fields.
It also shows a decrement on the array factor, as the inter-element spacing was increased, for instance, when the element spacing of 0.47λ was mathematically evaluated, the correspondent value of array factor is 0.9984 as the value of the list spacing considered in this evaluation. While 0.9924 is obtained as the array factor, when 1.0λ spacing is considered as the widest distance between the elements. On the other hand, the numerical value of directivity increases as the inter-element spacing increases.
Figure 4. Graph of array factor against (a) directivity (b).
Factor against Directivity and inter-element spacing of the Table 1 above. The two graphs obtained show that as the inter-element spacing of the matrix array is increased, the Array Factor decreases and the Directivity of the radiated pattern of the field increases.
7. Radiation Patterns of Varying Element Spacing
Figure 5 below are the simulation results obtained both in 2D plots and polar plots of Array Factor against the elevation (azimuthal). It is the various radiation patterns obtained by varying the inter-element spacing along and axis of the matrix array as and are varied, while the number of elements, N and progressive phase shift, theta/phi = 45, number of elements N = M = 4 are kept constant.
ax=ay=0.47λ ax=ay=0.50λ ax=ay=0.55λ ax=ay=0.60λ ax=ay=0.65λ ax=ay=0.70λ ax=ay=0.75λ ax=ay=0.80λ ax=ay=0.85λ ax=ay=0.90λ ax=ay=0.95λ ax=ay=1.0λ
Figure 5. Radiation patterns of varying element spacing.
8. Discussion of Results
Figure 2 above shows the various simulated radiation patterns obtained by varying the inter-element spacing of the 4 × 4 matrix array. While the values of progressive phase shift and the number of elements N, is kept constant throughout the variation.
The results have a great effect on the array factor as it is clearly shown by radiation patterns of Figure 2 above. Element spacing of 0.47λ to 0.55λ has broad radiation pattern and occupied large bandwidth, which decrease with increase in the spacing. But from element spacing of 0.60λ, the radiation patterns shows an emergence of grating lobe whose size and level increases with the increase in and respectively. When element spacing equals 0.70λ, the size of the grating lobe equals that of the main lobe in opposite direction. Further increase to 0.75λ shows a beginning in the split of the grating lobe. Complete splitting in the grating lobe is eminence when is increased to 0.85λ, 0.90λ, 0.95λ and 1.0λ. Increasing the element spacing towards l results in an increased directivity and grating lobe effect with a maximum grating lobe amplitude equal to the main lobe magnitude at an element spacing of 1.0λ, see Figure 2.
1) As inter-element spacing increases, the grating lobes increase in size, number and levels.
2) Increase in inter-element spacing also increases the directivity of the main lobe.
3) Increasing the spacing to 1.0λ, is a bad element spacing as more grating lobes are developed equal in size and level with the main lobe.
4) Therefore, the best directivity (radiation pattern) can only be obtained when the element spacing is within 0.1 - 0.5λ.
In this paper, a low profile planar array of monopole antenna has been presented. The purpose was to design a low level grating lobe with high directivity using planar array of monopole antenna by varying the inter-element spacing within the array matrix. A 4 × 4 matrix array of monopole antennas has been used in order to enhance the directivity and directivity as high as 39.67 dB has been achieved. The results of Table 1 and Figure 1 above show that directivity of the radiated field decreases as the inter-elements spacing is increased.
The idea of carrying out the present investigation originated from Prof. S. A. Adekole, FAS, who was the author’s M. Eng. Supervisor and Engr. S. Ojomu who is always ready with his mentorship.