# Modelling and Analysis of Constant False Alarm Rate Performance in Presence of Jamming Environments

Authors: Xiaozhou Chen; Mengzhong Hu; Shu Lu
Journal: Mathematical Problems in Engineering
 (2022)

---

## Abstract

A novel “Bernoulli experiment model” of constant false alarm rate (CFAR) is presented in order to analyze CFAR performance in jamming condition. In this model, target detection can be treated as a sequence of independent subevents, and the detection process equals to pick up target energy from an energy pool diluted with noise, interference, and jamming pulses. In a multiple-jammer case, each jammer contributes independently, and the content of the energy pool only relates to the corresponding jammer. Impacts of jammer factors, such as signal-to-interference ratio (SIR), jamming-to-interference ratio (JIR), the number of false targets in reference cells, and the number of false targets in a cell under test (CUT), are analyzed and compared in detail, and two jamming operation principles are deduced as an application. The deduction procedure and conclusions drawn from the model can also be utilized on similar occasions.


---

## Body

## 1. Introduction

The detection of an unknown signal needs to compare the signal plus interference power in the cell under test (CUT) to a predetermined threshold. The case that interference exceeds the threshold with no target echo presented in CUT is defined as a false alarm. In Neyman–Pearson criterion, we need to control the false alarm rate to an acceptable range while maximizing the detection probability [1]. From that perspective, estimating the interference level and dynamically adjusting the detection threshold accordingly is required. This gives birth to constant false alarm rate (CFAR) techniques.The core idea of CFAR is to utilize some range cells lagging or leading the CUT as reference cells, design an algorithm to estimate interference from the reference cells, and determine the threshold according to the estimation. In this way, thresholds will track the changes of interference automatically, and the false alarm rate will be maintained to a constant value. The traditional CFAR scheme is cell-average (CA-) CFAR, where the estimation is based on the mean value of all reference cells. CA-CFAR performs well in homogeneous environments but suffers from self-masking and mutual target masking in heterogeneous environments. Improvements are designed by selecting or discarding a certain number of reference cells to overcome the above drawbacks, yielding greatest of (GO) [2–4], smallest of (SO) [5], trimmed mean(TM) [6, 7], and order statistic (OS) CFARs [8, 9]. All of the above schemes concentrate on the cell selecting strategies, but the basic models and fundamental computations are like CA-CFAR.Nowadays, sophisticated and intelligent jammers bring about unprecedented challenges for CFAR. Traditionally, jammers are able to create large amounts of false target “clouds,” covering most areas on a radar screen, where real targets are entirely submerged. If the false targets are dense enough, they act like continuous wave noise, effectively raising the dynamic threshold and reducing the probability of detection. Moreover, advanced techniques and strategies such as cooperative jamming [10, 11], distributive jamming [12, 13], and smart noise [14, 15] make it more complicating to evaluate effects of every single technique. A clear and precise mathematical model revealing the relationship between jamming parameters and CFAR characteristics is urgently required.Unfortunately, compared to the deep research on CFAR techniques, CFAR performance in complex electromagnetic environment, especially in jamming conditions, is not frequently discussed. Some researchers created rough models to analyze CFAR in heterogeneous environments, by which various CFAR schemes were compared [4, 16, 17]. However, the presented model was not widely applicable and did not unveil its physical meaning. Yan-Juan et al. [18] created a mathematical model and calculated the impact of jammer-to-signal ratio and distance of false targets, but the model is not sufficiently simplified, and complicated expressions concealed the direct relationship between each parameter. Bachmann et al. [19] focused on dynamic interaction between radars and jammers, and the optimization method under certain constrains was illustrated. Liu and Li [20] discussed the power cooperative and time cooperative jamming case, and effects from different jamming modulations were analyzed in detail. Neither Bachmann et al. [19] nor Liu and Li [20] presented a generalized model. Kunpeng et al. [21] deduced the CFAR detection capability under a type of coherent jamming (namely, intermittent sampling and repeated forwarding) and compared the effects of corresponding parameters. But the conclusion is unavailable in complicating random electronic environments since the deduction was not based on statistical analysis.To overcome the above shortcomings, this paper mainly makes the following contributions:(1)
A “Bernoulli experiment model” (BEM) is constructed to compute the probability of detection of CFAR. Unlike in reference [18–20], this model provides an analytical, closed-form, and adequately simplified formation to describe the performance in jamming conditions and can be calculated directly since the parameters are determined. Similar ideas can be absorbed to utilize in the detection of slight changes in a sensor network [22].(2)
BEM is more generalized and widely applicable, and can be transformed into several special cases, where models in reference [4, 16, 17] are only one of them. Moreover, it has a clear physical meaning, where the detection process can be seen as a sequence of Bernoulli experiments to pick the desired energy from the energy pool.(3)
Influences of the number of jamming cells, jamming-to-interference ratio (JIR) and signal-to-interference ratio (SIR), are quantitatively analyzed and compared in detail with the help of the BEM. As an application, we obtain the “optimum power distribution principle” and “maximum jamming cells principle” in Section4.4 to better understand and predict jammers’ operation before finding an efficient counter measurement. Since the model is based on statistical analysis, these influences and principles are more applicable than those in reference [21].(4)
An analytic deduction is presented by means of rigorous direct integral, where the calculation techniques used in the derivation process will provide valuable inspiration for solving similar problems. Simplification using the moment generating function (MGF) simplifies and generalizes the deduction, which is also shown in Section3.


### 1.1. Manuscript Notation

Italic letters denote scalar parameters, with capital letterZ and lower-case letter z denoting the random variable and constant variable, respectively. P⋅ or f⋅ denotes the probability density function (PDF) of a random variable, while P⋅|⋅ or f⋅|⋅ is the corresponding conditional counterparts. ∗ denotes the convolution operator, ! denotes factorial, and finally, ⋅M denotes rising factorials in the Pochhammer symbol.


## 1.1. Manuscript Notation

