dipole.m

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

% Radar Cross Section of Dipole
% -------------------------------

clc;clf;clear;

l=.03/4;  			 % Half Wave Resonance in X-band
Sig=[];Sigf=[];

% Vary angle aspect

M=8; 
for Ang=pi*(1/M:1/M:1-1/M);
fGhz=[];

	% Inner loop varies frequency,kl
	% and converts frequency to Ghz

	for fMhz=5000:100:50000;
	wl=300/fMhz;
	k=2*pi/wl;
	kl=k*l;
	fGhz=[fGhz;fMhz/1000]; fMhz=[];
	
   % Evaluation of Integral as function of kl
   
	N=32;ang=pi*(1/N:1/N:1-1/N);I=0;
	   for n=1:N-1;
	   I=I+pi/N*((cos(kl*cos(ang(n)))-cos(kl)).^2)/sin(ang(n));
	   end;

	G=2*((cos(kl*cos(Ang))-cos(kl))/sin(Ang))^2/I;
	sig=G^2*wl^2/pi;
	Sig=[Sig;sig];sig=[];
        end;

   % Convert for angle aspect
        
	Sig=Sig';
	Sigf=[Sigf;Sig];Sig=[];
end;

% Correction for phase between I1 and I2
% and polarization losses

korr=.5*(cos(4*pi*fGhz*l/.3)).^2;

% Plot Mean RCS

plot(fGhz,korr'.*mean(Sigf));
xlabel('Frequency - Ghz');
ylabel('Radar Cross Section - sq.m');
title(['Dipole RCS for All Orientations and Polarizations']);
grid;