0% found this document useful (0 votes)
312 views

Radar Simulation in Matlab PDF

Radar simulation in matlab

Uploaded by

Muhammad Aljamal
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
312 views

Radar Simulation in Matlab PDF

Radar simulation in matlab

Uploaded by

Muhammad Aljamal
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 23

Radar Simulation Using Matlab

345 7849 7:;/ -' + :+

 
  

 !" #$%& '() : 

Radar Simulation Using Matlab

  
 

?'> 7> RADAR> E45>D > RAdio Detection And Ranging  IJ>K(>
Q;4 )O .7
>>4>>+ [KJ & 7>>'Z>> >>QY 4>>T ' 7>>XU4>>$W 'T 7>>U;>>4 >>4 )^(7>>]4
U) a) a$X) `$ a'b 7$Y$^4:cg ? a+ d(eO> d>$ . O> ] `>$ a'>b 7>Xd>(e
; mKX` ) a ki? 7]4j i ) 454 h'Y UK) .EY$^4:c q4X4 >5 '>7?eoY
?(5 4>)?X) h>?+ a>+>U >Xe>U v+>; ,)] v>X>` >$u> a>) s wW45>Q''`>$
?U.
' aK I:5& 7'Z :h
.1 .Ground based 7$b
.2 )(' 7 .Airborne
.3 )('.Space borne 4$W4{| 7
.4 (.Ship based radar 7
' a>>K >>>>+ I:5 >>4g$:5 >>D 4>>'+ wW45>>D h>>?+
 `(4  ,7& U4$W 4T ' .7)QY'>4
' aK: #Y I:5' 7'U  7g$ Eo  4U.
) ZX& 7'Z ['XY EW  $e' aK  | +4o  a>) >$c$ k4>%>4X >
 a>>K>> ">> I' >>[ >> ['XY4>>$W ">> I' >>[ 4>>) 4> e4u>> h'Y>> >>IW4
'X   a) IK) EW ) a$XJ) a .8 U4>$W "> I' >[ >) a>4X $>b
) v  vU$)j 4 7$K$&4K$) 7o K& .7$  v>$KT&>^ a>+ E> >(>j K 4>$
'ic  "4:X.I5
34>T
4) 4e4u I:5 #Y4T  E 4U)QY #Y >| .[$cJ 4&i>D
'| 7)QY  I:5:h
 4 )'Y ):(Continues Wave CW
E  E 7]4^ JU) 'Y) [KJ 7$Y$^4:c ['XY> [>;) EW >[eoY>) EW
) .a$?5g:'? aK i > ['XY4>T '>'Y  (Unmodulated CW) 7>X) >$c m$>o
]4X 7 ; 7+U ) 7?( )] vX`  3 a>K' a>K(>5 4>)?X) h>?+ >X a>+
U Q;4 4T  7>X) >$c>| E4>4 >i >: )> a> >( >o| ['XY> ;>7+
U IJ( | E v$io.IW
 4 &:(Pulse Radar) 7${e
E4j] ['XY ) a4{e: ' .7Xi : ) a& 7'Z ' a>K >4>'+34 I:5 h>?+
 K 4{e: ) (Pulse Repetition Frequency PRF:h
.A K &4{e ):(Low PRF) gQ:
E( ['XY  XU 7+Y U& 3 4):+U.

  !" #$%& '() :

+'+ : 7849 7:;/ -345

Radar Simulation Using Matlab

 
  

.(Medium PRF) ;) 4{e& K  .B


:(High PRF) E4+ 4{e& K  .C
.U 7+;  (E| ['XY
: ?7 (` a$e E4 %4| [$cJ  h?+ 4'+34 I$:5 7eY:4 4)

IJ>K $g>;:T 7>oe^ a'>b mKX:> E> 7$Y>$^4:c)U 4> 7$4X'[   XY
` 4':$ .(The Over Horizon TOH) (Z: D I?D) |O I?D o E $Xe O
>>> >>>$Xe 4|4Y>>>'>>> Ke' i>>>&T E>>>| (UHF) 7>>>oW4g( >>> VHF) >>>7>>>$4X
.(Early Warning System EWS)

345 7849 7:;/ -' + :+

 !" #$%& '() : 

Radar Simulation Using Matlab

  
 

( :Xe
' [8 [KJ E4' jQ:5]? E .E{e:

 ">: >( K`) E>:>4 >`) a`) E>:' a'>b 7>?j>i 4>Z: .'> 4> 4>T '7>X
