main content

modeling radar detectability factors -凯发k8网页登录

this example shows how to model antenna, transmitter, and receiver gains and losses for a detailed radar range equation analysis. we start by computing the available signal-to-noise ratio (snr) at the radar receiver using the snr form of the radar equation. next, we define the detectability factor as a threshold snr required to make a detection with specified probabilities of detection, pd, and false alarm, pfa. we then estimate the maximum range of the system as the range at which the available snr is equal to the detectability factor, that is the maximum range at which the detection with the specified pd and pfa is still possible. the example further explores the impact of losses introduced by different components of the radar system on the estimated maximum range. we first consider the effects of the sensitivity time control (stc) and eclipsing on the available snr. scanning and signal processing losses that require an increase in the radar detectability factor are considered next. the example concludes by computing the resulting pd at the detector output to show the impact of the losses on the detection performance of the radar system.

available snr

the radar equation combines the main parameters of a radar system and allows a radar engineer to compute a maximum detection range, a required peak transmit power, or a maximum available snr for a radar system. the radar equation is typically a family of relatively simple formulas each corresponding to one of these three key performance characteristics. a common form of the radar equation for computing the maximum available snr at a range r is:

snr=ptτgtgrλ2σ(4π)3ktsr4l

where

pt is the peak transmit power

τ is the transmitted pulse width

gt is the transmit antenna gain

gr is the receive antenna gain

λ is the radar wavelength

σ is the radar target cross section (rcs)

k is boltzmann's constant

ts is the system noise temperature

l is the general loss factor that combines losses along the transmitter-target-receiver path that reduce the received signal energy.

on the right-hand side all parameters except the target range and rcs are under control of the radar designer. this equation states that for a target of a given size located at a range r, the snr available at the receiver can be increased by transmitting more power, increasing the antenna size, using lower frequency, or having a more sensitive receiver.

consider an s-band airport surveillance radar operating at the frequency of 3 ghz. the peak transmit power is 0.2 mw, the transmit and receive antenna gain is 34 db, the pulse duration is 11 μs, and the noise figure is 4.1 db. assume the radar is required to detect a target with 1 m2 rcs at the maximum range rm of 100 km.

lambda = freq2wavelen(3e9);         % wavelength (m)
pt = 0.2e6;                         % peak power (w)
tau = 1.1e-5;                       % pulse width (s)
g = 34;                             % transmit and receive antenna gain (db)
ts = systemp(4.1);                  % system temperature (k) 
rcs = 1;                            % target radar cross section (m^2)
rm = 100e3;                         % required maximum range (m)

to start, assume no losses, that is l = 0 db. we use the radar equation to compute the available snr at the receiver as a function of target range.

l = 0;                              % combined transmission line and propagation losses (db)
r = (1:40:130e3).';                 % range samples (m)
snr = radareqsnr(lambda,r,pt,tau,'gain',g,'ts',ts,'rcs',rcs,'loss',l);

calculate the available snr at the required maximum range of 100 km.

snratrm = snr(find(r>=rm,1))
snratrm = 18.3169

plot the maximum range requirement together with the computed available snr.

radarmetricplot(r*1e-3,snr,'metricname','available snr','maxrangerequirement',rm*1e-3,...
    'rangeunit','km','metricunit','db','radarname',{'surveillance radar'});
legend('location','best');

figure contains an axes object. the axes object with title available snr vs range, xlabel target range (km), ylabel available snr (db) contains 2 objects of type line, constantline. these objects represent surveillance radar, max range.

required snr

is the computed available snr high enough to make a detection? since the signal processed by a radar receiver is a combination of the transmitted waveform and the random noise, the answer to this question depends on the desired probability of detection pd and the maximum acceptable probability of false alarm pfa. these probabilities define the required snr, also known as the detectability factor (or detectability). the detectability factor is the minimum snr required to declare a detection with the specified probabilities of detection and false alarm. it also depends on the rcs fluctuation and the type of the detector. compute the detectability factor for a single pulse received from a steady (swerling 0) target by a square-law detector assuming pd = 0.9 and pfa = 1e-6.

pd = 0.9;
pfa = 1e-6;
d0 = detectability(pd,pfa,1,'swerling0')
d0 = 13.1217

