<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta>
      <journal-title-group>
        <journal-title>International Journal of Cardiology</journal-title>
      </journal-title-group>
    </journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.1016/j.hrthm.2017.09.003</article-id>
      <title-group>
        <article-title>Simulation of Cyclic Signals (Generalized Approach)</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Serhii Lupenko</string-name>
          <email>lupenko.san@gmail.com</email>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Iaroslav Lytvynenko</string-name>
          <email>iaroslav.lytvynenko@gmail.com</email>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Volodymyr Hotovych</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Ternopil Ivan Puluj National Technical University</institution>
          ,
          <addr-line>Ruska str. 56, Ternopil, 46001</addr-line>
          ,
          <country country="UA">Ukraine</country>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2021</year>
      </pub-date>
      <volume>17</volume>
      <issue>20</issue>
      <fpage>1259</fpage>
      <lpage>1267</lpage>
      <abstract>
        <p>The efficiency of cyclic signal modeling algorithms is significantly determined by their mathematical models. Using a mathematical model in the form of a cyclic random process and taking into account the obtained statistical information in the form of statistical estimates of mathematical expectation and variance, as well as the estimated discrete rhythm function of the studied cyclic signal, a computer modeling method was developed. This method allows modeling a wide range of signals (of different physical nature), which have a repetitive (cyclic) structure. Based on the information about the estimated rhythmic structure (discrete rhythm function) of the known signal, or by setting (simulating) the rhythm of unfolding in time of the cyclic signal, which can be constant (periodic) and variable (cycle lengths are not equal to each other), the method allows modeling signals with different rhythms. Computer modeling, methods of statistical processing, cyclic signal, cyclic random process</p>
      </abstract>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>Computer simulation of various cyclic signals and their processing using digital systems is an
important task. The simulation allows determining the possibilities of known and created methods of
processing cyclic signals at different stages of signal analysis, testing, for example, them using
simulated implementations. In addition, the simulation of cyclic signals allows the training of newly
created decision-making systems (diagnostic or prognostic), as well as methods for recognizing
biomedical images through their training and testing.</p>
      <p>
        Stochastic periodic processes (periodic random processes), such as periodically correlated random
process and periodically distributed random process with different types of distributions, can be
distinguished as mathematical models used for modeling cyclic signals, taking into account the
stochastic approach. They are widely used in the stochastic modeling paradigm for mathematical
models of cyclic signals [
        <xref ref-type="bibr" rid="ref1 ref2 ref3 ref4 ref5 ref6">1-19</xref>
        ]. In [20, 21] the generalization of stochastic periodic processes in the
form of a new mathematical model of cyclic signals is made. It is presented as a cyclic random
process, which, as a special case, includes a periodic random process.
      </p>
      <p>The model of a cyclic random process is substantiated and applied as a mathematical model of
various signals with cyclic spatiotemporal structure in different diagnostics information systems
(medical diagnostics, technical diagnostics, and forecasting), as well as in systems of analysis and
forecasting of mechanical and economic cyclic processes [23, 24].</p>
    </sec>
    <sec id="sec-2">
      <title>2. Purpose of the research</title>
      <p>This work is devoted to the development of a method of computer modeling of cyclic signals of
different physical nature on the basis of information obtained by statistical processing methods on the
probabilistic characteristics of cyclic random processes, based on a mathematical model in the form of
cyclic random process in medicine, technology, and economics.</p>
      <p>2021 Copyright for this paper by its authors.</p>
    </sec>
    <sec id="sec-3">
      <title>3. Proposed method of cyclic signals simulation</title>
      <p>Modeling of cyclic signals allows to carry out training of the created systems of decision-making
(diagnostic or prognostic), methods of recognition of biomedical images by their training and testing.</p>
      <p>Often for the needs of testing and analysis of various algorithms in practice, it is necessary to have
a large base of representative implementations of cyclic signals. Therefore, for these purposes, it is
convenient to have tools (methods and algorithms) for computer modeling, which allow you to
simulate the desired implementation of cyclic signals. Thus in cyclic signals which are modeled, if the
mathematical model allows, those parameters and characteristics to which the investigated method
(algorithm) reacts and which need to be investigated, to carry out an estimation by the created
methods are put. In the formation of real cyclic signals, their implementation should be preprocessed
by an expert who, using his own experience, will determine in advance representative cycles and
certain parameters and characteristics, which will then be taken into account in cyclic signals
modeled, analyzed, and determined by research methods.</p>
      <p>In computer simulation, the input data are the characteristics and parameters of the morphological