;[  hU:" 7^4; EW 'X/';[  e[ >U a$> EW EX>bT;>4 4eo>;3 a>+
^ ;4 ;4eo )' 4Q;4 'Y & mgU EW;4 4eo;3 d$  v$4T !D
| T;4  hU| E| 4':$ EW 4eo;3 >  v>$4>T )>5 UK) (EY>$^4:c> h'>[eoY
i { Q>i >5 .7>%4X'? v>$U> (>5 h>?+'4>)?X >U4 7>"4Q ^ a>+> ">:
' d$ 7%4X( X #YU 4$o  $D`) t E>:>i > v]cY4{>e:  >jo'> a$> 7|4Y>
U 44 ) 44 a$X iDO4e+3 ; 7+&4J ' 7UK)9 C = 3 *108 m / sec 7$Y$^4:c
& ej: 7]!X
Ct
2

=R

4+ [KJ   E{e:;[ 4j] [eoY ) a4{e: ' 4 ) E| b [KJ:E4

  !" #$%& '() :

+'+ : 7849 7:;/ -345

Radar Simulation Using Matlab

 
  

4>) 4>e4u  7{e: +( Inter Pulse Period IPP) T E 4{e: ki a$ 7$:)` g
4{>e:> K > > v>?o)(Pulse Repetition Interval PRI) 4{>e:> K a)`> IPPh+>
(Pulse Repetition Frequency PRF)
1
fr =
T
U a) W4X 5 Z: o| a)` 7$Y$^4:c)U 7 )4;  o PRI a)`! D
:7eY: a) #Y( dt v `))[ 4X( ( 4;T )| 4;T .7$oe' g E|
dt =

:7]!X4 4UXJ E 7j;' 7+4j;3 hjX E44


Pav = Pt .dt

