Electronica

COMPUTACION NUMERICA EN TEORIA ELECTROMAGNETICA
DAVID ANTONIO FUENTES VARGAS 2100161 GRUPO H1

Se presenta el programa de simulación de un resonador metálico en dos dimensiones. Se simula elcomportamiento de una fuente de onda electromagnética concentrada en un solo punto contenida en una caja de material conductor, intentando recrear su comportamiento lo más real posible.

%El programasimula el comportamiento de una fuente de onda concentrada %en un punto al interior de una frontera conductora
%Resonador metálico 2D TMz
%David Antonio Fuentes Vargas 2100161
%UNIVERSIDAD INDUSTRIALDE SANTANDER E3T
%Definición de constantes
muo=4*pi*1e-7;
epso=8.85*1e-12;
c0=1/sqrt(muo*epso);
%Definicion de mallado
fc=2.0e9;
fmax=3.0e9;
lamdamin=c0/fmax;
dx=lamdamin/20; %paso espacialdy=lamdamin/20;
dt=0.99*(dx/(c0*1.4142));%criterio de estabilidad
%Estado inicial de los campos
Ez=zeros(101,81);
Hx=zeros(101,81);
Hy=zeros(101,81);
is=51;
js=41;
imagesc(Ez)
%Definicion delas constantes de actualización
c1=dt/(muo*dx);
c2=dt/(epso*dx);
%****************************
%Loop actualización de campos
%****************************
for n=1:1000
for i=2:100for j=1:80
Hx(i,j)=Hx(i,j)-c1*(Ez(i,j+1)-Ez(i,j));
end
end
for i=1:100
for j=2:80
Hy(i,j)=Hy(i,j)-c1*(Ez(i,j)-Ez(i+1,j));
end
end
for i=2:100for j=2:80
Ez(i,j)=Ez(i,j)+c2*(Hy(i,j)-Hy(i-1,j)-Hx(i,j)+Hx(i,j-1));
end
end
%Termino de fuente
source=0.001*sin(2*pi*fmax*n*dt);
Ez(is,js)=Ez(is,js)+source;
%Graficapause(0.1)
imagesc(Ez)
end



¿Qué representa la simulación?
La simulación puede representar una cavidad resonante que rebota ondas a una frecuencia específica provenientes de una fuente dealimentación. Cuanta más energía en forma de ondas entra a la cavidad desde la fuente, más se recombina y refuerza con las ondas que ya están dentro, incrementando su intensidad.

CONCLUSIONES
La...