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

end
% Plot RCS versus aspect angle
delta = 180 /index;
angle = 0.001:delta:180;
plot (angle,10*log10(rcs));
grid;
xlabel ('Aspect angle - degrees');
ylabel ('RCS - dBsm');
if(indicator ==1)
title ('Viewing from large end');
else
title ('Viewing from small end');
end
Listing 11.9. MATLAB Function Ðrcs_cylinder.mÑ
function [rcs] = rcs_cylinder(r1, r2, h, freq, phi, CylinderType)
% rcs_cylinder.m
% This program computes monostatic RCS for a finite length
% cylinder of either curricular or elliptical cross-section.
% Plot of RCS versus aspect angle theta is generated at a specified
% input angle phi
% Last modified on July 16, 2003
r = r1; % radius of the circular cylinder
eps =0.00001;
dtr = pi/180;
phir = phi*dtr;
freqGH = num2str(freq*1.e-9);
lambda = 3.0e+8 /freq; % wavelength
% CylinderType= 'Elliptic'; % 'Elliptic' or 'Circular'
switch CylinderType
case 'Circular'
% Compute RCS from 0 to (90-.5) degrees
index = 0;
for theta = 0.0:.1:90-.5
index = index +1;
thetar = theta * dtr;
rcs(index) = (lambda * r * sin(thetar) / ...
(8. * pi * (cos(thetar))^2)) + eps;
end
% Compute RCS for broadside specular at 90 degree
thetar = pi/2;
index = index +1;
rcs(index) = (2. * pi * h^2 * r / lambda )+ eps;