crosseye.m

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

% Field Pattern for Cross-Eye Jamming
% ------------------------------------

clf;clc;clear;

% System parameters

k=.8;

f0=9e+9;c=3e+8;t=0;

P1=400;G1=100;
P2=P1/k;G2=100;

% Antenna separation

a=20;

x1=-a/2;y1=0;
x2=a/2;y2=0;

% Field extent (dx,dy) about (x,y)

x=0;y=2000;
dx=8;dy=.1;  

% Number grid points along axis

n=50;

Ex=[];E=[];

% Define loops 

for yf=y-dy/2:dy/(n-1):y+dy/2+dy/(n-1)/2;
	for xf=x-dx/2:dx/(n-1):x+dx/2+dx/(n-1)/2;

% Compute magnitudes of sinal vectors

r1=sqrt((x1-xf)^2+(y1-yf)^2);
r2=sqrt((x2-xf)^2+(y2-yf)^2);

% Compute field density at (xf,yf);

p1=P1*G1/4/pi/(r1^2);
p2=P2*G2/4/pi/(r2^2);

% Compute field strengths at (xf,yf);

E1=sqrt(277*p1);
E2=sqrt(277*p2);

% Compute field strength with 180 deg phase

e1=E1*sin(2*pi*f0*(t-r1/c)+pi);

% Compute field strength at (xf,yf)

e2=E2*sin(2*pi*f0*(t-r2/c));

% Combine vector fields at (xf,yf)

e=e1+e2;

% Matrix in x-direction

Ex=[Ex;e];

end;

Ext=Ex';
Ex=[];


% Matrix in y-direction

E=[Ext;E];

end;

% Plot sum fields

kx=10*log10(k);
subplot(211),mesh(E,[-45 45],[1 1 1/2]);
title(['Jamming sep=20 m , rel eff=',num2str(kx),...
'dB ,wl=3 cm,  phase front at 2 km']);
 xlabel('x');ylabel('y');zlabel('z');

subplot(212),contour(E,1);
title(['Jamming sep=20 m , rel eff=',num2str(kx),...
'db ,wl=3 cm, phase front at 2 km']);
xlabel(['Total ',num2str(dx),' m']);
ylabel(['Total  ',num2str(dy*100),' cm']);