:7]!X4 hjX   7{e: 7+4j; %&  4UXJ E h'ZX 7+4j;3 E Pt d$


E p = Pt . = Pav .T =

Pav
fr

E>| 7>:$e' 74( 4e+3 a$X iDO :+ 4':$ ?7(b 7|4Y'4 X T a)`? 7?4o' 7|4Y'
7|4Y'[ 4o hO 7{e: a+ q4: O 5 | E4[ KJ
R1 =

Ct
2

:E44 7$&48 7{e:  a+ 7%4: 7|4Y'[ 8' aK' v&| E&48 5 4':$
R2 =

345 7849 7:;/ -' + :+

Ct
2

 !" #$%& '() : 

Radar Simulation Using Matlab

 
  

:E44 X  a+ hO 7{e:  a+ 7%4& 7|4Y) [8' v&


R2 =

C (T + t )
2

:4{e:

4j] [8' i! 4' 4

:E? 4') [ 4Y( 4 i o


.['X .1
. 4UXJ E 7j;' 7+4j;3 .2
.7{e: 7+4j; .3
.4{e: K  .4
:7$4! D' h 4 i 4(
.7{e: + .1
.4{e: K a) .2
.h'ZX 7+4j;3 .3
:v;4$] ( D [ D [ )( )74  4 i a'b 7)QY' 3(' E4 % a$e

:U| 4 4)

MATLAB Function pulse_train.m


function [dt, prf, pav, ep, ru] = pulse_train(tau, pri, p_peak)
c = 3.0e+8;
dt = tau / pri;
prf = 1. / pri;

345 7849 7:;/ -' + :+

 !" #$%& '() : 

Radar Simulation Using Matlab

 
  

pav = p_peak * dt;


ep = p_peak * tau;
ru = 1.0e-3 * c * pri / 2.0;
return

:7|4Y' 7]
''>" .4>U:$ `>$$' > 7>4o'> O IJ> h>?+ ] >I5 E R [)4X 4U:+ eX
E]:>>" 4>%) M h> Y>o 4>'U:$ 7|4Y>' Rmax h> Rmin a>) >4>%) a'>b [>'X >
:(7)
M=

Rmax Rmin
R

> 4>' 4>%' a'>b [>)4[ KJ> [>?(]>[ ;> O h>?+ R 7|4Y>' 4U{>X a>+ [5g: E O
:[KJ E| b)

a'>b ['XY>'( '; 7 )e+ 4%) !D a) [?(| E]:5 4%' mg& E| o E O 4)
.4T 7%4X) 4$:o
. a$Y( [ a) R a) [$?o h?+ '?X v$''5)  E)QY) |4+ [KJ
R =

345 7849 7:;/ -' + :+

Ct C
=
2 2B

 !" #$%& '() : 

Radar Simulation Using Matlab

 
  

:7|4Y' 7][ 8' i! 4' 4


:4 i 4U)QY E 3(' E4 % a$e

:U| 4 4)
function [delta_R] = range_resolution(bandwidth,indicator)
% This function computes radar range resolution in meters
% the bandwidth must be in Hz ==> indicator = Hz.
% Bandwidth may be equal to (1/pulse width)==> indicator = seconds
c = 3.e+8;
if(indicator == 'hz')
delta_R = c / (2.0 * bandwidth);
else
delta_R = c * bandwidth / 2.0;
end
return

:?
.748 O 7 '(O a$ `$$'? si U 7+; 4$o ?   QY
i>> K a>K' .>U (> 7>%$& 4>T h>?+ [>"4( > '>` 4 &`I5 ? 4
mKX:$>; >U 7> E>| >j> $c |> E4>4 U 7 k4% #Y s 4e4; 4e )4&`3

345 7849 7:;/ -' + :+

 !" #$%& '() : 

Radar Simulation Using Matlab

  
 

^ h?+ 'U4| 7 i o ) a ; >^ h?+ 9 '> 4>') 7> >>^ 5>] h
' 7| mKX4 U ' a+ Xe ;  h4 ^ ' 7' 4 ):[KJ4 a$e

hjX ? :7]!X4


C v
f
C+v r

= f 'r

C +v
f
C v t

= f 't

C v
2v
= f = f '0 f
= f f
f
d
0 C +v 0 0 C v 0

|4 E  a$ 7 Q"[  a$ U 4X ; 7+U ' E| 4:[KJ

  !" #$%& '() :

+'+ : 7849 7:;/ -345

Radar Simulation Using Matlab

 
  

:[KJ iD|  ?


2v
f =
f COS
d 0

:?[  8' i 4


:E4 % #Y 3(' K  ? 4Y( 4 i o$;

:4
function [fd, tdr] = doppler_freq(freq, ang, tv, indicator)
% This function computes Doppler frequency and time dilation factor ratio
% tau_prime / tau
format long
c = 3.0e+8;
ang_rad = ang * pi /180.;
lambda = c / freq;
if (indicator == 1)
fd = 2.0 * tv * cos(ang_rad) / lambda;
tdr = (c - tv) / (c + tv);
else
345 7849 7:;/ -' + :+

10

 !" #$%& '() : 

Radar Simulation Using Matlab

 
  

fd = -2.0 * c * tv * cos(and_rad) / lambda;


tdr = (c + tv) / (c -tv);
end
return

: 74X)
d>$ (Omni directional antenna ) 44%3 $' E| 74Y) 7]4^ J EW 4: aK$
:7$4 7]!X4 g a) 7jo& E| 7]4j 7|48 ( aK'  4X &' v 4$WU a) : i

D =

Peak transmitted power watt


area of a sphere
m2
:7]!X4 4jX) R 7|4Y' h?+ 7]4j 7|48 K E44

D =

P
t
4R 2
. 4U?; E h'ZX 7+4j;3 E Pt d$

4>$WU ki>>` $'> g a>) `>$ a'>b 7>]4j 7>|48 4> `7U )4$W  X ['XY
:E a$?)4X ai a$ 7]!X Ae EWU? 74Xg 7(g4 G )[ 4X
G2
Ae =
4
:7]!X4 hjX 4$WU a) : i'[ XY i  4UXJ E 7]4j 7|48K| E44

D =

PG
t
4R 2

