=Paper=
{{Paper
|id=Vol-2623/paper4
|storemode=property
|title=Modification of the Software System for the Automated Determination of Morphological and Rhythmic Diagnostic Signs by Electrocardio Signals
|pdfUrl=https://ceur-ws.org/Vol-2623/paper4.pdf
|volume=Vol-2623
|authors=Serhii Lupenko,Iaroslav Lytvynenko,Nataliia Stadnyk,Halyna Osukhivska,Natalia Kryvinska
|dblpUrl=https://dblp.org/rec/conf/intelitsis/LupenkoLSOK20
}}
==Modification of the Software System for the Automated Determination of Morphological and Rhythmic Diagnostic Signs by Electrocardio Signals==
Modification of the Software System for the Automated Determination of Morphological and Rhythmic Diagnostic Signs by Electrocardio Signals Serhii Lupenko1[0000-0002-6559-0721], Iaroslav Lytvynenko1[0000-0001-7311-4103], Nataliia Stadnyk1[0000-0002-7781-7663], Halyna Osukhivska 1[0000-0003-0132-1378], Natalia Kryvinska2[0000-0003-3678-9229] 1 Ternopil Ivan Puluj National Technical University, Ternopil, Ukraine 2 University of Vienna, Vienna, Austria; Comenius University, Bratislava, Slovakia lupenko.san@gmail.com iaroslav.lytvynenko@gmail.com natalya.stadnik15@gmail.com osukhivska@tstu.edu.ua Ukraine natalia.kryvinska@univie.ac.at Abstract. In this paper we consider a modernized of the software system for the automated determination of morphological and rhythmic diagnostic signs of electrocardio signal. The modification of the system is to create a new method and appropriate software that involves processing electrocardio signals by re- ducing the discrete cyclic random process, as a model of electrocardio signal to isomorphic random periodic sequence. The use of a new mathematical model of electrocardio signals in the form of conditional discrete cyclic random process allowed to take into account and carry out automatic determination of both morphological and rhythmic diagnostic sings of electrocardio signal within the same mathematical model. The use of a new method of statistical processing based on the new model, allowed to obtain statistical characteristics that are in- fomative diagnostic signs (morphological and signs of rhythm) of the electro- cardio signal. The application of the method of reducing a discrete cyclic ran- dom process to isomorphic random periodic sequence before the procedures of statistical processing of the electrocardio signal, in particular before obtaining morphological and rhythmic sings allowed to increase the speed of automated processing of the electrocardio signal in comparison with the previously devel- oped methods which were based on it model in the form of a cyclic random process and did not account for the double stochastic model.In the structure of the modified software system, after the evaluation of the rhythmic structure of the electrocardio signal and the procedure of reduction to a random periodic se- quence of a discrete cyclic random process, its processing is branched into two parallel stages. The first stage carries out the morphological analysis, which in- volves the statistical processing of the electrocardio signal, normalization of sta- tistical estimates and their distribution in the Chebyshev base and decision mak- ing on the obtained morphological sings. Copyright © 2020 for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). IntelITSIS-2020 Keywords: Software System, Electrocardio Signal, Analysis Of Heart Rhythm, Morphological Analysis Of Electrocardio Signals. 1 Introduction Known processes and phenomena of reality are those that reflect over a repeating structure in time. Electrocardio signals are known and well-studied among such pro- cesses, signals. The current state of information technology development allows to approach the processing of such signals in a new way and solve the advanced tasks of modern medicine efficiently in the construction of diagnostic apparatus by creating new effective mathematical means which make it possible to increase the accuracy and informativeness of the processing of cyclic signals, in particular electrocardio signals [1-12]. Therefore, the development of diagnostic software systems for auto- mated diagnosing the human heart condition by registered electrocardio signals is a relevant scientific and technical task, the solution of which will enable to improve the quality and efficiency of diagnosing of the functional state of the heart and cardiovas- cular system of the human body as a whole. 2 Related works New mathematical models and methods of processing cyclic electrocardio signals have been developed and substantiated in papers [13, 14]. The use of new mathemati- cal models has permitted to improve the accuracy and reliability of diagnosing the functional state of the heart by increasing the informativeness of morphological anal- ysis and heart rhythm analysis. The software system for modeling and conducting morphological analysis, electrocardio signal analysis of heart rhythm on the basis of the developed mathematical models is presented in paper [15]. This system has been modernized and embodied newly developed methods for the electrocardio signals processing. 3 Overview of the Research Increasing the speed of methods of statistical processing of electrocardio signals, along with increasing the informativeness of their automated analysis is achieved through the use of a new mathematical model that reflects the double stochasticity of the investigated electrocardio signals (stochastic morphological and rhythmic struc- tures) and applying a new method of reducing a discrete cyclic random process to an isomorphic random periodic sequence. This paper is devoted to the improvement of the software system, where in contrast to the previous development [13, 14], the system includes the new method that pro- vides the processing of electrocardio signals by reducing to an isomorphic random periodic sequence of a discrete cyclic random process, as a model of electrocardio signal, which allowed to increase the speed of its processing in comparison with the previously developed methods by reducing the computational complexity of known statistical methods of estimating the probabilistic characteristics of cyclic random processes of a discrete argument. 4 Proposed model The purpose of the research is to modernize the software system for morphological analysis and heart rhythm analysis with increased informativeness, on the account of using a new mathematical model of heart signals, in the form of a conditional discrete circular random process, and a new statistical processing method, namely the method of reducing a discrete cyclic random process to an isomorphic random periodic se- quence, that by reducing the computational complexity of the electrocardiographic signal processing method, made it possible to speed up their processing in computer cardiac diagnostic systems compared to previously developed methods. 5 Results & Discussion Here are presented the basic mathematical relations that underlie mathematical sup- port in the modified software system. Still, we will focus on the developed part of the complex, which relates to the reduction to an isomorphic random periodic sequence of a discrete cyclic random process, as a model of the electrocardio signal. The software system enables the automated analysis of the electrocardio signal, particularly its morphological and rhythmic signs. The mathematical model of electrocardio signal in the form of a cyclic random process and conditional cyclic random process. It is known from [13] that a discrete random process: {ξ (ω , tml ), ω ∈ Ω, tml ∈ D} is called a cyclic discrete random process if there is such a discrete function T (t ml , n ) (rhythm function) that satisfies the conditions: 1) T (t ml , n ) > 0 , if n > 0 ; 2) T (t ml , n ) = 0 , if n = 0 ; 3) T (t ml , n ) < 0 , if n < 0 ; 4) for any t m1l1 ∈ D and t m2l2 ∈ D , for which t m2l2 > t m1l1 , for function T (t ml , n ) inequality ( ) ( ) t m1l1 + T t m1l1 , n < t m2l2 + T t m2l2 , n , ∀n ∈ Z is applied; that finite-dimensional vectors ( ξ (ω , t m1l1 ) , ξ (ω , t m2l2 ) ,..., ξ (ω , t mk lk ) ) and ( ) ( ξ (ω , t m1l1 + T t m1l1 , n ) , ( ) ξ (ω , t m2l2 + T (t m2l2 , n )) ,..., ξ (ω , t mk lk + T t mk lk , n ), n ∈ Z , at all integer k ≥ 1 is sto- chastically equivalent in a broad sense. The method of reducing the statistical processing (estimation, analysis, prediction) of a cyclic random process of a discrete argument to the corresponding statistical processing of an isomorphic periodic random sequence consists in gradual execution of the following steps: 1) transformation of a ω -realization ξ1ω (t ml ), ω ∈ Ω, t ml ∈ D of a cyclic random process ξ1 (ω , t ml ), ω ∈ Ω, t ml ∈ D in a ω -realization ξ 2ω (i ), i ∈ Z of an isomorphic (for the process) in relation to order and values of a L -periodic sequence ξ 2 (ω , i ), ω ∈ Ω, i ∈ Z , by the action of a scale transformation operator G y (t ml ) {} ⋅ with a scale transformation function y (t ml ) = L ⋅ (m − 1) + l ; 2) application of known methods of processing periodic random sequences and ob- taining their results (statistical point and interval estimates of certain probabilistic characteristics); 3) obtaining statistical estimates of the probabilistic characteristics of a cyclic ran- dom process ξ1 (ω , t ml ), ω ∈ Ω, t ml ∈ D , by the application inversion operator of the scale convension to previously obtained appropriate statistical estimations for L - periodic random sequence. We will assume, that is recorded M cycles by L counts in each cycle of the in- vestigated cyclic signal whose mathematical model is a cyclic random process ____ ____ ξ1 (ω , t ml ), ω ∈ Ω, t ml ∈ R, m = 1, M , l = 1, L (for the simplified further marking ξ1 (ω , t ml ) will be written). Relatively, mathematical model of a cyclic signal registro- ____ ____ gram will be a ω -realization ξ1ω (t ml ), t ml ∈ R, m = 1, M , l = 1, L (for the simplifying of further marking ξ1ω (t ml ) will be written) of this cyclic random process of a discrete argument. The isomorphic for the investigated discrete process in relation to order and ________ values of a L -periodic sequence ξ 2 (ω , i ), ω ∈ Ω, i = 1, M ⋅ L (for the simplifying of further marking ξ 2 (ω , i ) will be written) is obtained, by the action of a scale trans- formation operator with a scale transformation function y (t ml ) = L ⋅ (m − 1) + l to initial random process ξ1 (ω , t ml ) , namely: ξ 2 (ω , i ) = G y (t ml ) {ξ1 (ω , t ml )} (1) which is equivalent to a such system of equations: _____ ____ i = y (t ml ) = L ⋅ (m − 1) + l , m = 1, M , l = 1, L, ________ (2) ξ 2 (ω , i ) = ξ1 (ω , t ml ), i = 1, M ⋅ L, t ml ∈ R. The same scale transformation operator is related M -cyclic ω -realization cyclic ________ random process ξ1 (ω , t ml ) and M -cyclic ω -realization ξ 2ω (i ), i = 1, M ⋅ L (for the simplifying of further marking ξ 2ω (i ) will be written). The analytical formula for calculating the value of statistical estimation of the ini- tial moment function of the first order (mathematical expectation) L -periodic se- quence ξ 2 (ω , i ) , which is isomorphic in relation to the order and values of the cyclic process ξ1 (ω , t ml ) , has the form of 1 M −1 ___ mˆ ξ 2 (l ) = ∑ξ 2ω (l + L ⋅ n ), l = 1, L (3) M n=0 The analytical formula for the calculating the value of statistical estimation of cen- tral moment function of the second order (dispersion) L -periodic sequence ξ 2 (ω , i ) , which is isomorphic in relation to the order and values of the cyclic process ξ1 (ω , t ml ) , has the form of ∑( ) 1 M −1 ___ dˆξ 2 (l ) = ξ 2ω (l + L ⋅ n ) − mˆ ξ 2 (l ) 2 , l = 1, L (4) M − 1 n=0 Due to the statistical procedure for the electrocardio signal processing normaliza- tion received statistical estimation and their reduction in Chebyshev basis, that is in- vestigated in paper [13]. A mathematical model of electrocardio signal will be considered below. A model takes into account their double stochasticity, namely, morphological structure stochas- ticity and stochasticity of the rhythmic structures of electrocardio signal. The condi- tional cyclic random process is called a process {ξ (ω , ω ′, t ), ω ∈ Ω, ω ′ ∈ Ω′, t ∈ R} , that is set on the Cartesian product of two stochasticly independent probabilistic spac- es with the sample sets Ω and Ω′ , and on the set of real numbers R , and for which such conditions are satisfied: 1. a such random function exists T (ω ′, t , n ), ω ′ ∈ Ω′, t ∈ R, n ∈ Z , that for each ω ′ , relatively ω ′ -realization Tω ′ (t , n ) of this function, satisfies the conditions of rhythm function; 2. for each ω ′ from Ω′ finite-dimensional vectors ( ξ ω ′ (ω , t1 ) , ξ ω ′ (ω , t 2 ) ,…, ξ ω ′ (ω , t k ) ) and ( ξ ω ′ (ω , t1 + Tω ′ (t1 , n )) , ξ ω ′ (ω , t 2 + Tω ′ (t 2 , n ))) ,..., ξ ω ′ (ω , t k + Tω ′ (t k , n ) ), n ∈ Z , where {t1 , t 2 ,..., t k } - separable set of the process ξ ω ′ (ω , t ), ω ′ ∈ Ω′, ω ∈ Ω, t ∈ R , for all integer k ∈ N is a stochastic equivalent in a broad sense; 3. for any different ω1′ ∈ Ω′ and ω 2′ ∈ Ω′ random processes ξ ω1′ (ω , t ) and ξ ω 2′ (ω , t ) are isomorphic in relation to the order and values of the cyclic random process. A mathematical model of a rhythm cardio signal with increased resolution, accord- ing to the paper [14], is a discrete random process ____ T (ω ′, t ml , n), ω ′ ∈ Ω′, t ml ∈ R, m ∈ Z, l = 1, L, L ≥ 2, n ∈ Z , which is embedded in a random rhythm function T (ω ′, t , n), ω ′ ∈ Ω′, t ∈ R, n ∈ Z of a conditional cyclic ran- dom process {ξ (ω , ω ′, t ), ω ∈ Ω, ω ′ ∈ Ω′, t ∈ R} . The first stage of a heart rhythm analysis on the basis of a rhythm cardio signal with increased resolution is a for- mation of a vector of random stationary and stationary connected sequences ____ ______ Ξ L (ω ′, m) = Tl (ω ′, m), ω ′ ∈ Ω′, l = 1, L, m = 1, M . Then the statistical processing of a vector component is conducted. At the same time a mathematical expectation, dispersion, the type of distribution (checking it for normality) is evaluated, by the building a histogram and the use of Pearson's agreement criterion χ 2 . We present the basic mathematical relations for estimating the probabilistic characteristics of the components of this vector of random sequences. An expression for calculating a realization of a statistical estimation of сˆ1T of the l relative vector component of the first-order initial moment с1T (mathematical expec- l tation) of a stationary random sequence Tl (ω ′, m) , namely: 1 M ____ сˆ1T = ∑ Tlω′ (k ), l ∈ 1, L (5) l M k =1 where M - the number of cycles of registered realization of an electrocardio signal, Tlω′ (k ) - l vector component rhythmic cardio signal. Statistic estimation of autocorrelation function will have the form of: rˆ2T T (u ) = rˆ2T T (m1 − m2 ) = l1 l2 l1 l2 ∑( )( ) M − M1 Tl (k ) − сˆ1T ⋅ Tl 2 ′ (k + u ) − сˆ1T , 1 M − M 1 + 1 k = 0 1ω′ l1 ω l2 __________ _____ ___ u = 0, M 1 − 1, m1 , m2 ∈ 1, M 1 , l1 , l 2 ∈ 1, L (6) where M 1 - number of counts of correlation function, depth of correlation. The decreasing of the number of diagnostic signs in an information system of heart rhythm analysis on the basis of the main vector of a rhythm cardio signal with increased informativeness is achieved by the use of a spectral decompositions of tri- ˆ = rˆ ___ ____ angular matrix elements R T Tl1Tl2 2 (u ), l1 = 1, L , l 2 = l1 , L , in particular, by the us- ing of a discrete Fourier transform for the estimation of the autocorrelation and inter- correlation functions of this matrix. Namely, instead of a triangular matrix ˆ = rˆ ___ ____ R T 2Tl1Tl2 (u ), l1 = 1, L, l 2 = l1 , L of correlation functions, a triangular matrix ___ ____ Sˆ T = Sˆ 2T T (ν ), l1 = 1, L, l 2 = l1 , L could be used, elements of which are the Fourier l1 l2 transform of the image of relative estimations of a correlation function from the ma- trix R̂ T . Namely, the Fourier transform of the image from the matrix Ŝ T are calcu- lated as: M 1 −1 − j 2πuν Sˆ 2T T (ν ) = ∑ rˆ2T T (u ) ⋅ e __________ M1 , ν = 0, M 1 − 1, l1 l 2 l1 l 2 u =0 (7) ___ ____ l1 = 1, L, l 2 = l1 , L, j = − 1 The structural and functional scheme of the modified software system is presented on the fig. 1. The system of programs is implemented in a programming language Object Pascal. EVALUATION OF EVALUATION OF THE SEGMENT THE RHYTHMIC REDUCTION OF A DISCRETE CYCLIC RANDOM STRUCTURE OF STRUCTURE OF PROCESS (ELECTROCARDIOGRAM) TO AN THE THE ISOMORPHIC PERIODIC RANDOM SEQUENCE ELECTROCARDIO- ELECTROCARDIO- SIGNAL SIGNAL MAKING NORMALIZATION OF STATISTICAL DECISIONS BASED STATISTICAL PROCESSING OF ON ESTIMATES AND THEIR SCHEDULES IN ELECTROCARDIO- MORPHOLOGICAL THE CHEBYSHEV BASIS SIGNAL FEATURES DIAGNOSIS RESULT THE BLOCK OF MAKING SPECTRAL ANALYSIS STATIC PROCESSING THE FORMATION OF STATISTICAL DECISIONS BASED OF A VECTOR OF VECTOR`S ESTIMATES OF THE ON RHYTHMIC VECTOR `S RHYTHMOCARDIO- RHYTMOCARDIO- FEATURES RHYTHMOCARDIO- SIGNAL SIGNAL SIGNAL Fig. 1. Structural-functional scheme of the software system for the analysis of morphological and rhythmic diagnostic of the electrocardio signal signs The processing procedure of the investigated electrocardio signal involves the es- timation of the segmental by means of the segmentation methods. The estimation of rhythm function is conducted by interpolation of rhythmic structure (discrete rhythm function). After the evaluation of a rhythmic structure and the procedure of reduction to a random sequence of a discrete cyclic random process of processing branches into two parallel stages. The first stage conducts a morphological analysis which according to the given structure provides a statistical processing of the electro cardio signal, nor- malization of statistical estimations and their reduction in Chebyshev basis and mak- ing decision due to the obtained morphological sings. The second stage conducts the rhythm analysis and consists in the formation of the vector rhythm cardio signal, the statistical vector processing and spectral analysis of the obtained statistical estima- tions [13]. As an example, the Fig. 2 shows the general view of the program interface is given for the statistical processing of the electrocardio signal, which provides the use of the method of reducing the cyclic random process of a discrete argument to isomorphic random sequence. 120 ξ ω (t i ), y.o. 340 Tω ′ (t i ,1) 100 320 80 300 60 280 40 260 20 240 0 220 -20 200 -40 ti ti 180 0 500 1000 1500 2000 2500 0 400 800 1200 1600 2000 2400 a) b) Fig. 2. Results of the processing: а) a few cycles of the examined electro cardio signal; b) rhythmic structure of the electro cardio signal, the red color defines the samples corresponding to R-R - intervals Figures 3-6 show the graphs for the explanation of the stages of the processing of electrocardio signals by the software system. 120 ξ 2ω (i ), y.o. 100 80 60 40 T (t ml ,1) = T 2500 20 0 2000 -20 i ⋅101 1500 -40 t ml 1000 0 400 800 1200 1600 2000 0 500 1000 1500 2000 2500 a) b) Fig. 3. Graphs: a) a few cycles of the realization of the L -periodic random sequence, which is obtained from the electro cardiogram, the action of a scale transformation operator; b) the estimations of rhythm function T (t ml , n ) (at n = 1 ) of a L -periodic random sequence 120 ˆ 2 (i ) m 100 80 dˆ2 (i ) 60 80 40 60 20 40 0 20 i ⋅101 -20 -40 i ⋅101 0 0 400 800 1200 1600 2000 0 400 800 1200 1600 2000 a) b) Fig. 4. Graphs: а) a few estimation cycles m () ˆ ξ 2 i of the initial moment function of the first order of the L - periodic random sequence ξ 2 (ω , i ) ; b) a few estimation cycles dˆξ 2 (i ) of the central moment function of the second order of the L - periodic random sequence ξ 2 (ω , i ) A separate vector component of the rhythm cardio signal and the results of its sta- tistical processing is depicted on the figures 4-6. 86 84 T3ω ′ ( m) 82 80 78 76 74 72 70 68 66 64 62 60 m 0 20 40 60 80 100 120 140 160 180 200 a) b) Fig. 5. Graphs of: а) a realization T3ω ′ (m) of the third vector component of rhythmic cardio signal T3 (ω ′, m) , which describes the duration of the T -intervals in the electro cardio signal; b) histogram for the third component T3 (ω ′, m) , which describes the duration of the T - intervals in the electro cardio signal In this paper, the mathematical support of the software system for increasing the speed of electrocardio signal processing in comparison with previously known meth- ods is substantiated;the structural and functional scheme of the modernized software system is developed; developed a program that implements a new method of reducing a discrete random process to an isomorphic random periodic sequence, which will achieve faster processing of electrocardio signals in computer cardiodiagnostic sys- tems; the modernization of the modernized software system on real electrocardio signals was carried out. rˆ3T T (u ) 20 15 3 3 Sˆ3T T (u ) 3 3 10 15 5 10 -0 -5 5 -10 u u 0 0 10 20 30 40 50 60 70 0 5 10 15 20 25 30 35 a) b) Fig. 6. A realization graph of: а) the rˆ3T T (u ) statistical estimation of the autocorrelation function r3T T (u ) ( l1 = 3, l 2 = 3 ) 3 3 3 3 of the third vector component Ξ3 (ω ′, m) ; b) Sˆ3T T (ν ) statistical estimations of a spectral power densities S 3T T (ν ) ( l1 = 3, l 2 = 3 ) of 3 3 3 3 the third vector component Ξ3 (ω ′, m) 6 Conclusion Modernized software system due to the extension of its mathematical software, which is based on a new approach to the processing of electrocardiograms based on a math- ematical model in the form of a conditional cyclic random process and a meth-od of reducing their mathematical model in the form of a discrete cyclic random process to an isomorphic random sequence morphological analysis and analysis of the rhythm of the cardio signals with increased informativeness, which made it possible to increase the speed of their processing and increase the accuracy and reliability of diagnosis of the cardiovascular system of the human body. The created program system can be used as a component of specialized software in automated diagnostic systems for system morphoanalysis and heart rhythm analysis References 1. Indu Saini, Dilbag Singh, Arun Khosla. (2013) QRS detection using K-Nearest Neighbor algorithm (KNN) and evaluation on standard ECG databases. Journal of Advanced Re- search, Volume 4, Issue 4. July 2013. -pp. 331-344. doi.org/10.1016/j.jare.2012.05.007 2. M.K. Bhaskar, S.S. Mehta, N.S. Lingayat (2013) Probabilistic Neural Network for the Au- tomatic Detection of QRS-complexes in ECG using Slope. International Journal of Emerg- ing Technology and Advance Engineering Volume 3, Issue 6, June 2013. –pp.255-261. ISSN 2250-2459. ISO 9001:2008. 3. M. Rahimpour, M. E. Asl and M. R. Merati (2016) ECG fiducial points extraction using QRS morphology and adaptive windowing for real-time ECG signal analysis, 2016 24th Iranian Conference on Electrical Engineering (ICEE), Shiraz. -pp. 1925-1930. doi: 10.1109/IranianCEE.2016.7585836. 4. Israa Shaker Tawfic, Sema Koc Kayhan. (2017) Improving recovery of ECG signal with deterministic guarantees using split signal for multiple supports of matching pursuit (SS- MSMP) algorithm, Computer Methods and Programs in Biomedicine, vol. 139, 2017. – pp. 39-50. doi.org/10.1016/j.cmpb.2016.10.014 5. Ciucurel C., Georgescu L., Iconaru E.I. (2018), ECG response to submaximal exercise from the perspective of Golden Ratio harmonic rhythm, Biomedical Signal Processing and Control, 40. –pp. 156-162. doi.org/10.1016/j.bspc.2017.09.018. 6. Fumagalli F., Silver A.E., Tan Q., Zaidi N., Ristagno G. (2018), Cardiac rhythm analysis during ongoing cardiopulmonary resuscitation using the Analysis During Compressions with Fast Reconfirmation technology, Heart Rhythm, 15(2).-pp.248-255. doi: 10.1016/j.hrthm.2017.09.003. 7. Napoli N.J., Demas M.W., Mendu S., Stephens C.L., Kennedy K.D, Harrivel A.R, Bailey R.E., Barnes L.E. (2018), Uncertainty in heart rate complexity metrics caused by R-peak perturbations, Computers in Biology and Medicine, 103.-pp. 198-207. doi: 10.1016/j.compbiomed.2018.10.009. 8. Ebadollah Kheirati Roonizi, Reza Sameni. (2013) Morphological modeling of cardiac sig- nals based on signal decomposition, Computers in Biology and Medicine, vol. 43 (10), 2013. –pp. 1453-1461. doi.org/10.1016/j.compbiomed.2013.06.017. 9. Sandeep Raj, Kailash Chandra Ray. (2018) Sparse representation of ECG signals for au- tomated recognition of cardiac arrhythmias, Expert Systems with Applications, vol. 105, 2018. –pp. 49-64. doi.org/10.1016/j.eswa.2018.03.038. 10. Santanu Sahoo, Prativa Biswal, Tejaswini Das, Sukanta Sabut. (2016) De-noising of ECG Signal and QRS Detection Using Hilbert Transform and Adaptive Thresholding, Procedia Technology, vol. 25, 2016. –pp. 68-75. doi.org/10.1016/j.protcy.2016.08.082. 11. Athanasios Christou Micheas. (2018) Theory of Stochastic Objects Probability, Stochastic Processes and Inference. Chapman and Hall/CRC. January 24, 2018. –p. 408. ISBN- 13: 978-1466515208. 12. Scott Miller, Donald Childers. (2012) Probability and Random Processes: With Applica- tions to Signal Processing and Communications. 2nd Edition. Academic Press is an im- print of Elsevier. USA 2012. –p.522. ISBN: 9780123869814. doi.org/10.1016/C2010-0- 67611-5. 13. Lupenko S. (2018) Modeling and signals processing using cyclic random functions/ S. Lu- penko, O. Orobchuk, N. Stadnik, A. Zozulya//13th IEEE International Scientific and Technical Conference on Computer Sciences and Information Technologies (CSIT), Sep- tember 11-14. — Lviv, Ukraine, 2018. —T. 1, pp. 360-363. ISBN 978-1-5386-6463-6. IEEE Catalog Number: CFP18D36-PRT. 14. Lupenko S. (2018) Statistical analysis of human heart with increased informativeness/ S. Lupenko, N. Lutsyk, O. Yasniy, Ł. Sobaszek// Аcta mechanica et automatica, vol.12 no.4.– Białystok, Poland 2018. –pp.311-315. doi.org/10.2478/ama-2018-0047. 15. Lutsyk N.S., Lytvynenko Ia.V., Lupenko S.A., Zozulya A.M. (2016) Software complex for morphological and heart rate analysis with high informativeness. Journal of Vinnitsa Na- tional Technical UniversityInformation Technologies and Computer Engineering. Vinni- tsa, 2016. №1 (35). P. 13–22. ISSN 1999-9941.