<!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 />
    <article-meta>
      <title-group>
        <article-title>technology of neurological movements considering cognitive feedback influences of cerebral cortex signals</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Mykhaylo Petryk</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Mykhaylo Bachynskyi</string-name>
          <email>m.bachynskyi@gmail.com</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Vitaly Brevus</string-name>
          <email>v_brevus@tntu.edu.ua</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Ivan Mudryk</string-name>
          <email>i1mudryk@ukr.net</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Dmytro</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>International Workshop on Information Technologies: Theoretical and Applied Problems</institution>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>M. (ESPCI Paris Tech, ICEM CNRS)</institution>
          ,
          <addr-line>Wang J.-S., Louis E., Haubenberger D., Kalowitz D. and others</addr-line>
        </aff>
        <aff id="aff2">
          <label>2</label>
          <institution>Ternopil Ivan Puluj National Technical University</institution>
          ,
          <addr-line>56 Ruska str., Ternopil 46001</addr-line>
          ,
          <country country="UA">Ukraine</country>
        </aff>
      </contrib-group>
      <abstract>
        <p>A high-performance information technology for evaluating human neurological movements has been developed based on a hybrid model of wave signal analysis, considering the cognitive feedback effects of cerebral cortex neuro-nodes. With the use of hybrid Fourier transforms, a high-speed analytical solution of the model in vector form was implemented, which allows determining the elements of movements on each segment of a complex spiral trajectory performed by the patient with an electronic pen on a digital tablet, and identified the parameters of the studied neuro-systems with feedback.</p>
      </abstract>
      <kwd-group>
        <kwd>Neurological movements</kwd>
        <kwd>cognitive neuro-feedback signals</kwd>
        <kwd>abnormal tremor</kwd>
        <kwd>numerical</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>New modeling methods are used to provide an approach to the design of digital diagnostic health
systems for patients with neurological diseases. The creation of new software and hardware solutions
for medicine and automated diagnostic systems to identify new phenomena of the body and human
health is an urgent task.</p>
      <p>The latest information technologies and modeling methods in the design of a computer diagnostic
system improve the solution to the problem of treating critical diseases in the world, especially people
affected by neurological diseases, such as abnormal neurological movements (ANM) or tremors and
their extreme forms in the form of Alzheimer's, Parkinson's diseases [1]. ANM - unwanted oscillating
movements of certain part of the body (hands, organs of speech, eyeballs), arising as a result of
involuntary contraction of human muscles [2]. Movement regulation disorders signs of human
movements are an increase in their amplitude, a change in the frequency and form of oscillations. The
analysis of these ANM parameters is crucial for understanding the role of feedback dysfunction in
cerebral cortex (CC) neural nodes in cognitive control processes of human movements and early
detection of neuromotor disorders. The difficulty of identifying ANM lies in the imperfection of
existing diagnosis methods, their low accuracy, and the lack of mathematical and software tools for
identifying the neurofeedback influences of CC nodes on their behavior [2].</p>
      <p>Studies of neuro-systems related to the analysis of the behavior of patients with tremor symptoms
(T-objects) were conducted by a number of researchers, such as Pullman S. L., Legrand A.-P., Vidailhet
1ITTAP’2022: 2nd</p>
      <p>2022 Copyright for this paper by its authors.
abnormal states with complex, hard-to-predict behavior, which is inherent in real T-objects with a high
degree of tremor. Due to the imperfection of such methods, there is loss of 60-80% of important
information from the description of patients' real conditions, which de facto determines the low quality
level of such analysis.</p>
      <p>The authors proposed a high-performance information technology for ANM research, built on the
basis of a hybrid model analysis of neuro-system's signals, which describes the state and behavior of
the 3D elements of ANM trajectories of the T-object, considering the cognitive neuro-feedback effects
of the identified CC nodes. Using the methods of hybrid Fourier transform, high-speed analytical
solutions of the model were built in the form of vector functions that determine the elements of the
trajectories on each ANM segment. On their basis, high-performance algorithms for the identification
of tremor-movements parameters are proposed for the component-wise evaluation of neuro-feedback
effects, which allow parallelization of calculations.</p>
    </sec>
    <sec id="sec-2">
      <title>2. A brief description of the hardware used in the research</title>
    </sec>
    <sec id="sec-3">
      <title>2.1.The method of collecting movement data using a graphic tablet</title>
      <p>In addition to qualitative characterization, quantitative improvement of this test requires knowledge
