Informe - Final Antenas - Completo

UNIVERSIDAD POLITÉCNICA SALESIANA ANTENAS
Fecha: 21/06/2019 TELECOMUNICACIONES
PATRÓN DE RADIACIÓN DE LA ANTENA DIPOLO HERTZIANO,

MAGNÉTICO Y DIPOLO DOBLADO EN MATLAB
Carrillo Lagla Katerin Adriana
e-mail: [email protected]
Cuyago Loyola Michelle Estefanía
González Flores Fausto Alejandro
RESUMEN: En el siguiente informe se presenta el El patrón de radiación de una antena se puede

análisis y simulación del patrón de radiación para las representar gráficamente, presentando la distribución
Antenas Dipolo Hertziano, Dipolo Magnético y Dipolo relativa de potencia o propiedades de radiación en
Doblado mediante la herramienta matemática de Matlab. función de coordenadas espaciales. En la mayoría de
Además, se realizó el cálculo de campos para zonas casos, el patrón de radiación es evaluado considerando
apartadas de cada antena. condiciones en la región de campo lejano, y
representado en coordenadas direccionales.
PALABRAS CLAVE: Dipolo Hertziano, magnético,
dipolo doblado, Patrón de radiación.
A la traza de la potencia recibida a un radio constante se
le denomina patrón de potencia. Un sistema coordenado
1. OBJETIVOS
esférico es empleado en muchas aplicaciones prácticas
para desplegar la traza del patrón en 3 dimensiones. En
1.1 OBJETIVO GENERAL la Figura 1 se presenta un bosquejo de la geometría
 Comprender el comportamiento de las ondas esférica asociada. El plano de elevación se determina en
electromagnéticas, su propagación y patrones función de la variable θ, y el plano azimutal está en
de radiación para las antenas de dipolo función de la variable φ. La representación gráfica de la
Hertziano, Magnético y Doblado, mediante una variación espacial de los campos eléctrico o magnético,
simulación realizada en Matlab. con radio constante, se le denomina patrón de campo.
[1]
1.2 OBJETIVO ESPECÍFICOS
 Implementar un código en Matlab para

representar el patrón de radiación de cada
antena.
 Obtener los parámetros principales de cada
antena mediante el uso de sus respectivas
fórmulas.
2. MARCO TEÓRICO
2.1 PATRÓN DE RADIACIÓN Figura 1. Sistema de coordenadas polares empleado

para presentar la traza tridimensional del patrón de
El patrón de radiación de una antena se puede potencia de una antena. [1]
representar como una gráfica tridimensional de la
energía radiada vista desde fuera de esta. Los patrones
El patrón de radiación para un dipolo ideal es
de radiación usualmente se representan de dos formas,
representado por la función senθ, asociada con la
el patrón de elevación y el patrón de azimut. El patrón de
variación del campo eléctrico de radiación e, y también
elevación es una gráfica de la energía radiada por la
presenta características de radiación para el campo
antena vista de perfil. El patrón de azimut es una gráfica
magnético H. El plano que contiene el vector campo
de la energía radiada vista directamente desde arriba. Al
eléctrico de radiación se le denomina Plano-E; el plano
combinar ambas gráficas se tiene una representación
tridimensional de cómo es realmente radiada la energía perpendicular a éste se denomina Plano-H, el cual
desde la antena. [1] contiene el vector de campo magnético. [1]
1
2.2 ANTENA DIPOLO HERTZIANO

El campo eléctrico radiado está polarizado linealmente,
Una antena simple para la cual se pueden calcular los ^ , si se encierra al
la orientación del vector es según θ
campos de una manera directa no complicada, es el dipolo en una esfera, suponiendo orientación, z el
dipolo Hertziano o dipolo eléctrico elemental. Esta campo seguiría los meridianos. El plano E es el definido
antena ideal consiste de un elemento infinitesimal de por la dirección de máxima radiación y el campo
corriente (infinitesimal respecto de la longitud eléctrica), eléctrico en dicha dirección. [2]
de longitud “dl”, que transporta un fasor de corriente Î el
mismo que se asume constante en magnitud y fase a lo El campo eléctrico en zonas apartadas para este tipo de
largo de toda la longitud del segmento como se muestra estructura es el siguiente:
en la figura 2. [2]
^
E FF= j
Idl
4 πr
no β o sinθ e
−j β r ^ v
θ
o
m [ ]
( Ec .2)
Por lo tanto, el plano E es cualquier plano que contiene

al eje z, por ejemplo, el XZ o el YZ. En este caso la
definición del plano E n o es única, ya que el problema
tiene simetría de revolución. Cualquier corte φ = cte se
podría definir como plano E. Los cortes en la esfera
serían planos que incluyesen a los dos polos. [2]
2.2.1.2 Intensidad de Campo magnético ^

H
Los campos magnéticos generados por las corrientes y
que se calculan por la ley de Ampere o la ley de Biot-
Figura 2. Antena dipolo eléctrico elemental [2] Savart, se caracterizan por el campo magnético B
medido en Teslas. Pero cuando los campos generados
Así entonces, debido a que se conoce la distribución de
pasan a través de materiales magnéticos que por sí
corriente, el vector Potencial magnético en un punto mismo contribuyen con sus campos magnéticos
ubicado a una distancia radial r del origen de internos, surgen ambigüedades sobre que parte del
coordenadas y del dipolo, puede determinarse como: campo proviene de las corrientes externas, y que parte
la proporciona el material en sí. [2]
μ0 − jβ r ❑
^
A=
4π
e 0
∫ Î dz k^ (Ec .1) El plano H es el definido por la dirección de máxima
l
radiación y el campo magnético en dicha dirección. Para
Donde: el dipolo elemental es el plano XY. En la esfera sería el
plano ecuatorial.
μ0 : Es la permeabilidad magnética
Por tanto, la intensidad de campo magnético se describe
β 0 : Es la constante de fase de la siguiente manera:
r:
Î :
Es el radio
Es la corriente
^
H FF= j
Idl
β sinθ e
4 πr o
− jβ r
φ^
o
A
m [ ]
( Ec .3)
2.2.1.3 Vector de densidad de media potencia
2.2.1 CAMPOS EN ZONAS APARTADAS Se denomina vector de Poynting al vector cuyo módulo

representa la intensidad instantánea de energía
Debido a que el interés del estudio de antenas radica electromagnética que fluye a través de una unidad de
principalmente en el conocimiento de los campos en área perpendicular a la dirección de propagación de la
zonas apartadas (transmisión de información a onda electromagnética, y cuyo sentido es el de
distancia), las expresiones para E y H pueden propagación. Recibe su nombre del físico inglés John
simplificarse notablemente si se desprecian las
Henry Poynting. Se expresa mediante el símbolo S AV .
contribuciones del campo de inducción (proporcional a
1/r2) y del campo electrostático (proporcional a 1/r3). [3] [4]
^ El vector de Poynting en esta región lejana es radial y en

2.2.1.1 Intensidad de Campo eléctrico E el sentido de la propagación de la onda, se representa
2
mediante la siguiente ecuación:
^
S AV =
no I 2 dl 2 sen 2 θ W
32 π 2 r 2 [ ]
γ^ 2 ( Ec .4)
m
2.2.1.4 Potencia Radiada
La potencia total radiada por el dipolo se puede obtener

integrando la densidad de potencia radiada en una
superficie que encierre al dipolo elemental. [4] Figura 3. Patrón de radiación tridimensional de la antena
( )
2
dl
2 2 dipolo Hertziano. [2]
Prad =40 π I [ W ] ( Ec .5)
λ0 Se observa entonces, que la máxima radiación ocurre
para un ángulo θ = 90°, mientras que radiación cero para
cualquier punto ubicado sobre el eje Z. Las gráficas del
2.2.1.5 Resistencia de radiación patrón de radiación en los planos E y H dan patrones
como los que se indica en la figura 4, donde el patrón en
Representa la potencia o energía que la antena radia el plano E es direccional con máximos en θ = 90°,
hacia el espacio libre. Se define Rr como el valor de la siendo el mismo para cualquier ángulo φ, y el patrón en
resistencia que disiparía la misma potencia que la el plano H omnidirecional.
radiada por la antena. [4]
La resistencia de radiación se puede calcular a partir de

la potencia total radiada y la corriente en la entrada,
obteniendo así la siguiente ecuación:
Rr=80 π
2
( )
dl 2
λO
[ Ω ] (Ec .6)
2.2.1.6 Ganancia Directiva
La relación entre la potencia de un radiador isotrópico Figura 4. (a) Patrón de radicación plano E. (b) Patrón de
que produce la misma intensidad de campo eléctrico en radicación plano H. [2]
una dada dirección y la potencia que irradia la antena en
cuestión para lograr el mismo objetivo. [4] 2.3 ANTENA DIPOLO MAGNÉTICO
Se representa mediante la siguiente ecuación: Dipolo magnético elemental es un lazo conductor de
radio “a”, en el que la longitud totaldel lazo es
4π ^
2
S AV r sumamente pequeña comparada con la longitud de
D= ( Ec .7) onda, y a través de la cualcircula una corriente I que se
Prad asume es igual en magnitud y fase a lo largo de todo el
lazo(ver Figura 5). [2]
2.2.2 PATRÓN DE RADIACIÓN – DIPOLO
HERTZIANO
La gráfica del módulo de la densidad media de potencia

para valores constantes de r se conoce como el patrón
de radiación. Para este caso se obtiene un patrón
tridimensional como se indica en la figura 3. [2]
Figura 5. Dipolo magnético elemental. [2]
3
[ ]
uesto que se conoce la corriente, el vector potencial no I 2 a 4 β0 4 sen2 θ W
magnético en un punto situado a unadistancia radial r ^
S AV = γ^ 2 ( Ec .10)
del centro del lazo será: 32r 2 m
μ 0 I^ e− j β R
0
^
A=
4π
∮ R dl (Ec . 8) 2.3.1.4 Potencia Radiada
La potencia total radiada se puede obtener integrando la

Donde:
densidad de potencia radiada en una superficie que
encierre al dipolo elemental, obteniendo la siguiente
μ0 : Es la permeabilidad magnética ecuación:
β 0 : Es la constante de fase
( )
2
2 πa
Prad =15585 I 2
[ W ] ( Ec .11)
R: Es el radio λ0
Î : Es la corriente 2.3.1.5 Resistencia de radiación
Aunque lo habitual es definir el dipolo magnético como La resistencia de radiación se puede calcular a partir de
una pequeña espira o distribución de corriente, la potencia total radiada y la corriente en la entrada,
realmente lo que lo define es el campo que produce. obteniendo así la siguiente ecuación:
( )
Una partícula elemental, como el electrón, produce un
πa 2
campo magnético dipolar y por tanto es un dipolo Rr=31170 2
[ Ω ] ( Ec .12)
magnético aunque no sea una corriente eléctrica. [3] λ0
2.3.1 CAMPOS EN ZONAS APARTADAS
2.3.1.1 Intensidad de Campo eléctrico ^

E 2.3.2 PATRÓN DE RADIACION – DIPOLO
MAGNÉTICO
El campo eléctrico en zonas apartadas para este tipo de
estructura es el siguiente: La gráfica para valores constantes de r da una forma del
patrón de radiación tridimensional, similar al caso del
n o Î a β 0 dl
[ ]
2 2 dipolo Hertziano, como se observa en la figura 6. [2]
^ v
E FF= sinθ e− j β r φ^
0
( Ec .9)
4r m
Por lo tanto, el plano E es cualquier plano que contiene
al eje z, por ejemplo, el XZ o el YZ. En este caso la
definición del plano E n o es única, ya que el problema
tiene simetría de revolución. Cualquier corte φ = cte se
podría definir como plano E. Los cortes en la esfera
serían planos que incluyesen a los dos polos. [4]
2.3.1.2 Intensidad de Campo magnético ^

