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

%('enter 1 to view frustum from large end, 0 otherwise')
large_small_end = indicator;
if(large_small_end == 1)
% Compute normal incidence, large end
normal_incedence = (180./pi) * ((pi /2) + alpha)
% Compute RCS from zero aspect to normal incidence
for theta = 0.001:.1:normal_incedence-.5
index = index +1;
theta = theta * pi /180.;
rcs(index) = (lambda * z1 * tan(alpha) *(tan(theta - alpha))^2) / ...
(8. * pi *sin(theta)) + eps;
end
%Compute broadside RCS
index = index +1;
rcs_normal = factor * sin(alpha) / ((cos(alpha))^4) + eps;
rcs(index) = rcs_normal;
% Compute RCS from broad side to 180 degrees
for theta = normal_incedence+.5:.1:180
index = index + 1;
theta = theta * pi / 180. ;
rcs(index) = (lambda * z2 * tan(alpha) *(tan(theta - alpha))^2) / ...
(8. * pi *sin(theta)) + eps;
end
else
% Compute normal incidence, small end
normal_incedence = (180./pi) * ((pi /2) - alpha)
% Compute RCS from zero aspect to normal incidence (large end of frustum)
for theta = 0.001:.1:normal_incedence-.5
index = index +1;
theta = theta * pi /180.;
rcs(index) = (lambda * z1 * tan(alpha) *(tan(theta + alpha))^2) / ...
(8. * pi *sin(theta)) + eps;
end
%Compute broadside RCS
index = index +1;
rcs_normal = factor * sin(alpha) / ((cos(alpha))^4) + eps;
rcs(index) = rcs_normal;
% Compute RCS from broad side to 180 degrees (small end of frustum)
for theta = normal_incedence+.5:.1:180
index = index + 1;
theta = theta * pi / 180. ;
rcs(index) = (lambda * z2 * tan(alpha) *(tan(theta + alpha))^2) / ...
(8. * pi *sin(theta)) + eps;
end