of the position and pressure as a function of the time the patient uses the handle throughout the
experiment. This requires the use of a graphics tablet with specially adapted software. Basically, this
common tool is usually used for artistic drawing, which is provided by specialized software. The goal
here is not the same as above, we need to collect the X and Y coordinates as a function of time, the
frequency and accuracy of which depends on the manufacturer.</p>
      <p>For these reasons, the WACOM Bamboo Fun Medium graphics tablet was chosen. Its active
area (corresponding to 217mm x 137mm) is compatible with the generally accepted patterns of
FahnTolosa-Marin Tremor Rating Scale (FTRS). The reader has a pen input resolution of 2540 dpi, an
accuracy of 0.25 mm and a recognition speed of 133 points per second (according to the manufacturer's
specifications). In addition, this tablet allows you to display pen pressure while drawing on a plane and
measure movements at a distance of up to 16 mm above the surface of the pen, which allows you to
visualize the patient's movements in space.</p>
    </sec>
    <sec id="sec-4">
      <title>2.2. NeuroCom computer electroencephalography system produced by KHAI</title>
    </sec>
    <sec id="sec-5">
      <title>MEDYKA</title>
      <p>The NEUROKOM computer electroencephalograph, the fifth generation of developed computer
electroencephalography complexes [6], was chosen for the study of electroencephalography (EEG)
signals of the brain. The encephalograph in the complex is designed for registration, in-depth analysis
and interpretation of EEG and evoked potentials, conducting various analyzes for scientific research.</p>
      <p>The helmet with the installed hardware and software platform of the manufacturer provides
16channel selection of encephalograms and their transmission to a personal computer using the
appropriate protocol. Software conditioning of EEG signals and post-processing takes place on a PC.
Data is stored both in raw text and in a visualized representation at each point in time. Data is read from
each lead channel at 2ms intervals (500 Hz frequency).</p>
      <p>The hybrid model analysis of ANM built on the basis of the concept of propagation of a wave signal
