2015

FIEE-UNMSM
 “Universidad Nacional Mayor de San
Marcos”
Radiación de la Antena Lineal

PROFESOR:
- Ing. Hernan Villafuerte

ALUMNA:
- Ofracio Leiva, Estefani Marylim

CODIGO:
- 12190252

CURSO:
- Métodos Numéricos

HORARIO:
- Lunes 6-8 pm
 FIEE-UNMSM 3

Resumen

En este trabajo se utiliza los fundamentos del Método de los Momentos,

una poderosa herramienta para la solución de problemas de campo
electromagnético, mostramos una aplicación a la antena lineal
obteniendo el patrón de radiación.

3
 FIEE-UNMSM 4

Índice general

1. Introducción 2
1.1. Objetivo ................................................................................................. 2

2. Método de los momentos (MoM). 3

2.1. Planteamiento....................................................................................... 3
2.2. Método de los momentos .................................................................... 3
2.3. Principio de equivalencia .................................................................... 5
2.4. Funciones de base y prueba ................................................................ 6
2.4.1. Funciones base ......................................................................... 6
2.4.2. Funciones Prueba .................................................................... 7

3. Ecuaciones Integrales 8
3.1. Ecuación Integral del Campo Eléctrico(EFIE) .................................... 8

4. Aplicación a la antena lineal 10

4.1. Antena lineal....................................................................................... 10
4.2. Ecuación integral de Pocklington ..................................................... 11
4.3. Aplicando el método de momentos ................................................ 12
4.3.1. Matriz de impedancia ........................................................... 13
4.4. Programación .................................................................................... 14

5. Conclusiones 16
Bibliografía 21

4
 FIEE-UNMSM 5

1. Introducción

El método de momentos, aplicado a problemas de

electromagnetismo, fue introducido por Roger F. Harrigton in 1967 en su
artículo, “Matrix Methods for Field Problems”. La implementación del
método de momentos en Lawrence Livermore National Labs durante los
70s, establecido esta técnica de solución para el diseño de antenas

Para determinar la distribución de corriente en una antena lineal resultado

de una excitación arbitraria puede ser establecido en términos de una
ecuación integral. Esta ecuación emplea una función de Green el cual
relaciona un campo eléctrico conocido de las condiciones de contorno con
una distribución desconocida de corriente en la antena. El método de
momentos (MoM) aplica expansiones para convertir la ecuación integral
en un sistema de ecuaciones lineales. Funciones de base son usados para la
expansión de la corriente y funciones de prueba para el campo eléctrico.
La distribución de corriente es luego construido de los coeficientes de la
expansión. Las características de la radiación de la antena son derivadas
luego del calculo de la distribución de corriente.

En este trabajo se implementara el método para el análisis de una antena

lineal. La teoría matemática en el cual es basada es presentada y se derivan
las ecuaciones integrales que describen la corriente en la antena. La
solución de esta ecuación integral es realizada por el método de momentos
el cual se basa en la expansión de la corriente en un conjunto de funciones
base.

1.1. Objetivo
Estimar la distribución de la corriente y el patrón de radiación en una
antena lineal, implementando un programa para el análisis de la antena
usando el método de momentos.

5
 FIEE-UNMSM 6

2. Método de los momentos (MoM).

2.1. Planteamiento
Numerosos problemas físicos conducen a ecuaciones integro-diferenciales que pueden
expresarse de la forma:
𝐿𝑢 = 𝑣 (2.1)

Donde “u” es un elemento desconocido de un espacio de funciones U, “y” es un

elemento conocido (prefijado) de un espacio de funciones V (que puede coincidir con
U) y L es un operador integro-diferencial de U en V. La ecuación (2.1) estará
completada con algún tipo de condición de contorno aplicable a “u”.
En general, “u” constituye la respuesta del sistema físico considerando una excitación
“v”, el operador L representa los fenómenos físicos que relacionan excitación y
respuesta junto a datos tales como geometría del sistema.
En problemas electromagnéticos, la función “v” representa magnitudes de tipo
corrientes, potenciales y campos tanto eléctricos como magnéticos con valores
impuestos al sistema (condiciones de contorno) y la función “u” suele representar
corrientes, densidades de carga, etc.

2.2. Método de los momentos

El método de los momentos es un procedimiento general para obtener soluciones
aproximadas de ecuaciones de la forma (2.1). El primer paso consiste en representar la
función incógnita “u” como combinación lineal de infinitas funciones que se
denominan funciones base:
∞

𝑢 = ∑ 𝐼𝑛 𝑓𝑛 (2.2)
𝑛=1
donde 𝑓𝑛 son las funciones base y In son coeficientes desconocidos. En la práctica es
imposible trabajar con sumas infinitas , por lo que reducimos el sumatorio a un número
finito de términos N.

Figura 2.1: Diagrama del Método de los momentos.

