antenna theory balanis
thankx. but
I need just matlab code of example 8.3 8.4.
It seems to easy to code.
but i can't it.
-------[my matlab code ]-------------------------------------------------------------
% Moment Method for Line Antennas
clear all;
%Initial Value
ramda=.001;
f=3*10^8/ramda;
w=2*pi*f;
N=21;
l=0.47*ramda;
a=0.005*ramda;
b=a*2.3; % Zc=50 ohm Example 8.3 from balanis pp.449
delta=l/N;
e0=1/(36*pi)*10^-9;
u0=4*pi*10^-7;
k=w*sqrt(e0*u0);
Vs=1;
% Magnetic frill
n=1;
z=-l/2:l/(N-1):l/2;
R1=sqrt(z.^2+a^2);
R2=sqrt(z.^2+b^2);
for n=1:N;
Eiz(n,1)=-Vs/(2*log(b/a))*(exp(-j*k*R1(1,n))/R1(1,n)-exp(-j*k*R2(1,n))/R2(1,n));
end
%Poklinton integral eq.
n1=1;
for z=-l/2:l/(N-1):l/2
m1=1;
for z1=-l/2:l/(N-1):l/2
R=sqrt(a^2+(z-z1)^2);
Kz(n1,m1)=-1/(j*w*e0)*(exp(-j*k*R)/(4*pi*R^5)*((1+j*k*R)*(2*R^2-3*a^2)+(k*a*R)^2))*l/N;
m1=m1+1;
end
n1=n1+1;
end
% Result = [Iz]=inv([Kz])[Eiz]
Iz=inv(Kz)*Eiz;
plot(-l/2:l/(N-1):l/2,abs(Iz));
-------------------------------------------------------------------------------------------
Magnetic frill source modeling 'Eiz' are correct?
and Why N is bigger than 200 , change result of 'Iz' ?
use Moment method