H
La intensidad de campo magnético para una antena
dipolo magnético elemental se describe de la siguiente
manera:
− Î a β 0 dl
[ ]
2 2
−j β r ^ A
^
H FF= sinθ e 0
θ ( Ec .9)
4r m
Figura 6. Patrón de radiación tridimensional del dipolo
magnético elemental. [2]
2.3.1.3 Vector de densidad de media potencia
Los cortes del patrón de radiación en los planos E y H,
El vector de Poynting en esta región lejana es radial y en dan un patrón de radiación omnidireccional en el plano E
el sentido de la propagación de la onda, se representa y direccional en el plano H como se observa en la figura
mediante la siguiente ecuación:
4
7, Esto es, contrario de lo que ocurría con el dipolo de operación útil, para esta antena es mayor que para la
eléctrico. antena dipolo convencional. [2]
2.4.2 FUNCIONAMIENTO
El dipolo doblado se puede descomponer en el modo

par o modo antena, con la misma alimentación en los
dos brazos, y el modo impar o modo línea de
transmisión, con dos generadores con signos opuestos,
como se muestra en la figura 9. [2]
Figura 7. (a) Patrón de radicación plano E. (b) Patrón de

radicación plano H. [2]
2.4 ANTENA DIPOLO DOBLADO
2.4.1 DEFINICIÓN
Un dipolo doblado es una estructura formada por dos

dipolos paralelos, cortocircuitados en su extremo. Uno
de ellos se alimenta en el centro con un generador,
Figura 9. Dipolo doblado; geometría y modos de
como se muestra en la figura 8. [2]
funcionamiento. [2]
Puesto, que los conductores están separados por una
2.4.2.1 Modo impar o línea de transmisión
pequeñísima fracción de longitud de onda, hay una
diferencia de fase despreciable en el campo radiado
Para el modo línea de transmisión de la figura 3(b), la
desde cada conductor. Consecuentemente:
impedancia de entrada vista entre los terminales a y b o
 Se dobla campo entre e y f corresponde a la de una línea de transmisión
 Se cuadruplica la resistencia de radiación ≈ 300 terminada en corto. [2]
Ω
El modo impar equivale a dos líneas de transmisión en
cortocircuito, alimentadas en serie. Aunque no
contribuye a la radiación, sí afecta la impedancia a la
entrada. [2]
Como la diferencia de potencial entre los puntos a y b es

de V/2, y se aplica a una línea de transmisión de
longitud l/2, la corriente que fluye por la línea es:
V 1
I t= ∙ ( Ec .12)
2 Z¿
Donde Yin es la admitancia de entrada del segmento de
línea de transmisión terminada en corto circuito y por lo
Figura 8. Antena de Dipolo Doblado. [2] tanto está dada por:
La antena dipolo doblado por construcción tiene una

línea de transmisión equivalente que actúa como stub de
Y ¿=
1
Z¿
=− jY 0 cot
2( )
βl
(Ec .13)
sintonía compensando variaciones de impedancia de la
antena con la frecuencia. Así la banda de frecuencias
5
Siendo Yo la admitancia característica de la línea de

transmisión formada por dos conductores paralelos, la
misma que determinarse como:
Z 0=120 cosh
−1
( d
)
2 √ a 1 ∙ a2
(Ec .14)
It quedará:
I t=− j
V
2
∙ Y 0 cot
βl
2 ( )
(Ec .15)
2.4.2.2 Modo par o antena
En el modo antena que se representa en la figura 3(c),

los puntos c − d y g − h están cada uno de ellos situados
al mismo potencial y pueden conectarse entre sí para
formar un dipolo. Cada uno de los brazos del dipolo está
formado por un par de hilos muy juntos (s << λ) que van
desde la alimentación (c − d o g − h) hasta el corto. De
Figura 10. Dipolo equivalente. [3]
esta manera, la corriente en el modo antena viene dada
por: Siendo Yd la admitancia de entrada de una antena dipolo
V /2 construida con dos conductores paralelos conectados en
I a= ( Ec .16) los extremos y en el centro, donde la estructura dipolo
Zd
equivalente tiene un radio que puede determinarse con
donde Zd es la impedancia de entrada de un dipolo de
la siguiente relación:
longitud l y diámetro d. Para la configuración de la figura
3(c), el radio que se debe usar para calcular la
1
impedancia Zd puede ser la semidistancia entre los hilos ln a eq=ln a1 + 2
( u2 lnu+ 2u ln v ) ( Ec .20)
s/2 o un radio equivalente aeq. [2] (1+u )
La corriente de entrada puede ser encontrada por la De donde:
adición de las corrientes en el conductor principal debido
a las dos formas de excitación. [3] a2 d
u= y u=
En la figura 10 se muestra la antena dipolo equivalente a1 a1
para la cual la corriente 2I1, está dada por: Donde:
V 1 a eq :Radio Equivalente
2 I 1= ∙ (Ec .17)
2 Zd a 1 , a2 : Radios
O d : Distancia de separaci ó n
V 2.4.2.3 Admitancia de la antena
2 I 1= ∙ Y ( Ec .18)
2 d
De donde: Las corrientes totales son la suma de las corrientes en
los dos modos, descrita en la siguiente ecuación:
V
I 1= ∙Y ( Ec .19) I =I par + I impar [ A ] ( Ec .21)
4 d
Y la admitancia vista en los terminales de entrada
de la antena estará dada por:
I I1 I2
Y ant = = + ( Ec .22)
V V V
Obteniendo:
6
mismo contribuyen con sus campos magnéticos

internos, surgen ambigüedades sobre que parte del
( )
Yd Y0 βl
Y ant = − j cot ( Ec .23) campo proviene de las corrientes externas, y que parte
4 2 2 la proporciona el material en sí. [2]
El plano H es el definido por la dirección de máxima

 Debe notarse que la admitancia de la antena radiación y el campo magnético en dicha dirección. Para
dipolo (simple) se reduce en un factor de 4. el dipolo elemental es el plano XY. En la esfera sería el
plano ecuatorial [2].
 Y una admitancia de compensación es
Por tanto, la intensidad de campo magnético al igual que
adicionada en paralelo. Cuando l=λ/2, la
el campo eléctrico es el doble del dipolo largo y se
admitancia de compensación se hace 0 puesto
describe de la siguiente manera:
que βl/2=π/2.[3]
Debido a que el interés del estudio de antenas radica

^
H FF= j
Im
πr
∙ F (θ ) ∙ e o
[ ]
− jβ r A
m
( Ec .26)
principalmente en el conocimiento de los campos en

zonas apartadas (transmisión de información a 2.4.3.3 Potencia Radiada
distancia), las expresiones para E y H pueden
simplificarse notablemente si se desprecian las Puesto que la corriente aplicada por la línea de
contribuciones del campo de inducción (proporcional a transmisión es únicamente Im, la resistencia de
1/r2) y del campo electrostático (proporcional a 1/r3). [3] radiación referida a los terminales de entrada de la
antena es incrementada en un factor de4 sobre aquella
de la antena dipolo convencional. Esto es la potencia
radiada será: [3]
2.4.3.1 Intensidad de Campo eléctrico ^
E
2
El campo eléctrico radiado está polarizado linealmente, Prad =4∗36,5|I m| [ W ] (Ec .27)
^ , si se encierra al
la orientación del vector es según θ
dipolo en una esfera, suponiendo orientación, z el
campo seguiría los meridianos. El plano E es el definido 2.4.3.4 Resistencia de radiación
por la dirección de máxima radiación y el campo
eléctrico en dicha dirección. [2] Representa la potencia o energía que la antena radia
hacia el espacio libre. Se define Rr como el valor de la
El campo eléctrico para este tipo de estructura es el resistencia que disiparía la misma potencia que la
doble del campo eléctrico del dipolo largo y se describe radiada por la antena. [4]
de la siguiente manera:
La resistencia de radiación se puede calcular a partir de
n .I
[ ]
la potencia total radiada y la corriente en la entrada,
^ v
E FF= j 0 m ∙ F (θ ) ∙ e− j β r
o
(Ec .24) obteniendo así la siguiente ecuación:
πr m
Donde:
Rrad =8∗36,5=292,5 [ Ω ] ( Ec .28)
F (θ)=
(l
2 ) ( )
cos β 0 . .cosθ −cos β 0 .
l
2
(Ec .25)
2.4.4 PATRÓN DE RADIACIÓN – DIPOLO
DOBLADO
senθ El patrón de radiación de la antena dipolo doblado es
exactamente igual que para el caso de la antena dipolo
^ simple, al igual que su polarización, mientras que el
2.4.3.2 Intensidad de Campo magnético H ancho de banda puede incrementarse hasta en un 5%
de su frecuencia central, como se indica en la figura 11.
Los campos magnéticos generados por las corrientes y [3]
que se calculan por la ley de Ampere o la ley de Biot-
Savart, se caracterizan por el campo magnético B
medido en Teslas. Pero cuando los campos generados
pasan a través de materiales magnéticos que por sí
7
^
E FF= j
Idl
4 πr
no B o sin e− j B R θ^
v
m
o
( )
 Campo Magnético
^
H FF= j
Idl
4 πr
sin e− j B R φ^
o
A
m ( )
Figura 11. Patrón de radiación tridimensional de la
antena dipolo doblado. [3]
4.1.2 VECTOR DENSIDAD MEDIA DE POTENCIA
Se observa entonces, que la máxima radiación ocurre
para un ángulo θ = 90°, mientras que radiación cero
 SAV
para cualquier punto ubicado sobre el eje Z. Las gráficas
del patrón de radiación en los planos E y H dan patrones
como los que se indica en la figura 12. [3] ^
SAV ❑ = j
Idl
4 πr
sin e
− jB R
φ^
v
m
o
( )
4.1.3 POTENCIA, RESISTENCIA Y GANANCIA
 Potencia radiada
no l 2 dl 2 B o2
Pr= (w )
12 π
Figura 12. (a) Patrón de radicación plano E. (b) Patrón

de radicación plano H. [2]  Resistencia de radiación
( )
2
2 dl ( )
3. MATERIALES Y EQUIPO Rr=80 π Ω
λO
 PC
 Software Matlab
 Materia Clase  Ganancia directiva
4. DESARROLLO Y PROCEDIMIENTO D=4 πSAV r

2
4.1 PARÁMETROS – ANTENA DIPOLO

HERTZIANO.
Para el desarrollo de los cálculos y por ende de la 4.2 PARÁMETROS – ANTENA DIPOLO
emisión de las gráficas respectivas es esencial conocer MAGNÉTICO
e investigar que formulas son necesarias para ser
aplicadas dependiendo de cada antena. FÓRMULAS:

FÓRMULAS:
4.1.1 CAMPOS EN ZONAS APARTADAS  Campo Eléctrico
 Campo Eléctrico ^
E FF= j
Idl
4 πr
no B o sin e
−jB R ^ v
θ
m
o
( )
8
 Campo Magnético
^ ℑ
H FF= j F(θ)e
πr
−j B r A
m
o
( )
^
H FF= j
Idl
4 πr
sin e− j B R φ^
o
A
m ( )  F (θ)
4.2.2 VECTOR DENSIDAD MEDIA DE POTENCIA
 SAV

^
SAV ❑ = j
Idl
4 πr
sin e− j B R φ^
v
m
o
( )  Yin
4.2.3 POTENCIA, RESISTENCIA Y GANANCIA
 Potencia radiada
 Yant
