# MATLAB Simulations for Radar Systems Design

### Bassem Mahafza (view profile)

• 1 file
• 4.36765

11 Sep 2003 (Updated )

MATLAB Simulations for Radar Systems Design

rcs_frustum (r1, r2, h, freq, indicator)
```function [rcs] = rcs_frustum (r1, r2, h, freq, indicator)
% This program computes the monostatic RCS for a frustum.
% Incident linear Polarization is assumed.
% To compute RCP or LCP RCS one must use Eq. (2.24)
% When viewing from the small end of the frustum
% normal incedence occurs at aspect pi/2 - half cone angle
% When viewing from the large end, normal incidence occur at
% pi/2 + half cone angle.
% RCS is computed using Eq. (2.43). This program assumes a geomtry
% similar top Fig. 2.13
format long
index = 0;
eps = 0.000001;
lambda = 3.0e+8 /freq;
% Enter frustum's small end radius
%r1 =.02057;
% Enter Frustum's large end radius
%r2 = .05753;
% Compute Frustum's length
%h = .20945;
% Comput half cone angle, alpha
alpha = atan(( r2 - r1)/h);
% Compute z1 and z2
z2 = r2 / tan(alpha);
z1 = r1 / tan(alpha);
delta = (z2^1.5 - z1^1.5)^2;
factor = (8. * pi * delta) / (9. * lambda);
%('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
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 inicedence (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
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
end
% Plot RCS versus aspect angle
delta = 180 /index;
angle = 0.001:delta:180;
plot (angle,10*log10(rcs),'k','linewidth',1);
grid;
xlabel ('Apsect angle - degrees');
ylabel ('RCS - dBsm');
if(indicator ==1)
title ('Viewing from large end');
else
title ('Viewing from small end');
end

```