buoyjam.m

此程序给出了不同情况的雷达距离测量仿真平台

% Program Calculates J/S for Jamming with
% Decoys in Naval Scenario
% --------------------------------------- 

clc;clf;clear;

% Jamming Decoy Parameters

Pj1=100;Gj1=1; % Erp
Pj2=100;Gj2=1; % Erp
x1=5000;y1=100;% meters
x2=3000;y2=50; % meters
wl=.03;        % meters

% Ship at origin

Rcs=400;       % m^2

% Missile with cos^2 Antenna and Go=30 db
% Launch range of 30 km

Rmax=30000;Rmin=2000;
Put=40000;

% N Element Range Vector

N=1024; step=(Rmax-Rmin)/(N-1);
R=Rmax:-step:Rmin;

%  Angle Vectors fi to Ship

fi=zeros(1,N);

% Angle Vectors to Jamming Buoy #1

fi1=atan(y1./(R-x1));

% Angle Vectors to Jamming Buoy #2

fi2=atan(y2./(R-x2));

% Range to Jamming Buoy #1

R1=((R-x1).^2+y1^2).^(1/2);

% Range to Jamming Buoy #2

R2=((R-x2).^2+y2^2).^(1/2);

% Received Jamming Power from Buoy #1

Prj1=(Pj1*Gj1*gcos2(30,fi1)*wl^2)./...
((4*pi)^2*(R1).^2);

% Received Jamming Power from Buoy #2

Prj2=(Pj2*Gj2*gcos2(30,fi2)*wl^2)./...
((4*pi)^2*(R2).^2);

% Jam Power Dilution for 10 Mhz Band Width
% in 2 Mhz Receiver Band Width

Prj1=Prj1*2/10;
Prj2=Prj2*2/10;

% Zero Prj1 and Prj2 for R<=x1 and R<=x2

for n=1:N;
  if R(n)<=x1;Prj1(n)=0;end;
end;

for n=1:N;
  if R(n)<=x2;Prj2(n)=0;end;
end;

% Total Received Jamming Power

Prj=Prj1+Prj2;
Prjdbm=30+10*log10(Prj+1e-12);

% Recived Target Power from Ship

Prtar=(Put*gcos2(30,fi).^2*Rcs*wl^2)./...
((4*pi)^3*R.^4);
Prtardbm=30+10*log10(Prtar+1e-12);

% Plot Received Jammer and Target Power

plot(R/1000,Prjdbm,R/1000,Prtardbm);
title('Jam and Target Power for Buoy Jamming');
ylabel('Received Signal - dbm');
xlabel('Target Range - km');grid;