no l 2 dl 2 B o2
Pr= (w )
12 π
 Resistencia de radiación
( )
2 4.3.3 RADIO, DISTANCIA, Aeq, Zo.
2 dl ( )
Rr=80 π Ω
λO  Radio del conductor
a2
u=
 Ganancia directiva a1
2
D=4 πSAV r
 Distancia
d
v=
a1
4.3 PARÁMETROS – ANTENA DIPOLO
DOBLADO
FÓRMULAS:
4.3.1 CAMPOS EN ZONAS APARTADAS  A equivalente
 Campo Eléctrico 1
ln aeq=ln a 1+( )(u2 ln u+2 ulnv
1+u 2
^
E FF= j
no ℑ
πr
F (θ)e
− jB r v
m
o
( ) )
 Zo
 Campo Magnético −1
Zo ¿120 cosh (
9
d
¿( Ω )
2 √ a 1+a 2
4.4 PARÁMETROS ADICIONALES GENERALES

PARA TODAS LAS ANTENAS
 Longitud de onda
Figura8. Software Matlab - New Script

C
λ=
f 4.5.3 Se realiza una caratula de presentación como
inicio del programa.
 Radio
r =√ x + y + z ( m o km )
2 2 2
 Ángulos
()
z F
θ=cos−1 igura9. Script – Caratula.
r
4.5.4 Para el ingreso de datos se programa un Menú
 no donde se permita ingresar manualmente los
datos por parte del usuario y además da la
n o=120 π opción de escoger entre que antena se desea
obtener los parámetros.
MENU:
 Bo 1. Antena dipolo Hertziano
2. Antena dipolo Magnético
2π 3. Antena dipolo Doblado
Bo=
λ
4.5 SOFTWARE MATLAB Figura10. Menú de opciones de Antena.
Es necesario recurrir a un software que permita el 4.5.5 Al inicio de ejecutar el programa se visualiza la
cálculo de varios parámetros matemáticos, físicos y presentación y el menú de opciones que
además que permita visualizar las gráficas respetivas en permite escoger una de las dos antenas.
2D Y 3D.
PASOS A SEGUIR:
4.5.1 Se ingresa al Software Matlab.
4.5.2 Buscar la pestaña de New Script para abrir un

nuevo proyecto y empezar a desarrollar el
código.
Figura11. Diseño del Menú de opciones de Antena.
4.5.6 El usuario deberá escoger entra la opción 1 u

opción 2.
10
-OPCIÓN 1: Si se selecciona la opción 1 se despliega el

ingreso de datos de todos los parámetros a calcular de
la antena Dipolo Hertziano. 4.6.3 CÁLCULO DE PARAMETROS
-OPCIÓN 2: Si se selecciona la opción 2 se despliega el Se debe seguir un orden de cálculo para la resolución de
ingreso de datos de todos los parámetros a calcular de las formulas dadas y llegar a encontrar todos los
la antena Dipolo Magnético. parámetros, entonces el código está desarrollado en
base a lo siguiente:
-OPCIÓN 3: Si se selecciona la opción 3 se despliega el
ingreso de datos de todos los parámetros a calcular de PASOS:
la antena Dipolo Doblado.
Primero se observa los datos con los que se cuenta para
4.5.7 Una vez que se selecciona una de las dos la resolución del ejercicio.
antenas se despliega el menú de preguntas,
pidiendo al usuario que ingrese los valores de 1. Para el cálculo de los parámetros principales
los datos para calcular sus parámetros. de las antenas es necesario calcular los
parámetros generales como son:
4.6 CODIGO EN MATLAB – ANTENA DIPOLO 1.1 El Bo se calcula con el valor de la frecuencia y la
HERTZIANO Y MAGNETICO constante C para el valor de lambda.
4.6.1 PRESENTACIÓN.
CÁLCULO DE Bo
CARÁTULA DE PRESENTACIÓN disp ( '--------- ------------- ------------ ')

clc disp('Cálculos de Parametros')
clear
close all c=3*10^8;
disp ( '----------------------------------------------------------- ') landa=c/f1;
disp ( 'UNIVERSIDAD POLITECNICA SALESIANA ' ) bo=2*pi/landa;
disp ( 'Integrantes: ' ) no=120*pi;
disp ( '*CARILLO KAYHERIN ' ) disp(['bo=',num2str(bo),' [rad/m]'])
disp ( '*CUYAGO MICHELLE ' )
disp ( '*GONZÁLEZ ALEJANDRO ' )
disp ( 'TEMA: CALCULOS PARA UNA ANTENA
DIPOLO HERTZIANO, DIPOLO ELECTRICO Y 1.2 El Radio se calculó con el valor del punto de
DIPOLO DOBLADO ' ) coordenadas dado.
CÁLCULO DEL RADIO

4.6.2 Se define las variables de entrada con el
comando disp e input que en este caso viene a r=sqrt(x^2+y^2+z^2);
ser los datos del ejercicio, estos datos los da el r1=r*1000;
usuario. disp(['r=',num2str(r),' [km]'])
INGRESO DE DATOS 1.3 Los ángulos se calculan con respecto al punto de

disp ( '----------------------------------------------------------- ') coordenadas dado y realizando una estimación en un
disp ( '--------- INGRESO DE DATOS ---------- ') gráfico de la posición de los ángulos.
dl=input ( 'Ingrese el valor dl en [cm]: ');
f=input( 'Ingrese el valor Frecuencia en [MHZ]: ');
I=input ( 'Ingrese el valor de la Corriente en [uA]: ');
dl1=dl/100; CÁLCULO DE ÁNGULOS
f1=f*10^6;
I1=I*10^-6; ar=z/r;
disp(['dl:',num2str(dl1),' [m]']) angulot=acos (ar);
disp(['f1:',num2str(f1),' [HZ]']) grad=(angulot*180)/pi;
disp(['I:',num2str(I1),' [A]']) sen=sin(angulot);
disp('Ingrese el valor del punto P(x,y,z):') sen2=sen^2;
x=input ( 'Ingrese el valor de x en [km]: ');
y=input ( 'Ingrese el valor de y en [km]: ');
z=input ( 'Ingrese el valor de z en [km]: ');
2. Cálculo de los parámetros principales de la

antena dipolo hertziano y magnético.
11
2.1 Una vez obtenido los valores de los parámetros

generales se reemplaza en la fórmula del campo disp('***** Resistencia de Radiacion*****')
eléctrico. Rrad=80*(pi^2)*(dl1/landa)^2;
disp(['Rrad:',num2str(Rrad),' [ohm]'])
CÁLCULO DE CAMPO ELÉCTRICO
ar=z/r;
angulot=acos (ar);
grad=(angulot*180)/pi; 2.6 Una vez obtenido los valores de los parámetros
sen=sin(angulot); generales se reemplaza en la fórmula de ganancia
sen2=sen^2; directiva.
CALCULO GANANCIA DIRECTIVA

generales se reemplaza en la fórmula del campo disp('***** Ganancia Directiva *****')
magnético. Di=1.5*sen2;
disp(['Di:',num2str(Di),' []'])
CÁLCULO DE CAMPO MAGNÉTICO
disp('***** Campo Eléctrico*****') 3. Para la elaboración de los planos h y plano e

disp(['Eff:',num2str(eff),'*e^(',num2str(E1),' ) [v/m]']) para la visualización del patrón de radiación en
disp('***** Campo Magnetico*****') 3D.
hff=j*((I1*dl1)/(4*pi*r1))*bo*sen;
disp(['Hff:',num2str(hff),'*e^(',num2str(E1),' ) [A/m]'])
CÁLCULO GRÁFICA PLANO E

generales se reemplaza en la fórmula de densidad de figure(4)
potencia. plot(Y,Z); % Gráfica plano E
title('PLANO E')
xlabel('Y');
CÁLCULO DE DENSIDAD DE POTENCIA ylabel('Z');
%figure(4) %FIGURA 2
%GRAFICAR CAMPO MAGNETICO EN
disp('***** Densidad de Potencia*****') COORDENADAS POLARES
sav=(no*(I1^2)*(dl1^2)*sen2*(bo^2))/(32*(pi^2)*(r1^2)); polar(TETA,E);
disp(['SAV:',num2str(sav),' [w/m]']) hold on
polar(TETA,-E);
hold on

generales se reemplaza en la fórmula de potencia de
radiación.
CÁLCULO GRÁFICA PLANO H
CÁLCULO POTENCIA DE RADIACIÓN

figure(3)
plot(X,Y);
disp('***** Potencia de Radiacion*****') title('PLANO H') % Gráfica plano H
Prad=10*(I1^2)*(dl1^2)*(bo^2); xlabel('X');
disp(['Prad:',num2str(Prad),' [w]']) ylabel('Y');
end
2.5 Una vez obtenido los valores de los parámetros %% dipolo magnetico elemental
generales se reemplaza en la fórmula de resistencia de %%ejercicio
radiación. %3. Una antena dipolo magnético tiene un radio de
1cm,
%trabaja a una frecuencia de 6Mhz,
%y lo atraviesa una corriente de 15mA.
CALCULO RESISTENCIA DE RADIACION
%Determinar para el punto (10,5,15) Km:
12
A1=input ( 'Ingrese el valor del diámetro del conductor

1 [m]: ');
VER (ANEXO 1)
while (A1<=0)
disp('')
disp('-----ingrese valores positivos----')
1 [m]: ');
if (A1>=1)

2 [m]: ');
4.7 CÓDIGO EN MATLAB – ANTENA DIPOLO while (A2<=0)

DOBLADO disp('')
4.7.1 PRESENTACIÓN. A2=input ( 'Ingrese el valor del diámetro del conductor
2 [m]: ');
if (A2>=1)
CARÁTULA DE PRESENTACIÓN D=input ( 'Ingrese el valor de la distancia de

clc separacion [m]: ');
clear
close all
disp ( '----------------------------------------------------------- ') while (D<=0)
disp ( 'UNIVERSIDAD POLITECNICA SALESIANA ' ) disp('')
disp ( 'Integrantes: ' ) disp('-----ingrese valores positivos----')
disp ( '*CARILLO KAYHERIN ' ) D=input ( 'Ingrese el valor de la distancia de
disp ( '*CUYAGO MICHELLE ' ) separacion [m]: ');
disp ( '*GONZÁLEZ ALEJANDRO ' ) if (D>=1)
disp ( 'TEMA: CALCULOS PARA UNA ANTENA
DIPOLO HERTZIANO, DIPOLO ELECTRICO Y disp('Ingrese el valor del punto P(x,y,z):')
DIPOLO DOBLADO ' ) x=input ( 'Ingrese el valor de x en [m]: ');
y=input ( 'Ingrese el valor de y en [m]: ');
z=input ( 'Ingrese el valor de z en [m]: ');
4.7.2 Se define las variables de entrada con el

comando disp e input que en este caso vienes
a ser los datos del ejercicio, estos datos los da 4.7.3 CÁLCULO DE PARÁMETROS
el usuario.
Se debe seguir un orden de cálculo para la resolución de
las formulas dadas y llegar a encontrar todos los
INGRESO DE DATOS parámetros, entonces el código está desarrollado en
base a lo siguiente:
disp('----INGRESO DE DATOS ----') PASOS:
l=input ( 'Ingrese el valor dE Longitud (l) en [m]: ');
while (l<=0) Primero se observa los datos con los que se cuenta para
disp('') la resolución del ejercicio.
l=input ( 'Ingrese el valor dE Longitud (l) en [m]: '); 1. Para el cálculo de los parámetros principales
if (l>=1)
de la antena es necesario calcular los
parámetros generales como son:
f1=input( 'Ingrese el valor Frecuencia en [MHZ]: ');
while (f1<=0)
1.1 El Bo se calcula con el valor de la frecuencia y la
disp('')
disp('-----ingrese valores positivos----') constante C para el valor de lambda.
if (f1>=1)
CÁLCULO DE Bo
I1=input ( 'Ingrese el valor de la Corriente en [mA]: ');
while (I1<=0)
disp('')
if (I1>=1)
13
disp ( '--------- ------------- ------------ ') %cálculo de F(o)

disp('Cálculos de Parametros') F=(arg1-arg2)/sin(angulot);
E1= ((Bo)*(r));
c=3*10^8; EFF=((no*I)/(pi*r))*F
landa=c/f1; HFF=(I/(pi*r))*F
bo=2*pi/landa;
no=120*pi;
disp(['bo=',num2str(bo),' [rad/m]']) 2.2 Una vez obtenido los valores de los parámetros
generales se reemplaza en la fórmula del campo
eléctrico.
1.2 El Radio se calculó con el valor del punto de
coordenadas dado.
CÁLCULO DE CAMPO ELÉCTRICO
CÁLCULO DEL RADIO ar=z/r;

angulot=acos (ar);
r=sqrt(x^2+y^2+z^2); grad=(angulot*180)/pi;
r1=r*1000; sen=sin(angulot);
disp(['r=',num2str(r),' [km]']) sen2=sen^2;
1.3 Los ángulos se calculan con respecto al punto de 2.3 Una vez obtenido los valores de los parámetros
coordenadas dado y realizando una estimación en un generales se reemplaza en la fórmula del campo
gráfico de la posición de los ángulos. magnético.
CÁLCULO DE ÁNGULOS
CÁLCULO DE CAMPO MAGNÉTICO
ar=z/r;
angulot=acos (ar);
grad=(angulot*180)/pi; disp('***** Campo Eléctrico*****')
sen=sin(angulot); disp(['Eff:',num2str(eff),'*e^(',num2str(E1),' ) [v/m]'])
sen2=sen^2; disp('***** Campo Magnetico*****')
1.4 El valor de u y v son esenciales para los cálculos

próximos de una antena dipolo doblado, se basan en los
2.4 Cuando ya se tiene los valores del radio a1 y a2 se
radios del conductor a1 y a2 como también en la
distancia d. puede reemplazar en la fórmula de Aeq.
CÁLCULO DE U Y V CÁLCULO DE Aeq
%%RADIO EQUIVALENTE Aeq=log(a1)+((1)/(1+u)^2)*((u)^2*log(u)+2*u*log(v));

aeq=exp(Aeq)
u=a2/a1;
v=d/a1;
2.5 Es necesario realizar la INTERPOLACION una vez
que se tenga los valores de Bo para obtener los
2. Cálculo de los parámetros principales de la
antena dipolo doblado. parámetros de R(θ ) y X(θ ).
2.1 Para el cálculo de los campos eléctrico y magnético

se requiere de un parámetro que debe ser calculado con CALCULO DE LA INTERPOLACION
anterioridad denominado F(θ ).
%%INTERPOLACION
B=(pi/landa)*long;
if(B>=0.1 && B<0.2)
CÁLCULO DE F(θ ) RB=0.1506+((RB-0.1)*(0.7980-0.1506))/(0.1) ;
XB=1.010+((RB-0.1)*(2.3002-1.010))/(0.1);
14
end CÁLCULO DE Zdip y Zo

if(B>=0.2 && B<0.3)
RB=0.7980+((RB-0.2)*(1.821-0.7980))/(0.1) ; Z1=log(long/(2*aeq))-1
XB=2.302+((RB-0.2)*(3.818-2.302))/(0.1); cotag=cot((B*long)/2)
end zdip=RB-j*(120*Z1*cotag-XB);
if(B>=0.3 && B<0.4 ) ydip=((1)/(zdip))
RB=1.821+(((B-0.3)*(3.264-1.821))/(0.1));
XB=3.818+(((B-0.3)*(5.584-3.818))/(0.1));
end zo=120*acosh(d/(2*sqrt(a1*a2)));
if(B>=0.4 && B<0.5 ) yo=1/zo
RB=3.264+(((B-0.4)*(5.171-3.264))/(0.1)); vc=cot((B*long)/2)
XB=5.584+(((B-0.4)*(7.141-5.584))/(0.1));
end
if(B>=0.5 && B<0.6 )
RB=5.171+(((B-0.5)*(7.563-5.171))/(0.1)); 2.7 Una vez calculados todos los parámetros se finaliza
XB=7.141+(((B-0.5)*(8.829-7.141))/(0.1));
end reemplazando en la fórmula de Yant para así obtener el
if(B>=0.6 && B<0.7 ) valor de Zant que es lo que se requiere.
RB=7.563+(((Bl2-0.6)*(10.48-7.563))/(0.1));
XB=8.829+(((Bl2-0.6)*(10.68-8.829))/(0.1));
end CÁLCULO DE Zant
if(B>=0.7 && B<0.8 )
RB=10.48+(((Bl2-0.7)*(13.99-10.48))/(0.1));
XB=10.68+(((Bl2-0.7)*(12.73-10.68))/(0.1)); Yant=(ydip/4)-i*((yo/2)*vc)
end Zant=1/Yant
if(B>=0.8 && B<0.9 )
RB=13.99+(((B-0.8)*(18.16-13.99))/(0.1));
XB=12.73+(((B-0.8)*(15.01-12.73))/(0.1));
end 3. Para la elaboración de los planos h y plano e
if(B>=0.9 && B<1 ) para la visualización del patrón de radiación en
RB=18.16+(((B-0.9)*(23.07-18.16))/(0.1)) 3D.
XB=15.01+(((B-0.9)*(17.59-15.01))/(0.1))
end
if(B>=1 && B<1.1 )
RB=23.07+(((B-1)*(28.83-23.07))/(0.1))
XB=17.59+(((B-1)*(20.54-17.59))/(0.1))
end
if(B>=1.1 && B<1.2 )
RB=28.83+(((B-1.1)*(35.60-28.83))/(0.1)) DISEÑO DE GRÁFICA
XB=20.54+(((B-1.1)*(23.93-20.54))/(0.1))
end %%grafica
if(B>=1.2 && B<1.3 ) bo=(2*pi)/landa %beta
RB=35.6+(((B-1.2)*(43.55-35.6))/(0.1));
XB=23.93+(((B-1.2)*(27.88-23.93))/(0.1));
end %gráfica
if(B>=1.3 && B<1.4 ) fi=(0:.01:1)*2*pi; %fi
RB=43.55+(((B-1.3)*(52.92-43.55))/(0.1)) teta=(0:.01:1)*2*pi; %theta
XB=27.88+(((B-1.3)*(32.2-27.88))/(0.1))
end %fórmulas
if(B>=1.4 && B<1.5 ) teta1=acos(z/r)
RB=52.92+(((B-1.4)*(64.01-52.92))/(0.1)) disp ( 'Graficas en 3D, PlanoH y Plano E del dipolo
XB=32.2+(((B-1.4)*(38-32.2))/(0.1)) largo' )disp ( 'El Campo Eléctrico, Magnético, la
end impedancia y la resistencia de pérdida' )
if(B>=1.5 && B<(pi/2) ) disp ( 'Presione 1:' )
RB=64.01+(((B-1.5)*(73.12-64.01))/((pi/2)-1.5)) g=input('=> ' );
XB=38+(((B-1.5)*(42.46-38))/(0.1)) if(g==1)
end
if(B==(pi/2) )
RB=73.12
XB=42.46
End
CALCULO GRÁFICA PLANO E
plot(Y,Z); % Gráfica plano E

title('PLANO E')
2.6 Con los valores de Aeq y los que vota la xlabel('Y');
interpolación se puede encontrar el valor de Zip para ylabel('Z');
realizar un reemplazo y obtener Zo. subplot(2,2,3)
title('PLANO E')
15
polar(TETA,E) % Gráfica plano E Campo Eléctrico (EFF) 4.0276e-

subplot(2,2,2) 11i*e^(0-
335920.3432i )
[v/m]
Campo Magnético (HFF) 1.0683e-
13i*e^(0-
CALCULO GRÁFICA EN PLANO H 335920.3432i )
[A/m]
subplot(2,2,2) SAV (SAV) 2.151e-24[w/m]
Potencia radiada (Pr) 9.935e-13 [w]
plot(X,Y);
Resistencia de Radiación (Rr) 1.9871 [Ω]
title('PLANO H') % Gráfica plano H
xlabel('X'); Ganancia directiva 1.4286
ylabel('Y');
end
ANÁLISIS 1: Cuando se trabaja con corrientes altas, el
dipolo Hertziano tiene un comportamiento de irradiar
VER (ANEXO 2) potencias bajas, es decir es inversamente proporcional
el valor de la corriente con el valor de la Potencia
5. ANÁLISIS Y RESULTADOS Radiada, para el ejercicio se implementó una corriente
de 1uA a una f=70 MHz el cual al momento de realizar
5.1 SISTEMA DE PRUEBAS - ANTENA DIPOLO los cálculos se observa que arroja una Potencia radiada
HERTZIANO de un valor bajo igual a 9.935e-13 [w] por lo que se
verifica que el cálculo es correcto, caso contrario si se
Se realiza la estimación de parámetros con una escala tuviera un valor de Potencia Radiada alto quiere decir
de valores: que está trabajando con un amperaje muy pequeño.
 Prueba 1 - Valores Medios en longitud y
frecuencia.
 Prueba 2 - Valores Medios en longitud y
frecuencia.
5.1.1 PRUEBA 1: Se ingresa los siguientes datos con

valores altos con el fin de obtener los parámetros que
posee una Antena Dipolo Hertziano.
DATOS 1:
 dl= 21.5 cm
 f=70 Mhz
 I=1 uA
Figura 12. Patrón de radicación plano E-Dipolo
 P (200,100,50) km
HERTZIANO 1.
RESULTADOS 1:
Tabla 1: Resultados 1 – Parámetros Básicos Antena Dipolo
Hertziano.
PARÁMETROS BASICOS RESULTADO

Bo 1.4661 [rad/m]
Radio (r) 229.1288 [km]
Angulo Theta Radianes 1.3508 [rad]
Angulo Theta Grados 77.3956 [°]
Tabla 2: Resultados 1 – Parámetros Antena Dipolo

Hertziano.
PARÁMETROS RESULTADO
16
Figura 14. Patrón de radicación plano H-Dipolo

HERTZIANO 1.
Figura 13. Patrón de radicación de 0 a 360°-Dipolo Figura 15. Patrón de radicación 3D - Dipolo
HERTZIANO 1. HERTZIANO 1.
ANÀLISIS 2: El tipo de polarización que presenta esta 5.1.2 PRUEBA 2: Se ingresa los siguientes datos con
antena, está determinado por la posición del vector E valores medios con el fin de obtener los parámetros que
(vector del campo eléctrico) con respecto a una posee una Antena Dipolo Hertziano.
superficie reflectora, la que generalmente es la
superficie terrestre, en este caso se tiene una antena DATOS 2:
con polarización vertical, para conseguir polarización
 dl= 1,5 cm
vertical se debe ubicar el dipolo en posición vertical con
 f=20 Mhz
respecto a la tierra, como se puede observar el circulo
 I=1 uA
amarillo y el circulo azul tienden a aproximarse a su
 P (200,100,50) km
valor de ejes de 90 y 270, es decir que el patrón de
radiación aumenta de forma considerable, este cambio RESULTADOS 2:
se da porque a mayor frecuencia en este caso 70Mhz,
con una longitud de 21,5 cm y una corriente de 1uA, el Tabla 3: Resultados 2 – Parámetros Básicos Antena Dipolo
patrón de radiación crece, en el plano E en el eje de las Hertziano.
Y se puede observar que el patrón de radiación está en
PARÁMETROS BASICOS RESULTADO
un rango de -1.5 a 1.5 y en el plano z está en un rango
de -1 a 1, tiende a tener un equilibrio en sus laterales Bo 0.41 [rad/m]
como en su parte superior e inferior sus valores de Radio (r) 229.1288 [km]
rango definen una antena dipolo hertziano con un ancho
Angulo Theta Radianes 1.3508 [rad]
similar a su amplitud o altura, se puede corroborar en el
plano H donde se define de mejor manera sus rangos. Angulo Theta Grados 77.3956 [°]

Hertziano.
Campo Eléctrico (EFF) 8.0284e-
13i*e^(0-
95977.2409i )
[v/m]
Campo Magnético (HFF) 2.1296e-
15i*e^(0-
95977.2409i )
[A/m]
SAV (SAV) 1.756e-25[w/m]
Potencia radiada (Pr) 8.1106e-14 [w]
17
Resistencia de Radiación (Rr) 0.16221 [Ω] Figura 17. Patrón de radicación de 0 a 360°-Dipolo
Ganancia directiva 1.4286 HERTZIANO 2.
ANALISIS 4: El tipo de polarización que presenta esta

antena, es una polarización vertical como en el caso de
ANALISIS 3: Cuando se trabaja con corrientes altas, el
dipolo Hertziano tiene un comportamiento de irradiar la Prueba 1, como se puede observar el círculo amarillo
potencias bajas, es decir es inversamente proporcional y el círculo azul tienden a alejarse del valor de los ejes
el valor de la corriente con el valor de la Potencia de 90 y 270, es decir que el patrón de radiación
Radiada, para el ejercicio se implementó una corriente disminuye de forma considerable, este cambio se da
de 1uA a una f=20 MHz el cual al momento de realizar porque a menor frecuencia en este caso 20Mhz, con una
los cálculos se observa que arroja una Potencia radiada longitud de 21,5 cm y una corriente de 1uA, el patrón de
de un valor bajo igual a 8.1106e-14 [w] por lo que se radiación decrece, en el plano E en el eje de las Y se
verifica que el cálculo es correcto, caso contrario si se puede observar que el patrón de radiación está en un
tuviera un valor de Potencia Radiada alto quiere decir rango de -5 a 5 y en el eje z está en un rango de -3 a 3,
que está trabajando con un amperaje muy pequeño, en comparando con los rangos del plano E de la Prueba 1
el caso de la Prueba 1 se trabajó con una frecuencia de se tiene que los valores son más extensos y alejados del
70 MHz lo que quiere decir que la frecuencia no influye cero porque la antena dipolo tiende a ensancharse en
el valor de la Potencia Radiada ya que se mantiene en
los laterales y aplanarse en la parte superior e inferior,
un valor bajo lo que es correcto para esa corriente.
da la apariencia como una dona aplastada.
Figura 18. Patrón de radicación plano H-Dipolo

Figura 16. Patrón de radicación plano E-Dipolo HERTZIANO 2.
HERTZIANO 2.
Figura 19. Patrón de radicación en 3D-Dipolo

HERTZIANO 2.
18
5.2 SISTEMAS DE PRUEBAS - ANTENA DIPOLO receptora ya que en la parte donde es implementada se
MAGNETICO trabaja con potencias muy altas y la antena dipolo
magnético necesariamente debe trabajar con corrientes
Se realiza la estimación de parámetros con una escala muy altas o con potencias muy altas, en este caso del
de valores: ejercicio se simula con una corriente alta de 25mA a una
frecuencia de 20Mhz y una longitud de 50 cm, con este
 Prueba 1 - Valores Medios en longitud y valor de corriente se lograr tener un nivel de potencia
frecuencia. radiada bajo igual a 0.00011869 [w], lo que es ideal
 Prueba 2 - Valores Bajos Medios en longitud y para este tipo de antenas.
frecuencia.

valores altos con el fin de obtener los parámetros que
posee una Antena Dipolo Magnético.
DATOS 1:
 dl= 50 cm
 f=20 Mhz
 I=25 mA
 P (45,200,100) km
RESULTADOS 1:
Tabla 5: Resultados 1 – Parámetros Básicos Antena Dipolo

Magnético.
PARÁMETROS BÁSICOS RESULTADO

Bo 0.4188 [rad/m]
Figura 20. Patrón de radicación plano E-Dipolo

MAGNÉTICO 1.

Magnético.
07*e^(0-
95542.0716i )
[v/m]
Campo Magnético (HFF) -1.0803e-
09*e^(0-
95542.0716i )
[A/m]
SAV (SAV) 2.199e-16[w/m]
Potencia radiada (Pr) 0.00011869 [w] Figura 21. Patrón de radicación plano H-Dipolo
Resistencia de Radiación (Rr) 0.3798 [Ω] MAGNÉTICO 1.
ANÀLISIS 6: El tipo de polarización que presenta esta

ANÀLISIS 5: El dipolo Magnético tiene un antena dipolo magnético, es una polarización vertical
comportamiento muy ineficiente como una antena como en el caso de la, como se puede observar el
19
circulo amarillo y el circulo azul mantienen una distancia 05*e^(0-

similar con respecto al valor de los ejes de 90 y 270 y 286626.2148i )
con respecto al punto 0 es decir que el patrón de [v/m]
radiación se encuentra como en un equilibrio ni muy Campo Magnético (HFF) -3.889e-
pequeño ni muy grande, este cambio se da porque a 08*e^(0-
menor frecuencia en este caso 20Mhz, con una longitud 286626.2148i )
de 50 cm y una corriente de 25mA, el patrón de [A/m]
radiación decrece en amplitud pero tiene un SAV (SAV) 2.850e-13[w/m]
decrecimiento moderado lo contrario con la Prueba 2 del
dipolo Hertziano que a pesar de tener la misma Potencia radiada (Pr) 0.15382 [w]
frecuencia tuvo un decrecimiento bastante notable. En el Resistencia de Radiación (Rr) 492.2192 [Ω]
plano E se puede observar que el patrón de radiación
tiene un valor de 0.5 y su corte tiende a estar por la
mitad en referencia a los ejes y en el plano H tiene un ANÀLISIS 7: El dipolo Magnético tiene un
valor de 0.001, la gráfica en 3D muestra como la antena comportamiento muy ineficiente como una antena
dipolo hertziano tiene un decrecimiento moderado no receptora ya que en la parte donde es implementada se
tiene mucho ensanchamiento en sus laterales. trabaja con potencias muy altas y la antena dipolo
magnético necesariamente debe trabajar con corrientes
muy altas o con potencias muy altas, en este caso del
ejercicio se simula con una corriente alta de 25mA a una
frecuencia de 60Mhz y una longitud de 100 cm, con este
valor de corriente se debería lograr tener un nivel de
potencia radiada bajo igual que la Prueba 1 de la antena
dipolo magnético, pero se tuvo un valor de 0.15382 [w],
lo que ya no es un valor tan bajo como se estimaba, esto
se debe a que aquí ya existe un factor muy importante
que influye y es el aumento de la frecuencia a 60 MHz .
Figura 22. Patrón de radicación 3D - Dipolo

MAGNETICO 1.

valores medios con el fin de obtener los parámetros que
posee una Antena Dipolo Hertziano.
DATOS:
 dl= 100 cm
 f=60 Mhz
 I=25 mA
 P (45,200,100) km
RESULTADOS 2:
Tabla 7: Resultados 2 – Parámetros Básicos Antena Dipolo Figura 23. Patrón de radicación plano E - Dipolo
Magnético. MAGNÉTICO 2.
PARÁMETROS BÁSICOS RESULTADO
Bo 1.2566 [rad/m]

Magnético.
PÁRAMETROS RESULTADO
20
Figura 25. Patrón de radicación 3D-Dipolo MAGNÉTICO

2.
5.3 SISTEMA DE PRUEBAS - ANTENA DIPOLO

Figura 24. Patrón de radicación de 0 a 360°-Dipolo DOBLADO.
MAGNÉTICO 2.
Se realiza la estimación de resultados con la variación
ANÁLISIS 8: El tipo de polarización que presenta esta de 3 parámetros en una escala de valores:
antena dipolo magnético, es una polarización vertical
como se puede observar el círculo amarillo y el círculo  Prueba 1 - Se varía el valor de longitud en 5
azul tienen una aproximación de un 95% hacia el valor estimaciones diferentes.
de los ejes de 90 y 270 , lo que quiere decir que al  Prueba 2 - Se varía el valor de corriente en 6
aumentar la frecuencia de 20 MHz de la Prueba 1 a estimaciones diferentes.
60MHz de la Prueba 2 el patrón de radiación va a crecer  Prueba 3 - Se varía el valor de frecuencia en 4
o aumentar en amplitud, también da a entender que la estimaciones diferentes.
longitud que se ingrese influye por ejemplo en el caso
anterior se tuvo un dl de 50 cm y el patrón mantuvo un 5.3.1 PRUEBA 1: Se ingresa los siguientes datos con
equilibrio, pero ahora con un dl de 100cm el patrón de valores constantes y con valores de longitud mostrados
radiación crece con una aproximación mayor lo que es en la tabla 9 con el fin de obtener los parámetros que
bastante notable. posee una Antena Dipolo Doblado.
En el plano E se puede observar que el patrón de DATOS 1:
radiación tiene un valor de 20 y su corte tiende a estar
cerca de los ejes de referencia y en el plano H tiene un  d=3.4 m
valor de 0.05, la gráfica en 3D muestra como la antena  a1=a2=0.0003 m
dipolo magnético tiene un crecimiento en la parte de sus  f=110 Mhz
laterales.  I=0.3 mA
 P (100,180,230) m
RESULTADOS 1:
Tabla 9: Resultados 1– Parámetros Básicos Antena Dipolo

Doblado variando la longitud.
LONG CAMPO ELECTRICO CAMPO MAGNÉTICO Zant

1 j5.31e-05*e^j(1939.6 ) j1.40e-07*e^j(1939.6 ) 6.22e-61-2.54e-13i
3 j0.0003*e^j(1939.6 ) j8.06e-07*e^j(1939.6 1.58e-59-7.62e-13i
6 j-0.0001*e^j(1939.6 ) j-4.20e-07*e^j(1939.6 ) 1.55e-58-1.52e-12i
8 j-1.35e-06*e^j(1939.6 j-3.59e-09*e^j(1939.6 4.14e-58-2.03e-12i
10 j-0.0002*e^j(1939.6 ) j-5.35e-07*e^j(1939.6 ) 8.94e-58-2.54e-12i
21
mismo, el patrón de radiación tuvo un aumento de 16

lóbulos como se puede observar en la en la Fig 27 y
gracias a esto se llega a definir que en una Antena
Dipolo Doblado el número de lóbulos será igual al doble
del valor de su longitud.
Ver (ANEXO 3)
Figura 26. Patrón de radicación de 0 a 360°-Dipolo

DOBLADO – 3m de longitud.
Figura 28. Patrón de radicación plano H - Dipolo

DOBLADO-3m de longitud.

DOBLADO – 8m de longitud.
Figura 29. Patrón de radicación plano E - Dipolo
ANÁLISIS 9: El tipo de polarización que presenta esta
antena dipolo doblado es una polarización vertical,
posee un comportamiento distinto al tener valores de
longitud diferentes, se establece un rango de valores
que empieza en 1m y termina en 10m, en la Fig 26 se
estima una prueba con una longitud de 3m y se observa
que tiene un comportamiento diferente en el patrón de
radiación ya que no cuenta con los 2 lóbulos
convencionales sino que tiende a incrementar los
lóbulos, en este caso al ingresar 3m sus lóbulos
aumentan el doble de su longitud es decir a 6 lóbulos,
para verificar si esta prueba es correcta se ingresa un
nuevo valor 8m y su resultado fue exactamente el
22

valores constantes y con valores de corriente mostrados
en la tabla 10 con el fin de obtener los parámetros que
posee una Antena Dipolo Doblado.
DATOS 1:
 Long=1 m
 d=3.4 m
 a1=a2=0.0003 m
 f=110 Mhz
 P (100,180,230) m
RESULTADOS 2:

Doblado variando la corriente.
Figura 30. Patrón de radicación 3D-Dipolo DOBLADO- ImA CAMPO ELECTRICO CAMPO MAGNÉTICO Zo
3m de longitud. 0.2 j2.87e-05*e^j(711.20 ) j7.63e-08*e^j(711.20) 2256.70+5425.46i
0.6 j8.63e-05*e^j(711.20 ) j2.29e-07*e^j(711.20 ) 2256.70+5425.46i
0.9 j3.43e-07*e^j(711.20) j0.00012*e^j(711.20 ) 2256.70+5425.46i
2 j0.00028*e^j(711.20 ) j7.63e-07*e^j(711.20) 2256.70+5425.46i
6 j0.00086*e^j(711.20 ) j2.29e-06*e^j(711.20 ) 2256.70+5425.46i
9 j0.00129*e^j(711.20 ) j3.43e-06*e^j(711.20 ) 2256.70+5425.46i
Corriente Vs Zant
2500
2000
1500
Zant
1000
500
0
Figura 31. Patrón de radicación plano H - Dipolo 0.2 0.6 0.9 2 mA 6 mA 9 mA
DOBLADO-8m de longitud. mA mA mA
Corriente
Figura 33. Gráfica Corriente vs. Zant.
ANÁLISIS 10: La antena dipolo doblado posee un

segundo comportamiento al tener valores de corriente
diferentes, se establece un rango de valores que
empieza en 0,2mA y termina en 9mA como se puede
observar en la Fig 33 el parámetro que se ve
notablemente afectado es el valor de Zant ya que
presenta una permanencia en sus valores, por lo que se
define que mientras va creciendo el valor de corriente, el
parámetro de Zant tiende a ser constante es decir que el
nuevo valor de Zant independientemente del valor de
corriente que se ingrese no va a cambiar.
Figura 32. Patrón de radicación plano E - Dipolo

23

DOBLADO – 0,2mA de corriente. Figura 36. Patrón de radicación plano E - Dipolo
DOBLADO-0,2mA de corriente.
ANÁLISIS 11: Gracias a esta variación el patrón de
radiación que presenta esta antena también sufre
cambios representativos por ejemplo el patrón de
radiación de la corriente de 0,2mA empieza con una
aproximación de un 95% hacia el valor de sus ejes de
90° y 270°, lo que quiere decir que el patrón aumenta su
amplitud, pero para las corrientes de 0,6mA y 0.9mA la
amplitud del patrón de radiación empieza a disminuir,
ahora a partir de la corriente de 2mA la amplitud del
patrón vuelve a aumentar en un 95% de aproximación
hacia sus ejes, pero, para las corrientes de 6mA y 9mA
la amplitud del patrón de radiación empieza a disminuir,
es decir que para corrientes bajas se tiene un patrón de
radiación con una amplitud que aumenta mientras que
con corrientes altas la amplitud va a disminuir.
Figura 37. Patrón de radicación 3D-Dipolo DOBLADO-

0,2mA de corriente.

valores constantes y con valores de frecuencia
mostrados en la tabla 11 con el fin de obtener los
parámetros que posee una Antena Dipolo Doblado.
DATOS 1:
 d=3.4 m
 Long=1m
 a1=a2=0.0003 m
 I=0.9 mA
 P (100,180,230) m
RESULTADOS 3:
DOBLADO – 0,2 mA de corriente.
24

Doblado variando la frecuencia.
FRE CAMPO MAGNETICO CAMPO ELECTRICO ZANT

C
110 j3.43e-07*e^j(711.20) j0.00012*e^j(711.20 ) 141.81-216.77i
150 j5.42e-07*e^j(969.83) j0.00020*e^j(969.83 ) 292.48+169.84i
210 j7.23e-07*e^j(1357.76 j0.00027*e^j(1357.76 486.00+792.97i
300 j4.22e-07*e^j(1939.6) j0.00015*e^j(1939.66) 129.37+89.19i

DOBLADO – 150 Mhz de frecuencia.

DOBLADO – 150 Mhz de frecuencia.
ANÁLISIS 12: La antena dipolo doblado posee un tercer

comportamiento al tener valores de frecuencia
diferentes, se establece un rango de valores que Figura 40. Patrón de radicación plano E - Dipolo
empieza en 110MHz y termina en 300MHz como se DOBLADO – 150 Mhz de frecuencia.
puede observar en la Tabla 11 el parámetro que se ve
notablemente afectado es el valor de Zant, ésta
impedancia tiene una tendencia de ir creciendo o
aumentado por lo que se define que mientras va
aumentando el valor de frecuencia, el parámetro de Zant
tiene la misma relación de crecimiento lo que quiere
decir son directamente proporcionales.
ANÁLISIS 13: Como se puede observar para el Patrón

de radiación con una frecuencia de 150MHz se tiene una
aproximación de un 95% hacia el valor de los ejes de
90° y 270°, como se observa en la Fig 38 , si
comparamos con una frecuencia de 210MHz resulta que
no tiene variación su aproximación es similar tiene un
95% hacia sus ejes lo que quiere decir que su amplitud
es constane.
Figura 41. Patrón de radicación 3D-Dipolo DOBLADO-

150 Mhz de frecuencia.
25
6. CONCLUSIONES de diagramas de radiación estos diagramas

se los puede obtener al variar longitud de la
 El software diseñado en este proyecto permite onda de la antena , conservando su
al usuario ingresar los distintos parámetros frecuencia, el punto de análisis y la corriente
tanto de la antena dipolo hertziano como la de se puede obtener las distintas gráficas de los
la antena de dipolo magnético para el cálculo campos y el patrón de radiación los cuales
de diagramas de radiación estos diagramas servirán para analizar el comportamiento de
se los puede obtener al variar longitud de la estas antenas.
onda de la antena , conservando su frecuencia,
el punto de análisis y la corriente se puede  La metodología de estudio en este trabajo me
obtener las distintas gráficas de los campos y permitió aprender a observar y comprender el
el patrón de radiación los cuales servirán para patrón de radiación en los campos de las ondas
analizar el comportamiento de estas antenas. electromagnéticas. En las 3 dimensiones, los
cortes del Azimut permiten ver la energía y el
 La metodología de estudio en este trabajo me aumento de esta, y de campo al igual sucede
permitió aprender a observar y comprender el con la vista del patrón de elevación. Las cuales
patrón de radiación en los campos de las ondas permitieron comprender como es o debería ser
electromagnéticas. En las 3 dimensiones, los el campo, la energía radiada y el
cortes del Azimut permiten ver la energía y el comportamiento real de ambos en una antena.
aumento de esta, y de campo al igual sucede
con la vista del patrón de elevación. Las cuales 7. RECOMENDACIONES
permitieron comprender como es o debería ser
el campo, la energía radiada y el Ingresar de las fórmulas a Matlab Verificando que cada
comportamiento real de ambos en una antena. una de ellas realice lo que se estima y comprobar con
los cálculos teóricos.
 Tanto el dipolo eléctrico como el magnético,
son antenas sumamente ineficientes esto es,
requieren de corrientes excesivamente altas 8. REFERENCIAS
para irradiar bajas potencias, o de muy altas
[1] UNIVERSIDAD NACIONAL DE LA PLATA,
potencias en el medio para inducir muy bajas
«Resumen de Fórmulas sobre Radiación y Antenas,»
corrientes, sin embargo, las expresiones
obtenidas para estas antenas, permiten [En línea]. Available:
simplificar de cierta manera el análisis de https://catedra.ing.unlp.edu.ar/electrotecnia/camposyo/G
estructuras más complejas. uia%20TP12_Parte1.pdf.
[2] J. Anguera y A. Pérez, «TEORÍA DE ANTENAS,»

 En conclusión, un dipolo doblado, equivale,
desde el punto de vista de radiación a un dipolo 2008. [En línea]. Available:
simple con corriente de valor doble, e http://www.salleurl.edu/semipresencial/ebooks/ebooks/e
impedancia 4 veces. book_teoria_antenas.pdf.
[3] J. J. M. Fuentes, «Radiación y Radiocomunicación,»

 Una comparación del dipolo doblado con el
dipolo simple es el hecho de que el segundo 2010. [En línea]. Available:
tiene un patrón de radiación omnidireccional https://personal.us.es/murillo/docente/radio/documentos/
totalmente, mientras que el patrón de radiación tema2.pdf.
del dipolo doblado tiene una forma más
parecida a la de un cardioide. Por último, el [4] L. D. A. Torres, «Propiedades de las antenas,» 2013.
dipolo doblado permite lograr un ancho de [En línea]. Available:
banda mayor al del dipolo simple, pudiendo https://es.slideshare.net/anderssoncruz/libro-1-antenas.
alcanzar un rango de operación del 5% a partir
de la frecuencia central.
 El dipolo doblado es, un dipolo cuyos brazos

han sido doblados por la mitad y plegados
sobre sí mismos. En esencia, una antena única
formada por dos elementos. Un elemento se
alimenta en forma directa, mientras que el otro
tiene acoplamiento inductivo en los extremos.
 El software diseñado en este proyecto permite

al usuario ingresar los distintos parámetros
tanto de la antena dipolo hertziano como la de
la antena de dipolo magnético para el cálculo
26
ANEXO 1: CÓDIGO EN MATLAB DEL DIPOLO MAGNETICO Y HERTZIANO
clc
clear
close all
disp ( '----------------------------------------------------------- ')
disp ( ' UNIVERSIDAD POLITECNICA SALESIANA ' )
disp ( 'Integrantes: ' )
disp ( '*CARILLO KAYHERIN ' )
disp ( '*CUYAGO MICHELLE ' )
disp ( '*GONZÁLEZ ALEJANDRO ' )
disp ( 'TEMA: CÁLCULOS PARA UNA ANTENA DIPOLO HERZZIANO Y DIPOLO ELECTRICO ' )
disp ( '----------------------------------------------------------- ')
disp ( 'ELIJA EL TIPO DE ANTENA ' )

disp ( 'PRESIONE 1 ANTENA DIPOLO HERTZIANO:' )
disp ( 'PRESIONE 2 ANTENA DIPOLO MAGNETICO2:' )
OPCION=input('=' );
if(OPCION==1)
disp ( '--------- ANTENA DIPOLO HERZTZIANO ---------- ')

disp ( '--------- INGRESO DE DATOS ---------- ')
dl=input ( 'Ingrese el valor dl en [cm]: ');
f=input( 'Ingrese el valor Frecuencia en [MHZ]: ');
I=input ( 'Ingrese el valor de la Corriente en [uA]: ');
dl1=dl/100;
f1=f*10^6;
I1=I*10^-6;
disp(['dl:',num2str(dl1),' [m]'])
disp(['f1:',num2str(f1),' [HZ]'])
disp(['I:',num2str(I1),' [A]'])
disp('Ingrese el valor del punto P(x,y,z):')
x=input ( 'Ingrese el valor de x en [km]: ');
y=input ( 'Ingrese el valor de y en [km]: ');
z=input ( 'Ingrese el valor de z en [km]: ');
disp(['x:',num2str(x),' y:',num2str(y),' z:',num2str(z),' [km]'])

disp ( '--------- ------------- ------------ ')
disp('Calculos de Parámetros')
c=3*10^8;
landa=c/f1;
bo=2*pi/landa;
no=120*pi;
disp(['bo=',num2str(bo),' [rad/m]'])
r=sqrt(x^2+y^2+z^2);
r1=r*1000;
disp(['r=',num2str(r),' [km]'])
ar=z/r;
27
angulot=acos (ar);
grad=(angulot*180)/pi;
sen=sin(angulot);
sen2=sen^2;
disp(['angulo Theta=',num2str(angulot),' [rad]'])

disp(['angulo Theta=',num2str(grad),' [°]'])
E1= (-j*(bo)*(r1));
eff=j*((I1*dl1)/(4*pi*r1))*no*bo*sen;
disp ( '--------- ------------- ------------ ')
disp('***** Antena Dipolo Hertziano *****')
disp('***** Campo Eléctrico*****')

disp(['Eff:',num2str(eff),'*e^(',num2str(E1),' ) [v/m]'])
disp('***** Campo Magnético*****')
disp('***** Densidad de Potencia*****')

sav=(no*(I1^2)*(dl1^2)*sen2*(bo^2))/(32*(pi^2)*(r1^2));
disp(['SAV:',num2str(sav),' [w/m]'])
disp('***** Potencia de Radiación*****')

Prad=10*(I1^2)*(dl1^2)*(bo^2);
disp(['Prad:',num2str(Prad),' [w]'])
disp('***** Resistencia de Radiación*****')

Rrad=80*(pi^2)*(dl1/landa)^2;
disp(['Rrad:',num2str(Rrad),' [ohm]'])
disp('***** Ganancia Directiva *****')

Di=1.5*sen2;
disp(['Di:',num2str(Di),' []'])
disp ( '***** TEMA: DIAGRAMA DE RADIACIÓN DE UN DIPOLO HERTZIANO***** ' )
CREACION DE LAS GRÁFICA
eo=8.854*10^(-12); %constante permitibidad en el vacio

u0=4*pi*10^-7; %constante permeabilidad en el vacio
w=2.*pi.*f1; %frecuencia angular
no=120*pi; % sqrt(uo/eo)
c=3.*10^8; %velocidad de la luz
fi=(0:.01:1)*2*pi; %fi
teta=(0:.01:1)*pi; %theta
disp ( '1:Gráfica en 3D, Plano H y Plano E' )

disp ( 'PRESIONE 1 :' )
28
g=input('=> ' );
if(g==1)
%EN 3D
[FI,TETA]=meshgrid(fi,teta); % FI y TETA asignar para variables gráfica con malla
E=abs((I1*dl.*(sin(TETA)).*(exp(-(i*bo*r1)))*(i*(bo^2/r1)+(bo/(r1^2))-i*(1/(r1^3))))/
(w*r1*eo*4*pi)); % Campo eléctrico del dipolo elemental campos cercanos
X=E.*sin(TETA).*cos(FI); % componente en X del campo eléctrico
Y=E.*sin(TETA).*sin(FI); % componente en Y del campo eléctrico
Z=E.*cos(TETA); % componente en Z del campo eléctrico
figure(1)
subplot(2,2,4)
S=surface(X,Y,Z,(abs(E))); % Gráfica del campo eléctrico en 3D

axis equal % ejes iguales
axis off
lighting gouraud;
shading interp
title('PATRÓN DE RADIACIÓN EN 3D')
view(0,40)
rotate3d
subplot(2,2,1)
title('PLANO E')
xlabel('Y');
ylabel('Z');
%%figure(4) %FIGURA 2
%GRAFICAR CAMPO MAGNÉTICO EN COORDENADAS POLARES
subplot(2,2,2)
title('PLANO E')
polar(TETA,E);
hold on
polar(TETA,-E);
hold on
subplot(2,2,3)
plot(X,Y);
xlabel('X');
ylabel('Y');
end
end
if(OPCION==2)
disp ( '********DIPOLO MAGNÉTICO ELEMENTAL ********* ')
disp ( '--------- INGRESO DE DATOS ---------- ')
ar=input ( 'Ingrese el valor del radio en [cm]: ');% radio
29
fre2=input( 'Ingrese el valor Frecuencia en [MHZ]: ');

I2=input( 'Ingrese el valor de la Corriente en [mA]: ');
xm=input ( 'Ingrese el valor de x en [km]: ');
ym=input ( 'Ingrese el valor de y en [km]: ');
zm=input ( 'Ingrese el valor de z en [km]: ');
%%calculo del campo magnetico

ar1=ar/100;%conversion del valor de radio de cm a m
fre1=fre2*(10^(6));% conversion de la frecuencia
I3=I2*10^-3; %corriente en [A]
c1=3*10^8; %constante c del espacio
landa1=c1/fre1; %landa
bo1=2*pi/landa1;
no1=120*pi;
disp(['bo=',num2str(bo1),' [rad/m]'])
r4=sqrt(xm^2+ym^2+zm^2);%coordenadas del punto

rt=r4*1000;
arg=zm/r4;%%argumento de la funciÓn trigonometrica

angulot2=acos (arg); %%Ángulo en radianes
grad1=(angulot2*180)/pi; %Ángulo en grados
sen2=sin(angulot2); %%calculo del angulo seno
sen4=sen2^2;
e4=(-j*(bo1)*(rt));
%%CALCULO DE CAMPO MAGNÉTICO

HFF=(-(I3*(ar1^(2))*bo1^2)/(4*rt))*sen2;
%%CAMPO ELÉCTRICO
EFF=(no1*I3*(ar1^(2)))*(bo1^(2))/(4*rt)*sen2;
%%CALCULO DE LA DENSIDAD MEDIA DE POTENCIA

A=3.1416*(ar1^2);
A1=(A/(landa1)^2)^2;
SAV=1860*(I3)^2*A1*(sen4/rt^2);
%%POTENCIA RADIADA
P=15585*((I3)^2)*(A/(landa1)^2)^2;
%%RESISTENCIA DE RADIACIÓN
Rrad1=31170*(A/(landa1)^2)^2;

disp(['Hff:',num2str(HFF),'*e^(',num2str(e4),' ) [A/m]'])
disp('***** Campo Eléctrico *****')

disp(['Eff:',num2str(EFF),'*e^(',num2str(e4),' ) [v/m]'])
disp('***** Densidad media de Potencia*****')
30
disp(['Prad:',num2str(SAV),' [w]'])
disp('***** Potencia de Radiación*****')

disp(['P:',num2str(P),' [w]'])
disp('***** Resistencia de Radiación*****')
disp(['Rrad:',num2str(Rrad1),' [ohm]'])
%%GRÁFICAS
disp ( '1:Gráfica en 3D, Plano H y Plano E' )

disp ( 'PRESIONE 1 :' )
opc=input('=> ' );
if(opc==1)
phi=(0:.01:1)*2*pi; %ÁNGULO PHI

teta=(0:.01:1)*pi; %ÁNGULO TETA
%%%CAMPO ELÉCTRICO
[PHI,TETA]=meshgrid(phi,teta); % CREACIÓN DE MESHGRID PARA GRAFICAR EN 3D
E=abs((((no1*I3*A^2*bo1^2)/(4*ar1)))*sin(TETA)*(exp(-imag(bo1*ar1)))); %ECUACION DEL CAMPO ELECTRICO
%%figure(1) %VENTANA 1
subplot(2,2,1)
%GRAFICAR CAMPO ELÉCTRICO EN COORDENADAS POLARES
polar(PHI,E);
title('PLANO E ')
hold on
polar(PHI,-E);
hold on
%CAMPO MAGNETICO
[PHI,TETA]=meshgrid(phi,teta); % CREACIÓN DE MESHGRID PARA GRAFICAR EN 3D
H=abs((-((I3*A^2*bo1^2)/(4*ar1)))*sin(TETA)*(exp(-imag(bo1*ar1))));%ECUACIÓN DEL CAMPO MAGNETICO
%CAMBIO DE COORDENADAS ESFERICAS A RECTANGULARES PARA GRAFICAR
X1=H.*sin(TETA).*cos(PHI);
Y1=H.*sin(TETA).*sin(PHI);
Z1=H.*cos(TETA);
R1=square(X1^2+Y1^2+Z1^2);%COMPONEN DEL RADIO
subplot(2,2,2)
%GRAFICAR CAMPO MAGNÉTICO EN COORDENADAS POLARES
polar(TETA,H);
title('PLANO H ')
hold on
polar(TETA,-H);
hold on
%PATRON DE RADIACIÓN EN 3D
[PHI,TETA]=meshgrid(phi,teta); % CREACIÓN DE MESHGRID PARA GRÁFICAR EN 3D
E=abs((((no1*I3*A^2*bo1^2)/(4*ar1)))*sin(TETA)*(exp(-imag(bo1*ar1))));%ECUACIÓN DEL CAMPO ELECTRICO
31
%CAMBIO DE COORDENADAS ESFERICAS A RECTANGULARES PARA GRÁFICAR

X=E.*sin(TETA).*cos(PHI);
Y=E.*sin(TETA).*sin(PHI);
Z=E.*cos(TETA);
R=square(X^2+Y^2+Z^2); %COMPONENTE DEL RADIO
subplot(2,2,3)
S=surface(X,Y,Z,(abs(E))); %GRÁFICA E 3D
mesh(X,Y,Z)
title('PATRÓN DE RADIACIÓN 3D')
%EJES SIMÉTRICOS, APAGADOS
axis equal
axis off
%ESTILO SOLIDO Y SOMBREADO
lighting gouraud
shading interp
view(0,30) %POSICIÓN DE LA GRÁFICA
end
end
ANEXO 2: CÓDIGO EN MATLAB DEL DIPOLO DOBLADO
disp('----INGRESO DE DATOS ----')

while (l<=0)
disp('')
if (l>=1)
end
end
while (f1<=0)
disp('')
if (f1>=1)
end
end
while (I1<=0)
disp('')
if (I1>=1)
end
end
32
A1=input ( 'Ingrese el valor del diametro del conductor 1 [m]: ');

while (A1<=0)
disp('')
if (A1>=1)
end
end
while (A2<=0)
disp('')
if (A2>=1)
end
end
D=input ( 'Ingrese el valor de la distancia de separacion [m]: ');
while (D<=0)
disp('')
D=input ( 'Ingrese el valor de la distancia de separacion [m]: ');
if (D>=1)
end
end

x=input ( 'Ingrese el valor de x en [m]: ');
y=input ( 'Ingrese el valor de y en [m]: ');
z=input ( 'Ingrese el valor de z en [m]: ');
%%constantes
no=120*pi;
c=3*10^(8);
long=l %longitud [m]
f=f1*10^(6);%Frecuencia en [H]
fre=f;
I=I1*10^(-3)%%[A] transformacion a Amperios
a1=A1
a2=A2
d=D
%%desarrollo
landa=(c/fre);
r=sqrt(x^2+y^2+z^2);%%[m]
ar=z/r ;%argumento angulo

angulot=acos (ar); %angulo [rad]
Bo=2*pi/landa;
arg1=cos(Bo*(long/2)*cos(angulot));
33
arg2=cos(Bo*(long/2));
F=(arg1-arg2)/sin(angulot);%calculo de F(o)
E1= ((Bo)*(r));
EFF=((no*I)/(pi*r))*F
HFF=(I/(pi*r))*F
%%impedancia
%%RADIO EQUIVALENTE
u=a2/a1;
v=d/a1;
Aeq=log(a1)+((1)/(1+u)^2)*((u)^2*log(u)+2*u*log(v));
aeq=exp(Aeq)
%%INTERPOLACION
B=(pi/landa)*long;
if(B>=0.1 && B<0.2)
RB=0.1506+((RB-0.1)*(0.7980-0.1506))/(0.1) ;
XB=1.010+((RB-0.1)*(2.3002-1.010))/(0.1);
end
if(B>=0.2 && B<0.3)
RB=0.7980+((RB-0.2)*(1.821-0.7980))/(0.1) ;
XB=2.302+((RB-0.2)*(3.818-2.302))/(0.1);
end
if(B>=0.3 && B<0.4 )

RB=1.821+(((B-0.3)*(3.264-1.821))/(0.1));
XB=3.818+(((B-0.3)*(5.584-3.818))/(0.1));
end
if(B>=0.4 && B<0.5 )

RB=3.264+(((B-0.4)*(5.171-3.264))/(0.1));
XB=5.584+(((B-0.4)*(7.141-5.584))/(0.1));
end
if(B>=0.5 && B<0.6 )

RB=5.171+(((B-0.5)*(7.563-5.171))/(0.1));
XB=7.141+(((B-0.5)*(8.829-7.141))/(0.1));
end
34
if(B>=0.6 && B<0.7 )

RB=7.563+(((Bl2-0.6)*(10.48-7.563))/(0.1));
XB=8.829+(((Bl2-0.6)*(10.68-8.829))/(0.1));
end
if(B>=0.7 && B<0.8 )

RB=10.48+(((Bl2-0.7)*(13.99-10.48))/(0.1));
XB=10.68+(((Bl2-0.7)*(12.73-10.68))/(0.1));
end
if(B>=0.8 && B<0.9 )

RB=13.99+(((B-0.8)*(18.16-13.99))/(0.1));
XB=12.73+(((B-0.8)*(15.01-12.73))/(0.1));
end
if(B>=0.9 && B<1 )

RB=18.16+(((B-0.9)*(23.07-18.16))/(0.1))
XB=15.01+(((B-0.9)*(17.59-15.01))/(0.1))
end
if(B>=1 && B<1.1 )

RB=23.07+(((B-1)*(28.83-23.07))/(0.1))
XB=17.59+(((B-1)*(20.54-17.59))/(0.1))
end
if(B>=1.1 && B<1.2 )

RB=28.83+(((B-1.1)*(35.60-28.83))/(0.1))
XB=20.54+(((B-1.1)*(23.93-20.54))/(0.1))
end
if(B>=1.2 && B<1.3 )

RB=35.6+(((B-1.2)*(43.55-35.6))/(0.1));
XB=23.93+(((B-1.2)*(27.88-23.93))/(0.1));
end
if(B>=1.3 && B<1.4 )

RB=43.55+(((B-1.3)*(52.92-43.55))/(0.1))
XB=27.88+(((B-1.3)*(32.2-27.88))/(0.1))
end
if(B>=1.4 && B<1.5 )

RB=52.92+(((B-1.4)*(64.01-52.92))/(0.1))
XB=32.2+(((B-1.4)*(38-32.2))/(0.1))
end
if(B>=1.5 && B<(pi/2) )

RB=64.01+(((B-1.5)*(73.12-64.01))/((pi/2)-1.5))
XB=38+(((B-1.5)*(42.46-38))/(0.1))
end
if(B==(pi/2) )
RB=73.12
35
XB=42.46
end
Z1=log(long/(2*aeq))-1
cotag=cot((B*long)/2)
zdip=RB-j*(120*Z1*cotag-XB);
ydip=((1)/(zdip))
zo=120*acosh(d/(2*sqrt(a1*a2)));
yo=1/zo
vc=cot((B*long)/2)
Yant=(ydip/4)-i*((yo/2)*vc)
Zant=1/Yant
IMPRESIONES EN PANTALLA
disp('****** CALCULOS Antena Dipolo DOBLADO ****** :')

disp(['Bo=',num2str(Bo),' [rad/m]'])
disp(['LANDA=',num2str(landa),' [m]'])
disp(['r=',num2str(r),' [m]'])
disp(['angulo=',num2str(angulot),' [rad]'])
disp(['F(Teta)=',num2str(F),' '])
disp(['u=',num2str(u),' '])
disp(['v=',num2str(v),' '])
disp(['radio equivalente [aeq]=',num2str(aeq),' [m] '])
disp(['interpolacion [B]=',num2str(B),' '])
disp([' [R]=',num2str(RB),' '])
disp([' [X]=',num2str(XB),' '])
disp([' [Zdip]=',num2str(zdip),' [OHM] '])
disp([' [ydip]=',num2str(ydip),' [1/OHM] '])
disp([' [Zo]=',num2str(zo),' [OHM] '])
disp([' [yo]=',num2str(yo),' [1/OHM] '])
disp('***** Campo Eléctrico*****')

disp(['EFF=','j',num2str(EFF),'*e^j(',num2str(E1),' ) [v/m]'])
disp(['HFF:','j',num2str(HFF),'*e^j(',num2str(E1),' ) [A/m]'])
disp ( '***** IMPEDANCIA DE LA ANTENA ***** ' )

disp([' [Zant]=',num2str(Zant),' [1/OHM] '])
%%grafica
bo=(2*pi)/landa %beta
%gráfica
36
fi=(0:.01:1)*2*pi; %fi
teta=(0:.01:1)*2*pi; %theta
%fórmulas
teta1=acos(z/r)
disp ( 'Graficas en 3D, PlanoH y Plano E del dipolo largo' )

disp ( 'El Campo Eléctrico, Magnético, la impedancia y la resistencia de pérdida' )
disp ( 'Presione 1:' )
g=input('=> ' );
if(g==1)
%EN 3D
subplot(2,2,1)
[FI,TETA]=meshgrid(fi,teta); % FI y TETA asignar para variables gráfica con malla
Fteta=(cos(bo*(l/2).*cos(TETA))-cos(bo*(l/2)))./(sin(TETA));
E=2*abs(((i*no*I).*Fteta).*(exp(-(i*bo*r)))/(2*pi*r)); % Campo eléctrico del dipolo elemental

campos lejanos
X=E.*sin(TETA).*cos(FI); % componente en X del campo eléctrico

Y=E.*sin(TETA).*sin(FI); % componente en Y del campo eléctrico
Z=E.*cos(TETA); % componente en Z del campo eléctrico
S=surface(X,Y,Z,(abs(E))); % Gráfica del campo eléctrico en 3D

axis equal % ejes iguales
axis off
lighting gouraud;
shading interp
title('Dipolo Largo en 3D')
view(0,40)
rotate3d
subplot(2,2,4)

title('PLANO E')
xlabel('Y');
ylabel('Z');
subplot(2,2,3)
title('PLANO E')
polar(TETA,E) % Gráfica plano E
subplot(2,2,2)
plot(X,Y);
xlabel('X');
37
ylabel('Y');
end
38
ANEXO 3: DIPOLO DOBLADO VARIACION DE LONGITUD
LONGITUD DE 1m
39
LONGITUD DE 6m
40
LONGITUD DE 10m
41
ANEXO 4: DIPOLO DOBLADO VARIACION DE CORRIENTE
DE 0,9mA
42
LONGITUD DE 6 mA
43
LONGITUD DE 9 mA
44
ANEXO 5: DIPOLO DOBLADO VARIACION DE FRECUENCIA
LONGITUD DE 300 MHZ
45
LONGITUD DE 210 MHZ
46

Informe - Final Antenas - Completo

Cargado por

Copyright:

Formatos disponibles

Informe - Final Antenas - Completo

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Informe - Final Antenas - Completo

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD POLITÉCNICA SALESIANA ANTENAS

Fecha: 21/06/2019 TELECOMUNICACIONES

PATRÓN DE RADIACIÓN DE LA ANTENA DIPOLO HERTZIANO,

RESUMEN: En el siguiente informe se presenta el El patrón de radiación de una antena se puede

 Implementar un código en Matlab para

2.1 PATRÓN DE RADIACIÓN Figura 1. Sistema de coordenadas polares empleado

2.2 ANTENA DIPOLO HERTZIANO

Por lo tanto, el plano E es cualquier plano que contiene

2.2.1.2 Intensidad de Campo magnético ^

2.2.1.3 Vector de densidad de media potencia

2.2.1 CAMPOS EN ZONAS APARTADAS Se denomina vector de Poynting al vector cuyo módulo

^ El vector de Poynting en esta región lejana es radial y en

mediante la siguiente ecuación:

2.2.1.4 Potencia Radiada

La potencia total radiada por el dipolo se puede obtener

La resistencia de radiación se puede calcular a partir de

2.2.1.6 Ganancia Directiva

La gráfica del módulo de la densidad media de potencia

Figura 5. Dipolo magnético elemental. [2]

La potencia total radiada se puede obtener integrando la

2.3.1.1 Intensidad de Campo eléctrico ^

2.3.1.2 Intensidad de Campo magnético ^

El dipolo doblado se puede descomponer en el modo

Figura 7. (a) Patrón de radicación plano E. (b) Patrón de

2.4 ANTENA DIPOLO DOBLADO

Un dipolo doblado es una estructura formada por dos

Como la diferencia de potencial entre los puntos a y b es

La antena dipolo doblado por construcción tiene una

Siendo Yo la admitancia característica de la línea de

2.4.2.2 Modo par o antena

En el modo antena que se representa en la figura 3(c),

mismo contribuyen con sus campos magnéticos

El plano H es el definido por la dirección de máxima

2.4.3 CAMPOS EN ZONAS APARTADAS

Debido a que el interés del estudio de antenas radica

principalmente en el conocimiento de los campos en

Figura 12. (a) Patrón de radicación plano E. (b) Patrón

4. DESARROLLO Y PROCEDIMIENTO D=4 πSAV r

4.1 PARÁMETROS – ANTENA DIPOLO

4.2.1 CAMPOS EN ZONAS APARTADAS

4.1.1 CAMPOS EN ZONAS APARTADAS  Campo Eléctrico

4.2.2 VECTOR DENSIDAD MEDIA DE POTENCIA

4.2.3 POTENCIA, RESISTENCIA Y GANANCIA

4.3.1 CAMPOS EN ZONAS APARTADAS  A equivalente

4.4 PARÁMETROS ADICIONALES GENERALES

Figura8. Software Matlab - New Script

4.5 SOFTWARE MATLAB Figura10. Menú de opciones de Antena.

4.5.1 Se ingresa al Software Matlab.

4.5.2 Buscar la pestaña de New Script para abrir un

Figura11. Diseño del Menú de opciones de Antena.

4.5.6 El usuario deberá escoger entra la opción 1 u

-OPCIÓN 1: Si se selecciona la opción 1 se despliega el

CARÁTULA DE PRESENTACIÓN disp ( '--------- ------------- ------------ ')

CÁLCULO DEL RADIO

INGRESO DE DATOS 1.3 Los ángulos se calculan con respecto al punto de

2. Cálculo de los parámetros principales de la

2.1 Una vez obtenido los valores de los parámetros

CÁLCULO DE CAMPO ELÉCTRICO

disp('* Campo Eléctrico*') 3. Para la elaboración de los planos h y plano e

%%RADIO EQUIVALENTE Aeq=log(a1)+((1)/(1+u)^2)((u)^2log(u)+2ulog(v));