Modul3_online_sft.asv

%unfinished block for Eplane%langsung gerak ke penugasan
%using spec B >>l=lambda*2.1:0.1:lambda*3
c=3e8;%agar bisa pakai Hz
f=3e9;%frekuensi kerjadalam Hz
lambda=c/f;%menentukan panjang gelombang
eta=377;%impedansi ruang bebas
I=1;%arus maksimum masuk antena
l=lambda*1.5;%panjang antena. bisa berubah
k=(2*pi)/lambda;%wavenuber
theta=(0:0.01:2*pi);%sudut elevasi
r=10*lambda;%radius spherical coordinate/ bisa dibilang magnitude.


%rumus Eplane. bisa diganti sin ke cos%rumus langsung ke penugasan
%ngerapihin menjadi 2 part hitungan Eplane
part1=(1j*eta*I*exp(-1j*k*r))/(2*pi*r);
part2=(cos((k*l/2).*cos(theta))-cos((k*l)/2))./sin(theta);
Eplane=part1.*part2;
%unfinished block for Eplane%

%unfinished array factor%
%spec B d=pi/2 and n=10,30,50
N=16;%jumlah array antena
d=lambda/2;%jarak antar elemen array/spacingnya lambda/2
beta=1*pi;%fase eksitasi
psi=k*d*cos(theta)+beta;
AF=(1/N)*((sin((N/2)*psi))./sin((1/2)*psi));
%unfinished array factor%

%ETotal radiasi kombinasi belum diconvert ke dB
Etotal=Eplane.*AF;

%Power pattern dan konversi magnitude ke dB
%rtheta sudah dikondisikan untuk entry plot 3D
rtheta=abs(Etotal);%bisa Eplane atau AF atau Etotal
Rtheta=10*log10(rtheta);

%Rnormalisasi
Rnorm=Rtheta-min(Rtheta);
%Rnorm perlu dinolkan di 1,1 sebelum plotting
Rnorm(1,1)=0;
%pembentukan matrix
phi_matrix=[];
for iterator=0:0.01:2*pi
    phi_matrix=[phi_matrix;0:0.01:2*pi];
end
theta_matrix=transpose(phi_matrix);
%bentuk matriks Rnorm
Rnorm_matrix=[];
for i_iterator=1:629
    for k_iterator=1:629
        Rnorm_matrix(k_iterator,i_iterator)=Rnorm(1,k_iterator);
    end
end
%pembentukan matrix


%konversi sphere-coord to cartesian
phi=(0:0.01:2*pi);
x=Rnorm_matrix .*cos(theta_matrix) .*cos(phi_matrix);
y=Rnorm_matrix .*cos(theta_matrix) .*sin(phi_matrix);
z=Rnorm_matrix .*sin(theta_matrix);

%plotting
%polarplot(Rtheta); ini sudah siap
figure(1);
polarplot(10*log10(abs(Eplane.*AF)));%tidak digabung konversinya ke dB dulu
figure(2);
surf(x,y,z,sqrt(x.^2+y.^2+z.^2),'EdgeColor','none','FaceColor','interp');
figure(3);