FACTA UNIVERSITATIS Series: Electronics and Energetics Vol. 33, N o 3, September 2020, pp. 459-476 https://doi.org/10.2298/FUEE2003459M © 2020 by University of Niš, Serbia | Creative Commons License: CC BY-NC-ND RELIABILITY ANALYSIS OF DIFFERENT RCIED ACTIVATION SIGNAL RESPONSIVE JAMMING TECHNIQUES AND THEIR COMPARISON TO ACTIVE JAMMING * Mladen Mileusnić, Predrag Petrović, Vladimir Kosjer, Aleksandar Lebl, Branislav Pavić IRITEL a.d., Belgrade, Serbia Abstract. In this paper we compared the time required for the successful jamming of remote controlled improvised explosive devices activation using active and responsive jamming methods. As a representative of active jamming method we analyzed jamming signal generation using frequency sweep. For the analysis of the possible activating signal presence based on responsive jamming procedures we first supposed Fast Fourier Transform (FFT) implementation and compared its analysis rate to the rate of sweep jamming. Taking into account the current technology state, it is proved that the time required to achieve the successful jamming relied on FFT analysis may be less than in the case of active sweep jamming. After that we considered pros and cons for energy detector and matched filter detector implementation in responsive jamming. For these two detector types it is shown how to determine the number of analysis blocks to achieve approximately the same number of collected samples as in the case of FFT implementation, starting from the probabilities of false detection and miss of detection. Key words: Active and Responsive jamming, RCIED - remote controlled improvised explosive devices, Frequency sweep, Fast Fourier Transform, Energy Detector, Matched Filter, Jamming Reliability 1. INTRODUCTION The common characteristic of all remote controlled improvised explosive devices (RCIED) is that they are activated by wirelessly transmitted messages. The results of RCIED activation message could be disastrous regarding people lives (VIP persons) and the equipments damages. All elements related to activation signal characteristics (signal power, frequency, implemented modulation method, message duration) are completely unknown. This fact produces great problems in the realization of RCIED activation jammers.  Received February 4, 2020; received in revised form March 3, 2020 Corresponding author: Aleksandar Lebl IRITEL a.d., 11080 Belgrade, Batajnički put 23, Serbia E-mail: lebl@iritel.com * The earlier version of this paper is awarded as the best one in the section Telecommunications at the 6 th IcETRAN Conference, Silver Lake, 3-6 June 2019, [1]. 460 M. MILEUSNIĆ, P. PETROVIĆ, V. KOSJER, A. LEBL, B. PAVIĆ Contributions [2] and [3] provide a general overview of jammers types, communications jamming requirements and their efficiency analysis. Modern communications jamming principles and techniques may be found in [4]. There are two basic approaches to the jammer implementation. The first one is active jamming, which consisted of continuous predefined jamming signals sending independently of the RCIED activating message characteristics. In this concept there are no „look through“ phases to detect the activation message existence and the jamming signal characteristics are selected in general using previous experience and expectations. The most important freely selected jamming signal parameter is the RF signal level. This level has to be as high as possible to successfully prevent activating message reception. Two key features which are not optimally chosen relate to continuous jamming regardless of RCIED activation message existence and the RF jamming signal level necessary for jamming successfulness due to the fact that the activation signal level is unknown. The alternative approach to jammer implementation is responsive jamming concept. In this case the jamming signal characteristics can be optimized using look through intervals to detect the activation message existence and its level. That’s why it is possible to send the jamming signal only during activation message presence and jamming signal level can be adjusted to the activation message level in order to successfully deny the threat. A wide range of active and responsive jammers may be found in [5]-[14]. It may be concluded from this short presentation of active and responsive jamming characteristics that active jamming is always successful, while responsive jamming efficiency depends on activation message detection reliability. The question is whether responsive jamming reliability may be higher than for active jamming. In this paper we compare the reliability of mostly implemented active jamming method – frequency sweep [15]-[19] to the reliability of a representative method for activating signal eventual presence detection in order to generate jamming signal according to the activation signal characteristics by implementation of Fast Fourier Transform (FFT) in the analysis [14]. A brief principle schematic of RCIED activation signal detection is explained in Section 2. After that the method for RCIED activation signal frequency spectrum estimation based on FFT analysis is presented in Section 3 with the emphasis on the required time for calculation. Sections 4 and 5 deal with the specificities of energy detector and matched filter detector implementation for RCIED activation message detection. The emphasis is on the determination of collected samples number. Section 6 is devoted to frequency sweep jamming and to determination of required time to realize one complete jamming cycle. In Section 7 jamming reliability on the basis of FFT analysis is compared to the frequency sweep jamming reliability, whereby two special purpose processors are considered for FFT calculation. Reliability estimation is based on the required time to allow successful jamming. Section 8 is focused on the presentation how to determine the necessary number of analysis blocks in energy detection or matched filter detection to achieve the comparable sample collection rate of these two detectors to FFT based detector. At the end, the paper conclusion is given in Section 9. 2. PRINCIPLES OF DETECTION PROCESS Main principles of RCIED activation signal detection may be explained using simplified block-schema presented in Fig. 1. Reliability Analysis of Different RCIED Activation Signal Responsive Jamming Techniques... 461 The first phase in detector function is signal samples collecting (block SCOL). After that follows processing of these samples (block PROC). The final step is making a decision about (eventual) presence of RCIED activation signal on the base of a set of comparison rules (block DECISION). These comparison rules are adjusted to the applied method of signal samples processing. This paper is mainly devoted to the block PROC. The analyzed methods are FFT, energy detector and matched filter. When the second or the third of these three methods is implemented, digital filter precedes the phase of processing. According to the available literature, there are also other methods which are less often applied for spectrum sensing, but they are possible candidates for RCIED activation signal detection. Some of them are waveform based detection, eigen-value based detection, wavelet based edge detection, ciclostationary feature detection [20] and so on. These methods, as generally less often applied ones, are beyond the scope of this paper. SCOL PROC DECISION Fig. 1 Block schema of RCIED activation signal detector 3. SIGNAL SPECTRUM ESTIMATION ON THE FFT BASE FFT is the calculation procedure, which allows relatively fast estimation of discretized signal frequency spectrum. Starting from n time samples of analyzed signal, this procedure gives a snapshot of signal frequency spectrum also in n points, i.e. n spectrum lines are obtained. FFT is the optimum method taking into account the required number of mathematical operations for signal spectrum determination. There are (n/2)·log2 (n) complex multiplications and n·log2 (n) complex additions [21]. The limitation for n is that the condition n=2 a must be satisfied, where a is the positive integer number. This is a significant saving in the number of mathematical operations and in the required calculation time comparing to the classical method of frequency spectrum estimation by Discrete Fourier Transform (DFT). Namely, it is necessary to perform n 2 complex multiplications and n 2 –n complex additions to obtain n frequency spectrum components by DFT on the base of n time samples. Let us suppose that fs is the frequency of analyzed signal sampling. The sample acquisition time is then: s n t f  (1) The frequency resolution on the base of sample acquisition time may be determined as: 1 sf df T n   (2) Therefore, frequency resolution is improved when acquisition time is increased, i.e. the space between frequency spectral components of the analyzed signal is lower. 462 M. MILEUSNIĆ, P. PETROVIĆ, V. KOSJER, A. LEBL, B. PAVIĆ Constant advancements in processor realization technology and mathematical algorithm improvements are visible in two aspects of FFT calculation progress. On the one side the number of points in which frequency spectrum is determined is constantly increased, and on the other side the required FFT calculation time for some exactly determined number of frequency spectrum components is constantly decreased, chronologically, successively according to presentations in [22], [23], [21], [24]. We selected two approaches referenced in [21] and [24] due to very fast processing algorithms. Data presented in [24] is related to the FFT calculation time as a function of the number of signal time samples implemented for FFT calculation, i.e. as a function of the obtained frequency components number in the analyzed spectrum. The presented data is for processor clock of 1GHz. It is further emphasized in [24] that improvement may be achieved by processor clock speed increase to 1.25GHz. Besides, it is stated in [25] that maximum processor clock frequency may be even 1.4GHz. On the base of these data, the FFT calculation time (Tcal in ms) is presented in Table 1 as a function of the number of points used in a calculation, for a processor clock of 1.25GHz and for 8 processor cores. The value of the constant K is 1024 in the first column of the Table 1. The time of FFT calculation (Tcal in μs) according to the data emphasized in [21] is presented in the Table 2. The processor clock in this case may be in the range between 60MHz and 150MHz [26]. That’s why data are presented for the mean processor frequency of 100MHz. FFT hardware accelerator (HWAFFT) is one of the parts in the processor implemented according to [21]. HWAFFT is intended for faster FFT calculation. Data in Table 2 are related to the case when HWAFFT is implemented. The number of points is relatively small (till 1024) where FFT is calculated comparing to the number of points, where FFT results are presented in Table 1. In accordance to Fig. 1, the total time, which is needed for signal analysis in a jammer (Tan) before (eventually) starting RCIED activation jamming signal emission, consists of three components: sample acquisition time (T), FFT calculation time (Tcal) and the time, which is necessary to compare obtained signal frequency components after FFT calculation (Tcomp) in order to determine whether it is necessary to start jamming. When considering the last component (Tcomp), there is not such a data in a literature, because calculation is very specific. For our analysis, we supposed that taking equal values of Tcal and Tcomp is a quite good approximation, i.e. 2an cal comp cal s n T T T T T f       (3) Table 1 The time of FFT calculation as a function of the number of calculation points for the processor presented in [24] Number of points for FFT calculation Calculation time Tcal [ms] (8 cores, 1.25GHz) 16K 0.1051 32K 0.1584 64K 0.2517 128K 0.5128 256K 0.9488 512K 2.4824 1024K 5.1226 Reliability Analysis of Different RCIED Activation Signal Responsive Jamming Techniques... 463 Table 2 The time of FFT calculation as a function of the number of calculation points for the processor presented in [21] Number of points for FFT calculation Calculation time Tcal [μs] (with HWAFFT, 100 MHz) 8 1.3 16 1.7 32 3.21 64 4.36 128 9.12 256 16.68 512 37.4 1024 73.15 4. RCIED ACTIVATION SIGNAL DETECTION BY ENERGY DETECTOR Energy detector is the simplest techniques for signal detection [27]. In the same time it is a very often applied technique. It is necessary first to measure signal energy in the pre-defined frequency band. The measured signal energy is then compared to the energy threshold according to the equation 2 1 ( ) ( ( )) N n E x x n     (4) where N is the number of samples implemented for signal energy estimation, x(n) is the amplitude of n th sample and γ is the threshold power. Although simple for implementation, energy detector performances are degraded due to noise uncertainty (noise level is variable during time) and background interference [28]. Noise uncertainty may be bounded or unbounded [29]. As a consequence of noise uncertainty, the detection by energy detector may become even impossible under relatively low value of signal to noise ratio (SNR) [30]. In other words, there exists a SNR wall: detection is possible only when signal power is higher than noise power uncertainty. For the analysis in this paper and for the comparison of energy detector characteristics with the characteristics of other methods for reactive jamming the most important parameter is the number of samples (N) to achieve necessary detection reliability. Our analysis is based on the formula for N from [27] [31]: 1 1 2 2 ( ( ) 1 2 ( )) 2 f dQ P SNR Q P N SNR         (5) where Pf is probability of false detection (detector announces signal presence although there is no signal), Pd is probability of successful detection and Q -1 is inverse Gaussian-Q function. In other words, Q -1 is the inverse of 2 1 21 ( ) 2 u x Q x e du        (6) 464 M. MILEUSNIĆ, P. PETROVIĆ, V. KOSJER, A. LEBL, B. PAVIĆ 1,E+01 1,E+02 1,E+03 1,E+04 1,E+05 1,E+06 1,E+07 -20 -19 -18 -17 -16 -15 -14 -13 -12 -11 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0 SNR(dB) N Pf=Pmd=0.1 Pf=Pmd=0.01 Pf=Pmd=0.001 Pf=Pmd=0.000001 Fig. 2 The necessary number of samples (N) when energy detector is applied as a function of signal-to-noise (SNR) ratio without noise uncertainty for different values Pmd =Pf. One additional important parameter in energy detection systems analysis is miss in detection Pmd (detector does not detect a signal although it exists). Probabilities Pd and Pmd are connected by the equation 1md dP P  (7) Fig. 2 presents the necessary number of samples (N) as a function of signal-to-noise ratio (SNR). The results are presented for equal values of Pf and Pmd. There is no noise uncertainty which means that optimum detector threshold value exists independently of SNR. For small SNR signal detection is always possible, but the value of N significantly increases. In our concrete implementation it is more important to achieve low value of Pmd than to achieve low value of Pf. In other words, consequences of miss in detection are more severe (RCIED is activated because there is no jamming) than if the detection is false (only jamming signal is waste generated). That’s why the results for probability values satisfying the condition Pmd1, is subtracted from the value SNR in the denominator of equation (5). This means that noise power, instead of having power equal to σ 2 when noise is completely defined, now has the value of power between (1/ρ)·σ 2 and ρ·σ 2 . SNR-wall is presented by the fourth graph in Fig. 3. Even for a very small noise uncertainty value ρ=0.25dB or ρ=1.059 when it is Pf=0.1 and Pmd=0.01 the value of N tends to infinity for SNR~-9.3dB and below this value -9.3dB it is not possible to detect a signal. As a conclusion it may be said that it is very important to constantly monitor the noise level and to adjust threshold value according to instantaneous noise level and in this way to avoid SNR-wall appearance. 5. RCIED ACTIVATION SIGNAL DETECTION BY MATCHED FILTER The second often implemented technique of spectrum analysis is based on the method of matched filters. The main property of such filters is that they are optimum linear filters applied for signal detection in white Gaussian noise, meaning that maximum SNR is achieved by their implementation [32]. Although this property contributes to easier and faster signal detection, the drawback of matched filter implementation is that it is necessary to precisely know time characteristics of the signal which has to be detected. Such knowledge is possible in some implementation areas, as for example, in cognitive radio [32], [33]. But, if considering RCIED activation signal jamming, there is a great variety of possible and, in the same time, unpredictable activation techniques. They usually depend on the devices which may be easily purchased in some country (region) and easily adapted for its malicious function. The number of applied solution types is not great in the analysis presented in [34] with the dominant implementation of one type, thus simplifying and limiting the necessary number of different matched filters. Nevertheless, application of matched filters is not quite suitable for RCIED activation signal detection and the analysis of this method has more theoretical than practical significance. Similarly to the analysis of energy detector, the necessary number of samples to achieve the desired probability of false alarm and probability of successful signal detection may be determined on the base of equation from [27]: 1 1 2 ( ( ) ( ))f dQ P Q P N SNR     (8) Fig. 4 presents the necessary number of samples (N) as a function of signal-to-noise ratio (SNR) when matched filter is implemented. The results are obtained by equation (8) and are presented for equal values of Pf and Pmd. After that, Fig. 5 presents the corresponding results if Pmd