and rhythmic nature of cyclic signals. Figure 1 shows a block diagram of the general algorithm of the
method of modeling cyclic signals using two approaches – deterministic and stochastic.</p>
      <sec id="sec-3-1">
        <title>Input of morphological characteristics and parameters</title>
      </sec>
      <sec id="sec-3-2">
        <title>Simulation of white noise with normal distribution law</title>
      </sec>
      <sec id="sec-3-3">
        <title>Enter the implementations of statistical estimates</title>
      </sec>
      <sec id="sec-3-4">
        <title>Enter the аоб implementation of a representative cycle</title>
        <p>  (t), t  W</p>
      </sec>
      <sec id="sec-3-5">
        <title>Enter the amount of cycles</title>
      </sec>
      <sec id="sec-3-6">
        <title>Enter the amount of zones</title>
        <p>C
Z
mˆ (t), t  W1
dˆ (t), t  W1</p>
      </sec>
      <sec id="sec-3-7">
        <title>The cyclic signal simulation</title>
        <p>f1(t), t  W1
fd (t), (t), Stest (t), t  W
Tˆ(t, n), t  W, n  Z</p>
      </sec>
      <sec id="sec-3-8">
        <title>Enter the characteristics and parameters of the rhythm</title>
      </sec>
      <sec id="sec-3-9">
        <title>The rhythmic structure estimation</title>
      </sec>
      <sec id="sec-3-10">
        <title>Enter the discrete</title>
        <p>rhythmic structure
r
o</p>
      </sec>
      <sec id="sec-3-11">
        <title>Formation of a discrete rhythmic structure</title>
        <p>Z
