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

CORSO DI LAUREA IN INGEGNERIA ELETTRONICA

ESAME DI ANTENNE

ANNO ACCADEMICO 1999 / 2000

 PROGRAMMI IN MATLAB

 Prof. Ing.  M. Bozzetti

Eseguito da :

MARINELLI VITO

Matr. 507566 U


DIAGRAMMA DI IRRADIAZIONE DI UN DIPOLO ELETTRICO

Listato del programma

Dipolo.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;

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

plot(x,yd)

title('grafico: 1.5*(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;

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

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

      end

   end

end

clear 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;

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

         [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;

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

         [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;

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

         [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

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

   end

end

 

figure

polar(TH,RD);

title('grafico: 1.5*(sin(x))^2')

figure

plot3(xd,yd,zd)

grid

xlabel('x')

ylabel('y')

zlabel('z')

axis equal

figure

plot(xxd,yyd)

xlabel('x')

ylabel('y')

axis equal

figure

plot(yyyyd,zzzzd)

xlabel('y')

ylabel('z')

axis equal

figure

plot(xxxd,zzzd)

xlabel('x')

ylabel('z')

axis equal


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

HyperCounter
BPath Contatore