Ingenieria
Páginas: 3 (628 palabras)
Publicado: 6 de noviembre de 2012
CODIGO DEL CAMPO ELECTRICO EN EL PLANO HORIZONTAL XY DE UN DIPOLO DE HERTZ
%plano horizontal(xy);k=2*pi;lambda=2*pi/k;epsilon0=8.854187817*10.^(-12);mu0=4*pi*10.^(-7);f=1/(epsilon0*mu0).^0.5;w=2*pi*f;nn=70;x1=linspace(-5,5*lambda,nn);y1=linspace(-5,5*lambda,nn);[X,Y]=meshgrid(x1,y1);axis([-5 5.*lambda -5 5.*lambda]);for t=0:400; title('Campo Eléctrico Horizontal de un dipolo deHertz');Etha=377;dl=lambda/2;I0=1;%Coeficiente del campoeléctrico;A=Etha*I0*dl/(4*pi*epsilon0*k);z=0;theta=pi/2;r=(X.^2+Y.^2+z.^2).^0.5;phi=atan(Y./X);Ey=A.*(((-r.^2.*sin(w.*t-k.*r).*k.*sin(theta)-sin(theta).*sin(w.*t))./r.^3).*cos(theta).*sin(phi)+(2.*sin(w.*t).*cos(theta)./r.^3).*sin(theta).*sin(phi));Ex=A.*(((-r.^2.*sin(w.*t-k.*r).*k.*sin(theta)-sin(theta).*sin(w.*t))./r.^3).*cos(theta).*cos(phi)+(2.*sin(w.*t).*cos(theta)./r.^3).*sin(theta).*cos(phi));E=((Ex).^2+(Ey).^2).^0.5;Eyn=Ey./E;Exn=Ex./E;quiver(X,Y,Exn,Eyn);pause(0.01);end |
CODIGO DEL CAMPO ELECTRICO EN EL PLANOVERTICAL YZ DE UN DIPOLO DE HERTZ
%plano horizontal(YZ);k=2*pi;lambda=2*pi/k;epsilon0=8.854187817*10.^(-12);mu0=4*pi*10.^(-7);f=1/(epsilon0*mu0).^0.5;w=2*pi*f;nn=90;y1=linspace(-5,5*lambda,nn);z1=linspace(-5,5*lambda,nn);[Y,Z]=meshgrid(y1,z1);axis([-5 5.*lambda -5 5.*lambda]);for t=0:300; title('Campo Eléctrico vertical de un dipolo de Hertz');Etha=377;dl=lambda/2;I0=1;%Coeficientedel campoeléctrico;A=Etha*I0*dl/(4*pi*epsilon0*k);X=0;theta=atan(Y./Z);r=(Y.^2+Z.^2+X.^2).^0.5;phi=pi./2;Ey=A.*(((-r.^2.*sin(w.*t-k.*r).*k.*sin(theta)-sin(theta).*sin(w.*t))./r.^3).*cos(theta).*sin(phi)+(2.*sin(w.*t).*cos(theta)./r.^3).*sin(theta).*sin(phi));Ez=A.*(((-r.^2.*sin(w.*t-k.*r).*k.*sin(theta)-sin(theta).*sin(w.*t))./r.^3).*-sin(theta)+(2.*sin(w.*t).*cos(theta)./r.^3).*cos(theta));E=((Ey).^2+(Ez).^2).^0.5;Eyn=Ey./E;Ezn=Ez./E;quiver(Y,Z,Eyn,Ezn);pause(0.00001);end |
CODIGO DEL CAMPO ELECTRICO EN EL PLANO VERTICAL XZ DE UN DIPOLO DE HERTZ
%plano horizontal...
%plano horizontal(xy);k=2*pi;lambda=2*pi/k;epsilon0=8.854187817*10.^(-12);mu0=4*pi*10.^(-7);f=1/(epsilon0*mu0).^0.5;w=2*pi*f;nn=70;x1=linspace(-5,5*lambda,nn);y1=linspace(-5,5*lambda,nn);[X,Y]=meshgrid(x1,y1);axis([-5 5.*lambda -5 5.*lambda]);for t=0:400; title('Campo Eléctrico Horizontal de un dipolo deHertz');Etha=377;dl=lambda/2;I0=1;%Coeficiente del campoeléctrico;A=Etha*I0*dl/(4*pi*epsilon0*k);z=0;theta=pi/2;r=(X.^2+Y.^2+z.^2).^0.5;phi=atan(Y./X);Ey=A.*(((-r.^2.*sin(w.*t-k.*r).*k.*sin(theta)-sin(theta).*sin(w.*t))./r.^3).*cos(theta).*sin(phi)+(2.*sin(w.*t).*cos(theta)./r.^3).*sin(theta).*sin(phi));Ex=A.*(((-r.^2.*sin(w.*t-k.*r).*k.*sin(theta)-sin(theta).*sin(w.*t))./r.^3).*cos(theta).*cos(phi)+(2.*sin(w.*t).*cos(theta)./r.^3).*sin(theta).*cos(phi));E=((Ex).^2+(Ey).^2).^0.5;Eyn=Ey./E;Exn=Ex./E;quiver(X,Y,Exn,Eyn);pause(0.01);end |
CODIGO DEL CAMPO ELECTRICO EN EL PLANOVERTICAL YZ DE UN DIPOLO DE HERTZ
%plano horizontal(YZ);k=2*pi;lambda=2*pi/k;epsilon0=8.854187817*10.^(-12);mu0=4*pi*10.^(-7);f=1/(epsilon0*mu0).^0.5;w=2*pi*f;nn=90;y1=linspace(-5,5*lambda,nn);z1=linspace(-5,5*lambda,nn);[Y,Z]=meshgrid(y1,z1);axis([-5 5.*lambda -5 5.*lambda]);for t=0:300; title('Campo Eléctrico vertical de un dipolo de Hertz');Etha=377;dl=lambda/2;I0=1;%Coeficientedel campoeléctrico;A=Etha*I0*dl/(4*pi*epsilon0*k);X=0;theta=atan(Y./Z);r=(Y.^2+Z.^2+X.^2).^0.5;phi=pi./2;Ey=A.*(((-r.^2.*sin(w.*t-k.*r).*k.*sin(theta)-sin(theta).*sin(w.*t))./r.^3).*cos(theta).*sin(phi)+(2.*sin(w.*t).*cos(theta)./r.^3).*sin(theta).*sin(phi));Ez=A.*(((-r.^2.*sin(w.*t-k.*r).*k.*sin(theta)-sin(theta).*sin(w.*t))./r.^3).*-sin(theta)+(2.*sin(w.*t).*cos(theta)./r.^3).*cos(theta));E=((Ey).^2+(Ez).^2).^0.5;Eyn=Ey./E;Ezn=Ez./E;quiver(Y,Z,Eyn,Ezn);pause(0.00001);end |
CODIGO DEL CAMPO ELECTRICO EN EL PLANO VERTICAL XZ DE UN DIPOLO DE HERTZ
%plano horizontal...
Leer documento completo
Regístrate para leer el documento completo.