calculate the detectability factor for a swerling 1 fluctuating target, which is a more accurate model for real-world targets. for the swerling 1 target the single-pulse detectability factor is significantly higher.

d1 = detectability(pd,pfa,1,'swerling1')
d1 = 21.1436

the resultant required snr is higher than the available snr, which means that a swerling 1 target will not be detected with a single pulse. a common way to lower the detectability factor is to perform pulse integration. calculate the detectability factor for n = 10 noncoherently integrated pulses.

n = 10;
dn = detectability(pd,pfa,n,'swerling1')
dn = 13.5033

this is lower than the available snr. thus, after a non-coherent integration of 10 pulses the radar system will be able to detect a 1 m2 target at the required maximum range of 100 km with the probability of detection 0.9 and a false alarm of 1e-6.

the detectability factor computed for a swerling 1 target and n pulses combines the effects of the integration gain and the fluctuation loss. the integration gain is a difference between the snr required to detect a steady target using a single pulse and the snr required to detect a steady target using n pulses.

gi = detectability(pd,pfa,1,'swerling0') - detectability(pd,pfa,n,'swerling0')
gi = 7.7881

the fluctuation loss is a difference between the snr required to detect a fluctuating target and the snr required to detect a steady target.

lf = detectability(pd,pfa,n,'swerling1') - detectability(pd,pfa,n,'swerling0')
lf = 8.1696

use the radarbudgetplot function to illustrate the components of the detectability factor.

radarbudgetplot([d0 -gi lf], {'single-pulse steady target','pulse integration gain','fluctuation loss'});
title('detectability factor')

figure contains an axes object. the axes object with title detectability factor, ylabel required snr (db) contains 10 objects of type line, patch, text. these objects represent detectability factor, losses, gains.

substitute the detectability factor into the range form of the radar equation as the minimum required snr to evaluate the actual maximum range of the system.

radareqrng(lambda,dn,pt,tau,'gain',g,'ts',ts,'rcs',rcs,'loss',l,'unitstr','km')
ans = 131.9308

to clearly indicate at which ranges the detection with the desired pd and the maximum acceptable pfa is possible, we add the computed detectability factor as a horizontal line to the snr vs range plot. we also use a stoplight chart to color code ranges and snr levels according to the computed detectability. at the ranges where the available snr curve passes through the green zone the radar meets the detection requirement, while at the range where it is in the red zone the detection with the specified pd and pfa is not possible.

radarmetricplot(r*1e-3,snr,dn, ...
    'metricname','available snr', ...
    'requirementname','detectability', ...
    'maxrangerequirement',rm*1e-3, ...
    'rangeunit','km','metricunit','db', ...
    'showstoplight',true, ...
    'radarname',{'surveillance radar'});
title([{'available snr vs range'}, {'(no losses)'}]);
legend('location','best');

figure contains an axes object. the axes object with title available snr vs range (no losses), xlabel target range (km), ylabel available snr (db) contains 5 objects of type patch, line, constantline. these objects represent pass, fail, surveillance radar, detectability, max range.

all ranges beyond the required maximum range are colored green and marked as pass.

this analysis assumes zero losses and therefore cannot adequately predict the range of the actual radar system. a real radar system with the specified parameters will have a shorter maximum range due to:

  • propagation effects caused by the earth's surface and atmosphere. these effects reduce the amount of available signal energy at the receiver.

  • various losses experienced throughout the radar system. some losses in this category reduce the available snr, while other losses result in increased detectability factor.

the following sections consider in more detail the impact of the losses belonging to the second category on the range performance of the radar system.

range-dependent factors

when designing a surveillance radar system, several factors must be included in the radar equation to account for the decrease in the available signal energy at the receiver.

eclipsing

pulse radar systems turn off their receivers during the pulse transmission. thus, the target echoes arriving from the ranges within one pulse length from the radar or within one pulse length around the unambiguous range will be eclipsed by the transmitted pulse resulting in only a fraction of the pulse being received and processed. the radar system considered in this example has a pulse width of 11μs. the closest range from which a full pulse can be received is the minimum range rmin.

rmin = time2range(tau)
rmin = 1.6489e 03

