Radar Cross Section (RCS) - MATLAB Simulations for Radar Systems Design - page 558

plot(ax6(ii),ay6(ii),'>r');
line([xc ax6(ii)],[yc ay6(ii)]);
plot(ax6,ay6,'g')
M6(ii) = getframe;
end
grid
xlabel('Ex')
ylabel('Ey')
title('Electric Field Locus; case 4')
rho = sqrt(ax6.^2 + ay6.^2);
major_axis = 2*max(rho);
minor_axis = 2*min(rho);
aspect6 = 10*log10(major_axis/minor_axis)
alpha6 = (180/pi) * atan2(ay6(1),ax6(1))
Listing 11.4. MATLAB Program Ðrcs_sphere.mÑ
% This program calculates the back-scattered RCS for a perfectly
% conducting sphere using Eq.(11.7), and produces plots similar to Fig.2.9
% Spherical Bessel functions are computed using series approximation and
recursion.
clear all
eps = 0.00001;
index = 0;
% kr limits are [0.05 - 15] ===> 300 points
for kr = 0.05:0.05:15
index = index + 1;
sphere_rcs = 0. + 0.*i;
f1 = 0. + 1.*i;
f2 = 1. + 0.*i;
m = 1.;
n = 0.;
q = -1.;
% initially set del to huge value
del =100000+100000*i;
while(abs(del) > eps)
q = -q;
n = n + 1;
m = m + 2;
del = (11.*n-1) * f2 / kr-f1;
f1 = f2;
f2 = del;
del = q * m /(f2 * (kr * f1 - n * f2));
sphere_rcs = sphere_rcs + del;