Italic letters denote scalar parameters, with capital letterZ and lower-case letter z denoting the random variable and constant variable, respectively. P⋅ or f⋅ denotes the probability density function (PDF) of a random variable, while P⋅|⋅ or f⋅|⋅ is the corresponding conditional counterparts. ∗ denotes the convolution operator, ! denotes factorial, and finally, ⋅M denotes rising factorials in the Pochhammer symbol.


## 2. CFAR Principles

Figure1 depicts the brief flow of CFAR. After an inclusion of a series of signal processing but not limited to A/D converter, doppler processing, pulse compressing, and array processing, digital samples are passed to a square law detector to convert complex values into their magnitude-squared values. During this noncoherent integration process, the phase information of complex data is discarded, leading to a loss in integration gain. A certain number of data from leading and lagging cells are stored and selected according to predetermined strategies before they are added up. At last, the threshold is obtained from the product of the summation and coefficient α, which represents our expected false alarm rate level.Figure 1
Basic principle of CFAR.The output of square law detector is a random variable. UnderH0 hypothesis (where only interference, yet no target return exists in the CUT), the PDF of output z can be expressed as(1)PZz=1σi2exp−zσi2,where σi2 is the average interference power. Interference originates from unwanted clutter, scatter, diffraction, and thermal noise (we do not exclude thermal noise from interference, since they are treated as a whole in the following analysis). Here, we assume that the interference is homogeneous; i.e., the interference is independent identically distributed in both leading, lagging reference cells, and CUT, which is restrictive especially in mountain or urban backgrounds. The nonhomogeneous scenario can be discussed according to the main idea from reference [23]. For distributed multiradar systems, the detection performance in nonstationary Weibull distributed clutter is discussed from reference [24].Assuming thatN reference cells are selected according to specific strategy and each has the same probability density as equation (1), the joint PDF of interference power is the product of each PDF, thereby(2)PZz=1σi2Nexp−∑n=1Nznσi2,where z=z1,…,zN. Then, the maximum likelihood estimate of average interference power in all the reference cells is(3)σ^i2=1N∑n=1Nzn.In CFAR, the thresholdUT is defined as the product of σ^i2 and the constant coefficient α; thus,(4)UT=ασ^i2,where α is determined by desired false-alarm probability PFA:(5)α=NPFA−1/N−1.Given the Swerling 1 or 2 target PDFf1x|H1 under H1 hypothesis (where CUT contains interference and target return) is still exponentially distributed similar to equation (1), and σi2 is replaced with total power of interference and target σi2+σt2; therefore,(6)fCUTxH1=1σi2+σt2exp−xσi2+σt2,x≥0,0,x<0.Without losing generosity, we induceγCUT as α times inverse of total power of interference σi2 and target σt2, namely,(7)γCUTα=σi2+σt2−1=σi21+η−1.Here,η≜σt2/σi2 is defined as SIR.The probability of detection can be calculated as the integral fromUT to infinity.(8)PD=∫UT∞fCUTx|H1dx.Substituting equation (6) into (8), we obtain that(9)PDγCUT|σ^i2=∫ασ^i2∞γCUTαexp−γCUTαxdx=exp−γCUT⋅σ^i2.Given the dictatedPFA, PD can be derived once σ^i2 is determined. Considering σ^i2 is a random variable, the expected value of PD is(10)PD¯γCUT=∫0∞PDγCUT|σ^i2fσ^i2σ^i2dσ^i2,where fσ^i2σ^i2 denotes PDF of total interference power in N reference cells.


## 3. Bernoulli Experiment Model

### 3.1. CFAR on False Targets

False targets will have great impacts onPD and PFA, and drastically degrade CFAR performance. On the one hand, false targets in reference cells will raise the estimated interference level as well as dynamic threshold, yielding the declination of PD; on the other hand, jamming pulses in CUT will cause false alarm and increase PFA.Assuming there is only one jammer source, radiatingM false points in reference cell (therefore M is defined as the number of jamming cells) and Q in CUT, respectively. Since the interference power in CUT has thus been lift up to Q times, γCUT (denoting as γ) in equation (7) can be rewritten as(11)γ≜αQσi2+σt2−1=ασi2Q+η−1.Accordingly,PD is transformed from equation (9) as(12)PDγ|x=exp−γx,lettingx to be the total power in reference cells.The expectation of detection probabilityPD¯ can be derived once γCUT and fσ^i2σ^i2 in equation (10) are known. Since the former is expressed in equation (11), the key step is the deduction of fσ^i2σ^i2. Here, we directly provide the final result, and corresponding derivation is listed in Appendix if readers are interested.(13)PD¯=1+αNQ+η−N1+1θαNQ+η−M≜P0NP1M,where 1/θ=JIR denotes the jammer-to-interference ratio.As shown in equation (13), PD¯ is composed of two factors, the former P0N=1+α/NQ+η−N is mainly determined by interference, and the latter P1M=1+1/θα/NQ+η−M reflects the jamming effects. The number of jamming cells M and JIR (denoted by 1/θ) are the main factors that influence CFAR performance. Moreover, M is more influential than JIR because of its exponential effects. Comparatively, Q impacts both factors, since jamming pulses in CUT become the “targets” to be detected.The corresponding false alarm ratePFA¯ can be derived by letting η=0, namely,(14)PFA¯=1+αNQ−N1+1θαNQ−M.The above equation illustrates that, under the jamming condition,PFA¯ is no more a constant value. Apparently, PFA¯ will decrease when JIR or M increases, inasmuch as rising power in reference cells also raises the threshold. Nevertheless, Q has an opposite effect, since the growing number of false targets in CUT are more likely to be detected mistakenly.


### 3.2. CFAR on Noise Jamming

