CC BY
5
Visnyk NTIJU KP1 Seriia Radiolekhnika tiadioaparat.obuduuannia, "2023, Iss. 92, pp. 54—59
UDC 621.3:519.2
Method of Complex Envelope Processing for Signal Edges Detection
Buhaiov M. V.
S. P. Korolov Military institute, Zhytomyr, Ukraine E-mai 1: karunen&ukr. ne I.
Problem statement. The need of information processing automation in modern radio monitoring systems stimulates development of flexible methods for signal detection and its parameters estimation in time domain. A priori uncertainty of signal time-frequency structure complicates the automatic determination of signals edges. Purpose. The purpose of the article is subsequent automation of radio frequency spectrum analysis process by developing and implementing a method for determining signals time edges under conditions of a known noise power and signal-t.o-noise ratio. Method. To determine signal time edges in given frequency channel, square of signals' complex envelope is first calculated, smoothed with moving average window and compared with threshold. Threshold is calculated as a quantile of gamma distribution using Wilson-Hilferty approximation of x2 distribution quantiles for a given probability of false alarm. An analytical expression is obtained for calculation length of moving average window depending on signal-t.o-noise ratio. An algorithm has been developed for determining signals' time parameters and filtering them by duration. Unknown noise power value in frequency channel can be replaced by its estimate under the assumption that frequency channel is not constantly occupied and noise level is estimated on signal-free time intervals. Conclusions. Proposed method makes it possible to automatically determine edges of signal with an arbitrary structure at signal-to-noise ratio values from -6 dB. Adjustable length of moving average window makes it possible to reduce the error in determining signal time parameters by 2-4 times with an increase in the signal-t.o-noise ratio compared to a fixed window length. Prospects for further research in this direction should be focused on development and implementation of methods for detection signal edges under conditions of an unknown noise level.
Keywords: signal edges: complex envelope: moving average: threshold: signal structure: signal-to-noise ratio DOI: 10.20535/RADAP. 2023.92.54-59
Introduction
Rapid development of wireless radio systems fl] requires automation of information processing in modern radio monitoring systems [2.3]. In particular, when conducting radio monitoring, it is necessary in real time to determine radio signals edges in order to establish the moment of start and end of recording in given frequency channel. Also, this problem arises when analyzing recordings of wide frequency bands, especially for processing pulse signals or signals transmitted by packets, when it is necessary to remove noise areas between signals to recognize the type of modulation and further processing. A priori uncertainty of signal time-frequency structure makes it difficult to automatically determine signal edges. Therefore, development of flexible methods for detection and estimation signals' time parameters is actual scientific and technical task.
1 Related works
The issue of detection pulse signals in time domain and estimation their parameters are considered in many works by researchers from all over the world. In particular, in [4] it is proposed to determine pulses duration and their repetition period using Haar filter and ratio of sliding sums, and in [5] signals edges are detected using continuous wavelet transform. Estimation of signals time parameters by smoothing its instantaneous energy using window functions is proposed in [6]. In [7] image processing methods were used to detect pulses edges. In [8] it is proposed to estimate signals time parameters by analyzing time-frequency transforms. The method works with a minimum signal-to-noise ratio (SNR) of -5 dB. To determine pulse signals parameters it is proposed in [9] to use a neural network. In [10] sequential analysis in form of a cumulative sum algorithm was used to detect pulse signals. In [11.12] rough detection of signal begin and end is carried out using a spectrogram with further refinement of its edges using a special filter [13]. To
Method of Complex Envelope Processing for Signal Edges Detection
55
determine signal edges it is proposed in [14] to calculate the smoothed complex envelope of signal.
As a result of latest research and publications analysis, it was established that the general drawback of considered methods is using algorithms with fixed parameters, which in case of known SNR value, does not make it possible to achieve potential accuracy of signal time parameters estimation. Also remains unsolved the task of automatically determining signals edges with an unknown time-frequency structure.
2 Problem statement
The purpose of the article is subsequent automation of radio frequency spectrum analyzing process by-developing and implementing a method for determining signals edges under conditions of a known noise power and signal-to-noise ratio.
3 System model
If the noise level is known, its value must be used to calculate threshold for signals edges detection. Noise power can be estimated using known methods [15], even with unknown occupancy of analyzed frequency band. Radio signals to be processed in receiving frequency channels can have coherent (with regular phase structure) or non-coherent (noise-like) structure. Vector of complex received signal in selected frequency channel may be written in following form:
x = s+e
£ = 6 + jÇq,
= \jah + °\Q = ^fc = ^«Q .
(4)
To determine signal edges in given frequency channel, we will calculate signal complex envelope as follows:
Ex= x2+xQ. (5)
This approach to the calculation of complex envelope will simplify further mathematical calculations, since E^ samples for noise will be subject to exponential probability density function (PDF).
Then threshold for probability of false alarm Pp is calculated according to the following equation:
7 = <7? ln(-PF) .
(6)
When smoothing complex envelope E^ with moving average window of L samples length, vector of smoothed samples C will be subject to gamma distribution [16] and will be described by the following
expression:
r,h-1 '
P C = -
- C Ll
tw
(7)
Threshold value 72 for smoothed complex envelope C can be calculated as a quantile of gamma distribution of p =1 — Pp level. To do this, we will use the qnantile approximation of gamma distribution using the chi-square distribution x2 [ ]:
lp = xl (2(L + 1)):
(8)
where (2 (L + 1)) is p level quantile for x2 distribution with K = 2 (L + 1) levels of freedom.
Using the Wilson-Hilferty approximation of x2 distribution quantiles [ ], threshold value can be calculated using the following expression:
72 =
1
2
9K
+ u
(9)
where
(1)
Vector of complex signal samples can be written as vectors sum of in-phase and quadrature signal components:
S = si + jsQ, (2)
where j = \f—T - imaginary unit.
By analogy, vector of complex Gaussian noise samples we will write in following form:
1, 24+0, 85H0'657 1+0, 0001F-3 + 2|8 :
H =-0, 96 ln
(1-iv).
(3)
where ^ - vector of in-phase noise component; -vector of quadrature noise component.
Vectors ^ mid £q are also normally distributed with zero mean and equal standard deviations (STD) a^j = . Then vector £ will also have zero mean and its STD calculates according to next equation:
(10)
Value of moving average window length L depends on SNR value and should be increased when SNR decreases. Therefore, signal processing in each frequency channel will depend on SNR. This approach will ensure optimal distribution of computing resources, since smoothing signal envelope requires additional calculations.
Consider case when signal has an incoherent structure, rectangular envelope and differs from noise only in power. With low values of SNR, L must be chosen large (hundreds). In this case, PDF of C can be considered normal and threshold can be calculated according to the following expression:
72 = mc* + aoc,,
where mC( = <t| - mean of C^; aC(
(ID
7L " STD of
C^; a - coefficient, which value depends on the false alarm probability.
Then L value will be determined on the condition that the minimum values of smoothed complex envelope of signal and noise mixture will be greater than threshold. Then for smoothed complex
u =
56
Buhaiov M. V.
envelope can be written similarly to expression (11):
72 = mCx + aaCx, (12)
where mc, = a2s + - mean of signal and noise mixture smoothed complex envelope:
as = \Ja"si + asQ ~ STD of complex signal;
ffc, = ^~ STD signal and noise smoothed complex envelope.
Equating expressions (11) and (12). we get:
_ 2
1=(13)
where q - SNR by power (not in dB).
Figure 1 shows graphs of moving average window length L dependences via SNR for different values of a parameter. From these curves, it can be seen that at high SNR (20 dB and more), value of window length L approaches to a2. Resulting L value can be used to smooth signals with any structure.
io3
102
101
\ - a = 3
\ \ V " \ v \ \ w \\ N. \ \ \ V \ — a = i -— a =
1 = 25 S. N \ V \ N N s.
L — 16
L=9
5 10
SNR, dB
15 20
Fig. 1. Moving window length dependence via SNR
begin
x, PF, q, W
Ex calculation (eq. 5)
I
L calculation (eq. 13)
Cx calculation 4 1
y2 calculation (eq. 11)
A = arg(Cx > y2) u 1
B = arg(a,+i -di = 1)
C = arg(£,+i - > 1)
dt = B(C, + 1), €j = B(c;)
begin = (bj, A(£))
end = (A(e,), bN)
T
10
11
Fig. 2. Block diagram of signal edges detection algorithm
Figure 2 shows block diagram of signal edges The following input data are required for algorithm detection algorithm in frequency channel. For a non- implementation: vector of received signal samples x, zero Pp value, it is possible to incorrectly determine noise STD Pp, SNR value ^d W (block 1). signal edges. As Pp decreases, probability of signal edge When calculating length of moving average window detection will also decrease. Therefore, it is advisable L according to expression ( ), value of a coefficient to choose exceeding specified threshold by a certain is recommended to be chosen depending on SNR: at number W of consecutive samples as a criterion for q > 12 dB - a = 3; at 12dB > q > 5dB - a = 4 and detection signal edges. at q < 5 dB - a = 5. In block 6 is formed a vector
4 Algorithm
Mot(w iin:iiiali(;iiiiH >iax:ouux mo>k paauociiriia-'iiu ihjihxom aiia.;ii:sy kdmujioktuoi o6ui,aiioi
57
of received signal samples numbers A that exceeded threshold In blocks 7-11 vector A is processed, as a result of which two new arrays of samples numbers are formed: begin - array of numbers of beginnings and end - endings of signal samples sequences. In blocks 12-14, values of begin and end arrays are filtered and only those signals whose duration is at least W samples remain.
If noise STD is unknown, then in order to detect signal edge in frequency channel, the unknown value of noise power can be replaced by its estimate, assuming that the frequency channel is not constantly occupied and noise power is estimated clt SI gnal-free time intervals. Moreover, method of noise power estimation should provide its' reliable estimate with a sample of small size. In this case, the signal may be present or absent in frequency channel at moment of observation. Then time of signal start (end) may be estimated by a sharp change in power level in given frequency channel.
However, only the known value of noise STD will not make it possible to determine SNR and calculate optimal L value. In this case, it is necessary to choose L and a values for the worst expected case, or based on such a minimum value of SNR, at which further
0 20 40 60 SO 100
t, ms
(a)
Fig. 3. Results of coherent signal
processing of signal (recognition, demodulation) will be impractical.
5 Simulation results and discussion
To study developed method were used signals with two types of structure: coherent a pnlse with a rectangular envelope and non-coherent (like OFDM) a fragment of white Gaussian noise also with a rectangular envelope. During simulation, sampling frequency was 10 MHz, carrier frequency of coherent signal was 1,6 MHz. Duration of studied signal and noise mixture fragment was 100 ms. Signal starts at 20 ms and ends at 40 ms. Probability of a false alarm during simulation was set at the level of 10-4.
Figure 3 shows results of signal processing with coherent structure in time domain for frequency channels with OdB (Fig. 3a) and 10 dB (Fig. 3b) SNR. For the first case, moving average window length was L = 146 at a = 5, and for the second - L = 21 at a = 4.
0 20 40 60 SO 100
t, ms (b)
ing for SNR 0 dB (a) and 10 dB (b)
Figure 4 shows results of signal processing with channels with SNR of OdB (Fig. 4a) and 10dB a noise-like structure in time domain for frequency (Fig. 4b).
(a) (b)
Fig. 4. Results of noisy signal processing for SNR 0 dB (a) and 10 dB (b)
58
Вугайов M. В.
Comparing Fig. 3 and Fig. 4 wo can conclude that for reliably signal edges determination, it is necessary to smooth complex envelope, regardless of signal structure and SNR. Moreover, for a coherent signal at high SNR (above 10dB). we can choose value of a parameter less than 2. At the initial stage of signal processing, as a rule, it is not known what structure signal has. therefore it is recommended to choose value of a parameter in accordance with recommendations from previous section.
Figure 5 shows the dependence of signal edges estimation error via SNR. which changed during the experiment from -6dB to 20 dB with a step of 1 dB.
a= 5 a = var
—
4
-5 0 5 10 15 20
SNR, dB
Fig. 5. Dependence of signal edges estimation error via SNR
It can be seen from this figure that the error of signal edges estimation decreases with the increase of SNR. This is explained by decrease of moving average window length L. In addition, at medium and high SNR values due to the adjusting a parameter, which determines L, error of signal edges estimates is reduced by 2-4 times, compared to fixed a.
Conclusions
Proposed method allows automatically determine signal edges with an arbitrary structure under the conditions of known SNR and noise power. Dependency between SNR and parameters of moving average window were established. It gives possibility to reduce by 2-4 times error of signal edges estimates in case of SNR increasing compared to fixed length of window with moving average. Developed algorithm can be implemented in existing and prospective radio monitoring systems.
Prospects for further research in this area should be focused on development and implementation of methods for determining signal edges under conditions of an unknown noise level and the development method for estimation noise level using samples of restricted volume.
References
[1] Liang Y.-C. ("2020). Dynamic Spectrum Management. From Cognitive Radio to Blockchain and Artificial Intelligence. Springer, 180 p. doi: 10.1007/978-981-15-0776-2.
[2] Elmasry F. G. (2021). Dynamic Spectrum Access Decisions. Local. Distributed. Centralized, and Hybrid Designs. Wiley, 728 p.
[3] Captain К. M., Joshi M. V. (2022). Spectrum Sensing for Cognitive Radio. Fundamentals and Applications. CRC Press, 256 p. doi: 10.1201/9781003088554.
[4] Ranney К.1., Tom K., Tadas D.. Tesny N.. Magill A., and Diehl W. (2022). Magnitude-based pulse width estimation via efficient edge detection. .Journal of Applied Remote Sensing, Vol. 16. Iss. 1. doi: 10.1117/1..IRS.16.016509.
[51 Bruno .I.M.. Mark B.L.. Tian Z. (2016). An Edge Detection Approach to Wideband Temporal Spectrum Sensing. IEEE Global Communications Conference, pp. 1-6. doi: 10.1109 / G L О С О M .2016.7841584.
[6] Adam A. A.. Adegboye В.Л.. Ademoh LA. (2016). InterPulse Analysis of Airborne Radar Signals Using Smoothed Instantaneous Energy. International .Journal of Signal Processing Systems, Vol. 4. No. 2. pp. 139 143. doi: 10.12720/ijsps.4.2.139-143.
[7] Li C.. et al. (2020). Pulse Train Detection Algorithm Based on Edge Enhancement. icWCSN '20: Proceedings of the 2020 International Conference on Wireless Communication and Sensor Networks, pp. 30 34. doi: 10.1145/3411201.3411202.
[8] Ahmad A. A.. Lawan. S.. Ajiya. M.. Yusuf. Z. Y.. & Bello. L. M. (2020). Extraction of the pulse width and pulse repetition period of linear FM radar signal using time-frequency analysis. .Journal of Advances in Science and Engineering, Vol. 3. No. 1. pp. 1 8. doi: 10.37121/jase.v3il.69.
[9] Zhu M. Zhang Z.. Li C.. Li Y. (2021). .IMRPE-Net: Joint modulation recognition and parameter estimation of cognitive radar signals with a deep multitask network. IET Radar, Sonar & Navigation, Vol. 15. Iss. 11. pp. 1508 1524. doi: 10.1049/rsn2.12142.
[10] Fridman P. A. (2010). A method of detecting radio transients. Monthly Notices of the Royal Astronomical Society, Vol. 409. lss.2. pp. 808 820. doi:10.1111/j.l365-2966.2010.17346.x.
[11] Fan X.. Li Т.. Su S.. Chen Z. (2015). A new method of pulse edge detection in low SNR. Proc. SPIE, Vol. 9812. id. 981213. 7 p. doi: 10.1117/12.2205763.
[12] Millioz F.. Martin N. (2006). Detection of the boundaries of time-frequency patterns with a reassignment-based method. Seventh International Conference on Mathematics in Signal Processing, UK. 4 p.
[13] Liu X.. Fan X.. Su S. (2019). Adaptive pulse edge detection algorithm based on short-time Fourier transforms and difference of box filter. .Journal of Applied Remote Sensing, Vol. 13. Iss. 2. 29 p. doi:10.1117/l..IRS.13.024502.
[14] Buhaiov M. V. (2022). Metod vyznachennia chastotno-chasovykh mezh zainiatykh dilianok radiochastotnoho spektra [Method for determining occupied time-frequency spectrum regions]. 2022: Proceedings of the Sixteenth International Scientific and Technical Conference "MODERN CHALLENGES IN TELECOMMUNICATIONS", Kyiv. Ukraine, pp. 109-111.
Метод ви:шачо1шя часових меж ралюсигнашв шляхом aiia.;ii:sy комплексно! обвщшн
59
[15] Buhaiov M.V. (2022). Iterative Method for Noise Power Estimating at Unknown Spectrum Occupancy. Visnyk NTUU KP1 Seriia - Radiotekhnika
, Vol. 88, pp. 42—49. doi: 10.20535/RADAP.2022.88.42-49.
[16] Kay S. M. (2013). Fundamentals of Statistical Signal Processing: Practical Algorithm Development. New York: Prentice Hall, 475 p.
[17] Kobzar Л. 1. (2006). Prikladnaya matematicheskaya st.at.i-stika. Dlya inzhenerou i nauchnyh rabotnikou [Applied mathematical statistics. For engineers and scientistsJ, Moscow: FIZMATLIT, 816 p.
Метод визначення часових меж радю-сигнал!в шляхом анал!зу комплексно!" обвщно!"
Бугайов М. В.
Постановка задач!. Необх1дшсть автоматизацп процеов оброблеппя шформацп в сучаспих системах радюмошторипгу стимулюе розроблешш гпучких мето-д!в ВИЯВЛС1ШЯ та оцшюваппя часових параметр!в сиг-пал!в. Anpiopua певизпачешсть щодо частотпо-часово! структури сигпалу ускладшое автоматично визпачеппя часових меж сигпал1в. Мета статт. Метою статт! е подалына автоматизагця процесу апал!зу радючастотпо-го спектра за рахупок розроблеппя та реал!зацГ! методу визпачеппя часових меж радюсигпал!в в умовах в!до-мого значения потужпост шуму та в1дпошеппя сигнал-шум. Виклад матер!алу до с jii дженни. Для визпа-чешш часових меж сигпал!в у кожному 1з частотпих
капал!в розраховуеться квадрат комплексно! обв1дпо! сигналу, згладжуеться 1з внкористашшм в!кпа ковзаю-чого середнього та пор1вшоеться 1з порогом. Зпачеппя порогу розраховуеться як кваптильгамма-розподшу 1з використаппям апроксимацп В1лсопа-Х1лферт1 квапти-л!в розподшу х2 Для задано! ймов!рност! хибно! тривоги. Отримаио апал1тичпий вираз для розрахупку довжшш вшпа ковзаючого середпього в залежпост! в!д в!дпо-шмшя сигпал-шум. Розроблепо алгоритм визпачмшя часових параметр!в сигпал1в та !х фшьтрацп за триваль стю. Нев1доме зпачмшя потужпост! шуму в частотному капал! можпа замшити його оцшкою в припущепш. що частотпий капал зайпятий по постишо 1 па вглышх в!д сигпал1в штервалах часу проводиться оцшюваппя рш-пя шуму. Висновки. Запропоповапий метод дозволяв автоматично визпачати часов! меж! сигпал!в з довгль-пою структурою при зпачешшх в1дпошеш1я сигпал-шум в!д -6 дБ. Змшпа довжипа в!кпа ковзаючого середпього дозволяв у 2-4 рази змепшити помилку визпачеппя часових параметр!в сигпалу при збглынепш в1дпошеп-пя сигпал-шум у пор1впяпш з ф1ксовапою довжипою в!кпа. Перспективи подальших досл!джепь у дапому папрямку догцлыю зосередити па розроблепш та реаль зацп метод!в визпачешш часових меж сигпал1в в умовах пев!домого р!впя шуму.
Клюноог слова: часов! меж! сигналу: комплексна обв!дпа: ковзаюче середпе: пор!г: структура сигналу: в1дношмшя сигпал-шум