the echoes from the targets that are closer than 1649 m will arrive before the pulse transmission has been completed. similar effect can be observed for the targets located at or near the multiples of the unambiguous range. assuming the pulse repetition frequency is 1350 hz (pulse repetition interval t0.75 ms), calculate the unambiguous range of the system.

prf = 1350;                         % pulse repetition frequency
rua = time2range(1/prf)
rua = 1.1103e 05

the echoes arriving from the ranges rua±rminwill be eclipsed by the next transmitted pulse. this diagram illustrates pulse eclipsing. the arrowheads indicate the front of the pulse.

eclipsing.png

due to eclipsing, the available snr will have deep notches at 0 range and the ranges equal to the multiples of rua. add the eclipsing factor to the radar equation in order to account for the loss in the available snr due to pulse eclipsing.

du = tau*prf;                       % duty cycle
fecl = eclipsingfactor(r,du,prf);   % eclipsing factor
snr = radareqsnr(lambda,r,pt,tau,'gain',g,'ts',ts,'rcs',rcs,'customfactor',fecl,'loss', l);
radarmetricplot(r*1e-3,snr,dn, ...
    'metricname','available snr', ...
    'requirementname','detectability', ...
    'maxrangerequirement',rm*1e-3, ...
    'rangeunit','km','metricunit','db', ...
    'showstoplight',true, ...
    'radarname',{'surveillance radar'});
title([{'available snr vs range'}, {'(with eclipsing)'}]);
legend('location','best');

figure contains an axes object. the axes object with title available snr vs range (with eclipsing), xlabel target range (km), ylabel available snr (db) contains 5 objects of type patch, line, constantline. these objects represent pass, fail, surveillance radar, detectability, max range.

the real-world radar systems utilize prf diversity in order to prevent the eclipsing loss and to extend the unambiguous range of the system.

sensitivity time control (stc)

a typical surveillance radar system must transmit a considerable amount of power to detect targets at long ranges. although the available energy rapidly decays with range, at very close ranges even small targets can have very strong returns due to high peak transmit power. such strong returns from small nuisance targets (birds, insects) might result in undesirable detections, while the regular size targets or nearby clutter can saturate the receiver. it is highly desirable for the surveillance radar system to avoid these kind of nuisance detections. to solve this problem the radar systems use stc. it scales the receiver gain up to a cutoff range rstc to maintain a constant signal strength as a target approaches the radar.

rstc = 60e3;                        % stc cutoff range (m)
xstc = 4;                           % stc exponent selected to maintain target detectability at ranges below rstc (since the signal power is inversely proportional to r^4)
fstc = stcfactor(r,rstc,xstc);      % stc factor
snr = radareqsnr(lambda,r,pt,tau,'gain',g,'ts',ts,'rcs',rcs,'customfactor',fecl fstc,'loss',l);
radarmetricplot(r*1e-3,snr,dn, ...
    'metricname','available snr', ...
    'requirementname','detectability', ...
    'maxrangerequirement',rm*1e-3, ...
    'rangeunit','km','metricunit','db', ...
    'showstoplight',true, ...
    'radarname',{'surveillance radar'});
title([{'available snr vs range'}, {'(with stc and eclipsing for 1 m^2 target)'}]);
legend('location','best');
ylim([-30 30])

figure contains an axes object. the axes object with title available snr vs range ( w i t h blank s t c blank a n d blank e c l i p s i n g blank f o r blank 1 blank m squared baseline blank t a r g e t ), xlabel target range (km), ylabel available snr (db) contains 5 objects of type patch, line, constantline. these objects represent pass, fail, surveillance radar, detectability, max range.

after adding the stc factor, the plot shows that the 1 m2 rcs target is still detected everywhere up to the maximum range rm, while a small target with rcs of 0.03 m2 will not be able to reach the required pd of 0.9 at any range and thus will be rejected.

snrsmallrcs = radareqsnr(lambda,r,pt,tau,'gain',g,'ts',ts,'rcs',0.03,'customfactor',fecl fstc,'loss',l);
radarmetricplot(r*1e-3,snrsmallrcs,dn, ...
    'metricname','available snr', ...
    'requirementname','detectability', ...
    'maxrangerequirement',rm*1e-3, ...
    'rangeunit','km','metricunit','db', ...
    'showstoplight',true, ...
    'radarname',{'surveillance radar'});