Noise jamming scenario can be treated as a special case of false target, where the number of jamming cellsM is the same as the total number of reference cells N. In fact, as a generalized model, equation (13) can be simplified to some expressions to adapt to concrete special occasion; for example, the steps are as follows:(i)
No jamming case. When settingM=0 and Q=1, we can get the theoretical PD of CFAR, namely,(15)PD0¯=1+αN1+η−N.(ii)
Multitarget case. The average of jamming power equals to that of interference; thusθ=1/η, and Q=1; therefore,(16)PD−MT¯=1+αN1+η−N1+αηN1+η−M.(iii)
Continuous wave jamming (or noise jamming) case. In this case, jamming signal “fills up” all the reference cells, soM=N and Q=1, namely,(17)PD−CWJ¯=PD0¯1+1θαN1+η−N.As for jamming waveform, a continuous wave requires much less JIR than discrete pulses. In Section4.3, we will see that even the slightest JIR would produce considerable effects on CFAR.


### 3.3. Physical Explanation

The framework of equation (13) can be interpreted from two levels: level I is from a subevent perspective (i.e., the explanation of P0 or P1), and level II is how the subevents are combined. The former is discussed in Section 3.3.1, where each subevent reflects the probability of picking the “effective target energy” from an “energy pool.” The latter is analyzed in Section 3.3.2, where the total probability is the product of this subprobability, which means each subevent is mutually independent, and the results can also be extended to multiple jammer scenarios. Deepening the research on the physical meaning may provide us with an intuitive understanding of CFAR scheme.


#### 3.3.1. Energy Pool Model of Subevents

Noting that factorsP0 and P1 in equation (13) can be transformed as(18)P0=1+PFA0−1/N−1Q+σt02/σi02−1=σt02+Qσi02σt02+Qσi02+PFA0−1/N−1σi02,(19)P1=1+σj02σi02⋅PFA0−1/N−1Q+σt02/σi02−1=σt02+Qσi02σt02+Qσi02+PFA0−1/N−1σj02,where σt02, σi02, σj02 refers to the average power in each cell of target return, interference, and jamming, respectively.Equation (18) can be interpreted as the following “energy pool model.” Considering there are two colors of balls uniformly merged in a pool. Red balls represent effective target power, and white balls represent effective interference power. The number of two colors of balls is σt02+Qσi02 and PFA0−1/N−1σi02, respectively. To detect targets is equivalent to picking up “red balls” from the pool, or more precisely, to select “effective target energy” from the “energy pool,” and the probability is apparently described as equation (18) (Figure 2(a)). Equation (19) has a similar interpretation, with interference energy replaced by jamming energy (see Figure 2(b)). From this perspective, the detection process can be seen as a sequence of subevents concatenating one another (as depicted in Figure 3), and the number of subevents equals the number of cells containing interference/jamming signals (namely, N or M). Only if all the subevents happen, targets can be detected.Figure 2
Energy pool model of subevents. (a) Target detection from interference. (b) Target detection from jamming.
(a)(b)Figure 3
Sequence of sub-event explanation of (13).


#### 3.3.2. Bernoulli Experiment Model

Equation (13) implicitly indicates that, the event that signal power in CUT exceeding detection threshold can be partitioned into 2 classes of independent events: N subevents with the probability mainly related to interference impact and M sub-events mainly affected by jamming source. From statistical perspective, this can be named as “Bernoulli experiment model,” since interference and jamming source act independently. Inspired by the intuition above, the Bernoulli experiment model can be an extent to a multijammer case.Theorem 1.
If there areS jamming sources, each radiating Mn false targets in reference cell and Qn in CUT, with JIR equal to 1/θnn=1,2,…, then the overall detection probability of CFAR is(20)PD=∏n=0S1+1θnαNQ+η−Mn,where n=0 denotes the interference without jamming, with the definition θ0≜1, M0≜N, and(21)Q=∑i=1SQi.Proof.
As for thenth jamming source, the overall power of Mn false targets has Erlang distribution similar to (A.1), by which the moment generating function (MGF) is(22)Fnt=∫−∞∞fnxe−txdt=∫0∞λnMnxMn−1Mn−1!e−λn+txdt=1+tλn−Mn,
where the parameterλn is defined as the inverse of average jamming power in Mn cells.(23)λn=Mnσn2.
As shown above, the total PDF of jamming plus interference power is the convolution of each PDF, where MGF (denoted asFt) is the product of them.(24)fx=f0x∗f1x…∗fSx,(25)Ft=∏n=0SFnt.
In the multijammer case, equations (11) and (12) still hold when Q is redefined as equation (21). By substituting equation (12) into (10), we attain that,(26)PD¯=∫0∞exp−γxfxdx=Ft|t=γ.
Since(27)Fnγ=1+γλn−Mn=1+σn2Mnασi2Q+η−Mn=1+σn2/Mnσi2/NαNQ+η−Mn=1+1θnαNQ+η−Mn.
The final solution of equation (20) can be attained by substituting equation (27) into (25). QED.
Equation (20) is not a coincidence, and there are two conditions leading to the multiplication form: (i) the prior probability of detection PDγ|x can be expressed in the form of exponential function such as equation (12); (ii) the jamming/interference signals out of square law detector are exponentially distributed, and the summation has Erlang distribution. Both conditions can easily be met in reality; therefore, the independence test model is widely applicable.
In the multijammer case, the framework in Figure3 can be extended as Figure 4. Equation (20) demonstrates that effects of different jamming sources work independently and can be analyzed, respectively. It should be noted that the number of jammer source is mathematical rather than the physical concept. Coherent jammers where pulses are from different transmitters are radiated at the same time, with energy synthesized via spatial power combining technique, can be treated as one jammer with doubled JIR; while multifunctional jammers which have different working mode (or function) are supposed to be divided into several unique monofunctional jammer sources.Figure 4
Framework of Bernoulli experiment model.


## 3.1. CFAR on False Targets