6
 FIEE-UNMSM 7

𝑢 = 𝑢𝑛 = ∑ 𝐼𝑛 𝑓𝑛 (2.3)
𝑛=1
Las funciones de la forma un estarán contenidas en el espacio funcional 𝑢𝑛 definido por
la base {𝑢1 … 𝑢𝑛 }. Como se representa en la figura. Si se sustituye el desarrollo de u
(2.3) en (2.1) y por la linealidad del operador.
𝑁

∑ 𝐼𝑛 𝐿𝑓𝑛 = 𝑣 (2.4)
𝑛=1

Esta expresión es válida si el operador L se puede aplicar sobre las funciones base, si su
elección es adecuada , puede obtenerse a partir del generador L un operador extendido
con las mismas propiedades de L que se pueda aplicarse sobre las funciones base.
El espacio generado por las N funciones 𝐼𝑛 𝐿𝑓𝑛 , en general, no contiene la funcion v. Por
tanto, al sustituir por su aproximación aparecerá un error.
𝑁

∑ 𝐼𝑛 𝐿𝑓𝑛 − 𝑣 = 𝑒𝑁 (2.5)
𝑛=1

Los coeficientes 𝐼𝑛 deberán escogerse de forma que minimicen la función error 𝑒𝑁.
En el método de los momentos este error se minimiza de la siguiente forma:
1. Se define un producto escalar valido tanto en V como en 𝐿𝑓𝑛 .
2. Se definen tantas funciones de peso o prueba, 𝑤𝑚 linealmente independientes como
funciones base N.
3. Se escogen los coeficientes 𝑤𝑚 de forma que los N productos escalares de la función
error 𝑒𝑁 por las N funciones de peso se cancelen:

〈𝑒𝑁, 𝑤𝑚 〉 = 0 (2.6)

Multiplicando escalarmente (2.5)

𝑁

∑ 𝐼𝑛 〈𝐿𝑓𝑛 , 𝑤𝑚 〉 = 〈𝑣, 𝑤𝑚 〉 (2.7)

𝑛=1

donde m = 1, 2, ..., N , que constituye un sistema de N ecuaciones, una por cada función
de peso, y N incógnitas, los coeficientes 𝐼𝑛 .
Sustituyendo los valores obtenidos al resolver (2.7) en (2.3) se obtiene la solución
aproximada buscada.
El conjunto de ecuaciones (2.7) se puede escribir de forma matricial como:

[Z][I] = [V] (2.8)

donde Z es la matriz del sistema (N × N ), denominada matriz de impedancias, con

𝑍𝑚,𝑛 〈𝑤𝑚 , 𝐿𝑓𝑛 〉, I es el vector de pesos incógnita (N × 1), con 𝐼𝑛 = 𝐼𝑛 y V es el
vector columna de valores conocidos (N × 1), con 𝑉𝑚 = 〈𝑣, 𝑤𝑚 〉
Despejando el vector de incógnitas:
[I] = [Z]-1 [V] (2.9) 7
 FIEE-UNMSM 8

En caso particular de que las funciones base y peso sean idénticas, al método de
los momentos se le suele denominar método de Galerkin.

2.3. Principio de equivalencia

El problema general que se pretende resolver es de la forma representada en la figura.

Figura 2.2: Problema general de dispersión

El campo electromagnético total en el medio 1 se puede descomponer en un campo

incidente o impreso (𝐸 ⃗⃗⃗𝑖 , ⃗⃗⃗⃗
𝐻𝑖 ) que sería generado por las fuentes suponiendo que no
existe obstáculo y un campo inducido o reflejado (𝐸 ⃗⃗⃗⃗𝑠 , ⃗⃗⃗⃗
𝐻𝑠 ) que es la perturbación debida
a la presencia del obstáculo. En el medio 2 no se realiza ninguna descomposición y el
campo recibe el nombre de transmitido (𝐸 ⃗⃗⃗𝑡 , ⃗⃗⃗⃗
𝐻𝑡 ).
El teorema de equivalencia permite anular el campo en el medio 2 sin modificarlo en el
medio 1, ello implica la introducción una distribución superficial de corrientes eléctricas
y magnéticas en la superficie de separación, cuyo valor puede obtenerse aplicando las
condiciones de frontera.

Figura 2.3: Problema equivalente.

8
 FIEE-UNMSM 9

⃗⃗𝐽𝑠 = 𝑛̂ × (𝐻
⃗⃗⃗⃗𝑖 + ⃗⃗⃗⃗
𝐻𝑠 ) (2.10)
⃗⃗⃗⃗
𝐾𝑠 = 𝑛̂ × (𝐸𝑖 + ⃗⃗⃗⃗
⃗⃗⃗ 𝐸𝑠 ) (2.11)

