This example calculates the performance of the two linearly-polarized rectangular patch antennas with results published in . The first antenna has a low-epsilon thin dielectric substrate while the second antenna has a thick high-epsilon dielectric substrate.
Create a rectangular patch antenna with a length of 40 mm and a width of 30 mm on a 80 mm x 60 mm ground plane. The lossless substrate has a dielectric constant of 2.2 and thickness of 1.57mm. The feed is offset by 5.5 mm from the origin along the x-axis.
p1 = patchMicrostrip; p1.Length = 40e-3; p1.Width = 30e-3; p1.Height = 1.57e-3; p1.GroundPlaneLength = 80e-3; p1.GroundPlaneWidth = 60e-3; p1.FeedOffset = [5.5e-3 0];
Antenna Toolbox™ has a list of substrates supported as part of it's dielectric catalog. To open the catalog use the following commands.
The substrate specified for this patch is not part of the dielectric catalog. You can add it to the catalog if desired, so that the material can be used next time without specifying its electrical properties. You can also specify the substrate property directly as well as shown below.
p1.Substrate = dielectric('Name','material1','EpsilonR',2.33); figure show(p1);
The impedance of the antenna shows a resonance at 2.37 GHz. This value is very close to the results published in the paper.
The antenna pattern at resonance also shows uniform illumination above the ground plane with minimum leakage below it.
figure; pattern(p1, 2.37e9);
Below is the mesh used to model the antenna. The triangles are used to discretize the metal regions of the patch, and terahedra are used to discrertize the volume of the dielectric substrate in the patch. These are indicated by the colors yellow and green respectively. The total number of unknowns is the sum of the unknowns for the metal plus the unknowns used for the dielectric. As a result, the solution slows down significantly as compared to pure metal antennas.
As seen in the figure above, only a single layer of tetrahedra are used to model the dielectric substrate. This gives good results when the substrate thickness is small as in this case. However for thicker substrates this might be a limitation. Manual meshing may be required to to get better accuracy. This is explained in the following section.
Thickness of a dielectric substrate is measured with respect to the wavelength. In the case above the wavelength in free-space is about 126 mm. The wavelength in dielectric is approximated by dividing the above number by the square root of the dielectric constant. This value is about 85mm. So the substrate thickness is sbout 1/50th of the wavelength in dielectric. This is a thin substrate.
The next case considers a patch antenna on a substrate whose thickness is 1/10th of the wavelength in dielectric. This is a thick substrate. A method of moments solver requires at least 10 elements per wavelength to give an accurate solution. So if the dielectric thickness becomes greater than this, the solution accuracy is affected. In Antenna Toolbox if the substrate thickness is greater than 1/10th of the wavelength in dielectric, manual meshing is recommended for solving the antenna. As the 10 element per wavelength criteria is not satsified, certain amount of error in the solution is expected.
Create a rectangular patch antenna with a length of 36 mm and a width of 48 mm on a 55 mm x 80 mm ground plane. The lossless substrate has a dielectric constant of 9.29 and thickness of 3.82mm. The feed is offset by 4 mm from the origin along the x-axis.
p2 = patchMicrostrip; p2.Length = 36e-3; p2.Width = 48e-3; p2.Height = 3.82e-3; p2.GroundPlaneLength = 55e-3; p2.GroundPlaneWidth = 80e-3; p2.FeedOffset = [4.0e-3 0]; p2.Substrate = dielectric('Name','material2','EpsilonR',9.29); figure; show(p2);
This is perhaps the most complicated case from the numerical point of view - a thick high-epsilon dielectric with significant fringing fields close to the patch edges. The figure below shows the impedance plot of the antenna with resonance close to 1.22 GHz. The paper shows a resonance close to 1.27 GHz.
Below is the mesh used for calculating the antenna performance.
Better results can be achieved by refining the mesh. The mesh can be refined using the maximum edge length criteria. In the case below the maximum edge length is set to 2mm. As can be seen below, 12000 tetrahedra are generated. As the mesh is very fine, the results take a longer time to compute and hence are saved in a .mat file.
Plotting the impedance shows the resonance close to 1.27 GHz as expected.
load thickpatch figure; plot(freq*1e-9, real(Z), 'r', freq*1e-9, imag(Z), 'b', 'LineWidth',2); legend('resistance', 'reactance'); title('Impedance'); ylabel('Impedance (ohms)'); xlabel('Frequency (GHz)'); grid on;
The directivity information for this patch is also precomputed and stored. This can be plotted using the patternCustom function as shown below. A significant back lobe is observed for the thick patch antenna.
figure; patternCustom(D.', 90-el, az); h = title('Directivity (dBi)'); h.Position = [-0.4179, -0.4179, 1.05];
 S. N. Makarov, S. D. Kulkarni, A. G. Marut and L.C. Kempel 'Method of Moments Solution for a Printed Patch/Slot Antenna on a Thin Finite Deielectric Sustrate Uisng the Volume Integral Equation', IEEE Trans. on Antenna and Propagation, vol. 54, No. 4, April 2006, pp 1174- 1184.