False targets will have great impacts onPD and PFA, and drastically degrade CFAR performance. On the one hand, false targets in reference cells will raise the estimated interference level as well as dynamic threshold, yielding the declination of PD; on the other hand, jamming pulses in CUT will cause false alarm and increase PFA.Assuming there is only one jammer source, radiatingM false points in reference cell (therefore M is defined as the number of jamming cells) and Q in CUT, respectively. Since the interference power in CUT has thus been lift up to Q times, γCUT (denoting as γ) in equation (7) can be rewritten as(11)γ≜αQσi2+σt2−1=ασi2Q+η−1.Accordingly,PD is transformed from equation (9) as(12)PDγ|x=exp−γx,lettingx to be the total power in reference cells.The expectation of detection probabilityPD¯ can be derived once γCUT and fσ^i2σ^i2 in equation (10) are known. Since the former is expressed in equation (11), the key step is the deduction of fσ^i2σ^i2. Here, we directly provide the final result, and corresponding derivation is listed in Appendix if readers are interested.(13)PD¯=1+αNQ+η−N1+1θαNQ+η−M≜P0NP1M,where 1/θ=JIR denotes the jammer-to-interference ratio.As shown in equation (13), PD¯ is composed of two factors, the former P0N=1+α/NQ+η−N is mainly determined by interference, and the latter P1M=1+1/θα/NQ+η−M reflects the jamming effects. The number of jamming cells M and JIR (denoted by 1/θ) are the main factors that influence CFAR performance. Moreover, M is more influential than JIR because of its exponential effects. Comparatively, Q impacts both factors, since jamming pulses in CUT become the “targets” to be detected.The corresponding false alarm ratePFA¯ can be derived by letting η=0, namely,(14)PFA¯=1+αNQ−N1+1θαNQ−M.The above equation illustrates that, under the jamming condition,PFA¯ is no more a constant value. Apparently, PFA¯ will decrease when JIR or M increases, inasmuch as rising power in reference cells also raises the threshold. Nevertheless, Q has an opposite effect, since the growing number of false targets in CUT are more likely to be detected mistakenly.


## 3.2. CFAR on Noise Jamming

Noise jamming scenario can be treated as a special case of false target, where the number of jamming cellsM is the same as the total number of reference cells N. In fact, as a generalized model, equation (13) can be simplified to some expressions to adapt to concrete special occasion; for example, the steps are as follows:(i)
No jamming case. When settingM=0 and Q=1, we can get the theoretical PD of CFAR, namely,(15)PD0¯=1+αN1+η−N.(ii)
Multitarget case. The average of jamming power equals to that of interference; thusθ=1/η, and Q=1; therefore,(16)PD−MT¯=1+αN1+η−N1+αηN1+η−M.(iii)
Continuous wave jamming (or noise jamming) case. In this case, jamming signal “fills up” all the reference cells, soM=N and Q=1, namely,(17)PD−CWJ¯=PD0¯1+1θαN1+η−N.As for jamming waveform, a continuous wave requires much less JIR than discrete pulses. In Section4.3, we will see that even the slightest JIR would produce considerable effects on CFAR.


## 3.3. Physical Explanation

The framework of equation (13) can be interpreted from two levels: level I is from a subevent perspective (i.e., the explanation of P0 or P1), and level II is how the subevents are combined. The former is discussed in Section 3.3.1, where each subevent reflects the probability of picking the “effective target energy” from an “energy pool.” The latter is analyzed in Section 3.3.2, where the total probability is the product of this subprobability, which means each subevent is mutually independent, and the results can also be extended to multiple jammer scenarios. Deepening the research on the physical meaning may provide us with an intuitive understanding of CFAR scheme.


### 3.3.1. Energy Pool Model of Subevents

Noting that factorsP0 and P1 in equation (13) can be transformed as(18)P0=1+PFA0−1/N−1Q+σt02/σi02−1=σt02+Qσi02σt02+Qσi02+PFA0−1/N−1σi02,(19)P1=1+σj02σi02⋅PFA0−1/N−1Q+σt02/σi02−1=σt02+Qσi02σt02+Qσi02+PFA0−1/N−1σj02,where σt02, σi02, σj02 refers to the average power in each cell of target return, interference, and jamming, respectively.Equation (18) can be interpreted as the following “energy pool model.” Considering there are two colors of balls uniformly merged in a pool. Red balls represent effective target power, and white balls represent effective interference power. The number of two colors of balls is σt02+Qσi02 and PFA0−1/N−1σi02, respectively. To detect targets is equivalent to picking up “red balls” from the pool, or more precisely, to select “effective target energy” from the “energy pool,” and the probability is apparently described as equation (18) (Figure 2(a)). Equation (19) has a similar interpretation, with interference energy replaced by jamming energy (see Figure 2(b)). From this perspective, the detection process can be seen as a sequence of subevents concatenating one another (as depicted in Figure 3), and the number of subevents equals the number of cells containing interference/jamming signals (namely, N or M). Only if all the subevents happen, targets can be detected.Figure 2
Energy pool model of subevents. (a) Target detection from interference. (b) Target detection from jamming.
(a)(b)Figure 3
Sequence of sub-event explanation of (13).


### 3.3.2. Bernoulli Experiment Model

