% Directed Energy Power for Seeker Burn-Out
% -----------------------------------------

clear;clf;clc

% Input Range

R=500:10000;		% Range - m

% Input Parameters

D=3;			   	% Dish Diameter - m
eff=.7;				% Antenna Efficiency
G=23;			   	% Seeker Antenna Gain - dB
f=10;			   	% Frequency - GHz
rx=[0 1 4];			% Rain Rate - mm/hr
atm=.024;			% Atm Attenuation - db/km
atr=.037;			% Rain Attenuation - db/km
Pr=12;				% Mixer Burn-Out - w 
Lx=2;			   	% Seeker Front End Loss - dB
Ly=3;			   	% DEW Polarization Loss - dB
Lx=10^(Lx/10);Ly=10^(Ly/10);

G=10^(G/10);
wl=3e8/(f*1e9);   % Wavelength - m

% Calculate Power Level

for i=1:3;
r=rx(i);
Gt=eff*(pi*D/wl).^2;
Lt=10^(-5*(atm+(r*atr)));
Pt=Pr*(4*pi*R).^2*Lt*Lx*Ly/(Gt*G*wl^2);
Ptx(:,i)=Pt';
end;

% Plot Power vs Range

plot(R,Ptx);grid;
xlabel('Range - m');
ylabel('Transmitter Power - w');
title('DEW Power for Seeker Burnout'); 
text(7500,3.25e7,'r=4 mm/hr');
text(7500,1.25e7,'r=1 mm/hr');
text(8000,.25e7,'clear');