аоб ,n
Z tW
,n __ ,Z
__ ,1
W =
,t ,j
ˆ=,,,it)(1nTСi ˆ=,,,it)(n1TСij
Ti</p>
      </sec>
      <sec id="sec-3-12">
        <title>Enter cycles</title>
        <p>durations
r Formation of
o cycles durations</p>
      </sec>
      <sec id="sec-3-13">
        <title>Determining the duration of cycles</title>
      </sec>
      <sec id="sec-3-14">
        <title>Enter the durations of the zones on the cycle</title>
      </sec>
      <sec id="sec-3-15">
        <title>Enter the ro value of the period</title>
        <p>ˆ
T</p>
      </sec>
      <sec id="sec-3-16">
        <title>Determination of zone durations</title>
        <p>Dˆc = tˆi , i = 1,С</p>
        <p>Tij</p>
        <p> ____ 
Dˆz = tˆi , i = 1, С, j = 1, Z 
 j </p>
        <p>Features of the form of cyclic signals are taken into account in the block of morphological
characteristics and parameters. Here, a representative cyclic signal cycle can be introduced
f1(t), t  W , using a deterministic approach to simulation of a cyclic signal. Implementations of
statistical estimates of mathematical expectation can be or are introduced mˆ (t), t  W and
dispersions dˆ (t), t  W , using a stochastic approach to simulating the cyclic signal. In addition, to
add additional stochasticity to the simulated data, if necessary, you can add white noise with a normal
distribution law, as an additive component of the simulation model. In addition to morphological
characteristics, the modeling takes into account the parameters of segmental structures, namely, the
number of cycles C and zones Z , as well as the characteristics and parameters of the rhythm in the
corresponding block on the block diagram.</p>
        <p>The following data can be entered in the rhythm generation unit (simulation): discrete rhythmic
structure, if known, for cases of identified segmental cyclic structure T (ti , n), ti  W, n  Z, i = 1, C
and the identified segmental zone structure T (ti , n), ti  W, n  Z , i = 1, C, j = 1, Z ; sets of
j j
samples of segmental cyclic structure Dc = ti , i = 1, C or multiple samples of the segmental zone
structure Dz = ti , i = 1, C, j = 1, Z .</p>
        <p>j</p>
        <p>In the case where certain durations of the respective segment segments are known Ti , i = 1, C or
zone segments Ti , i = 1, C, j = 1, Z or an estimate of the value of the period Tˆ (stable rhythm) then
j
the rhythm can be set taking into account these parameters. The input data for modeling the
implementation of the cyclic signal is a continuous rhythm function, which can be evaluated in the
block "estimation of rhythmic structure", taking into account one of the developed methods, or if it is
known in advance, can be entered in the block characteristics and rhythm parameters.</p>
        <p>For the stochastic case of computer simulation of cyclic signal implementations, stochastic,
independent cycles are modeled and sequentially combined:</p>
        <p>In the case of simulation taking into account the segments-zones:</p>
        <p>C
 ( , t) = U i ( , t),   Ω, t  W .</p>
        <p>i=1</p>
        <p>C Z
 ( , t) = UU i ( , t),   Ω, t  W .</p>
        <p>i=1 j=1 j</p>
        <p>Taking into account the relationship between segmental cyclic and segmental band structures of
cyclic signal implementations:</p>
        <p>Z
 i ( , t) = U i ( , t),   Ω, t  W .</p>
        <p>j=1 j
Computer simulation of a sequence of stochastic, equivalent, independent cycles:</p>
        <p>i (, t) = ( yi (t1), gi (, t1)), t1  W1, i = 1, C, t  W ,
where yi (t1) – is a scale transformation function that takes into account rhythmic structures;
gi ( , t1) – is a set of independent basic (representative) cycles:
(1)
(2)
(3)
(4)</p>
        <p>In such modeling, the attributes are the implementation of statistical estimates (mathematical
expectation and variance) or the implementation of a representative cycle of the cyclic signal (in the
case of a deterministic approach to modeling the cyclic signal):
mg (t1) = mˆ (t1 + T (t1, n)), d g (t1) = dˆ (t1 + T (t1, n)), n  Z;</p>
        <p>mg (t1) = f1(t1 + T (t1, n)), n  Z.</p>
        <p>In block diagram (see Fig. 1) the characteristics of rhythm (rhythmic structures) can be obtained by
taking into account the segmental structures of the set either or
Dz = ti , i = 1, C, j = 1, Z .</p>
        <p>j</p>
        <p>In addition, we can take into account the variable rhythm by specifying either the appropriate sets
of cycle durations Ti , i = 1,С  or sets of duration zones Tij , j = 1, Z , i = 1, С  on the cycle. Here
Dc = ti , i = 1, C
we take into account the stable rhythm by entering the value of the period Tˆ .</p>
      </sec>
    </sec>
    <sec id="sec-4">
      <title>4. Experiments and obtained results</title>
      <p>In this part of the study, we consider examples of the results of the application of the method of
computer modeling of cyclic signals, taking into account different rhythmic structures in different
areas (Figures 2 - 5).</p>
      <p>Mathematical expectation and variance estimates were used to simulate the implementation of the
electrocardiographic signal (Figure 2). Estimation of the mathematical expectation of the
electrocardiographic cycle is the input data for modeling – as a representative cycle that characterizes
a particular pathology. In this case it is hypertrophy of the right and left ventricles. The estimation of
the mathematical expectation determines the morphology of the cycle corresponding to the pathology,
and the variation of the values on the cycle when simulating the implementation of the
electrocardiographic signal is determined by the estimation of the variance. The rhythm of unfolding
the values of the implementation of the electrocardiographic signal in time is set by the rhythmic
structure (discrete rhythm function) given in Figure 3, a. The result of simulation the implementation
of the electrocardiographic signal is presented in Figure 3, b.</p>
      <p>t, s
1
0,06
0,05
0,04
0,03
0,02
0,01
0
dˆ (t), mV 2
0,3
0,2
0,1
0
-0,1
-0,2
-0,3
mˆ  (t ), mВ
0
0,2
0,4
0,6
0,8
0
0,2
0,4
0,6
0,8
a) b)
Figure 2: Input data for computer simulation of cyclic signals (electrocardiographic signals),
implementation of estimates of mathematical expectation and variance: a) Assessment of
mathematical expectation electrocardiographic signal (II lead), diagnosis – hypertrophy of the right
and left ventricles; b) Assessment of variance electrocardiographic signal
(5)
(6)
t, s
1</p>
      <p>T (t,1), s
1
dˆ (t), c. u.2</p>
    </sec>
    <sec id="sec-5">
      <title>5. Discussion</title>
      <p>Proposed by the authors model resembles the classical autoregressive model for the analysis of a
noisy signal. However, the new model differs from the autoregression one in that it allows taking into
account both the morphology of the simulated signals (by taking into account statistical estimates:
mathematical expectation, variance) of the cyclic process and information about the rhythm of its
deployment which can be constant (periodic) and variable (differ from each other among particular
cycles).</p>
      <p>Comparison of the mathematical apparatus offered in this work with others needs separate
