Antenne.m

 

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

% e di un dipolo

% confronti

 

clear all

close all

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

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

yd=1.5*(sin(x)).^2;

plot(x,yd,x,y,'r')

legend('dipolo elementare','antenna filiforme');

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;

         RD=1.5*(sin(th(i)))^2;

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

         [xd(h),yd(h),zd(h)]=sph2cart(phi(j),pi/2-th(i),RD);

      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;

         RRD=1.5*(sin(th(i)))^2;

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

         [xxd(h),yyd(h),zzd(h)]=sph2cart(phi(j),pi/2-th(i),RRD);

      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;

         RRRD=1.5*(sin(th(i)))^2;

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

         [xxxd(h),yyyd(h),zzzd(h)]=sph2cart(phi(j),pi/2-th(i),RRRD);

      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;

         RRRRD=1.5*(sin(th(i)))^2;

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

         [xxxxd(h),yyyyd(h),zzzzd(h)]=sph2cart(phi(j),pi/2-th(i),RRRRD);

      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;

      RD(i)=1.5*(sin(TH(i)))^2;

   end

end

 

figure

polar(TH,R,'r');

hold on

polar(TH,RD);

hold off

legend('antenna filiforme','dipolo elementare');

figure

plot3(xd,yd,zd,x,y,z,'r')

grid

legend('dipolo elementare','antenna filiforme');

xlabel('x')

ylabel('y')

zlabel('z')

axis equal

figure

plot(xxd,yyd,xx,yy,'-r')

xlabel('x')

ylabel('y')

legend('dipolo elementare','antenna filiforme');

axis equal

figure

plot(yyyyd,zzzzd,yyyy,zzzz,'r')

xlabel('y')

ylabel('z')

legend('dipolo elementare','antenna filiforme');

axis equal

figure

plot(xxxd,zzzd,xxx,zzz,'r')

xlabel('x')

ylabel('z')

legend('dipolo elementare','antenna filiforme');

axis equal