Software for statistical processing and modeling of a set of synchronously registered cardio signals of different physical nature Serhii Lupenkoa, Iaroslav Lytvynenkoa, Andrii Sverstiukb, Andrii Horkunenkob and Borys Shelestovskyia a Ternopil Ivan Puluj National Technical University, Ruska Street 56, Ternopil, 46000, Ukraine b I. Horbachevsky Ternopil National Medical University, Maidan Voli 1, Ternopil, 46000, Ukraine Abstract It has been developed the software complex, that allows to perform mutual statistical processing of synchronously registered cardiosignals on the basis of the vector model of the cyclic rhythmically connected stochastic processes. The developed software provides the ability to form a discrete rhythm function, depending on the number of cycles and zones of the cardio signal to be modeled. The developed set of programs allows to perform statistical processing of synchronously registered cardiosignals of different physical nature and to obtain estimates of mathematical expectation, dispersion, autocorrelation and mutual correlation functions taking into account the rhythm function and the period. The software also allows you to simulate cardio signals, using the loaded input data of the cardiosignal, rhythm functions, the number of cycles and zones of the signal to be modeled. Keywords 1 Synchronously registered cardiosignals, statistic processing methods, software complex 1. Introduction Cardiovascular diseases in the world are a serious problem. The medical-social burden of the circulatory system diseases, especially exacerbated by the rapid spread of COVID-19 and is that they significantly affect the duration and life quality, the loss economic potential indicators of countries. That is why the fight against diseases of the circulatory system at the present stage is a priority problem of modern medicine [1, 2]. Therefore, the development of modern technologies for instrumental diagnosis of the cardiovascular system (CVS) pathology in the early disease stages, monitoring the effectiveness of preventive and curative measures, as well as monitoring the vital functions of the body in critical conditions is an urgent problem of modern medicine and technology [3, 4]. A promising direction to overcome this problem is to improve existing methods of instrumental diagnosis of CVS by developing new diagnostic and prognostic features based on the choosing adequate mathematical analysis methods of synchronously registered signals. To date, a significant number of functional methods for studying the state of the CVS, which uses the synchronous registration of several identical or different in origin cardiosignals. One of the main ones is the synchronous registration of the electrocardiogram (ECG) in different leads 5 and in standardized formats [6]. CMIS-2021: The Fourth International Workshop on Computer Modeling and Intelligent Systems, April 27, 2021, Zaporizhzhia, Ukraine EMAIL: lupenko.san@gmail.com (S. Lupenko); iaroslav.lytvynenko@gmail.com (I. Lytvynenko); sverstyuk@tdmu.edu.ua (A. Sverstiuk); horkunenkoab@tdmu.edu.ua (A. Horkunenko); kaf_vm@tu.edu.te.ua (B. Shelestovskyi) ORCID: 0000-0002-6559-0721 (S. Lupenko); 0000-0001-7311-4103 (I. Lytvynenko); 0000-0001-8644-0776 (A. Sverstiuk); 0000-0002- 2021-006X (A. Horkunenko); 0000-0002-5606-6290 (B. Shelestovskyi) © 2021 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). CEUR Workshop Proceedings (CEUR-WS.org) When monitoring the vital functions of the body, which occurs during the provision of emergency and urgent medical care, synchronously recorded key indicators of the body: ECG, curves of arterial and central venous pressure, external respiration, oxygen and carbon dioxide in the blood and others, depending on specific clinical situation and capabilities of the applied equipment 7, 8. Mechanical devices to support blood circulation in advanced heart failure are considered in the work [9]. However, and in this situation the diagnostic analysis is conducted on each signal in particular which is compared with age norm or magnitude which is defined and corrected by the doctor. An urgent problem today is the improvement of existing methods of instrumental diagnosis of the cardiovascular system by developing a software package for the analysis of cardio signals. A smart system for monitoring and predicting heart disease, which is based on in-depth training, is presented in the work [10]. The article [11] presents cardio diagnostic systems that use artificial intelligence and machine learning, as it concerns the health of the cardiovascular system. Cyclic analysis between meteorological rhythms of respiratory and cardiovascular diseases is considered in [12]. Diseases of the circulatory system under the influence of ionizing radiation are considered in this work [13]. Algorithms and methods of cardiac signal processing are created on the basis of their adequate mathematical models. Thus, it can be stated that today the methods of compatible automated processing of the set of synchronously registered cardiosignals (SRCS) have not received significant development. This is due to the insufficient level of development of unifying ideas in the construction of mathematical models of different types of cardiosignals and diversity of methods for their processing. The development of other highly informative diagnostic methods also had a significant impact. However, it can be stated that the development of optimal ways to analyze SRCS or vital parameters of the organism with access to new integrated indicators and empirical confirmation of their higher diagnostic value, compared to existing ones, can raise to a new level the nature of diagnostic process and health management treatment. Despite significant advances in mathematical modeling of cyclic heart signals, their existing models do not take into account the following facts. For example, new diagnostic features for cardio diagnostic systems based on the model of electrocardiographic signals in the form of a vector of cyclic random processes are considered in [14]. The above arguments indicate the relevance of software for statistical processing and modeling of cardiosignals based on a mathematical model in the form of a vector of synchronously registered cardio signals for the needs of automated cardiodiagnostics. The effectiveness of modern cardiodiagnostic systems largely depends on the hardware and software components on which they are based. The use of various processing methods included in the software significantly expands the functionality and increases the reliability of human heart diagnosis [15]. While the methods are based on mathematical models that establish the possibility and specifics of processing. In particular, the simultaneous processing of cardiosignals of different physical nature can be performed only if their mathematical models are in some way correlated and have a similar structure. 2. Analysis of recent researches A stochastic approach to their simulation is mainly used in the direction of cardiac signal processing. Thus, for this purpose, use the mathematical models of cyclostationary signals [16-18], modeling diagnostic features [19-20], analytical approaches for myocardial fibrillation signals [21], modeling the pulse signal by wave-shape function [22]. An interesting approach to the description of cardio signals is the mathematical model presented in [23]. It allows to take into account both morphological features of cardiac signals and their rhythmic structure. It is important that this mathematical model takes into account the common rhythm of synchronously recorded cyclic cardio signals of different physical nature, which do not allow to implement other mathematical models using a stochastic approach. In order to be able to apply the proposed mathematical model, a necessary condition is the availability of information about the characteristic elements of cardio signals (cycles and zones, such as zones P, Q, R, S, T - electrocardiogram). It is possible to obtain tau information about the characteristic elements of cyclic cardiac signals by applying their segmentation methods, for example, the methods presented in [24], other methods of analysis of cardio signals are presented in [25-27]. 3. Objectives The work is devoted to the creation and application of software for statistical processing and modeling of synchronously registered cardiac signals as components of cardio diagnostic systems. The importance of these studies lies in the possibility of using a common rhythm in the modelling and static processing of different in nature cyclic signals but formed by a single source. 4. Method According to paper [22], we give a definition of the vector of cyclic rhythmically connected stochastic processes. Definition 1. If there is a function such as T t , n , which satisfies the conditions of the rhythm function that finitely measurable vectors {i1 ( , t1 ), i2 ( , t 2 ),..., ik ( , tk )} and {i1 ( , t1  T (t1 , n)), i2 ( , t2  T (t2 , n)),..., ik ( , tk  T (tk , n))} n  Z, i1 ,..., ik  1, N , where t1 ,..., tk - multiple separability of the vector Θ , t  , for all the integers k  1 are stochastic N equivalent in the broadest sense, we will call the vector Θ N  , t  of cyclic stochastic processes  ____    i  , t , i  1, N ,  ∈Ω, t ∈W  as the vector of strictly rhythmically connected stochastic   processes and the processes as strictly rhythmically connected. Area of definition W vector of cyclic rhythmically connected stochastic processes can be as ordered discretely W  R  {t ml  R , m  Z, l  1, L} or continuous W  R set of real numbers. In the case of discrete domain definition W  D for its elements if m2  m1 , or if m2  m1 , а l 2  l1 , in other cases t m1l1  t m2l2 ; m1 , m2  Z , l1 , l2  1, L there is a type of linear ordering: t m1l1  t m2l2 . Moreover 0  t m , l 1  t m , l   .   The rhythm function T t , n determines the law of changing the time intervals between the single-phase values of the vector of cyclic rhythmically connected stochastic processes. The function of the rhythm satisfies such conditions according to the theorem which is proved in the paper: а) T t , n  0 , if n  0 ( T t ,1  ∞ ); b) T t , n  0 , if n  0 ; с) T t , n  0 , if n  0 , t ∈W ; (1) for any t1 ∈W and t2 ∈W , for which t1  t2 , for function T t , n should be performed a strict inequality: T t1 , n  t1  T t2 , n  t2 , ∀ n ∈Z ; (2) function T t , n is the smallest modulo ( T t , n  ≤ T t , n  ) among all such functions T t, n,  ∈Γ, which satisfy (1) and (2).  In the partial case, if the rhythm function is T t , n  n T T  0, n ∈Z , we will call the vector Θ N  , t  as the vector T - periodically connected stochastic processes. Let us consider the properties of some probabilistic characteristics of the vector Θ N  , t  cyclic rhythmically connected stochastic processes. So, for its compatible k - measurable distribution function takes place equation: Fk ... ( x1 ,..., xk ; t1 ,..., t k )  i1 k1  Fki ... k ( x1 ,..., xk ; t1  T (t1 , n),..., t k  T (t k , n)), n  Z, i1 ,..., ik  1, N , t1 ,..., t k  W (3) 1 1 k Combined central moments function of order p  R : j 1 j rpi ...i (t1 ,..., t k )  M{( i1 ( , t1 )  mi1 (t1 )) R1  ...  ( i p ( , t k )  mik (t k )) Rk }  1 (4) k  rpi ...i (t1  T (t1 , n),..., t k  T (t k , n)), t1 , t 2 ,..., t k  W, i1 ,..., ik  1, N , n  Z. 1 k 5. Results and Discussion It was made the series of experiments on processing of the cardiaosignals of the same and different physical origin which were investigated for the purpose of approving of the greater effectiveness of the simultaneous processing of synchronously registered cardiosignals (SRCS) based on the model of the vector of cyclic rhythmically connected stochastic processes in comparison with the well-known method of their processing. As the example it is shown on the Figure 1 the realizations of SRCS electrocardiosignal (ECS) and phonocardiosignal (PCS) and on the figures 2-4 are represented the results of a comparative analysis.  (t ), mV  (t ), conditional units 0,3 0,5 0,2 0,4 0,1 0,3 0 0,2 0,1 -0,1 0 -0,2 -0,1 t , s -0,3 t, s 0 2 4 6 8 10 0 2 4 6 8 10 a) b) Figure 1: The realizations of SRCS: (a) electrocardiogram, (b) phonocardiogram. RˆT (t1 , t2 ), mV 2 RˆT (t ,1) (t1 , t2 ), mV 2 0,06 0,06 0 0 -0,06 -0,06 t2 , c t2 , c 0,8 0,8 0,6 0,8 t1 , s 0,6 0,8 t1 , s 0,4 0,6 0,4 0,6 0,4 0,4 0,2 0,2 0,2 0,2 0 0 а) b) Figure 2: The realization graphs of statistical estimations of the autocorrelation function of the ECS while its processing on the basis of: (a) the vector of periodically connected stochastic processes; (b) the vector of cyclic rhythmically connected stochastic processes. RˆT (t1 , t2 ), conditional units 2 RˆT (t,1) (t1 , t2 ), conditional units 2 0,012 0,012 0 0 -0,012 -0,012 t2 , c t 2 , 0,8 c 0,8 0,6 0,8 t1 , s 0,6 t1 , s 0,8 0,4 0,6 0,4 0,6 0,4 0,4 0,2 0,2 0,2 0,2 0 0 a) b) Figure 3: The realization graphs of statistical estimations of the autocorrelation function of PCS while its processing on the basis of: (a) the vector of periodically connected stochastic processes; (b) the vector of cyclic rhythmically connected stochastic processes. Rˆ , (t1 , t 2 ), m V  conditional un its Rˆ , (t1 , t2 ), mV  conditional units 0,09 0,012 0 0 -0,09 -0,012 t2 , c t2 , c 0,8 0,8 0,6 0,8 t1 , s 0,6 0,8 t1 , s 0,4 0,6 0,4 0,6 0,4 0,4 0,2 0,2 0,2 0,2 0 0 а) b) Figure 4: The realization graphs of statistical estimations of the mutual correlation function of the ECS and PCS while its processing on the basis of: (a) the vector of periodically connected stochastic processes; (b) the vector of cyclic rhythmically connected stochastic processes. It was set in the result of the comparative analysis of the statistical processing of sets of SRCS which was made that the method of statistical processing of the analyzed cardiosignals on the basis of the vector of cyclic rhythmically connected stochastic processes significantly reduces the negative effect of "blurring" of statistical estimations of mutual correlation functions as a part of set of synchronous cyclic signals of a heart that is strongly-pronounced in results of statistical processing of the cardiac signals which were investigated on the basis of the vector of periodically connected stochastic processes. It is proved by the fact that the new method of compatible statistical processing takes into account the variability of the SRCS rhythm, in contrast to the well-known methods. 5.1 Software operating principle It was developed the software complex which allows to perform SRCS statistical processing and its structural scheme is shown in Figure 5. In this complex was implemented: block of SRCS formation, SRCS filtering block, trend removal block, SRCS segmentation block, block for determining the SRCS rhythm function, statistical analysis block of the rhythm function, valuation block of statistical estimates of SRCS, statistical analysis block of SRCS, dimension minimization block of the vector of informative features, diagnostic block. The exported synchronously registered cardio signals were used from PhysioNet in accordance with the article [28] in the SRCS implementation generation unit. Pre-processing units are designed to filter and extract trends from SRCS. The segmentation unit is used to obtain the time-time structure of SRCS to determine their rhythmic functions. Blocks of statistical processing are intended for the normalization of statistical estimations, the analysis of investigated cardiac signals. The block of minimization of dimensionality of a vector of informative signs is necessary for minimization of their quantity representing norm or a certain pathology of the cardiovascular system. The diagnostic unit allows you to diagnose the cardiovascular system using SRCS (Fig. 5). BLOCKS OF SRCS PRE-PROCESSING BLOCK OF SRCS SRCS TREND SRCS FORMATION FILTERING REMOVAL SEGMENTATION BLOCK BLOCK BLOCK BLOCKS OF SRCS STATISTICAL PROCESSING BLOCK FOR VALUATION DETERMINING STATISTICAL BLOCK OF STATISTICAL THE SRCS ANALYSIS BLOCK OF STATISTICAL ANALYSIS BLOCK RHYTHM THE RHYTHM ESTIMATES OF OF SRCS FUNCTION FUNCTION SRCS DIMENSION MINIMIZATION BLOCK OF THE DIAGNOSTIC VECTOR OF BLOCK INFORMATIVE FEATURES Figure 5: Software complex for statistical processing of synchronously registered cardiosignals. For the compatible statistical processing of SRCS (electrical, magnetic, acoustic nature) the module “Compatible statistical processing of SRCS” is used, the structural-functional scheme of which is shown in fig. 6. Data entry Estimation of compatible correlation (two synchronously registered function taking into account cardiosignals and their rhythm functions mathematical expectations) Data storage Data visualization Visualization of the obtained results Figure 6: Structural‐functional module scheme of the program "Compatible statistical processing SRCS". The input data for compatible statistical processing are oversampled values of signals and their estimates of mathematical expectation. The result of joint statistical processing is a cross- correlation function obtained taking into account the rhythm function (Fig. 6). To simulate cardiosignals, a simulation unit SRCS was developed taking into account the rhythm function, the structural-functional scheme of which is shown in Fig. 7. In order to simulate cardio signals taking into account the rhythm function, it is first necessary to open the input data of the cardio signal and rhythm function, set the number of cycles and signal zones, load the input cycle (Fig. 7). Data entry Data entry of a new (cardio signal, rhythm rhythm function (zones function) and cycles) Enter quantity Forming a durations table cycles and zones of the rhythm function cardiosignal Cardio signal modeling taking Formation of a new into account the rhythm rhythm function function Saving of modeling results Visualization of results and new rhythm function Figure 7: Structural‐functional scheme of the SRCS simulation unit taking into account the rhythm function. 5.2 Software design In Fig. 8 the general view of program interfaces for data input (cardiosignals, rhythm functions) and their visualization is given. a) b) Figure 8: General view of the program interfaces for data entry (a) and their visualization (b). Two methods were used in the program for statistical processing of SRCS: the method of averaging the values of signal realization taking into account the period and the method of averaging the values of signal realization taking into account the rhythm function. An example of the program interface for statistical processing taking into account the period is shown in Fig. 9. The method of averaging the values of signal realization taking into account the rhythm function, in comparison with known methods, significantly reduces the negative effect of "blurring" of statistical estimates of probabilistic characteristics of SRCS, which allows to increase the accuracy, reliability and informative of complex automated diagnosis of the human cardiovascular system by synchronous cyclic heart signals (Fig. 9). a) b) Figure 9: General view of the program interfaces for statistical processing through period (a) and taking into account the rhythm function (b). During statistical processing, taking into account the rhythm function, it is necessary to specify additional parameters, namely: the sampling step of the signal and rhythm function (required during resampling of signal and rhythm function), as well as the number of heartbeat cycles and the zones number (per cycle) of registered signal implementation. Performing consistent steps: rhythm function interpolation; signal interpolation; rhythm function resampling; signal oversampling we obtain the necessary data for statistical processing, taking into account the rhythm function. In Fig. 9 shows that as a result of statistical processing we can obtain estimates of the following probabilistic characteristics: mathematical expectation, variance and autocorrelation function. All received data can be saved in text files if necessary. The general view of the program interfaces for visualization of the estimation of the mathematical expectation of the phonocardiogram (a) and the integrated rheogram (b), by the method of statistical processing taking into account the rhythm function, is given in Fig. 10. a) b) Figure 10: General view of the program interfaces for estimation visualization of the mathematical expectation of the phonocardiogram (a) and the integrated rheogram (b). The input data for compatible statistical processing are oversampled values of signals and their estimates of mathematical expectation. The result of compatible statistical processing is a mutual correlation function obtained taking into account the rhythm function. An example of the program interface for entering synchronously registered data (a) and their compatible processing (b) is shown in Fig. 11. a) b) Figure 11: General program interface view for entering synchronously registered data (a) and their compatible processing (b). In order to simulate cardiosignals taking into account the rhythm function, it is first necessary to open the input data of the cardiosignal and rhythm function, set the number of cycles and signal zones, load the input cycle. The program interfaces are used for input data (Fig. 12 (a)) and for visualization of the input cycle (Fig. 12 (b)). а) b) Figure 12: General view of the interface of the program of input data input (a) and visualization of the input cycle (b). The developed software provides the ability to form a discrete rhythm function, depending on the number of cycles and zones of the cardiosignal to be modeled. The program interfaces for forming the rhythm function and visualization of the modeled cardio signal are shown, respectively, in Fig. 13 (a, b). а) b) Figure 13: General view of the program interface for forming the rhythm function (a) and the modeled cardiosignal generated by the rhythm function (b). The developed set of programs makes it possible to perform statistical processing of SRCS of different physical nature and obtain estimates of mathematical expectation, dispersion, autocorrelation and mutual correlation functions (taking into account the rhythm function and taking into account the period). The software also allows you to simulate cardio signals. The developed software provides the ability to form a discrete rhythm function, depending on the number of cycles and zones of the cardio signal to be modeled. 6. Conclusions A set of programs has been developed that allows for statistical processing of cardiac signals of different physical origin and to obtain estimates of mathematical expectation, variance, autocorrelation and cross-correlation (taking into account the function of rhythm and period). The software allows you to simulate the implementation of cardio signals of different physical nature, taking into account their constant rhythm (taking into account the period), variable (taking into account the rhythm function) and the common rhythm. The modeling method includes information obtained at the stage of statistical processing of cardiac signals, in particular, morphological features are taken into account on the basis of estimates of mathematical expectation and variance of cardiac cycle cycles. The developed software provides the possibility of forming a discrete rhythm function depending on the number of cycles and zones of the simulated cardio signal. This software can be used as a component of diagnostic systems of the human heart. 7. References [1] Anastasia P. Nesterova, Eugene A. Klimov, Maria Zharkova, Sergey Sozin, Vladimir Sobolev, Natalia V. Ivanikova, Maria Shkrob, Anton Yuryev, Diseases of the circulatory system, in: Disease Pathways, 1st. ed., Elsevier, 2020, pp. 327-390. doi:10.1016/C2018-0-00586-1. [2] Jer-Hwa Hsu, I-Chia Chien, Ching-Heng Lin, Increased risk of ischemic heart disease in patients with bipolar disorder: A population-based study, Journal of Affective Disorders 281 (2021), 721- 726, doi:10.1016/j.jad.2020.11.083. [3] Jun Yang, Wenjing Xiao, Huimin Lu, Ahmed Barnawi, Wireless high-frequency NLOS monitoring system for heart disease combined with hospital and home, Future Generation Computer Systems 110 (2020) 772-780. doi:10.1016/j.future.2019.11.001. [4] M.T. Donofrio, A.J. Moon-Grady, L.K. Hornberger, J.A. Copel, M.S. Sklansky, A. Abuhamad, et al. Diagnosis and treatment of fetal cardiac disease: a scientific statement from the American heart association Circulation 129 (2014) 2183-2242. doi:10.1161/01.cir.0000437597.44550.5d. [5] Margherita Padeletti, Giuseppe Bagliani, Roberto De Ponti, Fabio M. Leonelli, Emanuela T. Locati, Surface Electrocardiogram Recording: Baseline 12-lead and Ambulatory Electrocardiogram Monitoring, Cardiac Electrophysiology Clinics 11 (2019) 189-201, doi:10.1016/j.ccep.2019.01.004. [6] Raymond R. Bond, Dewar D. Finlay, Chris D. Nugent, George Moore, A review of ECG storage formats, International Journal of Medical Informatics 80 (2011) 681-697, doi:10.1016/j.ijmedinf.2011.06.008. [7] Gail Baura, Electrocardiograph design lab, in Medical Device Technologies, 2d. ed., Academic Press, 2020, pp. 545-548. doi:10.1016/B978-0-12-811984-6.00022-2. [8] Abhinav Nair, Alec Vishnevsky, Daniel R. Frisch, Atrial fibrillation or unstable angina? Utilization of a mobile electrocardiographic device to diagnose acute coronary syndrome, Cardiovascular Digital Health Journal 1 (2020) 52-54. doi:10.1016/j.cvdhj.2020.06.001. [9] Jefferson L. Vieira, Hector O. Ventura, Mandeep R. Mehra, Mechanical circulatory support devices in advanced heart failure: 2020 and beyond, Progress in Cardiovascular Diseases 63 (2020) 630-639. doi:10.1016/j.pcad.2020.09.003. [10] Farman Ali, Shaker El-Sappagh, S.M. Riazul Islam, Daehan Kwak, Amjad Ali, Muhammad Imran, Kyung-Sup Kwak, A smart healthcare monitoring system for heart disease prediction based on ensemble deep learning and feature fusion, Information Fusion 63 (2020) 208-222, doi:10.1016/j.inffus.2020.06.008. [11] Arman Kilic, Artificial Intelligence and Machine Learning in Cardiovascular Health Care, The Annals of Thoracic Surgery 109 (2020) 1323-1329. doi: 10.1016/j.athoracsur.2019.09.042. [12] Pan Ma, Shigong Wang, Ji Zhou, Tanshi Li, Xingang Fan, Jin Fan, Siyi Wang, Meteorological rhythms of respiratory and circulatory diseases revealed by Harmonic Analysis, Heliyon 6 (2020). doi:10.1016/j.heliyon.2020.e04034. [13] Soile Tapio, Mark P. Little, Jan Christian Kaiser, Nathalie Impens, Nobuyuki Hamada, Alexandros G. Georgakilas, David Simar, Sisko Salomaa, Ionizing radiation-induced circulatory and metabolic diseases, Environment International 146 (2021). doi:10.1016/j.envint.2020.106235. [14] Martsenyuk V.P., Sverstiuk A.S., Klos-Witkowska A., Horkunenko A.B., Rajba S. Vector of Diagnostic Features in the Form of Decomposition Coefficients of Statistical Estimates Using a Cyclic Random Process Model of Cardiosignal, The 10th IEEE International Conference on Intelligent Data Acquisition and Advanced Computing Systems: Technology and Applications, volume 1 of IDAACS. Metz, France, 2019, pp. 298–303. doi: 10.1109/IDAACS.2019.8924398 [15] Edward B. Panganiban, Arnold C. Paglinawan, Wen Yaw Chung, Gilbert Lance S. Paa, ECG diagnostic support system (EDSS): A deep learning neural network based classification system for detecting ECG abnormal rhythms from a low-powered wearable biosensors, Sensing and Bio- Sensing Research 31 (2021), 100398, https://doi.org/10.1016/j.sbsr.2021.100398. [16] W. A. Gardner, Cyclostationarity: Half a century of research, Signal Processing 86 (2006) 639–697. doi:10.1016/j.sigpro.2005.06.016. [17] W. A. Gardner Exploitation of cyclostationarity for identifying the Volterra kernels of non– linear systems, IEEE Transactions on Information Theory 39 (1993) 535–542. doi:10.1109/18.212283. [18] W. A. Gardner Fraction–of–time probability for time–series that exhibit cyclostationarity, Signal Processing 23 (1991) 273–292. doi:10.1016/0165-1684(91)90005-4. [19] W. A. Gardner, C. M. Spooner Higher–order cyclostationarity, in: International Symposium on Information Theory and Applications, ISITA '90, Honolulu, HI, 1990, pp. 355–358. [20] W.A. Gardner, T.L. Archer Exploitation of cyclostationarity for identifying the Volterra kernels of non-linear systems, IEEE Transactions on Information Theory 39 (1993) 535-542. doi:10.1109/18.212283. [21] W.A. Gardner, C.M. Spooner The cumulant theory of cyclostationary time-series. II. Development and applications, IEEE Transactions on Signal Processing 42 (1994) 3409-3429. doi:10.1109/78.340776. [22] S. Lupenko, N. Lutsyk, Y. Lapusta, Cyclic linear random process as a mathematical model of cyclic signals, Acta mechanica et automatica 9 (2015) 219–224. doi:10.1515/ama-2015-0035. [23] A. S. Sverstyuk, Argumentation and verification of the mathematical model of synchronously registered cardiosignals using the vector of cyclic rhythmically connected stochastic processes, Measuring and Computing Devices in Technological Processes 1 (2009) 143–147. doi:10.1515/ama-2015-0035. [24] I. Lytvynenko, A. Horkunenko, O. Kuchvara, Y. Palaniza, Methods of processing cyclic signals in automated cardiodiagnostic complexes, in: Proceedings of the 1st International Workshop on Information-Communication Technologies and Embedded Systems, Mykolaiv, Ukraine, 2019, pp. 116-127. [25] S.-W. Chen, H.-C. Chen, H.-L. Chan, A real-time QRS detection method based on moving- averaging incorporating with wavelet denoising, Computer Methods and Programs in Biomedicine 82 (2006) 187–195. doi.org/10.1016/j.cmpb.2005.11.012. [26] V. S. Chouhan, S. S. Mehta, N. S. Lingayat, Delineation of QRS-complex, P and T-wave in 12-lead ECG, IJCSNS International Journal of Computer Science and Network Security 8 (2008) 185–190. [27] Bemmel J.H. van, Past and future research goals for computerized ECG processing, in: Comput. ECG Anal.: Towards Stand. Proc. IFIP - IMIA Work Cont., Amsterdam e.a., 1986, pp. 367-381. [28] A.L. Goldberger, L.A.Amaral, L. Glass, J.M. Hausdorff, P.Ch. Ivanov, R.G. Mark, J.E. Mietus, G.B. Moody, C.-K. Peng, H.E. Stanley. PhysioBank, PhysioToolkit, and PhysioNet: Components of a New Research Resource for Complex Physiologic Signals. Circulation 101(23):e215-e220 2000 (June 13). ULR : http://circ.ahajournals.org/content/101/23/e215.full.