next up previous   
Next: Antenna filiforme: esecuzione del programma Up: Introduzione Previous: Dipolo elettrico: esecuzione del programma
Intro Gen: Introduzione Generale  Home: Home page

DIAGRAMMA DI IRRADIAZIONE DI UNA ANTENNA FILIFORME A l/2

Listato del programma

Filiforme.m

% disegno del grafico della direttivitā di una antenna filiforme a lambda/2

 

clear all

close all

x=0.01:0.01:2*pi-0.01;

y=1.64*(cos((pi/2).*cos(x))./sin(x)).^2;

plot(x,y)

title('grafico: 1.64*(cos((\pi/2)*cos(x)/sin(x))^2')

zoom

clear all

phi=0:.1:2*pi; % č il TH del MATLAB

th=0:.1:pi;  % č pi/2-PHI del MATLAB

h=0;

for j=1:length(phi)

   for i=1:length(th)

      if sin(th(i))~=0

         h=h+1;

         R=1.64*(cos((pi/2)*cos(th(i)))/sin(th(i)))^2;

         [x(h),y(h),z(h)]=sph2cart(phi(j),pi/2-th(i),R);

      end

   end

end

clear R TH phi th RD

phi=-0.01:0.01:2*pi; % č il TH del MATLAB

th=pi/2;   % č pi/2-PHI del MATLAB

h=0;

for j=1:length(phi)

   for i=1:length(th)

      if sin(th(i))~=0

         h=h+1;

         RR=1.64*(cos((pi/2)*cos(th(i)))/sin(th(i)))^2;

         [xx(h),yy(h),zz(h)]=sph2cart(phi(j),pi/2-th(i),RR);

      end

   end

end

clear R TH phi th

phi=0;      % č il TH del MATLAB

th=-0.01:0.01:2*pi;   % č pi/2-PHI del MATLAB

h=0;

for j=1:length(phi)

   for i=1:length(th)

      if sin(th(i))~=0

         h=h+1;

         RRR=1.64*(cos((pi/2)*cos(th(i)))/sin(th(i)))^2;

         [xxx(h),yyy(h),zzz(h)]=sph2cart(phi(j),pi/2-th(i),RRR);

      end

   end

end

clear R TH phi th

phi=pi/2; % č il TH del MATLAB

th=-0.01:0.01:2*pi;   % č pi/2-PHI del MATLAB

h=0;

for j=1:length(phi)

   for i=1:length(th)

      if sin(th(i))~=0

         h=h+1;

         RRRR=1.64*(cos((pi/2)*cos(th(i)))/sin(th(i)))^2;

         [xxxx(h),yyyy(h),zzzz(h)]=sph2cart(phi(j),pi/2-th(i),RRRR);

      end

   end

end

 

clear R TH

TH=0:0.01:2*pi;

for i=1:length(TH)

   if sin(TH(i))~=0

      R(i)=1.64*(cos((pi/2)*cos(TH(i)))/sin(TH(i)))^2;

   end

end

 

figure

polar(TH,R);

title('grafico: 1.64*(cos((\pi/2)*cos(\theta))/sin(\theta))^2')

figure

plot3(x,y,z)

grid

xlabel('x')

ylabel('y')

zlabel('z')

axis equal

figure

plot(xx,yy)

xlabel('x')

ylabel('y')

axis equal

figure

plot(yyyy,zzzz')

xlabel('y')

ylabel('z')

axis equal

figure

plot(xxx,zzz)

xlabel('x')

ylabel('z')

axis equal


next up previous   
Next: Antenna filiforme: esecuzione del programma Up: Introduzione Previous: Dipolo elettrico: esecuzione del programma
Intro Gen: Introduzione Generale  Home: Home page
Vito Marinelli
13-8-2000

HyperCounter
BPath Contatore