=Paper= {{Paper |id=Vol-2864/paper17 |storemode=property |title=Software for Statistical Processing and Modeling of a Set of Synchronously Registered Cardio Signals of Different Physical Nature |pdfUrl=https://ceur-ws.org/Vol-2864/paper17.pdf |volume=Vol-2864 |authors=Serhii Lupenko,Iaroslav Lytvynenko,Andrii Sverstiuk,Borys Shelestovskyi,Andrii Horkunenko |dblpUrl=https://dblp.org/rec/conf/cmis/LupenkoLSSH21 }} ==Software for Statistical Processing and Modeling of a Set of Synchronously Registered Cardio Signals of Different Physical Nature== https://ceur-ws.org/Vol-2864/paper17.pdf
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.