Con esta suposición el campo impreso estará generado por las fuentes originales, y el
campo inducido se deberá a las corrientes superficiales denominadas corrientes
inducidas. Para ambos casos se considera espacio infinito, lineal, homogéneo e isótropo
con las características del medio 1.

2.4. Funciones de base y prueba

2.4.1. Funciones base

Las funciones de base se pueden clasificar en dos categorías [3]:

- Funciones definidas en todo el dominio del operador.

Este tipo de funciones se caracteriza por anularse en un número finito de puntos
del dominio. Sobre geometrías particulares estas funciones permiten obtener un
número reducido de incógnitas.

- Funciones base de subdominios, es subdividir la antena en pequeños segmentos

y modelar la distribución de corriente en cada segmento por una construcción
geométrica que puede ser rectangular, triangular o sinusoidal. La amplitud de
estas construcciones representa los coeficientes de la función expandida.

Varios tipos de función base definidas en un subdominio:

Función base pulso

Función base triangulo lineal

Función base triangulo sinusoidal

donde β es la constante de variación de fase de la función a representar.

9
 FIEE-UNMSM 10

2.4.2. Funciones Prueba

Análogamente a las funciones base, las funciones de prueba se pueden clasificar en [3]:

- Funciones de prueba definidas en el dominio del operador

Estas funciones no suelen utilizarse en la práctica debido a los largos cálculos que
originan.

- Funciones de prueba definidas en un intervalo del dominio.

Dentro de esta categoría se incluyen las funciones presentadas anteriormente y

añadimos algunas.
Funciones prueba delta de Dirac δ(x − 𝑥𝑛 )
La elección de este tipo de funciones han sido utilizadas en el análisis de antenas
sencillas, modeladas con subdominios de dimensiones similares, pero cuando la
geometría de la antena se complica o los subdominios tienen dimensiones diferentes los
resultados se vuelven inestables.
Funciones de prueba pulso de excitación
Este tipo de funciones intentan mejorar los resultados que se obtiene con las funciones
pulso.

Presentan la ventaja de permitir una representación exacta de la excitación, salvo en los

extremos del intervalo, y el inconveniente de requerir cálculos más complicados que la
función pulso.

10
 FIEE-UNMSM 11

3. Ecuaciones Integrales
3.1. Ecuación Integral del Campo Eléctrico
(EFIE)
Para un sistema de cargas y corrientes que varían con el tiempo, podemos efectuar un
análisis de Fourier de la dependencia temporal y tratar de forma separada cada una de
las componentes. Por tanto, no perdemos generalidad si consideramos que los
potenciales, los campos y la radiación debidos a un sistema localizado de cargas varían
sinusoidalmente con el tiempo.[2]

Para obtener las magnitudes físicas tomaremos la parte real de las expresiones. Los
potenciales y campos electromagnéticos presentan el mismo tipo de dependencia con el
tiempo.

El campo eléctrico en función de los potenciales esta dado por:

,de la condición de Lorentz:

11
FIEE-UNMSM 12

12
 FIEE-UNMSM 13

4. Aplicación a la antena lineal

4.1. Antena lineal
Un dipolo eléctrico radiante es una antena lineal, que puede ser vista como un conductor
perfecto cilíndrico con radio a y longitud l en posición a lo largo del eje z alimentada
por su centro como se muestra en la figura. La variable R representa la distancia entre la
fuente de corriente y el punto

Figura 4.1: Antena lineal [5]

de observación del campo. La distribución de corriente 𝐼𝑍 (𝑧 , ) es definido a lo largo de

la longitud de la antena desde 𝑧 , = −𝐿/2 hasta 𝑧 , = 𝐿/2 Luego nosotros podemos
asumir que la corriente en el dipolo existe solamente como una corriente superficial ⃗⃗𝐽𝑠 .

Asumimos que a << λ y l >> a, esto es referido como la aproximación de hilo delgado
[1] . El campo eléctrico total 𝐸⃗ puede ser separado en dos componentes.

𝐸⃗ = 𝐸⃗ 𝑖 +𝐸⃗ 𝑠 (4.2)

donde, el campo eléctrico impreso 𝐸⃗ 𝑠 es debido a la excitación, que es distinto de cero

solo en el gap de alimentación

13
 FIEE-UNMSM 14

y el campo eléctrico dispersado 𝐸⃗ 𝑠 debido a la corriente inducida en la superficie de la

antena.

4.2. Ecuación integral de Pocklington

El campo 𝐸⃗ 𝑠 se relaciona con la ecuación (3.4)

esta formulación para la antena lineal es conocida como la ecuación integral de

Pocklington. Las características de la radiación son determinadas del conocimiento de la
14
 FIEE-UNMSM 15

distribución de corriente en la antena de las diversas técnicas disponibles para resolver

esta ecuación integral, el método de momentos es una de las mas populares en la
industria.

4.3. Aplicando el método de momentos

El procedimiento de solución se inicia definiendo la desconocida distribución de

