República Bolivariana de Venezuela
Universidad de Carabobo
Facultad de Ingeniería
Antenas y Radiación
Antena Yagi Optimizada
Profesor: Ing. Rafael Panza Estudiante:
Br. Sequera Leynner
CIV: 21.531.203
Bárbula, 06 de diciembre de 2020
Bases Teóricas
Antena Tipo Yagi Uda
Básicamente la antena tipo Yagi Uda es un arreglo de elementos formado por: un
elemento activo (excitado), por un reflector y uno o más elementos directores (elementos
parásitos) que le permite ser direccional, tal como se muestra en la siguiente figura:
Figura 1. Antena Tipo Yagi convencional para TV.
Fue diseñada en 1926 por Hidetsugu Yagi y por Shintaro Uda. Esta antena fue un
avance para la época. Gracias a que su estructura le permite ser direccional la lleva a ser la
más utilizada en recepción de TV. Su estructura básica se forma por los siguientes
elementos:
El elemento excitado es un dipolo que tiene una longitud de 0.5λ.
El reflector tiene una longitud un poco mayor (en un 5%) 0.55λ
Los directores son ligeramente más cortos (en un 5%).
La separación entre los elementos es de 0.15λ a 0.25λ, aunque depende del fabricante.
La corriente que circula por el dipolo alimentado irradia un campo electromagnético,
el cual induce corrientes en los elementos parásitos (reflector y director) de la antena.
Provocando a su vez que estos últimos generen inducción electromagnética en los demás
elementos. Por lo que, la corriente que circula en cada uno de los elementos es el resultado
de la interacción entre todos los elementos.
La fase de la corriente que circula en el elemento parásito dependerá de la distancia
con el dipolo alimentado, así como de la longitud y diámetro de este último. La amplitud
también dependerá de los mismos factores, pero será de la misma magnitud que la corriente
del elemento alimentado.
Polarización
Si la antena Yagi es paralela al plano de la tierra, la componente eléctrica de la onda es
paralela al plano de la tierra: se dice que tiene polarización horizontal.
Si la antena Yagi es perpendicular al plano de la tierra, la componente eléctrica de la onda
es perpendicular al plano de la tierra: se dice que tiene polarización vertical.
En HF, y en VHF en clase de emisión banda lateral única se prefiere la polarización
horizontal, y en VHF en clase de emisión frecuencia modulada, la polarización vertical.
Por razones de dimensiones de las antenas, las redes de antenas Yagi se utilizan mucho en
VHF y UHF.
En nuestro caso de estudio se trabajará a la frecuencia de diseño de 600 MHz.
Por lo que estaremos en el rango de UHF.
A diferencia de la antena dipolo, es sumamente difícil modelizar con ecuaciones
matemáticas una antena Yagi. Por lo tanto, existen programas de simulación numérica de
antenas que permiten obtener una aproximación. En nuestro caso nos apoyaremos con los
softwares Matlab y EZNEC+.
El diseño más básico de una antena Yagi-Uda es aquella formada por un arreglo de
dipolos.
Diseño de Antena Yagi-Uda Optimizada:
Para el diseño e implementación de la Antena Yagi-Uda se basó en un arreglo de 6
elementos, conformado por: Un (1) dipolo reflector, un (1) dipolo excitado o alimentado y
cuatro (4) dipolos directores.
Cabe destacar que al agregar más dipolos directores la directividad de la antena
incrementa.
Los parámetros a tomar en cuenta para poder lograr optimizar un arreglo de antenas
Yagi pueden ser: longitud de separación entre los dipolos directores, longitud de cada
elemento, así como el diámetro de estos.
Como se adelantó para obtener los valores de diseño de nuestro arreglo nos apoyamos
del software de cómputo numérico Matlab, la cual nos brinda las herramientas necesarias
para modelar la antena deseada, más en específico con su “Antenna Toolbox” el cual realiza
sus cálculos basado en el método de los momentos (MoM).
Por otro lado, haremos uso del programa EZNEC+, un software basado en el Código
Electromagnético Numérico (NEC) con el cual podremos hacer el estudio de las
características de la antena diseñada en un rango de frecuencias específico.
Simulación de la optimización
Las condiciones planteadas para el diseño de la antena Yagi son los siguientes:
Frecuencia de Operación: 600Mhz.
Cantidad de Directores: 4.
Cantidad de Dipolos (energizado): 1.
Reflector: 1.
Frecuencia de estudio: 650 Mhz – 750 Mhz.
Diámetro de los dipolos: 19 cm.
Velocidad de la luz o espacio libre: 299.792458x106 m⁄s
Para el caso del modelo de antena Yagi-Uda (sin optimizar) se tomaron las siguientes
características:
Longitud del dipolo activo (m) 0.4963
Longitud del dipolo reflector (m) 0.5000
Tamaño de los directores (m) 0.2498
Distancia entre dipolos (m) 0.1249
Distancia entre reflector y dipolo 0.1499
energizado (m)
Estableciendo los parámetros mencionados obtenemos la representación de la antena
Yagi.
Figura 2. Diseño de Antena Yagi (sin optimizar) en Matlab.
Teniendo en cuenta esto, procedemos a simular la antena para obtener las gráficas de
diagrama de radiación tanto en 3D como en 2D.
Parámetros de diseño:
Comenzamos el proceso de diseño de la antena Yagi-Uda en Matlab usando la
técnica optimización de Antenna Toolbox.
Se Elige la frecuencia de diseño en el centro de la banda UHF:
freq = 600e6; 600 Mhz
wirediameter = 19e-3; 19 mm
c = 2.99792458e8; 299.79x𝟏𝟎𝟔 𝐦⁄𝐬
lambda = c/freq; longitud de onda.
Crear antena Yagi-Uda
El elemento accionado para la antena Yagi-Uda es un dipolo plegado. Esto es un
excitador estándar para dicha antena debido a la mayor entrada requerida de impedancia.
Se debe ajustar los parámetros de longitud y anchura del dipolo plegado dado que modelamos
estructuras cilíndricas como elementos metálicos, el ancho se calcula utilizando una función
de utilidad disponible en la Antena Toolbox (TM).
La longitud se elige para ser 0.5λ en el diseño de frecuencia.
d = dipoleFolded;
[Link] = lambda/2;
[Link] = cylinder2strip(wirediameter/2);
[Link] = 0.005;
[Link] = [Link]/60;
Se crea una antena Yagi-Uda con un dipolo doblado como elemento excitado. La
longitud del reflector y del director será 0.5λ.
La distancia del reflector respecto al director es de 0.3λ y el espaciado entre los
directores es de 0.25λ.
Estas opciones son una suposición inicial y servirán como punto de partida para
del procedimiento de optimización.
Numdirs = 4;
refLength = 0.5;
dirLength = 0.5*ones(1,Numdirs);
refSpacing = 0.3;
dirSpacing = 0.25*ones(1,Numdirs);
initialdesign = [dirLength refSpacing dirSpacing].*lambda;
yagidesign = yagiUda;
[Link] = d;
[Link] = Numdirs;
[Link] = refLength*lambda;
[Link] = dirLength.*lambda;
[Link] = refSpacing*lambda;
[Link] = dirSpacing*lambda;
show(yagidesign)
Se traza el patrón de radiación a la frecuencia de diseño antes de ejecutar el proceso
de optimización en 3D.
fig1 = figure;
pattern(yagidesign,freq);
fig7 = figure;
pattern(yagidesign,freq,0,[Link],'Type','efield','Normalize',true);
fig8 = figure;
pattern(yagidesign,freq,90,[Link],'Type','efield','Normalize',true);
Como se esperaba, esta antena no tiene una mayor directividad en la dirección que
deseamos, es decir, en el cenit (elevación = 90 grados). En general el diseño inicial de la
antena Yagi-Uda es un radiador mal diseñado como se muestra a continuación:
Figura 3. Diagrama de Radiación 3D en Matlab.
Figura 4. Diagrama de Radiación 2D (Elevación)en Matlab.
Figura 5. Diagrama de Radiación 2D (Azimutal)en Matlab.
Donde se obtiene que la ganancia de la Antena Yagi- Uda sin optimizar es de 4.49 dBi.
Configuración de los parámetros de optimización:
Los criterios para la optimización serán las longitudes de los elementos parásitos y el
de espacios entre elementos.
Establecemos límites inferior y superior de la longitud del reflector y longitudes del
director, así como los espacios entre ellos en términos de lambda (λ).
Elegimos un sector de 60 grados alrededor del cenit para calcular la directividad, con
el objetivo final de maximizar lo dirección de radiación en este sector en comparación con
los lóbulos laterales y el lóbulo posterior.
dirLengthBounds = [0.4 0.4 0.4 0.4; 0.495 0.495 0.495 0.495];
% límite inferior en la longitud del director
% límite superior en la longitud del director
refSpacingBounds = [0.05;0.3]
% límite inferior en el espacio del reflector
% límite superior en el espacio del reflector
dirSpacingBounds = [0.05 0.05 0.05 0.05; 0.23 0.23 0.23 0.23] ;
% límite inferior en el espacio de los directores
% límite superior en el espacio de los directores
LB = [dirLengthBounds(1,:) refSpacingBounds(1) dirSpacingBounds(1,:) ].*lambda;
UB = [dirLengthBounds(2,:) refSpacingBounds(2) dirSpacingBounds(2,:) ].*lambda;
parasitic_values = [ [Link], [Link], [Link]];
elang = [60 120]; % ángulos de ancho de haz de elevación en az = 90
Optimización de búsqueda directa
Global Optimization Toolbox (TM) proporciona una búsqueda directa basada en la
Función de optimización llamada: patternsearch.
Usamos esta función con dos opciones especificadas las cuales son:
psoptimset function.
En cada iteración, se traza el mejor valor de la función objetivo y se limita el número
total de iteraciones a 30.
La función objetivo tiene como fin maximizar el de diferencia entre la directividad
máxima del haz principal y el máximo de lóbulos laterales y el lóbulo posterior.
Se pasa la función objetivo a función patternsearch mediante el uso de una función
anónima junto con el límite y la estructura de opciones. La función objetivo utilizada durante
el proceso de optimización por búsqueda de patrones está disponible en el archivo
yagi_objective_function de Matlab.
Se evalúa la directividad en diferentes direcciones correspondientes a la región
angular definida para radiación máxima, así como el máximo del lóbulo lateral y el nivel
del lóbulo trasero a través de la función calcular_objetivos disponibles dentro de
yagi_objective function.
patternoptions = psoptimset(@patternsearch);
[Link] = @psplotbestf;
[Link] = 35;
optimdesign = patternsearch(@(x)
yagi_objective_function(yagidesign,x,freq,elang),parasitic_values,[],[],[],[],LB,UB,[],patternoptions);
Traza patrón optimizado. El patrón de antena optimizado ahora se traza a la frecuencia de
diseño.
[Link] = optimdesign(1:4);
[Link] = optimdesign(5);
[Link] = optimdesign(6:9);
fig6 = figure;
show(yagidesign)
fig2 = figure;
pattern(yagidesign,freq)
Cortes del patrón del plano E y H.
Para obtener una mejor visión del comportamiento en los dos planos ortogtonales,
graficamos la magnitud normalizada del campo eléctrico en los planos E y H, es decir:
azimut = 0 y 90 grados respectivamente.
fig3 = figure;
pattern(yagidesign,freq,0,[Link],'Type','efield','Normalize',true);
fig4 = figure;
pattern(yagidesign,freq,90,[Link],'Type','efield','Normalize',true);
El diseño optimizado muestra una mejora significativa en el patrón de radiación con
una mayor directividad lograda en la dirección deseada hacia el cenit. El lóbulo posterior es
pequeño, lo que da como resultado un buen frente hacia atrás de relación para esta antena.
Calculemos la directividad en el cenit, la relación de frente-atrás y el ancho del haz en los
planos E y H, respectivamente.
D_max = pattern(yagidesign,freq,0,90);
D_back = pattern(yagidesign,freq,0,-90);
F_B_ratio = D_max - D_back;
Eplane_beamwidth = beamwidth(yagidesign,freq,0,[Link]);
Hplane_beamwidth = beamwidth(yagidesign,freq,90,[Link]);
Tabulación de diseño inicial y optimizado
Las conjeturas de diseño iniciales y los valores de diseño optimizados finales son que se
muestra en la tabla a continuación:
yagiparam= {'Director Length - 1'; 'Director Length - 2';
'Director Length - 3'; 'Director Length - 4';
'Reflector Spacing'; 'Director Spacing - 1';
'Director Spacing - 2';'Director Spacing - 3';
'Director Spacing - 4'};
initialdesign = initialdesign';
optimdesign = optimdesign';
T = table(initialdesign,optimdesign,'RowNames',yagiparam)
Realizando la simulación obtendremos los valores de la antena optimizada,
redimensionando así las longitudes de los dipolos y las distancias entre los mismos. La
representación de la antena optimizada se muestra en la siguiente imagen.
Figura 6. Diseño de Antena Yagi (Optimizada) en Matlab.
Figura 7. Diagrama de Radiación 3D en Matlab.
Figura 8. Diagrama de Radiación (Azimutal) Optimizado 2D en Matlab.
Figura 9. Diagrama de Radiación (elevación) Optimizado 2D en Matlab.
Donde se obtiene que la ganancia de la Antena Yagi- Uda optimizada de 6.04 dBi
mejorando respecto a los 4.49 dBi calculado anteriomente, lo que representa un 35% de
mejoría en ganancia.
A su vez Matlab nos arroja los valores de los parámetros definidos a modificar en la
optimización siendo estos:
Director Length - 1 0.24983 0.20827
Director Length - 2 0.24983 0.24733
Director Length - 3 0.24983 0.23952
Director Length - 4 0.24983 0.21608
Reflector Spacing 0.1499 0.13427
Director Spacing - 1 0.12491 0.11492
Director Spacing - 2 0.12491 0.11492
Director Spacing - 3 0.12491 0.11492
Director Spacing - 4 0.12491 0.068045
Las medidas de los dipolos y distancias entre ellos se muestran en la siguiente tabla.
Reflector (m) 0,50000
Dipolo Energizado (m) 0.4963
Director 1 (m) 0.20827
Director 2 (m) 0.24733
Director 3 (m) 0.23952
Director 4 (m) 0.21608
Distancia reflector (m) 0.13427
Distancia director 1 (m) 0.11492
Distancia director 2 (m) 0.11492
Distancia director 3 (m) 0.11492
Distancia director 4 (m) 0.068045
De esta forma nos queda por representar es el estudio en el rango de frecuencias
especificado anteriormente (de 650 a 750 Mhz) a pasos de 10 Mhz para lo cual se usará el
software EZNEC+, replicando así las dimensiones de la antena descritas en la tabla anterior
Figura 10. Dimensiones de los cables de la antena optimizada en EZNEC+.
Figura 11. Antena Yagi-Uda Optimizada modelada en EZNEC.
Obtenemos los respectivos diagramas de radiación y SWR para el rango de
frecuencias, esto con lectura de datos cada 10 Mhz.
Figura 12. SWR de la Antena Yagi-Uda Optimizada en EZNEC+.
De este estudio obtenemos que el máximo SWR es de 15,20 a la frecuencia de 580 MHz
con una pérdida de retorno de 1,10 dB y un coeficiente de reflexión de 0,875.
Figura 13. SWR de la Antena Yagi-Uda Optimizada en EZNEC+.
A la frecuencia de diseño de 600 MHz se obtiene que el SWR es de 12,30 con una pérdida
de retorno de 1,40 dB y un coeficiente de reflexión de 0,849.
El grafico del SWR nos permite medir el rendimiento de la antena en el rango de 550 MHz
a 650 MHz donde se concluye que en el rango de frecuencia de 570 Mhz a 620 Mhz la
Antena Yagi-Uda optimizada obtiene un buen funcionamiento.
Para comprobar la eficiencia de la optimización se realiza un estudio del SWR en la misma
ventana de frecuencia de la Antena Yagi-Uda sin optimizar donde se obtuvo lo siguiente:
Figura 14. SWR de la Antena Yagi-Uda sin optimizar en EZNEC+.
En ella se observa un comportamiento constante del SWR en el rango correspondiente, sin
embargo, podemos percatar que el valor del SWR varía entre los 6,80 a 5,60 siendo estos
valores muy por debajo del máximo obtenido al optimizar.
Así como también una mayor eficiencia de la antena optimizada en los rangos de frecuencia
próximos a la frecuencia de diseño.
En este mismo orden de ideas, se comprueba los diagramas de radiación tanto en elevación
como azimut de la Antena tanto optimizada como sin optimizar para comparar el
rendimiento de la antena obteniendo:
Figura 15. Diagrama de Radiación (elevación) sin optimizar 2D en EZNEC+.
Figura 16. Diagrama de Radiación (elevación) optimizado 2D en EZNEC+.
Figura 17. Diagrama de Radiación (azimut) sin optimizar 2D en EZNEC+.
Figura 18. Diagrama de Radiación (azimut) optimizado 2D en EZNEC+.
Podemos realizar en EZNEC+ la comparación donde se obtiene que la ganancia máxima de
la antena Yagi-Uda sin optimizar es de 4.45 dBi y en la optimizada de 6.22 dBi y obtener
una mejora del 39% siendo semejante a lo obtenido en Matlab. A su vez se puede observar
la relación frente-atrás en ambos casos siendo estos valores de 3.31 dB y 3.25 dB donde de
nuevo a una leve mejora. La energía irradiada en los lóbulos laterales disminuye y se
aprovecha en gran medida hacia la dirección de interés.
Con los resultados obtenidos se concluye que al optimizar la antena Yagi-Uda se obtiene una
antena más directiva aumentando también de forma discreta la relación frente atrás de la
antena lo cual permitirá ser más resistente a señales que emanen de la parte trasera.
Bibliografía
Antenna Theory Analysis and Design Second Edition, Constatine A. Balanis.
Trabajo Practico del Laboratorio de Antenas, Rafael J. Albornoz A.
Diseño de Practicas de laboratorio de Antenas utilizando NEC, César Ruiz.
Programming and Engineering Computing with Matlab 2019, Huei-Huang Lee.