title([{'available snr vs range'}, {'(with stc and eclipsing for 0.03 m^2 target)'}]);
legend('location','best');
ylim([-30 20])

figure contains an axes object. the axes object with title available snr vs range ( w i t h blank s t c blank a n d blank e c l i p s i n g blank f o r blank 0 . 0 3 blank m squared baseline blank t a r g e t ), xlabel target range (km), ylabel available snr (db) contains 5 objects of type patch, line, constantline. these objects represent pass, fail, surveillance radar, detectability, max range.

it is clear from these plots that the stc only scales the available snr up to the specified cutoff range and does not affect the available snr at the maximum range of interest.

scanning

a radar system can scan a search volume either by mechanically rotating the antenna or by using the phased array antenna and performing electronic scanning. an imperfect shape of the antenna beam and the process of sweeping the beam across a search volume introduces additional losses to the system.

beam shape loss

the radar equation uses a peak value of the antenna gain assuming each received pulse has the maximum amplitude. in reality, as the beam passes the target, the received pulses are modulated by the two-way pattern of the scanning antenna resulting in beam shape loss. computing the exact value of this loss would require knowing the exact antenna pattern. this information might not be available in the early stages of the radar system design when this type of analysis is usually performed. instead, the shape of the main lobe of a typical practical antenna can be well approximated by a gaussian shape. assuming that the radar system performs dense sampling in the spatial domain (the beam moves by less than 0.71 of the half power beamwidth), calculate the beam shape loss for one dimensional scanning.

lb = beamloss
lb = 1.2338

the beam shape loss is doubled if the radar system scans in both azimuth and elevation.

beamloss(true)
ans = 2.4677

scan sector loss

in this example we assume the radar system employs an electronically steered phased array to perform scanning. using the phased array antenna will cause an increase in the required snr due to two effects: 1) beam broadening due to the reduced projected array area in the beam direction, and 2) reduction of the effective aperture area of the individual array elements at off-broadside angles. to account for these effects, add the scan sector loss to the detectability factor. assume that the system in the example scans only in the azimuth dimension and the scan sector spans from –60 to 60 degrees. compute the resultant loss.

theta = [-60 60];
larray = arrayscanloss(pd,pfa,n,theta,'swerling1')
larray = 2.7745

signal processing

prior to detection the received radar echoes must pass through the radar signal processing chain. the purpose of different components in the signal processing chain is to guarantee the required probabilities of detection and false alarm, reject unwanted echoes from clutter, and account for variable or non-gaussian noise. we further consider several components of the signal processing loss that must be accounted for in a surveillance radar system.

mti

moving target indicator (mti) is a process of rejecting fixed or slowly moving clutter while passing echoes from targets moving with significant velocities. typical mti uses 2, 3, or 4-pulse canceller that implements a high-pass filter to reject echoes with low doppler shifts. passing the received signal through the mti pulse canceller introduces correlation between the noise samples. this in turn reduces the total number of independent noise samples available for integration resulting in mti noise correlation loss. additionally, the mti canceller significantly suppresses targets with velocities close to the nulls of its frequency response causing an additional mti velocity response loss. assuming a 2-pulse canceller is used, calculate these two components of the mti loss.

m = 2;
[lmti_a, lmti_b] = mtiloss(pd,pfa,n,m,'swerling1')
lmti_a = 1.4468
lmti_b = 8.1562

in a system that uses a single prf, the mti velocity response loss can be very high for the high required probability of detection. to eliminate this loss, prf diversity is almost always used in real radar systems.

binary integration

binary integration is a suboptimal noncoherent integration technique also known as the m-of-n integration. if m out of n received pulses exceed a predetermined threshold, a target is declared to be present. the binary integrator is a relatively simple automatic detector and is less sensitive to the effects of a single large interference pulse that might exist along with the target echoes. therefore, the binary integrator is more robust when the background noise or clutter is non-gaussian. since the binary integration is a suboptimal technique, it results in a binary integration loss compared to optimal noncoherent integration. the optimal value of m is not a sensitive selection and it can be quite different from the optimum without significant penalty resulting in the binary integration loss being lower than 1.4 db. calculate the binary integration loss when n is 10 and m is set to 6.

