I.
PROCEDIMIENTO
alfa=2*pi*F*sqrt(((E*U)/2)*((sqr
t(1+(S/(2*pi*F*E))^2))-1));
1. para comprender y analizar el
comportamiento de las ondas betha=2*pi*F*sqrt(((E*U)/2)*((sq
electromagnéticas es necesario conocer rt(1+(S/(2*pi*F*E))^2))+1));
sus características principales razón por la
cual desarrollamos una función en velocidad=2*pi*F/betha;
MATLAB que a partir de las longitud=(2*pi)/betha;
características del medio (ε, μ, σ) y la end
frecuencia calcule la constante de
if ( S~=0 && E==10e-9/(36*pi) &&
atenuación, fase, velocidad y longitud de
U~=12.6e-7)
onda). la función debe indica que tipo de medio=4;
medio es. (Espacio libre, dieléctrico, disp('Buen Conductor');
dieléctrico disipativo o buen conductor). alfa = sqrt((2*pi*F*U*S)/2);
betha = alfa;
Acontinuacion se muestra el código o velocidad=2*pi*F/betha;
longitud=(2*pi)/betha;
la función realizada para calcular los
end
parámetros descritos en el inciso 1 de este end
laboratorio.
En el código se muestra la función la
Código #1 cual por medio de compraciones muestra
que tipo de medios es . ademas
function[medio,alfa,betha,veloci
dad,longitud ] =
dependiendo de que medio es esta
calcmedio(E,U,S,F) función calcula los parámetros requidores
en el inciso uno utilizando las ecuaciones
if (S==0 && E==10e-9/(36*pi) && características para cada medio.
U==12.6e-7) Un aspecto importante que cabe resaltar
medio=1; es que para ingresar los valores se
disp('Espacio Libre'); insertan por medio de scrib tal y como se
alfa=0;
betha=(2*pi*F)/3e8;
explicara en los procedimientos
velocidad=3e8; posteriores.
longitud=(2*pi)/betha;
end
2. Utilizando la función que realizamos
if ( S==0 && E~=10e-9/(36*pi) && en el inciso anterior desarrollamos otra
U~=12.6e-7) función en MATLAB que dibuje la
medio=2;
disp('Dielectrico sin
intensidad del campo eléctrico en función
perdidas'); de la distancia, utilice t como constante,
S=0; t=0. Expresando la distancia en
longitudes de onda, y el campo en V/m
alfa=0; (dibuje hasta 3λ). Una característica de
betha=2*pi*F*sqrt(U*E); esta función es que debía pedir al
usuario que ingresara la amplitud y las
velocidad=2*pi*F/betha; características de medio.
longitud=(2*pi)/betha;
end
Código #2
function [I] =
if (S~=0 && E~=10e-9/(36*pi) &&
punto2(medio,alfa,betha,longitud
U~=12.6e-7)
,tiempo,frecuencia,amplitud,dist
medio=3;
ancia)
disp('Dielectrico con
if (medio==1 )
perdidas');
que medios se utilizan las mismas
I=amplitud.*cos((2*pi*frecuencia comparaciones pero anexando una
*tiempo)-(betha.*distancia));
end
adiociional que me identificara de que tipo
if (medio==2) de medio nos estamos referiendo y a partir
se inserta la formula de la intensidad de
I=amplitud.*cos((2*pi*frecuencia campo para cada medio en especifico.
*tiempo)-(betha.*distancia));
end
3. en esta oportunidad Desarrollamos otra
if (medio==3)
I=amplitud.*exp(-
función en MATLAB que dibujaba la
alfa.*distancia).*cos((2*pi*frec intensidad del campo eléctrico en función
uencia*tiempo)- del tiempo, utilizamos z como constante,
(betha.*distancia)); z=[Link] el tiempo en periodos T,
end y el campo en V/m (dibujamos hasta 3T)
if(medio==4) debimaos solicitar al usuario E0 y la
I=amplitud.*exp(-
alfa.*distancia).*cos((2*pi*frec
frecuencia.
uencia*tiempo)-
(betha.*distancia)); Código #3
end
figure; function [tiempo,periodo] =
title('Intensidad del campo punto3(frecuencia,amplitud)
electrico en funcion de la periodo=1/frecuencia;
distanci') tiempo=0:0.01*periodo:3*periodo;
plot(distancia/longitud,I); I=amplitud.*cos((2*pi*frecuencia
end .*tiempo));
figure;
acontinuacion se muestran los plot(tiempo/periodo,I);
resultados obtenidos al title('Intensidad del campo
ejecutar el codigo 2.ç electrico en funcion del
tiempo');
end
Fig1: grafica del campo eléctrico en función de la
distancia. Fig2 : grafica del campo eléctrico en función del
timepo.
en el código 2 se pretende mostrar la
forma de la grafica de la intensidad del Pricipalmente en la grafica se puede
campo eléctrico en función de la distancia observar el campo eléctrico graficado
utilizando como base principal la hasta su tercer periodo que es lo que nos
función descrita en el punto1 , esto hace piden el insico 3 de este laboratorio.
con la finalidad de simplicar las
En cuanto a la función podemos decir que
operaciones arealizar. Cabe resaltar que
en esencia es la misma que la resializada
para no realizar las comparaciones de
en punto 2 ya que en esta oportunidad
solomente se modifican los parámetros para
que este pueda ser graficada hasta el tercer
periodo. Ademas cabe resaltar que esta
señal tiene diferente frecuencia.
4. a continuación Utilizando las funciones
anteriores dibujamos en 3D la intensidad
del campo eléctrico en función de z y t. y
expresamos la distancia en longitudes de
onda (hasta 3λ), expresamos el tiempo en
periodos T (hasta 3T) y el campo en V/m. Fig3: grafico en 3D del campo eléctrico.
Código #4 En el código 4 se muestra esencialmente
una combinación de todas las funciones
function [Intensidad] = anteriores solo que en esta ocasión se
punto4(medio,x,t,longitud,period grafica en tercera dimensión y en función
o,frecuencia,betha,alfa,amplitud
)
del tiempo y la distancia para comprender
mejor el comportamiento ondulatorio de
[X,T]=meshgrid(x,t); las ondas electromagnéticos.
if (medio==1 )
A continuación se hace la explicación de
Intensidad=amplitud.*cos((2*pi*f la forma como se insertan los parámetros
recuencia*T)-(betha.*X)); en el script que se diseñó para mejorar la
end
if (medio==2)
sencillez de los temas realizados.
Intensidad=amplitud.*cos((2*pi*f %%punto 1 (aqui se ingresan las
recuencia*T)-(betha.*X)); caracteristicas del medio
end E= 1* 10e-9/(36*pi);
U= 1* 12.6e-7 ;
S= 0;
if (medio==3) F= 1000;
Intensidad=amplitud.*exp(-
[medio,alfa,betha,velocidad,long
alfa.*X).*cos((2*pi*frecuencia*T
itud ] = calcmedio(E,U,S,F);
)-(betha.*X));
%%punto 2
end
amplitud=1;
if(medio==4) distancia=0:0.01*longitud:3*long
Intensidad=amplitud.*exp(-
itud;
alfa.*X).*cos((2*pi*frecuencia*T
tiempo=0;
)-(betha.*X));
E= 1* 10e-9/(36*pi);
end
U= 1* 12.6e-7 ;
figure; S= 0;
mesh(X/longitud,T/periodo,Intens
F= 1000;
idad);
[medio,alfa,betha,velocidad,long
title('Grafica de la intensidad
itud ] = calcmedio(E,U,S,F);
del campo electrico en funcion
[I] =
del tiempo y la distancia'); punto2(medio,alfa,betha,longitud
zlabel('V/m');
,tiempo,F,amplitud,distancia);
xlabel('longitud de onda');
%%punto 3
ylabel('Periodo');
amplitud=1;
end
frecuencia=1;
[tiempo,periodo] =
punto3(frecuencia,amplitud);
%%punto 4
E= 1* 10e-9/(36*pi);
U= 1* 12.6e-7 ;
S= 0;
F= 1000;
amplitud=1;
[medio,alfa,betha,velocidad,long
itud ] = calcmedio(E,U,S,F);
distancia=0:0.01*longitud:3*long
itud;
[Intesidad] =
punto4(medio,distancia,tiempo,lo
ngitud,periodo,F,betha,alfa,ampl
itud);
Basicamente en el scrib se insertaran
todos los datos necesrios para poder
correr todas las funciones. En el scrip
porm memeio de comentarios se explica
como insertar cada parametro