corriente 𝐼𝑍 (𝑍 , ) en términos de un conjunto ortogonal de funciones base. En la
figura se muestra algunas construcciones de las funciones base,

Figura 4.2: Funciones base en subdominios.[5]

donde se mantiene la continuidad de la distribución de corriente a lo largo de la antena.

Discretizamos el dominio físico de las fuentes en un números N de tramos tomando
N+1 puntos constante h = L/N entonces 𝑍 , = 𝑛ℎ con n = 0, 1, 2,...N . Luego elegimos
las funciones base de tal forma que se asemejen a la distribución de corriente,
elegimos las del tipo triángulo sinusoidal:

15
 FIEE-UNMSM 16

donde n = 1, 2, ..., N − 1. La amplitud de estas funciones representan los coeficientes de

la función expandida.

Definimos las funciones de prueba en términos del delta de Dirac.

𝑊𝑛 = δ(z − 𝑧𝑚 ) (4.7)

donde 𝑧𝑚 son los puntos especificos en la antena en el cual las condiciones de contorno
se cumplen, corresponden al punto medio de de cada función de base, es decir 𝑧𝑚 = 𝑚ℎ
con m = 1, 2, ..., N

4.3.1. Matriz de impedancia

La matriz de impedancia esta dado por

𝑍𝑚,𝑛 = 〈𝑤𝑚 , 𝐿𝑓𝑚 〉

reemplazando, tenemos

16
 FIEE-UNMSM 17

2
La integral en 𝑧 , se puede resolver asumiendo que (𝑘 2 + 𝜕𝑧
𝜕 ,
2 )𝑔(𝑟, 𝑟 ) se mantiene

uniforme en el subdominio fuente [4]. de donde obtenemos:

17
 FIEE-UNMSM 18

