Instituto Politécnico Nacional
Escuela Superior de Ingeniería Mecánica y Eléctrica
“Zacatenco”
Teoría de Radiadores Electromagnéticos.
Grupo: 5CV9
Integrantes:
Duran Sanabria Hugo
Marquez Rojas A. Alexis
Introducción.
Este trabajo tiene como propósito visualizar en gráficas los cálculos realizados de forma teórica
utilizando como herramienta el software MatLab. De esta forma se podrá tener una mayor
perspectiva de los resultados obtenidos y mejorar el entendimiento de la forma en la que radian las
antenas, en este caso un dipolo.
La radiación de un dipolo en el espacio libre en un plano perpendicular a la dirección del hilo del
dipolo es exactamente igual en todas direcciones, mientras que en el plano del dipolo radia con un
máximo en la dirección perpendicular al hilo y un mínimo en la dirección del hilo.
En el programa el usuario puede elegir la gráfica que quiere visualizar, ya sea una tridimensional o
una en el plano. Esto con la finalidad de que se pueda interpretar de mejor manera la forma en la
que la onda es radiada.
Gráfica en 3D Gráfica en 2D
Código.
clc;
clear;
disp ( 'Instituto Politécnico Nacional ' )
disp ( 'Integrantes: ** Duran Sanabria Hugo ' )
disp ( ' ** Marquez Rojas A. Alexis ' )
disp ( 'TEMA: DIAGRAMA DE RADIACIÓN DE UN DIPOLO' )
disp ( 'OPCIÓN DE GRAFICAS:' )
disp ( 'Para lambda presione 1' )
disp ( 'Para lambda/2 presione 2' )
disp ( 'Para lambda/4 presione 3' )
disp ( 'Para lambda/8 presione 4' )
disp ( 'Para 2*lambda presione 5' )
disp ( 'INGRESE OPCION :' )
op=input('=> ' );
if(op==1)
% Introducción de variables.
la=1; % lambda.
b=2*pi/la; % beta.
l=la; % longitud del dipolo.
fi=(0:.01:1)*2*pi; % fi.
teta=(0:.01:1)*pi; % theta.
disp ( '1:Grafica en 2D' )
disp ( '2:Grafica en 3D' )
disp ( 'INGRESE OPCION :' )
g=input('=> ' );
if(g==1)
% Gráfica en 2D.
E=abs((cos(b.*l./2.*cos(teta))-cos(b.*l./2))./sin(teta));
polar(teta,E)
title('PARA LAMBDA')
end
if(g==2)
% Gráfica en 3D.
[FI,TETA]=meshgrid(fi,teta);
E=abs((cos(b.*l./2.*cos(TETA))-cos(b.*l./2))./sin(TETA));
X=E.*sin(TETA).*cos(FI);
Y=E.*sin(TETA).*sin(FI);
Z=E.*cos(TETA);
S=surface(X,Y,Z,(abs(E)));
axis equal
axis off
lighting gouraud;
shading interp
title('PARA LAMBDA')
view(0,40)
rotate3d
end
end
if(op==2)
% Introducción de variables.
la=1; % lambda.
b=2*pi/la; % beta.
l=la/2; % longitud del dipolo.
fi=(0:.01:1)*2*pi; % fi.
teta=(0:.01:1)*pi; % theta.
disp ( '1:Grafica en 2D' )
disp ( '2:Grafica en 3D' )
disp ( 'INGRESE OPCION :' )
g=input('=> ' );
if(g==1)
% Gráfica en 2D.
E=abs((cos(b.*l./2.*cos(teta))-cos(b.*l./2))./sin(teta));
polar(teta,E)
title('PARA LAMBDA/2')
end
if(g==2)
% Gráfica en 3D.
[FI,TETA]=meshgrid(fi,teta);
% Creación del campo.
E=abs((cos(b.*l./2.*cos(TETA))-cos(b.*l./2))./sin(TETA));
X=E.*sin(TETA).*cos(FI);
Y=E.*sin(TETA).*sin(FI);
Z=E.*cos(TETA);
S=surface(X,Y,Z,(abs(E)));
axis equal
axis off
lighting gouraud;
shading interp
title('PARA LAMBDA/2')
view(0,40)
rotate3d
end
end
if(op==3)
% Introducción de variables.
la=1; % lambda.
b=2*pi/la; % beta.
l=la/4; % longitud del dipolo.
fi=(0:.01:1)*2*pi; % fi.
teta=(0:.01:1)*pi; % theta.
disp ( '1:Grafica en 2D' )
disp ( '2:Grafica en 3D' )
disp ( 'INGRESE OPCION :' )
g=input('=> ' );
if(g==1)
% Gráfica en 2D.
E=abs((cos(b.*l./2.*cos(teta))-cos(b.*l./2))./sin(teta));
polar(teta,E)
title('PARA LAMBDA/4')
end
if(g==2)
% Gráfica en 3D.
[FI,TETA]=meshgrid(fi,teta);
% Creación del campo.
E=abs((cos(b.*l./2.*cos(TETA))-cos(b.*l./2))./sin(TETA));
X=E.*sin(TETA).*cos(FI);
Y=E.*sin(TETA).*sin(FI);
Z=E.*cos(TETA);
S=surface(X,Y,Z,(abs(E)));
axis equal
axis off
lighting gouraud;
shading interp
title('PARA LAMBDA/4')
view(0,40)
rotate3d
end
end
if(op==4)
% Introducción de variables.
la=1; % lambda
b=2*pi/la; % beta.
l=la/8; % longitud del dipolo.
fi=(0:.01:1)*2*pi; % fi.
teta=(0:.01:1)*pi; % theta.
disp ( '1:Grafica en 2D' )
disp ( '2:Grafica en 3D' )
disp ( 'INGRESE OPCION :' )
g=input('=> ' );
if(g==1)
% Gráfica en 2D.
E=abs((cos(b.*l./2.*cos(teta))-cos(b.*l./2))./sin(teta));
polar(teta,E)
title('PARA LAMBDA/8')
end
if(g==2)
% Gráfica en 3D.
[FI,TETA]=meshgrid(fi,teta);
E=abs((cos(b.*l./2.*cos(TETA))-cos(b.*l./2))./sin(TETA));
X=E.*sin(TETA).*cos(FI);
Y=E.*sin(TETA).*sin(FI);
Z=E.*cos(TETA);
S=surface(X,Y,Z,(abs(E)));
hold on
axis equal
axis off
lighting gouraud;
shading interp
title('PARA LAMBDA/8')
view(0,40)
rotate3d
end
end
if(op==5)
% Introducción de variables.
la=1; % lambda.
b=2*pi/la; % beta.
l=2*la; % longitud del dipolo.
fi=(0:.01:1)*2*pi; % fi.
teta=(0:.01:1)*pi; % theta.
disp ( '1:Grafica en 2D' )
disp ( '2:Grafica en 3D' )
disp ( 'INGRESE OPCION :' )
g=input('=> ' );
if(g==1)
% Gráfica en 2D.
E=abs((cos(b.*l./2.*cos(teta))-cos(b.*l./2))./sin(teta));
polar(teta,E)
title('PARA 2*LAMBDA')
end
if(g==2)
% Gráfica en 3D.
[FI,TETA]=meshgrid(fi,teta);
E=abs((cos(b.*l./2.*cos(TETA))-cos(b.*l./2))./sin(TETA));
X=E.*sin(TETA).*cos(FI);
Y=E.*sin(TETA).*sin(FI);
Z=E.*cos(TETA);
S=surface(X,Y,Z,(abs(E)));
axis equal
axis off
lighting gouraud;
shading interp
title('PARA 2*LAMBDA')
view(0,40)
rotate3d
end
end