MATLAB Implementation of Yagi-Uda Antenna Using Circuit Theory Approach
MATLAB Implementation of Yagi-Uda Antenna Using Circuit Theory Approach
MATLAB Implementation of Yagi-Uda Antenna Using Circuit Theory Approach
T.Anushalalitha
Assistant Professor, Department of Telecommunication Engineering
BMS College of Engineering, Bangalore, Karnataka, India
Abstract - In this paper, the radiation pattern of yagi-uda antenna is simulated in MATLAB using circuit theory
approach. Our basic approach was to simulate the radiation pattern for a symmetrically shaped antenna and then
maximizing the output parameters by using various techniques such as using reflector surfaces wherever the loss in
antenna was due to side lobes. We have designed this antenna using circuit theory and adopting Pocklingtons integral
equation to calculate the performance given a set of parameters and adjust them to optimize the gain. and have made
improvements in the previous designs to have better electric field intensity and directivity.
Keywords Yagi-Uda, Pocklington Integral equation, MATLAB
I. INTRODUCTION
Yagi-Uda antennas are directional along the axis perpendicular to the dipole in the plane of the elements, from the
reflector toward the driven element and the director(s). Typical spacing's between elements vary from about 1/10 to
1/4 of a wavelength, depending on the specific design. The lengths of the directors are smaller than that of the driven
element, which is smaller than that of the reflector(s). These elements are usually parallel in one plane, supported on
a single crossbar known as a boom.
II. DESIGN
Kirchhoffs laws establish relations between current and voltage in electric networks where electrical parameters are
considered concentrated in certain points in meshes. They lead to Circuit Theory, which is an engineering area that
models electric networks accurately enough while the mesh dimensions are less than wavelength. However, as the
dimensions become greater than or equal to the wavelength, the more electrical parameters are distributed through
the network. In this way, Circuit Theory should be modified in order to model accurately distributed phenomena and
current-voltage relations.
2.1 Finite Difference Time Domain (FDTD) method
The Finite Difference Time Domain (FDTD) method is a general method to solve Maxwells partial differential
equations numerically in the time domain. It has been used extensively to model all kinds of electromagnetic
problems, such as radiation, scattering, and circuit problems. However, due to the limitation of computer resources,
it is difficult to apply the FDTD method in the modeling of electrically large objects. At the same time, the
implementation of the FDTD method requires absorbing boundary conditions to terminate the spatial grid, which
adds an additional burden on computational resources. Moreover, the recursive nature of the classic Yee algorithm
for the FDTD method implies that all the cells contained in a given volume must be computed. Therefore, when
applied to model scattering problems, the computation of free space nodes between scatterers is required even
though this is actually unnecessary. As a result, the FDTD method is characteristically time consuming and memory
demanding.
64
(2)
Figure 1. (a) Design steps (b) Flow chart showing the design and measurements
A weighted method - Gaussian method which states In numerical analysis, a quadrature rule is an approximation of
the definite integral of a function, usually stated as a weighted sum of function values at specified points within the
domain of integration is used for integration.
65
PARAMETERS
3-dB BEAMWIDTH IN THE E-PLANE PATTERN
3-dB BEAMWIDTH IN THE H-PLANE PATTERN
FRONT-TO-BACK RATIO IN THE E-PLANE
FRONT-TO-BACK RATIO IN THE H-PLANE
DIRECTIVITY
CASE 1
CASE 2
CASE 3
60.01
DEGREES
78.05
DEGREES
13.8937 dB
61.77
DEGREES
82.26
DEGREES
32.30
DEGREES
34.16
DEGREES
15.4882 dB
8.1523 dB
13.8876 dB
8.994 dB
15.4827 dB
8.599 dB
8.1366 dB
12.782 dB
PARAMETERS
NUMBER OF MODES
PER ELEMENT
NUMBER OF
ELEMENTS
LENGTH OF DIRECTOR
LENGTH OF
REFLECTOR
LENGTH OF DRIVEN
ELEMENT
SEPARATION
BETWEEN DRIVEN
ELEMENT & 1ST
DIRECTOR
SEPARATION
BETWEEN REFLECTOR
& DRIVEN ELEMENT
RADIUS FOR ALL
ELEMENTS USED
CASE 1
CASE 2
CASE 3
10
#
#
#
#
#
#
#
#
1
2
3
4
5
6
7
8
0.36000
0.34500
0.33000
0.31500
0.30000
0.28500
0.27000
0.37500
0.43900
WAVELENGTHS
0.43900
WAVELENGTHS
0.50200
WAVELENGTHS
0.50200
WAVELENGTHS
0.56250 WAVELENGTHS
0.45200
WAVELENGTHS
0.45200
WAVELENGTHS
0.41250 WAVELENGTHS
0.23800
WAVELENGTHS
0.23800
WAVELENGTHS
0.24000 WAVELENGTHS
0.13800
WAVELENGTHS
0.13800
WAVELENGTHS
2.00000e-003
WAVELENGTHS
1.00000e-003
WAVELENGTHS
0.24000 WAVELENGTHS=
SEPARATION BETWEEN DIRECTORS
3.75000e-002 WAVELENGTHS
Consequently, these antennas are often empirical designs using an element of trial and error, often starting with an
existing design modified according to one's hunch. The result might be checked by direct measurement or by
66
computer simulation. A well-known reference employed in the latter approach is a report published by the National
Bureau of Standards (NBS) (now the National Institute of Standards and Technology (NIST)) that provides six basic
designs derived from measurements conducted at 400 MHz and procedures for adapting these designs to other
frequencies. These designs, and those derived from them, are sometimes referred to as "NBS yagis." There are so
many ways to optimize the directivity and other antenna parameters such as the most primitive trial and error
method and the simple genetic algorithm method.
Table -3 Program output plots for the Yagi-Uda Array
CASE 1
PARAMETERS
CASE 2
CASE 3
Linear plot
Current distribution
Polar plot
Approx.
anticipated
gain(dB over
dipole)
5
7.5
8.5
9.5
IV.CONCLUSION
The gain of a Yagi antenna is governed mainly by the number of elements in the particular RF antenna. However the
spacing between the elements also has an effect. As the overall performance of the RF antenna has so many inter-
67
related variables, many early designs were not able to realize their full performance. Today computer programs are
used to optimize RF antenna designs before they are even manufactured and as a result the performance of antennas
has been improved.
V. FUTURE ENHANCEMENTS
The front to back ratio is important in circumstances where interference or coverage in the reverse direction needs to
be minimized. Unfortunately the conditions within the antenna mean that optimization has to be undertaken for
either Front to back ratio, or maximum forward gain. Conditions for both features do not coincide, but the front to
back ratio can normally be maximized for a small degradation of the forward gain. For Yagi-uda antenna, the results
of Matlab implementation and theoretically calculated results are within the reasonable error limits. Here also, the
final radiation possesses minimum radiation intensity in the back lobe (at angle of 180 degrees) and a maximum of
lobe power is concentrated in the main lobe along the axis of the antenna. The directivity of the antenna first
increases with frequency and then decreases. The simulations, that of Matlab, proved the above mentioned result.
These MATLAB implementations can also be expanded for other types of antennas and results are obtained.
Calculating the parameters of an antenna required for the work becomes cumbersome and prone to error if done
manually or in any other method so this MATLAB procedure finds its application.
VI. ACKNOWLEDGMENT
We thank LRDE,DRDO, govt of india for providing us an opportunity to perform this project.
REFERENCES
[1]
[2]
[3]
[4]
[5]
[6]
[7]
[8]
[9]
[10]
[11]
68