Equation (13) implicitly indicates that, the event that signal power in CUT exceeding detection threshold can be partitioned into 2 classes of independent events: N subevents with the probability mainly related to interference impact and M sub-events mainly affected by jamming source. From statistical perspective, this can be named as “Bernoulli experiment model,” since interference and jamming source act independently. Inspired by the intuition above, the Bernoulli experiment model can be an extent to a multijammer case.Theorem 1.
If there areS jamming sources, each radiating Mn false targets in reference cell and Qn in CUT, with JIR equal to 1/θnn=1,2,…, then the overall detection probability of CFAR is(20)PD=∏n=0S1+1θnαNQ+η−Mn,where n=0 denotes the interference without jamming, with the definition θ0≜1, M0≜N, and(21)Q=∑i=1SQi.Proof.
As for thenth jamming source, the overall power of Mn false targets has Erlang distribution similar to (A.1), by which the moment generating function (MGF) is(22)Fnt=∫−∞∞fnxe−txdt=∫0∞λnMnxMn−1Mn−1!e−λn+txdt=1+tλn−Mn,
where the parameterλn is defined as the inverse of average jamming power in Mn cells.(23)λn=Mnσn2.
As shown above, the total PDF of jamming plus interference power is the convolution of each PDF, where MGF (denoted asFt) is the product of them.(24)fx=f0x∗f1x…∗fSx,(25)Ft=∏n=0SFnt.
In the multijammer case, equations (11) and (12) still hold when Q is redefined as equation (21). By substituting equation (12) into (10), we attain that,(26)PD¯=∫0∞exp−γxfxdx=Ft|t=γ.
Since(27)Fnγ=1+γλn−Mn=1+σn2Mnασi2Q+η−Mn=1+σn2/Mnσi2/NαNQ+η−Mn=1+1θnαNQ+η−Mn.
The final solution of equation (20) can be attained by substituting equation (27) into (25). QED.
Equation (20) is not a coincidence, and there are two conditions leading to the multiplication form: (i) the prior probability of detection PDγ|x can be expressed in the form of exponential function such as equation (12); (ii) the jamming/interference signals out of square law detector are exponentially distributed, and the summation has Erlang distribution. Both conditions can easily be met in reality; therefore, the independence test model is widely applicable.
In the multijammer case, the framework in Figure3 can be extended as Figure 4. Equation (20) demonstrates that effects of different jamming sources work independently and can be analyzed, respectively. It should be noted that the number of jammer source is mathematical rather than the physical concept. Coherent jammers where pulses are from different transmitters are radiated at the same time, with energy synthesized via spatial power combining technique, can be treated as one jammer with doubled JIR; while multifunctional jammers which have different working mode (or function) are supposed to be divided into several unique monofunctional jammer sources.Figure 4
Framework of Bernoulli experiment model.


## 3.3.1. Energy Pool Model of Subevents

Noting that factorsP0 and P1 in equation (13) can be transformed as(18)P0=1+PFA0−1/N−1Q+σt02/σi02−1=σt02+Qσi02σt02+Qσi02+PFA0−1/N−1σi02,(19)P1=1+σj02σi02⋅PFA0−1/N−1Q+σt02/σi02−1=σt02+Qσi02σt02+Qσi02+PFA0−1/N−1σj02,where σt02, σi02, σj02 refers to the average power in each cell of target return, interference, and jamming, respectively.Equation (18) can be interpreted as the following “energy pool model.” Considering there are two colors of balls uniformly merged in a pool. Red balls represent effective target power, and white balls represent effective interference power. The number of two colors of balls is σt02+Qσi02 and PFA0−1/N−1σi02, respectively. To detect targets is equivalent to picking up “red balls” from the pool, or more precisely, to select “effective target energy” from the “energy pool,” and the probability is apparently described as equation (18) (Figure 2(a)). Equation (19) has a similar interpretation, with interference energy replaced by jamming energy (see Figure 2(b)). From this perspective, the detection process can be seen as a sequence of subevents concatenating one another (as depicted in Figure 3), and the number of subevents equals the number of cells containing interference/jamming signals (namely, N or M). Only if all the subevents happen, targets can be detected.Figure 2
Energy pool model of subevents. (a) Target detection from interference. (b) Target detection from jamming.
(a)(b)Figure 3
Sequence of sub-event explanation of (13).


## 3.3.2. Bernoulli Experiment Model

Equation (13) implicitly indicates that, the event that signal power in CUT exceeding detection threshold can be partitioned into 2 classes of independent events: N subevents with the probability mainly related to interference impact and M sub-events mainly affected by jamming source. From statistical perspective, this can be named as “Bernoulli experiment model,” since interference and jamming source act independently. Inspired by the intuition above, the Bernoulli experiment model can be an extent to a multijammer case.Theorem 1.
If there areS jamming sources, each radiating Mn false targets in reference cell and Qn in CUT, with JIR equal to 1/θnn=1,2,…, then the overall detection probability of CFAR is(20)PD=∏n=0S1+1θnαNQ+η−Mn,where n=0 denotes the interference without jamming, with the definition θ0≜1, M0≜N, and(21)Q=∑i=1SQi.Proof.
As for thenth jamming source, the overall power of Mn false targets has Erlang distribution similar to (A.1), by which the moment generating function (MGF) is(22)Fnt=∫−∞∞fnxe−txdt=∫0∞λnMnxMn−1Mn−1!e−λn+txdt=1+tλn−Mn,
where the parameterλn is defined as the inverse of average jamming power in Mn cells.(23)λn=Mnσn2.
As shown above, the total PDF of jamming plus interference power is the convolution of each PDF, where MGF (denoted asFt) is the product of them.(24)fx=f0x∗f1x…∗fSx,(25)Ft=∏n=0SFnt.
In the multijammer case, equations (11) and (12) still hold when Q is redefined as equation (21). By substituting equation (12) into (10), we attain that,(26)PD¯=∫0∞exp−γxfxdx=Ft|t=γ.
Since(27)Fnγ=1+γλn−Mn=1+σn2Mnασi2Q+η−Mn=1+σn2/Mnσi2/NαNQ+η−Mn=1+1θnαNQ+η−Mn.
The final solution of equation (20) can be attained by substituting equation (27) into (25). QED.
Equation (20) is not a coincidence, and there are two conditions leading to the multiplication form: (i) the prior probability of detection PDγ|x can be expressed in the form of exponential function such as equation (12); (ii) the jamming/interference signals out of square law detector are exponentially distributed, and the summation has Erlang distribution. Both conditions can easily be met in reality; therefore, the independence test model is widely applicable.
In the multijammer case, the framework in Figure3 can be extended as Figure 4. Equation (20) demonstrates that effects of different jamming sources work independently and can be analyzed, respectively. It should be noted that the number of jammer source is mathematical rather than the physical concept. Coherent jammers where pulses are from different transmitters are radiated at the same time, with energy synthesized via spatial power combining technique, can be treated as one jammer with doubled JIR; while multifunctional jammers which have different working mode (or function) are supposed to be divided into several unique monofunctional jammer sources.Figure 4
Framework of Bernoulli experiment model.


