Graphics Programs Reference
In-Depth Information
alpha = alphad*dtr;
beta = betad *dtr;
switch reference
case 'top face'
alpha = 0;
vi = pi/(2*pi-beta);
case 'bisector'
beta = alpha;
vi = pi/(2*pi-2*beta);
case 'on x-axis'
vi = pi/(2*pi-alpha-beta);
end
phip = phipd*dtr;
etar = sqrt(mur/epsr);
mu = mu0*mur;
eps = eps0*epsr;
lambda = 2.99e8/freq;
k = 2*pi/lambda; % free space wavenumber
ka = k*ar;
k1 = k*sqrt(mur*epsr); % wavenumber inside dielectric
k1a = k1*ar;
krhop = k*rhop;
omega =2*pi*freq;
% <<< Far field Calculations of Ez component >>>
% === Line source excitation ===
Nc =round(1+2*k*rhop); % number of terms for series summation
Term = pi*omega*mu0/(2*pi-alpha-beta);
Term0D = img*4*pi/(2*pi-alpha-beta);
Term0C = -img*4*pi/(2*pi-alpha-beta);
Term0 = 4*pi/(2*pi-alpha-beta);
for ip = 1:360
phii = (ip -1)*dtr;
xphi(ip) = ip-1;
if phii > alpha & phii < 2*pi-beta % outside the wedge region
EzFLs(ip) = 0;
for m = 1:Nc
v = m*vi;
ssterm = (img^v)*sin(v*(phip-alpha))*sin(v*(phii-alpha));
switch CapType
case 'Diel'
Aterm = k * besselj(v,k1a)*(dbesselj(v,ka)*bessely(v,krhop)...
-dbessely(v,ka)*besselj(v,krhop)) ...
+k1*dbesselj(v,k1a)*( bessely(v,ka)*besselj(v,krhop)...
Search WWH ::
Custom Search