rcsflta.m

来自「此程序给出了不同情况的雷达距离测量仿真平台」· M 代码 · 共 64 行

% RCS of a Flat Plate
% -------------------

clear;clc;clf;

% Flat Plate Parameters

l=10;		% Length - m
w=.15;		% Width - m
f=[10 3 1 .1];	% Frequency - GHz

wlx=3e8./(f*1e9);	% Wavelength - m

% Angle from Perpendicular

phi=0:.01:90;
phir=phi*pi/180+1e-12;

% Find RCS of Flat Plate

A=l*w;
for i=1:4;
wl=wlx(i);
sigo=4*pi*A^2./wl.^2;
k=2*pi./wl;
fac=k*l*sin(phir);
z=sin(fac)./fac;
RCS=sigo*(z.^2).*(cos(phir).^2);
RCSX(:,i)=RCS';
end;

% Plot RCS of Flat Plate

subplot(221);
semilogy(phi,RCSX(:,1));grid;
axis([0 30 1e-4 1e6]);
xlabel('Angle - deg');
ylabel('RCS - m^2');
title('Flat Plate RCS')
text(15,300,['f = ',num2str(f(1)),' GHz']);	

subplot(222);
semilogy(phi,RCSX(:,2));grid;
axis([0 30 1e-4 1e4]);
xlabel('Angle - deg');
ylabel('RCS - m^2');
title('Flat Plate RCS');
text(15,10,['f=  ',num2str(f(2)),' GHz']);	

subplot(223)
semilogy(phi,RCSX(:,3));grid;
axis([0 30 1e-4 1e3]);
xlabel('Angle - deg');
ylabel('RCS - m^2');
title('Flat Plate RCS');
text(15,10,['f=  ',num2str(f(3)),' GHz']); 

subplot(224)
semilogy(phi,RCSX(:,4));grid;
axis([0 30 1e-4 1e2]);
xlabel('Angle - deg');
ylabel('RCS - m^2');
title('Flat Plate RCS');
text(15,10,['f=  ',num2str(f(4)),' GHz']);