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

end
rcs(index) = abs(sphere_rcs);
sphere_rcsdb(index) = 10. * log10(rcs(index));
end
figure(1);
n=0.05:.05:15;
plot (n,rcs,'k');
set (gca,'xtick',[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15]);
%xlabel ('Sphere circumference in wavelengths');
%ylabel ('Normalized sphere RCS');
grid;
figure (2);
plot (n,sphere_rcsdb,'k');
set (gca,'xtick',[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15]);
xlabel ('Sphere circumference in wavelengths');
ylabel ('Normalized sphere RCS - dB');
grid;
figure (3);
semilogx (n,sphere_rcsdb,'k');
xlabel ('Sphere circumference in wavelengths');
ylabel ('Normalized sphere RCS - dB');
Listing 11.5. MATLAB Function Ðrcs_ellipsoid.mÑ
function [rcs] = rcs_ellipsoid (a, b, c, phi)
% This function computes and plots the ellipsoid RCS versus aspect angle.
% The roll angle phi is fixed,
eps = 0.00001;
sin_phi_s = sin(phi)^2;
cos_phi_s = cos(phi)^2;
% Generate aspect angle vector
theta = 0.:.05:180.0;
theta = (theta .* pi) ./ 180.;
if(a ~= b & a ~= c)
rcs = (pi * a^2 * b^2 * c^2) ./ (a^2 * cos_phi_s .* (sin(theta).^2) + ...
b^2 * sin_phi_s .* (sin(theta).^2) + ...
c^2 .* (cos(theta).^2)).^2 ;
else
if(a == b & a ~= c)
rcs = (pi * b^4 * c^2) ./ ( b^2 .* (sin(theta).^2) + ...
c^2 .* (cos(theta).^2)).^2 ;
else
if (a == b & a ==c)
rcs = pi * c^2;