## 4. Numerical Analysis

In this section, several numerical experiments are presented to examine the quantitative impacts of jamming parameters on CFAR performance. There are mainly three factors which include thepower of target echo (ascribing toSIR), jamming power in CUT (ascribing to Q), and jamming power in reference cells (ascribing to M and JIR). PD or PFA is computed through equation (13) or (14) with parameters listed in Table 1, respectively (here, we take SIR and JIR instead of η and 1/θ for convenience).Table 1
Parameter setting.
SectionSIRdBMJIRdBQ4.10∼300/1/5/151514.2200∼10151/3/5/74.3(a)201/5/10/20-5∼2014.3(b)200∼205/10/15/201General settingPFA0=10−6, α=19.91, N=20


### 4.1. Impact of SIR

We depicted the relationship ofPD and SIR (or η in the above equations) in Figure 5. The increase in SIR corresponding to an increase in PD is intuitive, since stronger target return has more possibility to exceed threshold.Figure 5
Impact of SIR onPD.Nevertheless, once jamming power gets into the reference cell, a loss inPD is produced. The more cells are plunged by jamming power, or equivalently more JIRs are used as the worse performance in detection. We need to take measures to raise SIR to “combat,” or to some extent, to “compensate” effects of jamming cells. As shown in the curves of M=0 and M=1, we need to improve SIR to about 4 dB to achieve the corresponding PD as the no-jamming case. And as M increases to 5, the compensation would be as high as 10 dB.


### 4.2. Impact ofQ

Q is the number of those points, originated by jammer, just positioned in CUT, and mistakenly identified as targets. Physically, those points tend to raise the interference level in CUT, so PD and PFA will rise simultaneously as Q increases. The impact is described in Figure 6.Figure 6
Impact of Q. (a)PD vs Q; (b) PFA vs Q.
(a)(b)Although bothPD and PFA are positively influenced by the value of Q, the extents are quite different. Compared to the very limited improvement on PD (PD’s improvement is no more than 0.003 for every increase in Q, and curves in Figure 6(a) almost overlap), it has notable effects on PFA (tens of dozens of dB’s exacerbation for every increase in Q). That is to say, once jamming pulses enter CUT, PFA would be raised to a great extent.


### 4.3. Impact ofM and JIR

The impact ofM and JIR is depicted in Figure 7. Evidently, the rising of both parameters will exacerbate CFAR’s detect capability, and either parameter may determine the PD’ s reduction speed with respect to the other parameter. Moreover, M has relatively greater impact on PD than JIR since lines in Figure 7(b) are more concave than those in Figure 7(a). When M=20 (implying all the reference cells are padded by jamming signals), merely 5 dB JIR may cause PD to reduce from 0.82 to 0.44, almost by half. From this perspective, continuous wave jammer (usually in the form of smart noise or coherent noise) would be especially destructive without much JIR.Figure 7
Impact of jammer onPD. (a) PD vs JIR; (b) PD vs M.
(a)(b)


### 4.4. Jamming Principle

As radars confronting with overheating challenges nowadays, sophisticated technologies enable jammers to produce a complicated waveform and comprehensive electronic attacks. Before finding an efficient way to counter them, we need to shift our eyes from radar to jammers’ operation on the top of which lies in the optimum power distribution problems.


#### 4.4.1. Optimum Jamming Power Distribution Principle

Advanced jammers usually produce hybrid signals combined with a dense false target (mode I) and narrow-band continuous wave noise (mode II) in order to raisePFA and decrease PD simultaneously. A common problem for an electronic warfare engineer is given as total jamming power σΣ2of how to distribute it to different jammer modes in order to minimize PD.We can distinguish respective effects from each mode with the help of equation (20). Assuming power distributed to dense false target mode is rσΣ20<r<1, and to continuous noise mode is 1−rσΣ2. Considering false targets occupying M reference cells, the JIR for each mode is rσΣ2/Mσi2 and 1−rσΣ2/Nσi2, respectively. Then, PD under this complex jammer condition is(28)PD=PD01+rσΣ2Mσi2β−M1+1−rσΣ2Nσi2β−N,where β≜k−1=α/NQ+η. Equation (28) can be written in a logarithmic way as(29)lnPD=lnPD0−Mln1+rσΣ2Mσi2β−Nln1+1−rσΣ2Nσi2β.Letting the derivation oflnPD with respect to r equal to 0, we can get(30)dlnPDdr=−MσΣ2/Mσi2β1+rσΣ2/Mσi2β+NσΣ2/Nσi2β1+1−rσΣ2/Nσi2β=0.namely,(31)r=MM+N.In order to reducePD, the optimal power ratio of noise jamming and false target jamming should be M/N.


#### 4.4.2. Maximum Jamming Cell Principle