m = 6;
lbint = binaryintloss(pd,pfa,n,m)
lbint = 1.0549

the binaryintloss function computes the loss assuming a steady (swerling 0) target. since the fluctuation loss is included in the detectability factor, the same binary integration loss computation can be used in the case of a fluctuating target.

cfar

constant false alarm rate (cfar) detector is used to maintain an approximately constant rate of false target detections when the noise or the interference levels vary. since cfar averages a finite number of reference cells to estimate the noise level, the estimates are subject to an error which leads to a cfar loss. cfar loss is an increase in the snr required to achieve a desired detection performance using cfar when the noise levels are unknown compared to a fixed threshold with a known noise level. calculate the cfar loss assuming that total 120 cells are used for cell-averaging cfar.

nrc = 120;
lcfar = cfarloss(pfa,nrc)
lcfar = 0.2500

effective detectability factor

the scanning and the signal processing losses increase the detectability factor, which means that more energy is required to make a detection. the resulting detectability factor that includes effects of all these losses is called the effective detectability factor. use the radarbudgetplot function to create a waterfall chart showing the combined effect of the computed scanning and signal processing losses on the detectability factor.

d = [d0 -gi lf lmti_a lmti_b lbint lcfar larray lb];
radarbudgetplot(d, {'single-pulse steady target','pulse integration gain','fluctuation loss'...
    'mti loss', 'binary integration loss', 'cfar loss', 'scan sector loss', 'beam shape loss'});
title('detectability factor')

figure contains an axes object. the axes object with title detectability factor, ylabel required snr (db) contains 25 objects of type line, patch, text. these objects represent detectability factor, losses, gains.

the resulting effective detectability factor equals 28.42 db. by taking the scanning and the signal processing losses into account, the required snr increases by almost 15 db. the analysis shows that the system cannot actually meet the stated requirement of detecting a 1 m2 rcs target at 100 km with pd = 0.9 and pfa = 1e-6.

radarmetricplot(r*1e-3,snr,sum(d), ...
    'metricname','available snr', ...
    'requirementname','detectability', ...
    'maxrangerequirement',rm*1e-3, ...
    'rangeunit','km','metricunit','db', ...
    'showstoplight',true, ...
    'radarname',{'surveillance radar'});
title([{'available snr vs range'}, {'(with stc, eclipsing, scanning and signal processing losses)'}])
legend('location','best')
ylim([-10 30]);

figure contains an axes object. the axes object with title available snr vs range (with stc, eclipsing, scanning and signal processing losses), xlabel target range (km), ylabel available snr (db) contains 5 objects of type patch, line, constantline. these objects represent pass, fail, surveillance radar, detectability, max range.

this problem can be addressed either by increasing the available snr or decreasing the required snr. transmitting more power or increasing the antenna gains brings up the available snr, while increasing the integration time lowers the required snr. however, in some applications, a subset of the system parameters can be constrained by other requirements and thus cannot be changed. for example, if the analysis is performed for an existing system, increasing the available snr might not be an option. in that case, making adjustments to the signal processing chain to lower the detectability factor could be an acceptable solution. to lower the required snr, in the following sections we assume that the number of pulses n is increased from 10 to 40.

in addition, we can change the requirements on the maximum range and the detection probabilities. instead of a single number specifying the desired probability of detection or the maximum range, a pair of objective and threshold values can be defined. the objective requirement describes the desired performance level of the system that would be needed to fully satisfy the mission needs. the threshold requirement describes a minimum acceptable performance level of the system. using a pair of values to define a requirement instead of a single value provides more flexibility to the design and creates a trade-space for selecting the system parameters. in this example, we assume that the objective requirement for pd is 0.9 and set the threshold value to 0.8. similarly, the objective maximum range requirement remains 100 km, while the threshold value is set to 90 km. the detectability factor is now computed both for the objective and the threshold pd.

n = 40;
m = 18;
pd = [0.9 0.8];
[lmti_a, lmti_b] = mtiloss(pd,pfa,n,m,'swerling1');
dx = detectability(pd,pfa,n,'swerling1')   cfarloss(pfa,nrc)   beamloss ...
      lmti_a   lmti_b   binaryintloss(pd,pfa,n,m)   arrayscanloss(pd,pfa,n,theta,'swerling1')