4.4. Programación
%***************************************************************
%********************* ANTENA ***************************
%* Este es un programa desarrollado para resolver antenas Yagi-Uda y *
%* Reflector de esquina. Esta basado en la ecuacion integral de *
%* Pocklington, mediante el metodo de momentos.
%* OFRACIO LEIVA Estefani Marylim
%*********************************************************************
function []=YagiReflect
clear all
global M N D L LM YP radio N1 N2 M1 M2 GAMMA Inm
%***************************Numero de
MODOS********************************
M = input (' Ingrese el numero de modos: ');
%************************ Numero de
ELEMENTOS******************************
N = input (' \n Ingrese num de elementos: ');
%**********************Reflectores de Esquina por
lado*********************
NN = input('\n Ingrese el numero de reflectores de esquina por lado:
');
%****************Angulo de apertura del reflector
esquina******************
psi = input('\n Ingrese el angulo del reflector de esquina: ');
%*********************************************************************
*****
% INGRESO DE LONGITUD DE LOS ELEMENTOS EN EL ARREGLO Yagi - Uda
%*********************************************************************
*****
%*******************Longitud de los
DIRECTORES*****************************
if (N > 3)
fprintf (1, ' Todos los directores poseen la misma longitud?\n');
RESP = input (' RESPUESTA: (Y o N) ...... ', 's');
else
RESP = 'N';
end
fprintf (1, '\n');
if (RESP == 'Y') | (RESP == 'y');
LDIR = input ('Longitud uniforme de directores (en longitudes de
onda)= ');
L = LDIR * ones (1, N-2);
else (RESP == 'N') | (RESP == 'n');
a = 1;
while a <= (N-2);
fprintf(1,' Longitud del director # %2d (en longitudes de onda) =',
a);
b = input ('');
L (a) = b;
a = a + 1;
end
end
%*************************Longitud del
REFLECTOR***************************
LREFL = input (' Ingrese la longitud de los reflectores del arreglo
(enlongitudes de onda): ');
LCor = LREFL.*(ones(1,2*NN));

18
 FIEE-UNMSM 19

%***********************Longitud del ELEMENTO

ACTIVO***********************
LDRIVEN = input (' Ingrese la longitud del elemento activo (en
longitudesde onda): ');

L = [L LCor LREFL LDRIVEN];

NN=2*NN;
%*********************************************************************
*****
% INGRESO DE DISTANCIAS DE SEPARACION ENTRE ELEMENTOS
%*********************************************************************
*****
%*******Distancia de separacion entre 1er Director y elemento
activo*******
b = input ('\n Separacion entre elemento activo y 1er director = ',
's');
b = str2num (b);
D_1 = b;
%**********************Separacion entre
Directores*************************
if (N > 3)
fprintf (1, '\n Es la separacion entre directores uniforme?\n');
RESP = input (' RESPUESTA: (Y OR N) ...... ', 's');
else
RESP = 'N';
end
if (RESP == 'Y') | (RESP == 'y');
DDIR = input (' La separacion uniforme entre directores (en
longitudesde onda) = ', 's');
DDIR = str2num (DDIR);
D = DDIR * ones (1, N-2);
else (RESP == 'N') | (RESP == 'n');
a = 2;
while a <= (N-2)
fprintf (1, ' Separacion entre director# %2d y # %2d (en longitudesde
onda) =', a-1, a);
b = input (' ', 's');
b = str2num (b);
D (a) = b;
a = a + 1;
end
end
D (1) = D_1;
%*****************Distancia entre elemento activo y
reflector**************
b = input ('\n Separacion entre reflector y elemento activo (en
longitudesde onda) = ', 's');
b = str2num (b);
D (N-1) = b;
%*********************************************************************
*****
% INGRESO DE LONGITUD DE LOS ELEMENTOS EN EL ARREGLO Corner Reflector
%*********************************************************************
*****
DD = input('\n Ingrese la distancia de separacion entre reflectores
(en longitudes de onda): ');
%*********************************************************************
*****
% Radio de los elementos

19
 FIEE-UNMSM 20

%*********************************************************************
*****
b = input ('\n Radio para todos los elementos usados (en longitudes de
onda) = ', 's');
b = str2num (b);
radio = b;
%*********************************************************************
*****
%* Vector YPr = Vector de distancias de los reflectores de esquina *
%* sobre el eje y. *
%*********************************************************************
*****
SS=zeros(1, (NN/2) );
SS=1:NN/2;
SS=SS.*DD;
YPr=SS.*cos(psi);
YPr=[(YPr-D(N-1)),(YPr-D(N-1))];
%*********************************************************************
*****
%* Vector YP = Vector de distancias de todos los elementos *
%* sobre el eje y. *
%*********************************************************************
*****
YP=zeros(1,(N+NN));
for m=1:(N-2)
if m==1
YP(m)=D(m);
else
YP(m)=D(m)+YP(m-1);
end
end
YP(N+NN)=0;
YP(N+NN-1)=-D(N-1);
YP(1,(N-1):(N+NN-2))=YPr;
%*********************************************************************
*****
%* Vector XPr = Vector de distancias de los reflectores de esquina *
%* sobre el eje x. *
%* Vector XP = Vector de distancias de todos los elementos *
%* sobre el eje x. *
%*********************************************************************
*****
XP=zeros(1,(N+NN));
XPr=SS.*sin(psi);
XPr=[XPr,-XPr];
XP(1,(N-1):(N+NN-2))=XPr;
%*********************************************************************
*****
%* Matriz fracLM = Distancia Z de cada segmento sobre cada elemento *
%* sobre el eje z. *
%*********************************************************************
*****
LM=zeros((N+NN),M);
fracLM = L/(2*M-1);
fracLM = fracLM';
for m=1:M
LM(:,m) = fracLM*(m-1);
end
%*********************************************************************
*****
20
 FIEE-UNMSM 21

%* Matriz Anm = Matriz con los coeficientes de las corrientes, *

%* es la solucion a la ecuacion (2.66) *
%*********************************************************************
*****
Nel=N;
N=Nel+NN;
Anm = zeros (M*N,M*N);
for N1 = 1:N
for M1 =1:M
for N2=1:N
for M2=1:M
if N1==N2;
alpha=radio;
else
alpha=0;
end
GAMMA= alpha^2 + (YP(N1)-YP(N2))^2 + (XP(N1) - XP(N2))^2;
UL=L(N2)/2;
Parte1 = quadl(@integralG2,0,UL);
G2 = sumatoriaG2(UL);
Anm( ((N1-1)*M + M1), ((N2-1)*M + M2) ) = (((2*M2-1) *pi...
/ L(N2)) * (-1)^(M2 + 1) * G2 + ...
((2*pi)^ 2 - (2*M2-1)^2 * pi^2 / L(N2)^2) * Parte1);

end
end
end
end
%*********************************************************************
*****
%* fracLM2 = Distancia entre z=0 y primer segmento superior en el *
%* dipolo activo. *
%*********************************************************************
*****
fracLM2 = L(N)/(2*M-1);
%*********************************************************************
*****
%* Vector s = Vector de los campos incidentes en cada punto de *
%* adaptacion de cada dipolo. *
%*********************************************************************
*****
s = [zeros(1,(M*N-M)), -i/(fracLM2*2), zeros(1,M-1) ]';
%*********************************************************************
*****
%* Vector Inm = Vector de las corrientes complejas de los dipolos en *
%* cada modo. *
%*********************************************************************
*****
Inm = Anm\s;
%*********************************************************************
*****
%* Vector Etheta = Vector del patron de radiacion del campo en el *
%* plano E. Cada valor del campo es realizado cada grado *
%*********************************************************************
*****
%* EXPLICACION: Para graficar el patron de radiacion en el plano *
%* paralelo a la tierra. Para este calculo se debe notar que: *
%* 0 < theta < pi; y Phi = pi/2 y Phi = 3pi/2. *
%* De esta forma cubrimos todo el plano theta. *
%* de theta se van a tomar cada grado. *
21
 FIEE-UNMSM 22

%* Solucion de la Ecuacion 2.67 *

%*********************************************************************
*****
phi=1;
while phi <= 3 ;
for theta=0:1:180
SumEN=0;
for n=1:N
AAA = exp(i*2*pi*YP(n)*sin(theta*pi/180)*sin(phi*pi/2) + ...
i*2*pi*XP(n)*sin(theta*pi/180)*cos(phi*pi/2));
SumEM=0;
for m=1:M
SumEM = Inm( (M*(n-1)+m) ,1 )*( Zmas(m,(theta*pi/180),n)...
+ Zmenos(m,(theta*pi/180),n) ) + SumEM;
end

SumEN = SumEM * AAA * L(n) + SumEN;

end
if phi==1
Etheta( 1, theta+1 ) = j*(3*10^8)*(4*pi*10^-7)/4 * ...
sin(theta*pi/180) * SumEN;
else
Etheta( 1, (theta+180) ) = j*(3*10^8)*(4*pi*10^-7)/4 *...
sin(theta*pi/180) * SumEN;
end
end
phi = phi*3;
end
%*********************************************************************
*****
%* Vector EthetaDB = Vector logaritmico en dB de Etheta *
%*********************************************************************
*****
AbsEtheta = abs(Etheta);
MaxAbsTheta=max(AbsEtheta);
for theta=1:360
if (AbsEtheta(theta)/MaxAbsTheta) > (10^-6)
EthetaDB(theta) = 20 * log10(AbsEtheta(theta)/MaxAbsTheta );
else
EthetaDB(theta) = -120;
end
end
%*********************************************************************
*****
%* Figura 1 representa el patron de radiacion en el plano E, *
%* en coordenadas polares y dB *
%*********************************************************************
*****
theta=0:1:359;
figure(1);
polar_dB(theta,EthetaDB,-40,0,10,'-')
%*********************************************************************
*****
%* FTBeplane = Valor de la ganancia frente espalda en el plano E *
%*********************************************************************
*****
FTBeplane=-(EthetaDB(271))
%*********************************************************************
*****
%* DB3Etheta = Ancho del haz en el plano E *
22
 FIEE-UNMSM 23

%* EXPLICACION: Se asume un patron simetrico, por lo tanto *

%* si se conoce el angulo de apertura de un lado se multiplica *
%* por 2 y se encuentra el ancho del haz total *
%*********************************************************************
*****
EthetaDB1=EthetaDB(1,92:181);
V=find(EthetaDB1 >= -3);
V=fliplr(V);
if EthetaDB1(V(1)) == -3
DB3Etheta=2*V(1);
else
DB3Etheta=2* ( V(1)- (EthetaDB1( V(1) ) + 3) / (EthetaDB1(V(1)) - ...
EthetaDB( V(1) + 92 )) );
end
DB3Etheta
%*********************************************************************
*****
%* Vector Etheta = Vector del patron de radiacion del campo en el *
%* plano H. Cada valor del campo es realizado cada grado *
%*********************************************************************
*****
%* EXPLICACION: Para graficar el patron de radiacion en el plano *
%* perpendicular a la tierra. Para este calculo se debe notar: *
%* theta = pi/2; y 0 < Phi < 2pi *
%* De esta forma cubrimos todo el plano phi. *
%* Se toma valores cada grado. *
%* Solucion de la Ecuacion 2.67 *
%*********************************************************************
*****
theta1=pi/2;
for phi1=0:1:359
SumEN=0;
for n=1:N
AAA = exp(i*2*pi*YP(n)*sin(theta1)*sin(phi1*pi/180) + ...
i*2*pi*XP(n)*sin(theta1)*cos(phi1*pi/180));
SumEM=0;
for m=1:M
SumEM = Inm( (M*(n-1)+m) ,1 )*( Zmas(m,(phi1*pi/180),n)...
+ Zmenos(m,(phi1*pi/180),n) ) + SumEM;
end

SumEN = SumEM * AAA * L(n) + SumEN;

end
Hphi( 1, phi1+1 ) = j*(3*10^8)*(4*pi*10^-7)/4 *...
sin(theta1) * SumEN;
end
%*********************************************************************
*****
%* Vector HphiDB = Vector logaritmico en dB de Hphi *
%*********************************************************************
*****
AbsHphi = abs(Hphi);
MaxAbsHphi=max(AbsHphi);
for phi1=1:360
if (AbsHphi(phi1)/MaxAbsHphi) > (10^-6)
HphiDB(phi1) = 20 * log10(AbsHphi(phi1)/MaxAbsHphi );
else
HphiDB(phi1) = -120;
end
end
23
 FIEE-UNMSM 24

%*********************************************************************
*****
%* Figura 2 representa el patron de radiacion en el plano H, *
%* en coordenadas polares y dB *
%*********************************************************************
*****
phi1=0:1:359;
figure(2);
polar_dB(phi1,HphiDB,-40,0,10,'-')
%*********************************************************************
*****
%* FTBhplane = Valor de la ganancia frente espalda en el plano H *
%*********************************************************************
*****
FTBhplane=-(HphiDB(271))
%*********************************************************************
*****
%* DB3Hphi = Ancho del haz en el plano H *
%* EXPLICACION: Se asume un patron simetrico, por lo tanto *
%* si se conoce el angulo de apertura de un lado se multiplica *
%* por 2 y se encuentra el ancho del haz total *
%*********************************************************************
*****
HphiDB1=HphiDB(1,92:181);
V=find(HphiDB1 >= -3);
V=fliplr(V);
if HphiDB1(V(1)) == -3
DB3Hphi=2*V(1);
else
DB3Hphi=2* ( V(1) + (HphiDB1( V(1) ) + 3) / (HphiDB1(V(1)) -...
HphiDB( V(1) + 92 ) ) );
end
DB3Hphi
%*********************************************************************
*****
%* Do = Valor de la directividad *
%* EXPLICACION: Se resuelve la Ecuacion 2.26, *
%* El Prad se integra sobre toda el area de una esfera imaginaria*
%* El Umax es la energia de radiacion en la direccion maxima, *
%* en nuestro caso en: theta = pi/2 y phi = pi/2 *
%*********************************************************************
*****
Umax=Umax(pi/2,pi/2);
Prad= dblquad(@energiaU,0,pi,0,2*pi,1e-3,@quadl);
Do=4*pi*Umax/abs(Prad);
Do=10*log10(Do)
%*********************************************************************
*****
%* Iactivo = repreenta el valor de lo corriente sobre todo *
%* el dipolo activo. *
%* Figura 3 = Grafico de la distribucion de corriente del Dipolo
activo*
%*********************************************************************
*****
dZ=L(N)/100;
Iactivo = zeros(1,51);
for O=1:51
for m=1:M
Iactivo(1,O)= Inm((M*(n-1)+m) , 1)*cos((2*m-1) *...
pi*dZ*(O-1)/L(N))+ Iactivo(1,O);
24
 FIEE-UNMSM 25

end
Iactivo2(1,O)=Iactivo(1,O);
Iactivo(1,O)=abs(Iactivo(1,O));
end
Iactivo=[ fliplr(Iactivo(1,2:O)), Iactivo];
figure(3);
plot( linspace(-L(N)/2 , L(N)/2 , 101),Iactivo);
xlabel('Distancia del Dipolo Activo en lambda');
ylabel('Magnitud de la corriente del Dipolo Activo');
title('Distribucion de Corriente del Dipolo Activo');
%*********************************************************************
*****
%* Iact = Suma de todas las corrientes de todos los modos en el *
%* dipolo activo *
%* Zin = Impedancia de entrada de la antena en sus terminales *
%*********************************************************************
*****
Iact=0;
for m=1:M
Iact= Inm((M*(N-1)+m) , 1) + Iact;
end
Iact
Zin =- 2 / (Iact)
%*********************************************************************
*****
%* Figura 4 = representacion grafica de la antena Yagi-Uda o *
%* reflector de esquina *
%*********************************************************************
*****
a=[0, 0];
figure(4)
plot3(a,[YP(N-1) YP(N)],a,'linewidth',2.2,'color',[1 0 0]); hold on;
plot3([XP(N-1) XP(Nel+NN/2-2) ],[YP(N-1) YP(Nel+NN/2-
2)],a,'linewidth',...
2.2,'color',[0 1 0]); hold on;
plot3([XP(N-1) XP(Nel+NN-2)] ,[YP(N-1) YP(Nel+NN-2)],a,'linewidth',...
2.2,'color',[0 1 0]); hold on;
for n=1:N
a=[XP(n) , XP(n)];
b=[YP(n) , YP(n)];
c=[L(n)/2 , -L(n)/2 ];
plot3(a,b,c,'linewidth',2);hold on;
end
grid on;
%********************************************************************
%********************* FIN PROGRAMA ********************************
%*********************************************************************

%********************************************************************
%***************** Funcion: integralG2 ****************************
%*********************************************************************
function y=integralG2(ZZ)
global N D L LM N1 N2 M1 M2 GAMMA
Rplus = sqrt( GAMMA + (LM(N1,M1)-ZZ).^2);
Rminus = sqrt( GAMMA + (LM(N1,M1)+ZZ).^2);
y = ((exp(-i*2*pi.*Rminus)./(Rminus)) + ((exp(-
i*2*pi.*Rplus)./(Rplus))))...
25
 FIEE-UNMSM 26

.*cos((2*M2-1)*pi.*ZZ./L(N2));
%*********************************************************************
%*********************************************************************

%*********************************************************************
%*************** Funcion: sumatoriaG2 ******************************
%*********************************************************************
function G2=sumatoriaG2(UL)
global GAMMA N1 M1 LM
Rminus = sqrt (GAMMA + (LM(N1,M1) - UL) ^ 2);
Rplus = sqrt (GAMMA + (LM(N1,M1) + UL) ^ 2);
G2 = exp (-j*2*pi* Rplus) / (Rplus) + exp (-j*2*pi* Rminus) / (
Rminus);
%*********************************************************************
%*********************************************************************

%*********************************************************************
%******************* Funcion: Zmas *********************************
%*********************************************************************
function y=Zmas( m , theta , n)
global L
y1= (((2*m -1)*pi)/ L(n) + 2*pi.*cos(theta)).*L(n)/2;
[W1 W2] = size(y1);
if W2 == 1
if y1 == 0
y=1;
else
y=sin(y1)./y1;
end
else
[W1,W2] = size(find(y1==0));
if W2 == 0
y=sin(y1)./y1;
else
W1=find(y1==0);
y1(1,W1)=1;
y=sin(y1)./y1;
y(1,W1)=1;
end

end
%*********************************************************************
%*********************************************************************

%*********************************************************************
%******************* Funcion: Zmenos ********************************
%*********************************************************************
function y=Zmenos( m , theta , n)
global L
y1= (((2*m -1)*pi)/ L(n) - 2*pi.*cos(theta)).*L(n)/2;
[W1 W2] = size(y1);
if W2 == 1
if y1 == 0
y=1; 26
 FIEE-UNMSM 27

else
y=sin(y1)./y1;
end
else
[W1,W2] = size(find(y1==0));
if W2 == 0
y=sin(y1)./y1;
else
W1=find(y1==0);
y1(1,W1)=1;
y=sin(y1)./y1;
y(1,W1)=1;
end

end
%*********************************************************************
%*********************************************************************

%*********************************************************************
%***************** Funcion: energiaU *******************************
%*********************************************************************
%Esta funcion plantea la funcion que se integra para determinar el
Prad----
function y =energiaU(theta,phi)
global M N YP L Inm
SumEN=0;
for n=1:N
AAA = exp(i*2*pi*YP(n).*sin(theta).*sin(phi));
SumEM=0;
for m=1:M
SumEM = Inm( (M*(n-1)+m) ,1 ).*( Zmas(m,theta,n)+ Zmenos(m,theta,n) )
+ SumEM;
end
SumEN = SumEM.*AAA.*L(n) + SumEN;
end
y =3.75*pi.*(abs(SumEN).^2).*(sin(theta).^3);
%*********************************************************************
%*********************************************************************

%*********************************************************************
%***************** Funcion: Umax *************************************
%*********************************************************************
function y=Umax(theta,phi)
global M N YP L Inm
SumEN=0;
for n=1:N
AAA = exp(i*2*pi*YP(n)*sin(theta)*sin(phi));
SumEM=0;
for m=1:M
SumEM = Inm( (M*(n-1)+m) ,1 )*( Zmas(m,(theta),n)...
+ Zmenos(m,(theta),n) ) + SumEM;
end

SumEN = SumEM * AAA * L(n) + SumEN;

end
y=3.75 * pi * (abs(SumEN).^2) * (sin(theta).^2);
%*********************************************************************
%********************************************************************* 27
 FIEE-UNMSM 28

5. Conclusiones
Usamos el método de los momentos aplicado a la antena lineal para poder determinar de
forma aproximada la distribución de corriente y el patrón de radiación. La simulación de
la antena se realizó con el programa de Matlab, el patrón de radiación calculado fue de
acuerdo a lo esperado Al mismo tiempo, se pueden realizar nuevas implementaciones y
mejoras al código cambiando las funciones base y las funciones de prueba.

28
 FIEE-UNMSM 29

6. Bibliografía
 J.D. Jackson, Electrodinámica cl´asica,2da edición, Jhon Wiley-Sons, España,
1980.
 V.V. Nikolski, Electrodin´amica y propagaci´on de ondas de radio, MIR, Moscú,
1980.
 J. L. Fernández, Contribución al estudio de antenas en las cercanías de cuerpos
conductores aplicando el método de los momentos y modelado por hilos.
Universidad Politécnica de Madrid, 1985.
 A. Zozaya. Caracterización de antenas lineales usando el método de los
momentos. Laboratorio de electromagnetismo aplicado. Universidad de
Carabobo.
 The Method of Moments: A Numerical Technique for Wire Antenna Design. By
W.D. Rawle. Smiths Aerospace
o http://www.highfrequencyelectronics.com/Archives/Feb06/HFE0206_Ra
wle.pdf
 http://home.ict.nl/~arivoors/

29