Jammers usually have tunable false target density. If this parameter increases,M and Q would both rise, which will lead to a remarkable drop in PD and increase in PFA.On the one hand, from the detection perspective, jamming power will produce stronger effects, when it is evenly distributed to as many reference cells as possible, rather than just concentrating in one or a few reference cells. That can be referred to as “maximum jamming cells principle.”This interesting fact can be observed from Figure7. Noting PDJIR=20dB,M=2=0.21 is much larger than PDJIR=10dB,M=20=0.13, and the discrepancy can be as high as 61.5%. Apparently, the more dispersive the jamming energy is, the more obvious deterioration happens in PD.It is not difficult to understand via the “Bernoulli experiment model.” The distribution of jamming power toM cells would create M subevents and M “compound reduction” in PD, which usually descends faster than the “single test.”On the other hand, the increase inQ will raise the false alarm rate, but more M will counteract this effect (see Figure 6(b)). A compromise should be made between those two factors. In order to fully quantify jamming effects, a comprehensive metric Jamming effect factor (JEF) can be defined as(32)JEF≜PD¯PD0¯+μκPFA0¯PFA¯,where the subscript “0” means parameters without jamming, and κ is a scaling factor to balance the scale of improvements on PD¯ and PFA¯. The higher JEF value, the more influences jamming have on CFAR. μ is a weight coefficient to judge the significance of PD¯ and PFA¯. Generally, μ=1, but in some protective weapons where PD¯ is primarily advocated, μ is supposed to be less than 1, whereas in remote early warning devices where PFA¯ should be rigorously controlled to a restricted range, we can set μ>1. To summarize, here lists steps to find optimum parameters for jamming CFAR.Step 1: To give the initial value of false target densityStep 2: To estimateM and Q according to false target densityStep 3: To computePD¯ and PFA¯ through equations (13) and (14)Step 4: To compute JEF from equation (32) and adjust false target densityStep 5: Change the value of false target density and redo Step 1∼4 until JEF convergents to its maximum value


## 4.1. Impact of SIR

We depicted the relationship ofPD and SIR (or η in the above equations) in Figure 5. The increase in SIR corresponding to an increase in PD is intuitive, since stronger target return has more possibility to exceed threshold.Figure 5
Impact of SIR onPD.Nevertheless, once jamming power gets into the reference cell, a loss inPD is produced. The more cells are plunged by jamming power, or equivalently more JIRs are used as the worse performance in detection. We need to take measures to raise SIR to “combat,” or to some extent, to “compensate” effects of jamming cells. As shown in the curves of M=0 and M=1, we need to improve SIR to about 4 dB to achieve the corresponding PD as the no-jamming case. And as M increases to 5, the compensation would be as high as 10 dB.


## 4.2. Impact ofQ

Q is the number of those points, originated by jammer, just positioned in CUT, and mistakenly identified as targets. Physically, those points tend to raise the interference level in CUT, so PD and PFA will rise simultaneously as Q increases. The impact is described in Figure 6.Figure 6
Impact of Q. (a)PD vs Q; (b) PFA vs Q.
(a)(b)Although bothPD and PFA are positively influenced by the value of Q, the extents are quite different. Compared to the very limited improvement on PD (PD’s improvement is no more than 0.003 for every increase in Q, and curves in Figure 6(a) almost overlap), it has notable effects on PFA (tens of dozens of dB’s exacerbation for every increase in Q). That is to say, once jamming pulses enter CUT, PFA would be raised to a great extent.


## 4.3. Impact ofM and JIR

The impact ofM and JIR is depicted in Figure 7. Evidently, the rising of both parameters will exacerbate CFAR’s detect capability, and either parameter may determine the PD’ s reduction speed with respect to the other parameter. Moreover, M has relatively greater impact on PD than JIR since lines in Figure 7(b) are more concave than those in Figure 7(a). When M=20 (implying all the reference cells are padded by jamming signals), merely 5 dB JIR may cause PD to reduce from 0.82 to 0.44, almost by half. From this perspective, continuous wave jammer (usually in the form of smart noise or coherent noise) would be especially destructive without much JIR.Figure 7
Impact of jammer onPD. (a) PD vs JIR; (b) PD vs M.
(a)(b)


## 4.4. Jamming Principle

As radars confronting with overheating challenges nowadays, sophisticated technologies enable jammers to produce a complicated waveform and comprehensive electronic attacks. Before finding an efficient way to counter them, we need to shift our eyes from radar to jammers’ operation on the top of which lies in the optimum power distribution problems.


### 4.4.1. Optimum Jamming Power Distribution Principle

Advanced jammers usually produce hybrid signals combined with a dense false target (mode I) and narrow-band continuous wave noise (mode II) in order to raisePFA and decrease PD simultaneously. A common problem for an electronic warfare engineer is given as total jamming power σΣ2of how to distribute it to different jammer modes in order to minimize PD.We can distinguish respective effects from each mode with the help of equation (20). Assuming power distributed to dense false target mode is rσΣ20<r<1, and to continuous noise mode is 1−rσΣ2. Considering false targets occupying M reference cells, the JIR for each mode is rσΣ2/Mσi2 and 1−rσΣ2/Nσi2, respectively. Then, PD under this complex jammer condition is(28)PD=PD01+rσΣ2Mσi2β−M1+1−rσΣ2Nσi2β−N,where β≜k−1=α/NQ+η. Equation (28) can be written in a logarithmic way as(29)lnPD=lnPD0−Mln1+rσΣ2Mσi2β−Nln1+1−rσΣ2Nσi2β.Letting the derivation oflnPD with respect to r equal to 0, we can get(30)dlnPDdr=−MσΣ2/Mσi2β1+rσΣ2/Mσi2β+NσΣ2/Nσi2β1+1−rσΣ2/Nσi2β=0.namely,(31)r=MM+N.In order to reducePD, the optimal power ratio of noise jamming and false target jamming should be M/N.


### 4.4.2. Maximum Jamming Cell Principle