dx = 2×1
   24.2522
   18.0494
rm = [100e3 90e3];
radarmetricplot(r*1e-3,snr,dx(1),dx(2), ...
    'metricname','available snr', ...
    'requirementname','detectability', ...
    'maxrangerequirement',rm*1e-3, ...
    'rangeunit','km','metricunit','db', ...
    'showstoplight',true, ...
    'radarname',{'surveillance radar'});
title([{'available snr vs range'}, {'(n=40)'}])
legend('location','best')
ylim([-10 30]);

figure contains an axes object. the axes object with title available snr vs range (n=40), xlabel target range (km), ylabel available snr (db) contains 8 objects of type patch, line, constantline. these objects represent pass, warn, fail, surveillance radar, threshold detectability, objective detectability, objective max range, threshold max range.

the snr vs range plot now has a yellow warn zone indicating the snr values and target ranges where the performance of the system is between the objective and the threshold requirements. we can see that up to approximately 70 km the system meets the objective requirement for pd. from 70 km to 100 km the objective requirement for pd is violated while the threshold requirement is still satisfied.

effective probability of detection

the snr vs range plot above shows that the detection performance of the radar system varies with range. a 1 m2 target at ranges below 70 km will be detected with the probability of detection greater or equal to 0.9, while between 70 km and 100 km it will be detected with pd of at least 0.8. since some of the considered losses depend on the probability of detection, the actual pd at the detector output varies with range. we can use the roc curve to compute pd as a function of range.

% generate a vector of probability values at which to compute the roc curve
p = probgrid(0.1,0.9999,100);
% compute the required snr at these probabilities
[lmti_a, lmti_b] = mtiloss(p,pfa,n,m,'swerling1');
dx = detectability(p,pfa,n,'swerling1')   cfarloss(pfa,nrc)   beamloss ...
      lmti_a   lmti_b   binaryintloss(p,pfa,n,m)   arrayscanloss(p,pfa,n,theta,'swerling1');
% plot the roc curve
helperradarpdvssnrplot(dx,p,[0.1 0.9999]);

figure contains an axes object. the axes object with xlabel effective detectability factor (db), ylabel probability of detection contains an object of type line.

the effective probability of detection at the detector output can now be computed by interpolating this roc curve at the available snr values.

% interpolate the roc curve at the available snr
pdeff = rocinterp(dx,p,snr,'snr-pd');
% plot the effective pd as a function of range
radarmetricplot(r*1e-3,pdeff,pd(1),pd(2), ...
    'metricname','effective p_d', ...
    'requirementname','p_d', ...
    'maxrangerequirement',rm*1e-3, ...
    'rangeunit','km', ...
    'showstoplight',true, ...
    'radarname',{'surveillance radar'});
legend('location','best')
ylim([0.5 1.0])

figure contains an axes object. the axes object with title effective p_d vs r a n g e effective p_d vs range, xlabel target range (km), ylabel effective p indexof d baseline effective p_d contains 8 objects of type patch, line, constantline. these objects represent pass, warn, fail, surveillance radar, threshold p_d, objective p_d, objective max range, threshold max range.

this result shows that due to application of the stc the probability of detection is almost constant for ranges from 2 km to 60 km. for a target with 1 m2 rcs it is higher than 0.92. in the range between 70 km and 87 km the effective pd is above 0.85. at the threshold value of the maximum range requirement the probability of detection is approximately 0.84, and at the objective range of 100 km it is slightly above 0.8.

summary

this example shows how various losses impact the detection performance of a radar system. it starts with a radar equation and introduces the concepts of the available snr and the detectability factor. considering an example surveillance radar system, it shows how the available snr is reduced by stc and eclipsing, while the detectability factor is increased by scanning and signal processing losses. finally, the example demonstrates how to compute the effective probability of detection at the receiver output for different target ranges.

references

  1. barton, d. k. radar equations for modern radar. artech house, 2013.

  2. richards, m. a., scheer, j. a. & holm, w. a. principle of modern radar: basic principles. scitech publishing, 2010.

网站地图