7$(j>; 4>$ h> > 4>') >4>%) a'>b >] U jY j  ;4UXJ E 7]4j
7>]4j ki> 7>$'' X 44%3 $' E| 7$Y$^4:c)U 7]4^ 4X 4$ ki o U j; h?+
>X> $ )>[ 4X [>)X ki> >[ a>+ >eX d>$ 4U:) 7+:5' 4' EW4`$g v?K ,U % h?+
:7]!X4 hjX v `)( Radar Cross Section RCS) jo'
P
= r m2

P
D

345 7849 7:;/ -' + :+

11

 !" #$%& '() : 

Radar Simulation Using Matlab

 
  

: 4U?eoY E 7$?K 7]4j K E44 .U a+ 7YKX:' 7]4j Pr d$

P G
Ae
P = t
Dr
(4R 2 ) 2
:7$4 7]!X h?+ [5(& Ae 7'$] X
P G 22
P = t
Dr
(4 )3 R 4
>U 4i> #eY> i> q$%{> Xe> 7>?'() >[ >e] a>) 7?eoY' 4T K E?'X ] E|
7>|48 4 vg" aK' EWJ+ K q$%{ i . 4U?'XY E  $' E| 7u) $u
7> 7>)`(> X N J> 7+4j>;> e .(Power Spectral Density PSD) 7>]4j I>$^
: B?
N = NoisePSD * B
:7]!X4 hjX EWU[ D J 7+4j;
N i = kTe B

k = 1.38 *10 23 joul / kelvin

:F J 5 X)[ 4X J a+ eX


F=

SNR
S /N
i = i i
SNR
S /N
o
o o

X4>.[eoY>'> D[ D h?+ s( Signal Noise Ratio) q$%{? 4T 7eY&  SNR d$
:7o4Y 74X' E|
Si = kTe BF ( SNR)o

:7]!X4 4U:+ eX 4U4eo; aK' 4 c" | E44


S

min

= kTe BF ( SNR)

o min

: 4UgJK E h'ZX 7|4Y' %&

345 7849 7:;/ -' + :+

12

 !" #$%& '() : 

Radar Simulation Using Matlab

 
  

P G 22

t
Rmax =
3
(4 ) KT B( SNR)
e
o min

1
4

:|4K E
P G 22
t
(SNR) =
0
3
(4 ) KT BFR 4
e
:7]!X e5 L 4U `) E  4+4$b 4e+3 a$X iDO :+
P G 22
t
(SNR) =
0
3
(4 ) KT BFLR 4
e

: 74X) a+ eX i 4


:E4 % #Y 4Q'! D' K SNR 4Y( 4 i o$;

345 7849 7:;/ -' + :+

13

 !" #$%& '() : 

Radar Simulation Using Matlab

 
  

:U| 4 4)

MATLAB Function radar_eq.m


function [out_par] = radar_eq(pt, freq, g, sigma, te, b, nf, loss, input_par, option,
rcs_delta1, rcs_delta2, pt_percent1, pt_percent2)
c = 3.0e+8;
lambda = c / freq;
p_peak = base10_to_dB(pt);
lambda_sq = lambda^2;
lambda_sqdb = base10_to_dB(lambda_sq);
sigmadb = base10_to_dB(sigma);
for_pi_cub = base10_to_dB((4.0 * pi)^3);
k_db = base10_to_dB(1.38e-23);
te_db = base10_to_dB(te)
b_db = base10_to_dB(b);
if (option == 1)
temp = p_peak + 2. * g + lambda_sqdb + sigmadb - ...
for_pi_cub - k_db - te_db - b_db - nf - loss - input_par;
out_par = dB_to_base10(temp)^(1/4)
% calculate sigma(+-)10dB (rcs +- rcs_delta1,2)
sigmap = rcs_delta1 + sigmadb;
sigmam = sigmadb - rcs_delta2.;
% calculate.pt_percent1 * pt and pt_percent2% * pt
pt05 = p_peak + base10_to_dB(pt_percent1);
pt200 = p_peak + base10_to_dB(pt_percent2);
index = 0;
% vary snr from.5 to 1.5 of default value
for snrvar = input_par*.5: 1: input_par*1.5
index = index + 1;
range1(index) = dB_to_base10(p_peak + 2. * g + lambda_sqdb + ...
sigmam - for_pi_cub - k_db - te_db - b_db - nf - loss - snrvar) ...
^(1/4) / 1000.0;
range2(index) = dB_to_base10(p_peak + 2. * g + lambda_sqdb + ....

345 7849 7:;/ -' + :+

14

 !" #$%& '() : 

Radar Simulation Using Matlab

 
  

sigmadb - for_pi_cub - k_db - te_db - b_db - nf - loss - snrvar) ...