determines the segment-by-segment description of the elements of ANM trajectories (hands of the
patient) considering the matrix of cognitive effects of groups of neuro-nodes of the CC on the movement
segments [5]. The implementation is based on the method of determining the position of the patient's
hand with an electronic pen, which reproduces the trajectory of the template (Archimedes' spiral) on
the screen of the interactive tablet [7, 8]. The deviation of the trajectory of the pen movement from the
template has a complex shape (Figure 4) and provides digitized information for determining the patient's
neurological condition. The pen movement trace is broken down into simpler elements for the purpose
of decomposing complex ANM-movements within the schematization of the model. The number of
divisions depending on the complexity of the ANM image can be chosen arbitrarily. The hybrid model
provides quantitative amplitude and frequency characteristics of ANM.</p>
      <p>For considering the cognitive neurofeedback effects of the system, the obtained digital sets of
indicators of EEG signals are used, which synchronously with the movement of the electronic pen in
the patient's limb, come from a certain set of CC neuro-nodes. EEG signals in general determine the
dynamics of ANM for each  -th segment of the trace,  = 1,  1 + 1 where  1 is the number of points
of breakdown of the ANM track (Figure 4). In the model, the partition can be set automatically in an
arbitrary way, with any finite number of segments. Their lengths can be different depending on the level
of detail of the traffic sections.</p>
      <p>Model description. In order to present a mathematical solution in the form of implementing the
procedure for the functional identification of the amplitude components and the phase velocity of the
ANM wave propagation  2,  = 1,  1 + 1 as functions of time in the framework of the model
decomposition, taking into account the conditions that the traces (observation data in the form of digital
data of ANM-movements of the patient) of the solution for each  -th segment are known,  =
1,  1 + 1, a system of initial-boundary value problems (micromodels) is obtained for successive
segments of ANM [6, 9]:
 2  2
  2   ( ,  ) =  2   2   +   ∗( ,  ),</p>
      <p>with initial conditions:
  ( ,  )| =0 = 0, | = 0,  = 1,  1 + 1,</p>
      <p>=0
Boundary conditions on each of the ANM segments along z:</p>
      <p>
        −1( ,  )| =  −1 =     −1 ,   ( ,  )| =  =    ,  = 1,  1 + 1, (
        <xref ref-type="bibr" rid="ref1">3</xref>
        )
Selection of the residual functional. We assume that the components of the phase velocity of the
ANM wave propagation  ,  = 1,  1 + 1 of the boundary value problem (1) - (
        <xref ref-type="bibr" rid="ref1">3</xref>
        ) are unknown
functions of time. With known values of the position of the pen   ( ,  ) at observation points on
ANMsegments  ⊂   ,  = 1,  1 + 1
and
  ( ,  )|  =    ( ,  )| , (4)
 
problem (1) - (4) can be considered for each point z for each  1-th segment of the ANM-trajectory
will consist in finding the functions   ∈  , where
 = { ( ,  ):  |  1 ∈  (  1 ) ,  &gt; 0,  = 1,  1 + 1}.
(1)
(2)
Functional-incoherence  1 ∈   1, according to [6, 10] will be written in the form
      </p>
      <p>1  (5)
  (   ) = ∫ (‖  ( ,  ,   ) −   ∗‖2 )  ,</p>
      <p>2 0</p>
      <p>The construction and mathematical substantiation of the solution of the hybrid model is included in
[9].</p>
    </sec>
    <sec id="sec-6">
      <title>4. Modeling and identification of parameters of neurological movements signals of a person under the influence of cognitive neuro-feedback effects of cerebral cortex neuro-nodes</title>
      <p>To configure the model of identification, fragments of the ANM-trajectory of the spiral type, which
was performed by the patient on a digital tablet (Figure 1) were used.</p>
      <p>Modeling and identification of the parameters of ANM movements was carried out within the
framework of the task of identifying parameters of cognitive neuro-feedback effects of EEG on ANM
trajectories developed using the ANM hybrid model considering the feedback effects of EEG-signals.
To set up the identification model, we used a fragment of the ANM-trace, made by the patient using an
electronic pen on an interactive digital tablet according to Figure 5.</p>
      <p>The corrected fragment of the trace of this spiral example of the test pattern (Archimedes spiral) by
the patient with an electronic pen on a digital tablet in the number of discretized 2400 points - positions
of the track is presented in Figure 6. Here, the abscissa k is the number of positions of the deviation of
the pen from equilibrium during the passage of the spiral sample.</p>
      <p>The analysis used a test set of EEG signals of the Fp1- tap sensors of NeuroCom (according to the
CC node classification map), which was measured synchronously with the movement of the electronic
pen when the patient traced the test pattern (Figure 7).</p>
      <p>Setting up the model ANM track and their step-by-step and segment-by-segment identification
(amplitude and frequency parameters for each segment taking into account the integrity of the system)
to a specific sample of the trace performed by the patient (observation curve or experimental curve)
was performed according to the feedback scheme and the analytical solution of the hybrid ANM model
tracks (Figure 8).</p>
      <p>The degree of tremor in this part of the recording ETG (essential tremors graphic), which is based
on the matrix of interactions of brain signals (EEG)  .
Математичний опис елементів траєкторії АНТР на множині In1 у:
 lk ln+1 
 u1 (.t.., z )   l0 11(t.−.. , z, ) ...... ln 1,n1+1(..t. − , z, )   S1 (..., ) </p>
      <p>t  lk ln+1
 u j (..t., z )  = 0  l0 j1(t.−.. , z, ) ...... ln ... ...
un+1 (t, z ) lk n+1,1(t − , z, ) ... ln+1 n1+11,n1+1(t − , z, ) Sn ( , )
l0 ln 
j,n+1(t − , z, )    S ( , )  d d
  i 
  
ETG</p>
      <p>EEG
feedback</p>
      <p>The results of simulation of ANM of the limb with an electronic pen of the patient, their verification
according МtoихthайeлdоaПtЕaТРoИfКthe n(ТeНuТrУo)-experimкаeфnедtрaаcпcрoоrгdрiаnмgноtoїінFжigенuеrрeії4, w12h1i–chSofwtwearreeEancgicneoemrinpganied8b/3y65the
synchronous cognitive neuro-feedback effect of the Fp1 neuro-node (a set of EEG signals according to
Figure 7) and the use of a hybrid model of ANM analysis (based on the hybrid Fourier transform) are
presented in Figure 9-Figure 13. At first, we took a relatively small number of points in an effort to
reproduce the profile of the observation curve (the profile of the ANM-track performed by the patient
taking into account the picture of the EEG feedback curve (Figure 7).</p>
      <p>As can be seen from Figure 9, the accuracy of matching the model track and the patient's real trace
is very high (up to 1.5-2%) for 60 observation points. Amplitude and frequency characteristics due to
the hybrid spectral function built by us, built systematically for all segments of the breakdown (taking
into account their connectivity, not each separately), made it possible to obtain an almost complete
coincidence of the model track with the real patient trace. Then we gradually increased the number of
points on the track. For the number of points 200, the results turned out to be practically the same
(Figure 10).</p>
      <p>In the future, we again gradually increased the number of points to 600, 1200, 2400 and 4 and studied
the behavior of the model curves, evaluating their possible deviations from the patient's experimental
ANM-traces.</p>
      <p>In Figure 11, for a segment with 600 points, we observe a slight deviation in the saddle zone around
the 60th point of the track, about 3-5%. However, this problem can be technically solved by making the
track segmentation in this area smaller. By the way, the model itself allows for arbitrary partitioning
with arbitrary sizes of each segment and making them as small as necessary.</p>
      <p>A positive point is that on all graphs we see a complete reproduction of the frequency characteristics
of the track (the periodicity of the model curves almost completely corresponds to the periodicity of the
curves performed by the patient). As the number of examined points in individual saddle or ridge points
increases, it decreases. But this can be corrected by choosing a smaller division in these zones.</p>
      <p>As can be seen from the presented graphs, the developed model reproduces the patient's behavior at
a high level, displaying an ANM track that practically coincides with the one drawn by him on the
tablet. The most important thing is that the model includes the possibility of displaying the mechanisms
of its feedback effects of CC signal in the form of a matrix of EEG signals that determine the behavior
of these movements. Further research may include a change in this behavior, apparently for the better,
depending on the change in the magnitude of these EEG feedback effects after certain therapeutic
procedures and expand the limits of application.</p>
    </sec>
    <sec id="sec-7">
      <title>5. Conclusions</title>
      <p>Information technology for evaluating the neurological movements of a person is proposed, the one
is based on a hybrid model of wave signal analysis, considering the reverse cognitive effects of cerebral
cortex neuro-nodes. Using hybrid Fourier transform, a high-speed analytical solution of the model in
vector form was obtained and implemented. It determines the elements of trajectories on each segment
of a complex spiral-drawing performed by the patient with an electronic pen on a digital tablet with the
identification of the parameters of the researched feedback systems.</p>
      <p>The proposed hybrid model provides deep decomposition of the system without affecting its
integrity and connections, which is not possible for by classical methods of signal processing, that lead
to the loss of 60-80% of information about the real condition of the object. Weighting coefficients which
characterize the influence of the digital recordings sets of the cognitive influence signals produced by
neuro-nodes of the patient's cerebral cortex are pre-specified by machine learning methods during of
the patient-executed test examples of drawing spiral trajectories on a Wacom digital tablet.</p>
      <p>Compared to the paper method of drawing a spiral, the following indicators were obtained:
•reduction of test duration by 2.5-3 times (about 1 minute instead of 2.5-3 minutes of manual input);
•efficiency and accuracy of the assessment - a significant increase in productivity due to the high
speed of data analysis, the quality of the process of obtaining input data of the 3D drawing of the
spiral, the unambiguous interpretation of the results.</p>
      <p>This approach makes it possible to more qualitatively describe the complex mechanisms of human
neurological conditions under the influence of cognitive connections of the nervous system, which
determine abnormal behavior deviations caused by various man-made and other factors providing a
high degree of data completeness, clarifying diagnoses for treatment.
6. References
[1] A.N. Khimich, M.R. Petryk, D.N. Mykhalyk, I.V. Boyko, A.V. Popov, V.A. Sydoruk, Methods
for mathematical modeling and identification of complex processes and systems based on
visoproductive computing (neuro- and nanoporous cyber-physical systems with feedback,
models with sparse structure data, parallel computing). Monograph, Kiev: National Academy
of Sciences of Ukraine. Glushkov Institute of Cybernetics. 2019. - 176 p. ISBN:
978-966-029188-1.
[2] R. Bhidayasiri, Z. Mari, Digital phenotyping in Parkinson's disease: Empowering neurologists
for measurement-based care. Parkinsonism Relat Disord. 2020 Nov; 80:35-40. DOI:
10.1016/j.parkreldis.2020.08.038.
[3] D. Haubenberger, D. Kalowitz, F.B. Nahab, C. Toro, D. Ippolito, D.A. Luckenbaugh, L.</p>
      <p>Wittevrongel, M. Hallett, Validation of Digital Spiral Analysis as Outcome Parameter for
Clinical Trials in Essential Tremor. Movement Disorders 26 (11), 2073-2080, (2011)
[4] V. Rajaraman, D. Jack, S.V. Adamovich, W. Hening, J. Sage, H. Poizner, A novel quantitative
method for 3D measurement of Parkinsonian tremor. Clinical neurophysiology, 11(2), 187-369
(2000)
[5] H. Xie, Z. Wang, Mean frequency derived via Huang-Hilbert transform with application to
fatigue EMG signal analysis. Comput Meth Progr Biomed 2006; 82:114–20
[6] Electroencephalography complex NEUROKOM, NEUROLAB. Instructions for medical
application AINC.941311.001 I1 U 33.1-02066769-001-2002
[7] A.P. Legrand, I. Rivals, A. Richard, E. Apartis, E. Roze, M. Vidailhet, S. Meunier, E. Hainque,
New insight in spiral drawing analysis methods – Application to action tremor quantification.</p>
      <p>J Clinical Neurophysiology, 128 (10), 1823–1834. (2017)
[8] J.-S. Wang, F.-C Chuang. An Accelerometer-Based Digital Pen with a Trajectory Recognition
Algorithm for Handwritten Digit and Gesture Recognition. IEEE Transactions on Industrial
Electronics, 59(7), 2998-3007 (2012)
[9] E.D. Louis, A. Gillman, S. Böschung, C.W. Hess, Q. Yu, &amp; S.L. Pullman, High width
variability during spiral drawing: Further evidence of cerebellar dysfunction in essential tremor.</p>
      <p>Cerebellum, 11, 872-879 (2012).
[10]I. Mudryk, M. Petryk, Hybrid artificial intelligence systems for complex neural network
analysis of abnormal neurological movements with multiple cognitive signal nodes. 2020 IEEE
Third International Conference on Data Stream Mining &amp; Processing (DSMP) : Conference,
Lviv, 21-25 August 2020. P. 108–111. DOI: 10.1109/DSMP47368.2020.9204341
[11]P. Viviani, P.R. Burkhard, S.C. Chiuvé, C.C. dell’Acqua, P. Vindras, Velocity control in
Parkinson’s disease: a quantitative analysis of isochrony in scribbling movements. Exp Brain
2009; 194:259–83.
[12]A. Salarian, H. Russmann, C. Wider, P.R. Burkhard, F.J. Vingerhoets, and K. Aminian,
“Quantification of tremor and bradykinesia in Parkinson's disease using a novel ambulatory
monitoring system,” Biomedical Engineering, IEEE Transactions on, 54. Jg., Nr. 2, pp.
313322, 2007.
[13]G. LO, A.R. Suresh, L. Stocco, S. González-Valenzuela, and V. C. Leung, “A wireless sensor
system for motion analysis of Parkinson's disease patients,” (PERCOM Workshops), 2011
IEEE International Conference on. IEEE, pp. 372-375, 2011
[14]I. Mudryk, M. Petryk, High-performance modeling, identification and analysis of
heterogeneous abnormal neurological movement’s parameters based on cognitive neuro
feedback-influences. Innovative Solutions in Modern Science. 2021. New York. TK Meganom
LLC. V1(45). P. 235-249 doi: 10.26886/2414-634X.1(45)2021.16.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          <article-title>3. A hybrid mathematical model of ANM-analysis considering the cognitive neuro-feedback-influences of СС nodes</article-title>
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>