Dipolo finito
Páginas: 3 (619 palabras)
Publicado: 13 de marzo de 2012
--DIPOLO FINITO-----------
Modelo de corriente sinusoidal obtenido a partir de "abrir" una LdT
en circuito abierto y es vàlido para alambres infinitamente delgados.
>restart:t:=theta:p:=phi:Digits:=10:pi:=evalf(Pi):dBlim:=-30:
C:=299792458:eta:=C*4e-7*pi:
Longitud:=1.5: # Longitud Elèctrica del Dipolo.
Integral Analìtica para Obtener el Vector de Radiaciòn
Como lacorriente està dirigida en ^z, Ntrans tiene sòlo componente theta.
> Nth:=-sin(t)*simplify(int(Io*sin(k*(L/2-zz))*exp(I*k*zz*cos(t)),zz=0..L/2)+int(Io*sin(k*(L/2+zz))*exp(I*k*zz*cos(t)),zz=-L/2..0));
Definiciòn del Patròn de Campo
Se eliminan las costantes multiplicativas al Nth, la L se pasa a medidas elèctricas.
Se normaliza automàticamente usando una utilerìa, sin embargo si lalongitud elèctrica
es muy grande (aprox. Longitud>6.2 lambdas) habrìa que acotar màs el intervalo.
> Delta:=evalf(1*10^(-Digits+1)):Fs:=abs(eval(-(k/(2*Io))*Nth,[k=2*Pi/lambda,L=L_e*lambda]));
Opt:=Optimization[Maximize](eval(Fs,[L_e=Longitud,Pi=pi]),t=0..pi/2):
Lambda:=1/Opt[1];
F:=unapply(Lambda*eval(Fs,[L_e=Longitud,Pi=pi,t=t+Delta]),t):FdB:=(t)->piecewise(20*log10(F(t)+Delta)<dBlim,dBlim,20*log10(F(t)+Delta));
Los Requerimientos de Graficaciòn
> st1:='orientation=[-23,80],numpoints=1200,linestyle=DOT,axes=boxed':st2:='coords=spherical,orientation=[-140,75],numpoints=2000,linestyle=DOT,axes=boxed,
scaling=constrained':
st3:='thickness=3,color=black,axes=boxed,numpoints=200':
st4:='thickness=3,color=brown,axes=boxed,numpoints=200':po:='coords=polar,scaling=constrained':
gr:=proc(a,b)
plottools[rotate](plots[display](b,plottools[reflect](a,[[0,0],[1,0]])),pi/2) end proc:
deu:=proc(a,b) plots[display](a,b) end proc:radiales:=proc(FF,umin,umax)
plots[display]([seq(plottools[line]([0,0],
[cos(umin+h*(umax-umin)/12)*eval(FF,u=umin+h*(umax-umin)/12),
sin(umin+h*(umax-umin)/12)*eval(FF,u=umin+h*(umax-umin)/12)]),h=0..12)],...
Modelo de corriente sinusoidal obtenido a partir de "abrir" una LdT
en circuito abierto y es vàlido para alambres infinitamente delgados.
>restart:t:=theta:p:=phi:Digits:=10:pi:=evalf(Pi):dBlim:=-30:
C:=299792458:eta:=C*4e-7*pi:
Longitud:=1.5: # Longitud Elèctrica del Dipolo.
Integral Analìtica para Obtener el Vector de Radiaciòn
Como lacorriente està dirigida en ^z, Ntrans tiene sòlo componente theta.
> Nth:=-sin(t)*simplify(int(Io*sin(k*(L/2-zz))*exp(I*k*zz*cos(t)),zz=0..L/2)+int(Io*sin(k*(L/2+zz))*exp(I*k*zz*cos(t)),zz=-L/2..0));
Definiciòn del Patròn de Campo
Se eliminan las costantes multiplicativas al Nth, la L se pasa a medidas elèctricas.
Se normaliza automàticamente usando una utilerìa, sin embargo si lalongitud elèctrica
es muy grande (aprox. Longitud>6.2 lambdas) habrìa que acotar màs el intervalo.
> Delta:=evalf(1*10^(-Digits+1)):Fs:=abs(eval(-(k/(2*Io))*Nth,[k=2*Pi/lambda,L=L_e*lambda]));
Opt:=Optimization[Maximize](eval(Fs,[L_e=Longitud,Pi=pi]),t=0..pi/2):
Lambda:=1/Opt[1];
F:=unapply(Lambda*eval(Fs,[L_e=Longitud,Pi=pi,t=t+Delta]),t):FdB:=(t)->piecewise(20*log10(F(t)+Delta)<dBlim,dBlim,20*log10(F(t)+Delta));
Los Requerimientos de Graficaciòn
> st1:='orientation=[-23,80],numpoints=1200,linestyle=DOT,axes=boxed':st2:='coords=spherical,orientation=[-140,75],numpoints=2000,linestyle=DOT,axes=boxed,
scaling=constrained':
st3:='thickness=3,color=black,axes=boxed,numpoints=200':
st4:='thickness=3,color=brown,axes=boxed,numpoints=200':po:='coords=polar,scaling=constrained':
gr:=proc(a,b)
plottools[rotate](plots[display](b,plottools[reflect](a,[[0,0],[1,0]])),pi/2) end proc:
deu:=proc(a,b) plots[display](a,b) end proc:radiales:=proc(FF,umin,umax)
plots[display]([seq(plottools[line]([0,0],
[cos(umin+h*(umax-umin)/12)*eval(FF,u=umin+h*(umax-umin)/12),
sin(umin+h*(umax-umin)/12)*eval(FF,u=umin+h*(umax-umin)/12)]),h=0..12)],...
Leer documento completo
Regístrate para leer el documento completo.