^(1/4) / 1000.0;
range3(index) = dB_to_base10(p_peak + 2. * g + lambda_sqdb + ...
sigmap - for_pi_cub - k_db - te_db - b_db - nf - loss - snrvar) ...
^(1/4) / 1000.0;
end
index = 0;
for snrvar = input_par*.5: 1: input_par*1.5;
index = index + 1;
rangp1(index) = dB_to_base10(pt05 + 2. * g + lambda_sqdb + ...
sigmadb - for_pi_cub - k_db - te_db - b_db - nf - loss - snrvar) ...
^(1/4) / 1000.0;
rangp2(index) = dB_to_base10(p_peak + 2. * g + lambda_sqdb + ...
sigmadb - for_pi_cub - k_db - te_db - b_db - nf - loss - snrvar) ...
^(1/4) / 1000.0;
rangp3(index) = dB_to_base10(pt200 + 2. * g + lambda_sqdb + ...
sigmadb - for_pi_cub - k_db - te_db - b_db - nf - loss - snrvar) ...
^(1/4) / 1000.0;
end
snrvar = input_par*.5: 1: input_par*1.5;
figure (1)
subplot (2,1,1)
plot (snrvar,range2,snrvar,range1,snrvar,range3)
legend ('default RCS','RCS-rcs_delta1','RCS+rcs_delta2')
xlabel ('Minimum SNR required for detection - dB');
ylabel ('Detection range - Km');
%title ('Plots correspond to input parameters from example 1.4');
subplot (2,1,2)
plot (snrvar,rangp2,snrvar,rangp1,snrvar,rangp3)
legend ('default power','.pt_percent1*pt', 'pt_percent2*pt')
xlabel ('Minimum SNR required for detection - dB');
ylabel ('Detection range - Km')
else
range_db = base10_to_dB(input_par * 1000.0);
out_par = p_peak + 2. * g + lambda_sqdb + sigmadb - ...
for_pi_cub - k_db - te_db - b_db - nf - loss - 4.0 * range_db
% calculate sigma -- rcs_delta1,2 dB
sigma5 = sigmadb - rcs_delta1;
sigma10 = sigmadb - rcs_delta2;
% calculate pt_percent1% * pt and pt_percent2*pt
pt05 = p_peak + base10_to_dB(pt_percent1);
pt200 = p_peak + base10_to_dB(pt_percent2);
index = 0;
% vary snr from .5 to 1.5 of default value

345 7849 7:;/ -' + :+

15

 !" #$%& '() : 

Radar Simulation Using Matlab

 
  

for rangvar = input_par*.5 : 1 : input_par*1.5


index = index + 1;
var = 4.0 * base10_to_dB(rangvar * 1000.0);
snr1(index) = p_peak + 2. * g + lambda_sqdb + sigmadb - ...
for_pi_cub - k_db - te_db - b_db - nf - loss - var;
snr2(index) = p_peak + 2. * g + lambda_sqdb + sigma5 - ...
for_pi_cub - k_db - te_db - b_db - nf - loss - var;
snr3(index) = p_peak + 2. * g + lambda_sqdb + sigma10 - ...
for_pi_cub - k_db - te_db - b_db - nf - loss - var;
end
index = 0;
for rangvar = input_par*.5 : 1 : input_par*1.5;
index = index + 1;
var = 4.0 * base10_to_dB(rangvar * 1000.0);
snrp1(index) = pt05 + 2. * g + lambda_sqdb + sigmadb - ...
for_pi_cub - k_db - te_db - b_db - nf - loss - var;
snrp2(index) = p_peak + 2. * g + lambda_sqdb + sigmadb - ...
for_pi_cub - k_db - te_db - b_db - nf - loss - var;
snrp3(index) = pt200 + 2. * g + lambda_sqdb + sigmadb - ...
for_pi_cub - k_db - te_db - b_db - nf - loss - var;
end
end
rangvar = input_par*.5 : 1 : input_par*1.5;
2000 by Chapman & Hall/CRC
figure (2)
subplot (2,1,1)
plot (rangvar,snr1,rangvar,snr2,rangvar,snr3)
legend ('default RCS','RCS-rcs_delta1','RCS-rcs_delta2')
xlabel ('Detection range - Km');
ylabel ('SNR - dB');
%title ('Plots correspond to input parameters from example 1.4');
subplot (2,1,2)
plot (rangvar,snrp2,rangvar,snrp1,rangvar,snrp3)
legend ('default power','.pt_percent1*pt','pt_percent2*pt')
xlabel ('Detection range - Km');
ylabel ('SNR - dB');
Input file radar_reqi.m
% Use this input file to reproduce
clear all
pt = 1.5e+6; % peak power in Watts
freq = 5.6e+9; % radar operating frequency in Hz
g = 45.0; % antenna gain in dB
sigma = 0.1; % radar cross section in m square
te = 290.0; % effective noise temperature in Kelvins

