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

for theta = 0.0:.1:90-.1
index = index +1;
theta = theta * pi /180.;
rcs(index) = (lambda * a1 * sin(theta) / ...
(8 * pi * (cos(theta))^2)) + eps;
end
theta*180/pi;
theta = pi/2;
index = index +1;
rcs(index) = (2 * pi * h^2 * a1 / lambda )+ eps;
for theta = 90+.1:.1:180.
index = index + 1;
theta = theta * pi / 180.;
rcs(index) = ( lambda * a1 * sin(theta) / ...
(8 * pi * (cos(theta))^2)) + eps;
end
r = a1;
index = 0;
for aspect_deg = 0.:.1:180
index = index +1;
aspect = (pi /180.) * aspect_deg;
% Compute RCS using Eq. (11.37)
if (aspect == 0 | aspect == pi)
rcs_po(index) = (4.0 * pi^3 * r^4 / lambda^2) + eps;
rcs_mu(index) = rcs_po(1);
else
x = (4. * pi * r / lambda) * sin(aspect);
val1 = 4. * pi^3 * r^4 / lambda^2;
val2 = 2. * besselj(1,x) / x;
rcs_po(index) = val1 * (val2 * cos(aspect))^2 + eps;
end
end
rcs_t =(rcs_po + rcs);
angle = 0:.1:180;
plot(angle,10*log10(rcs_t(1:1801)),'k');
grid;
xlabel ('Aspect angle -degrees');
ylabel ('RCS -dBsm');
Listing 11.20. MATLAB Program ÐSwerling_models.mÑ
% This program computes and plots Swerling statistical models
% sigma_bar = 1.5;
clear all