research. For example, when comparing with mathematical apparatus of wavelets, the result of
modeling using wavelets will depend on the correctness of the choice of "mother wavelet". That is, it
should be close to the shape of the simulated cycle. There is no such shortcoming in the proposed
method of modeling, because the shape (morphology) of the signal cycle is given arbitrarily and is
determined by estimating the mathematical expectation and variance.</p>
    </sec>
    <sec id="sec-6">
      <title>6. Conclusions</title>
      <p>The developed method allows modeling a wide range of signals of different physical nature, which
have a repetitive (cyclic) structure. Based on the information about the estimated rhythmic structure
(that is discrete rhythm function) of the known signal, or by setting (simulating) the rhythm of time
development of the cyclic signal, which can be constant (periodic signal) and variable (cycle lengths
are not equal to each other), the proposed method allows simulation of signals with a different
rhythm.</p>
      <p>This method allows considering both a rhythm of signal development in time and features of
morphology of a signal representative cycle.</p>
      <p>The developed method based on a mathematical model in the form of a cyclic random process
allows computer modeling of cyclic signals that can be used in intelligent information systems in
economics, mechanics, and medicine. Also, this method is especially important at the testing stage of
new systems for processing signals of such nature.</p>
      <p>This paper concerns the method of modeling cyclic signals, the simulation results of which are
sufficiently interpreted in this paper. To confirm the results obtained in this work, the authors have
used mathematical statistics. However, it is not possible to place all the obtained results only in one
work. In further research, it is planned to conduct computer modeling of cyclic signals in the field of
energy consumption and to investigate the absolute and relative errors of modeling, in particular gas
consumption and energy consumption signals.</p>
    </sec>
    <sec id="sec-7">
      <title>7. References</title>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>W. A.</given-names>
            <surname>Gardner</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Napolitano</surname>
          </string-name>
          , L. Paura, Cyclostationarity:
          <article-title>Half a century of research</article-title>
          ,
          <source>Signal Processing</source>
          ,
          <year>2005</year>
          , №
          <volume>86</volume>
          (
          <year>2006</year>
          ), pp.
          <fpage>639</fpage>
          -
          <lpage>697</lpage>
          . https://doi.org/10.1016/j.sigpro.
          <year>2005</year>
          .
          <volume>06</volume>
          .016.
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>H. L.</given-names>
            <surname>Hurd</surname>
          </string-name>
          , Periodically Correlated Random Sequences :
          <source>Spectral Theory and Practice</source>
          , The University of North Carolina at Chapel Hill Hampton University.
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>P.</given-names>
            <surname>Kochel</surname>
          </string-name>
          ,
          <article-title>Periodically stationary Markovian decision models, Elektron</article-title>
          . Informationsverarb. Kybernet,
          <year>1980</year>
          , no.
          <issue>16</issue>
          , pp.
          <fpage>553</fpage>
          -
          <lpage>567</lpage>
          (in German).
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>N.</given-names>
            <surname>Singh</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J. M.</given-names>
            <surname>Kegan</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J. C.</given-names>
            <surname>Wilcox</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D.</given-names>
            <surname>Hadley</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D.</given-names>
            <surname>Plews</surname>
          </string-name>
          , and
          <string-name>
            <given-names>V.</given-names>
            <surname>Froelocher</surname>
          </string-name>
          , “
          <article-title>Heart Rate Variability: An Old Metric with New Meaning in the Era of using mHealth Technologies for Health and Exercise Training Guidance</article-title>
          .
          <source>Part One: Physiology and Methods”</source>
          ,
          <source>Arrhythmia &amp; electrophysiology Review</source>
          , vol.
          <volume>7</volume>
          ,
          <issue>2018</issue>
          , pp
          <fpage>193</fpage>
          -
          <lpage>198</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>C.</given-names>
            <surname>Hoang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>N.</given-names>
            <surname>Kien</surname>
          </string-name>
          , N. Dung, “
          <article-title>A Review of Heart Rate Variability and its Applications”</article-title>
          ,
          <source>APCBEE Procedia</source>
          , vol.
          <volume>7</volume>
          ,
          <issue>2013</issue>
          , pp
          <fpage>80</fpage>
          -
          <lpage>85</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>R.</given-names>
            <surname>Sassi</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S.</given-names>
            <surname>Cerutti</surname>
          </string-name>
          ,
          <string-name>
            <given-names>F.</given-names>
            <surname>Lombardi</surname>
          </string-name>
          , and etc., “
          <article-title>Advances in heart rate variability signal analysis: joint position statement by the e-Cardiology ESC Working Group and the European Heart</article-title>
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>