345 7849 7:;/ -' + :+

16

 !" #$%& '() : 

Radar Simulation Using Matlab

 
  

b = 5.0e+6; % radar operating bandwidth in Hz


nf = 3.0; % noise figure in dB
loss = 0.0; % radar losses in dB
option = 1; % 1 ===> input_par = SNR in dB
% 2 ===> input_par = Range in Km
input_par = 20;
rcs_delta1 = 5.0; % rcs variation choice 1
rcs_delta2 =10.0; % rcs variation choice2
pt_percent1 = 0.5; % peak power variation choice 1
pt_percent2 =2.0; % peak power variation choice 2
:Q 4K K

345 7849 7:;/ -' + :+

17

 !" #$%& '() : 

Radar Simulation Using Matlab

 
  

:(Low PRF Radar) gQ:' 4{e: K   74X)


Pt E>>> h>>>'ZX 7+4j>>>;3 T 4{>>>e:>>> K E>>> >>> 4UXJ>>> E>>> 7{>>>e:>>> + aK$>>>
:4eo;3'[ + 4Y aK' dt= /T d$ Pav=Pt*dt E 7j;' 7+4j;3
dr =

T
= 1 * f r
T

[>'X |> E4>4( T>> ) ( |>Low PRF) gQ:>' 4{>e:> K > > 74 E|
:e5; SNR 74X' %& dt1
P G 22 * n
t
p
(SNR) np =
(4 )3 R 4kT BFL
e
where
n p = Ti f r
P G 22 *T f
i r
(SNR) np = t
(4 )3 R 4kT BFL
e
>K > fr 4{>e:? >U> X a>) > Ti > U h>?+ 7j]4Y> 4{>e:> + > np d$
.4{e:

345 7849 7:;/ -' + :+

18

 !" #$%& '() : 

Radar Simulation Using Matlab

 
  

:(Low PRF) 74X) a+ eX i 4


4>Q'! D>'> K gQ:>' 4{>e:> K  > 74X) ;[ $8' 4 i o
:E4 % #Y

: 4

MATLAB Function lprf_req.m


function [snr_out] = lprf_req (pt, freq, g, sigma, te, b, nf, loss, range, prf, np,
rcs_delta,
pt_percent, np1, np2)
% This program implements the LOW PRF radar equation.
c = 3.0e+8;
lambda = c / freq;
p_peak = base10_to_dB(pt);
lambda_sq = lambda^2;
345 7849 7:;/ -' + :+

19

 !" #$%& '() : 

Radar Simulation Using Matlab

 
  

lambda_sqdb = base10_to_dB(lambda_sq);
sigmadb = base10_to_dB(sigma);
for_pi_cub = base10_to_dB((4.0 * pi)^3);
k_db = base10_to_dB(1.38e-23);
te_db = base10_to_dB(te)
b_db = base10_to_dB(b);
np_db = base10_to_dB(np);
range_db = base10_to_dB(range * 1000.0);
% Implement Eq. (1.65)
snr_out = p_peak + 2. * g + lambda_sqdb + sigmadb + np_db - ...
for_pi_cub - k_db - te_db - b_db - nf - loss - 4.0 * range_db
% Generate plots
index = 0;
n1 = np_db;
n2 = base10_to_dB(np1);
n3 = base10_to_dB(np2)
for range_var = 25:5:400 % 25 - 400 Km
index = index + 1;
rangevar_db = base10_to_dB(range_var * 1000.0);
snr1(index) = p_peak + 2. * g + lambda_sqdb + sigmadb + n1 - ...
for_pi_cub - k_db - te_db - b_db - nf - loss - 4.0 * rangevar_db;
snr2(index) = p_peak + 2. * g + lambda_sqdb + sigmadb + n2 - ...
for_pi_cub - k_db - te_db - b_db - nf - loss - 4.0 * rangevar_db;
snr3(index) = p_peak + 2. * g + lambda_sqdb + sigmadb + n3 - ...
for_pi_cub - k_db - te_db - b_db - nf - loss - 4.0 * rangevar_db;
end
figure(1)
var = 25:5:400;
plot(var,snr1,'k',var,snr2,'k--',var,snr3,'k--.')
legend('np = 1','np1','np2')
xlabel ('Range - Km');
ylabel ('SNR - dB');
%title ('np = 1, np1 = 10, np2 =100');
% Generate plots
sigma5 = sigmadb - rcs_delta.;
pt05 = p_peak + base10_to_dB(pt_percent);
index = 0;
for nvar =1:10:500 % 500 pulses
index = index + 1;
ndb = base10_to_dB(nvar);
snrs(index) = p_peak + 2. * g + lambda_sqdb + sigmadb + ndb - ...
for_pi_cub - k_db - te_db - b_db - nf - loss - 4.0 * range_db;
snrs5(index) = p_peak + 2. * g + lambda_sqdb + sigma5 + ndb - ...
for_pi_cub - k_db - te_db - b_db - nf - loss - 4.0 * range_db;

345 7849 7:;/ -' + :+

20

 !" #$%& '() : 

Radar Simulation Using Matlab

 
  

end
index = 0;
for nvar =1:10:500 % 500 pulses
index = index + 1;
ndb = base10_to_dB(nvar);
snrp(index) = p_peak + 2. * g + lambda_sqdb + sigmadb + ndb - ...
for_pi_cub - k_db - te_db - b_db - nf - loss - 4.0 * range_db;
snrp5(index) = pt05 + 2. * g + lambda_sqdb + sigmadb + ndb - ...
for_pi_cub - k_db - te_db - b_db - nf - loss - 4.0 * range_db;
end
nvar =1:10:500;
figure (2)
subplot (2,1,1)
plot (nvar,snrs,'k',nvar,snrs5,'k --')
legend ('default RCS','RCS-delta')
xlabel ('Number of coherently integrated pulses');
ylabel ('SNR - dB');
%title ('delta = 10, percent = 2');
subplot (2,1,2)
plot (nvar,snrp,'k',nvar,snrp5,'k --')
legend ('default power','pt * percent')
xlabel ('Number of coherently integrated pulses');
ylabel ('SNR - dB');

Input file lprf_reqi.m


pt = 1.5e+6; % peak power in Watts
freq = 5.6e+9; % radar operating frequency in Hz
g = 45.0; % antenna gain in dB
sigma = 0.1; % radar cross section in m square
te = 290.0; % effective noise temperature in Kelvins
b = 5.0e+6; % radar operating bandwidth in Hz
nf = 3.0; % noise figure in dB
loss = 0.0; % radar losses in dB
np = 1; % 1 number of coherently integrated pulses
prf = 100 ; % PRF in Hz
range = 250.0; % target range in Km
np1 = 10; % choice 1 of np
np2 = 100; % choice 2 of np
rcs_delta = 10.0; % rcs variation
pt_percent = 2.0; % pt variation

345 7849 7:;/ -' + :+

21

 !" #$%& '() : 

Radar Simulation Using Matlab

 
  

:4KO K

345 7849 7:;/ -' + :+

22

 !" #$%& '() : 

Radar Simulation Using Matlab

 
  

References:
1. Radar Analyses and Design Using Matlab, Bassem R. Mahafza, Ph.D.
2. WWW.RADARTUTORIAL>.COM

345 7849 7:;/ -' + :+

23

 !" #$%& '() : 

You might also like