Lab 1
Lab 1
Lab 1
AREQUIPA
AREQUIPA PERU
2018
1 Escriba una función en MATLAB que determine la impedancia característica de un dipolo de
longitud H ,longitud de onda lambda y diámetro de conductor a dado por la formula
Ho
2 Empleando la función generada en el problema anterior obtenga la gráfica de ZO vs
lambda
plot(ejex,R)
title('Impedancia Caracteriztica vs Ho/lambda')
xlabel('H0/lambda')
ylabel('Zo')
grid on
Impedancia Caracteriztica vs Ho/lambda
0
-50
-100
-150
Zo
-200
-250
-300
-350
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
H0/lambda
3 Obtenga la impedancia de una antena que opera a 75 MHz y que tiene una longitud de
H=3/4*lambda. Considerando que se emplea una línea de transmisión de impedancia
característica de 50 ohms grafique en la carta de Smith la curva que representa esta carga
desacoplada.
iantena =
105.4040
100
80
Ganancia Directiva
60
40
20
-20
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
H0/lambda
7.7 Grafique los patrones de radiación de las antenas lambda/2, lambda/4, 3lambda/4 en unplano
cartesiano en lugar de un plano polar
Para resolver este ejercicio usaremos la función radres ya definida en la pregunta 3, además
definimos la fucnion de ganacia directiva máxima “Gdmax” y la función intensidad de campo
“campo”
FUNCION GDMAX
function[G]=Gdmax(Rv,bet,H)
G=(120/Rv)*(1-cos(bet*H))^2;
end
FUNCION CAMPO
function [E,theta]=campo(R,beta,H)
theta=0:0.1:2*pi;
E=(60/R)*(cos((beta*H)/2)-cos(beta*H*cos(theta)/2))./sin(theta);
end
Calculo de patrones de radiación para las antenas requeridas.
f=100e6;
lambda=300e6/f;
beta=(2*pi)/lambda;
thet=0:0.001:2*pi;
%lambda/4
H=lambda/4;
R=radres(beta,lambda,H)
G=Gdmax(R,beta,H);
[E,thet]=campo(R,beta,H);
subplot(1,3,1);
plot(thet,E)
title('lambda/4')
grid on
%lambda/2
H2=lambda/2;
R2=radres(beta,lambda,H2)
G2=Gdmax(R2,beta,H2);
[E2,thet2]=campo(R2,beta,H2);
subplot(1,3,2);
plot(thet,E2)
title('lambda/2')
grid on
%3*lambda/4
H3=lambda*3/4;
R3=radres(beta,lambda,H3)
G3=Gdmax(R3,beta,H3);
[E3,thet3]=campo(R3,beta,H3);
subplot(1,3,3);
plot(thet,E3)
grid on
title('3*lambda/4')
0.2 0.8
0.3
0.15 0.6
0.2
0.1 0.4
0.1
0.05 0.2
0 0 0
-0.05 -0.2
-0.1
-0.1 -0.4
-0.2
-0.15 -0.6
-0.3
-0.2 -0.8
-0.25 -0.4 -1
0 2 4 6 8 0 2 4 6 8 0 2 4 6 8
7 Genera el patrón de radiación de 3D de una antena isotrópica.
disp ( 'OPCIN 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)
%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
g=2;
%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
if(op==2)
%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
g=2
%EN 3D
[FI,TETA]=meshgrid(fi,teta);
%creacion 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
if(op==3)
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
g=2;
%EN 3D
[FI,TETA]=meshgrid(fi,teta);
%creacion 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
if(op==4)
%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
g=2;
%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
if(op==5)
%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
g=2;
%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