Jammers usually have tunable false target density. If this parameter increases,M and Q would both rise, which will lead to a remarkable drop in PD and increase in PFA.On the one hand, from the detection perspective, jamming power will produce stronger effects, when it is evenly distributed to as many reference cells as possible, rather than just concentrating in one or a few reference cells. That can be referred to as “maximum jamming cells principle.”This interesting fact can be observed from Figure7. Noting PDJIR=20dB,M=2=0.21 is much larger than PDJIR=10dB,M=20=0.13, and the discrepancy can be as high as 61.5%. Apparently, the more dispersive the jamming energy is, the more obvious deterioration happens in PD.It is not difficult to understand via the “Bernoulli experiment model.” The distribution of jamming power toM cells would create M subevents and M “compound reduction” in PD, which usually descends faster than the “single test.”On the other hand, the increase inQ will raise the false alarm rate, but more M will counteract this effect (see Figure 6(b)). A compromise should be made between those two factors. In order to fully quantify jamming effects, a comprehensive metric Jamming effect factor (JEF) can be defined as(32)JEF≜PD¯PD0¯+μκPFA0¯PFA¯,where the subscript “0” means parameters without jamming, and κ is a scaling factor to balance the scale of improvements on PD¯ and PFA¯. The higher JEF value, the more influences jamming have on CFAR. μ is a weight coefficient to judge the significance of PD¯ and PFA¯. Generally, μ=1, but in some protective weapons where PD¯ is primarily advocated, μ is supposed to be less than 1, whereas in remote early warning devices where PFA¯ should be rigorously controlled to a restricted range, we can set μ>1. To summarize, here lists steps to find optimum parameters for jamming CFAR.Step 1: To give the initial value of false target densityStep 2: To estimateM and Q according to false target densityStep 3: To computePD¯ and PFA¯ through equations (13) and (14)Step 4: To compute JEF from equation (32) and adjust false target densityStep 5: Change the value of false target density and redo Step 1∼4 until JEF convergents to its maximum value


## 4.4.1. Optimum Jamming Power Distribution Principle

Advanced jammers usually produce hybrid signals combined with a dense false target (mode I) and narrow-band continuous wave noise (mode II) in order to raisePFA and decrease PD simultaneously. A common problem for an electronic warfare engineer is given as total jamming power σΣ2of how to distribute it to different jammer modes in order to minimize PD.We can distinguish respective effects from each mode with the help of equation (20). Assuming power distributed to dense false target mode is rσΣ20<r<1, and to continuous noise mode is 1−rσΣ2. Considering false targets occupying M reference cells, the JIR for each mode is rσΣ2/Mσi2 and 1−rσΣ2/Nσi2, respectively. Then, PD under this complex jammer condition is(28)PD=PD01+rσΣ2Mσi2β−M1+1−rσΣ2Nσi2β−N,where β≜k−1=α/NQ+η. Equation (28) can be written in a logarithmic way as(29)lnPD=lnPD0−Mln1+rσΣ2Mσi2β−Nln1+1−rσΣ2Nσi2β.Letting the derivation oflnPD with respect to r equal to 0, we can get(30)dlnPDdr=−MσΣ2/Mσi2β1+rσΣ2/Mσi2β+NσΣ2/Nσi2β1+1−rσΣ2/Nσi2β=0.namely,(31)r=MM+N.In order to reducePD, the optimal power ratio of noise jamming and false target jamming should be M/N.


## 4.4.2. Maximum Jamming Cell Principle

Jammers usually have tunable false target density. If this parameter increases,M and Q would both rise, which will lead to a remarkable drop in PD and increase in PFA.On the one hand, from the detection perspective, jamming power will produce stronger effects, when it is evenly distributed to as many reference cells as possible, rather than just concentrating in one or a few reference cells. That can be referred to as “maximum jamming cells principle.”This interesting fact can be observed from Figure7. Noting PDJIR=20dB,M=2=0.21 is much larger than PDJIR=10dB,M=20=0.13, and the discrepancy can be as high as 61.5%. Apparently, the more dispersive the jamming energy is, the more obvious deterioration happens in PD.It is not difficult to understand via the “Bernoulli experiment model.” The distribution of jamming power toM cells would create M subevents and M “compound reduction” in PD, which usually descends faster than the “single test.”On the other hand, the increase inQ will raise the false alarm rate, but more M will counteract this effect (see Figure 6(b)). A compromise should be made between those two factors. In order to fully quantify jamming effects, a comprehensive metric Jamming effect factor (JEF) can be defined as(32)JEF≜PD¯PD0¯+μκPFA0¯PFA¯,where the subscript “0” means parameters without jamming, and κ is a scaling factor to balance the scale of improvements on PD¯ and PFA¯. The higher JEF value, the more influences jamming have on CFAR. μ is a weight coefficient to judge the significance of PD¯ and PFA¯. Generally, μ=1, but in some protective weapons where PD¯ is primarily advocated, μ is supposed to be less than 1, whereas in remote early warning devices where PFA¯ should be rigorously controlled to a restricted range, we can set μ>1. To summarize, here lists steps to find optimum parameters for jamming CFAR.Step 1: To give the initial value of false target densityStep 2: To estimateM and Q according to false target densityStep 3: To computePD¯ and PFA¯ through equations (13) and (14)Step 4: To compute JEF from equation (32) and adjust false target densityStep 5: Change the value of false target density and redo Step 1∼4 until JEF convergents to its maximum value


## 5. Conclusion

This paper presents a Bernoulli experiment model to quantitatively describe CFAR performance under the jamming condition. The detection process is modeled as a Bernoulli experiment with a success probability that has a clear physical essence from energy perspective. Different jammers act independently, and effects of multifunctional jammers can be seen as overall effects of several single jammers. BEM is a generalized model and can be transformed into specific occasions.There exist four factors:SIR, Q, M, and JIR that greatly impact CFAR performance. Once jamming pulses enter the reference cells, an improvement in SIR is required to compensate the deterioration. Jamming pulses in CUT will greatly improve the probability of the false alarm to dozens of decibels. Continuous noise jamming is destructive to CFAR with relatively lower JIR.As applications of BEM, two jamming principles are obtained in terms of power distribution of comprehensive multimode jammers and parameters setting principles, where limited power should be rationally distributed to working modes and reference cells.It is recommended that further research be undertaken in the following areas: (1) the joint impacts of jamming and heterogeneous clutter to CFAR performance; (2) detection capabilities and comparisons of different CFAR techniques in jamming environment; (3) the quantitative effectiveness of electronic counter-counter measurements from detection perspective, etc.



---

*Source: 3780252-2022-01-31.xml*