The analysis of these physiological signals on the plane of complex frequencies with use of the Proni procedure Текст научной статьи по специальности «Медицинские технологии»

Инфокоммуникационные технологии и системы

DOI: 10.14529/ctcr180108

THE ANALYSIS OF THESE PHYSIOLOGICAL SIGNALS ON THE PLANE OF COMPLEX FREQUENCIES WITH USE OF THE PRONI PROCEDURE

A.N. Ragozin, [email protected],

V.F. Telezhkin, [email protected]

South Ural State University, Chelyabinsk, Russian Federation

The method of the spectral analysis on the plane of complex frequencies with use of the Proni procedure and also a graphic form of its representation is offered. In physiological studies and in assessing the functional state is very important to data analysis of physiological signals as oscillatory processes. In this widely accepted model of the oscillatory process is the representation of the process under study through the superposition (sum) of sinusoidal periodic processes, each of which is constant in time amplitude, frequency and initial phase. The calculated set of raw signal parameters harmonic sets the range of the investigated process. It should be noted that physiological signals in the General case, represent highly organized in time oscillations. Therefore, in the literature, physiological processes are mostly treated as random processes. In this case, as a measure of the distribution of the frequency of oscillatory activity of FS use a statistical evaluation of the power spectral density (PSD) of a signal, which reflects the dependence of the distribution of the average power (oscillatory activity) of the signal frequency

Keywords: time row, Proni method, complex frequencies, spectral density. physiological signals, oscillatory processes, random processes, oscillatory activity.

Introduction

In physiological researches and diagnostics of functional states the important place is taken by the analysis of these physiological signals (PS) as oscillatory processes. At the same time widespread model of oscillatory process is representation of the studied process through superposition (sum) of sinusoidal periodic processes, each of which is characterized by constants in time amplitude, frequency and an initial phase. A set of the parameters of harmonicas (sinusoidal periodic processes) calculated from an initial signal sets a range of the studied process (differently, distribution of amplitudes and phases of harmonicas depending on frequency). It should be noted that physiological signals generally represent the oscillatory processes which are difficult organized in time. Therefore in literature physiological processes generally are considered as casual processes. In this case as a distribution measure for the frequency of oscillatory activity of PS use statistical assessment of the spectral density of power (SPM) of a signal which reflects dependence of distribution on average of the power (oscillatory activity) of a signal from frequency.

1. Framework for the analysis of SPM

As a basic method for calculation of SPM of the studied signal the periodogrammny method of Welch, with the procedure of the fast transformation of Fourier (FTF) which is its cornerstone usually is used [1, 2]. For improvement of quality of the counted dependence of SPM of the studied process on frequency use also more difficult methods of calculation of SPM with use of procedures of parametrical modeling [1]. Initially concept SPM is entered as the characteristic of spectral structure of stationary casual process, that is process on average uniform in time (the invariance of dispersion, an average, etc.).

At the same time the dependence of SPM on frequency calculated for physiological non-stationary process on some final interval of time of process needs to be understood as an average spectral charac-

teristic of this process on this interval of time. It is necessary to refer various transition processes arising in a human body to non-stationary processes (for example, transition processes of PS) [3, 4]. Temporary changes of the parameters causing not stationarity of process aren't reflected dependence of SPM.

Possible method of the analysis of a spectral characteristic of the studied non-stationary process is calculation of set of the dependences of SPM counted on short time intervals which are consistently displaced in time. However reduction of an interval of time for calculation of dependence of SPM leads to its washing out on frequency (to deterioration in resolution on frequency) and according to its quality [1]. For the analysis of a spectral characteristic of both stationary, and non-stationary FS the method of the spectral analysis on the plane of complex frequencies is offered. This method is realized in the computer program "Spectral Analysis of Physiological Signals" and underwent approbation at the solution of diagnostic and research tasks [4-7]. Also the possibilities of the analysis of SPM PS are considered by parametrical methods [6].

2. Mathematical Model

We will define a concept of complex frequency. Sinusoidal periodic harmonious process is expressed by the known formula:

X (t) = A sin (2%f0t + 9) = Im [ A cos (2nf0t + 0) + jAsin (2%f0t + 9)] =

= Im

Ae< 2f+0)'

= Im

X (t )

(1)

where A - amplitude;f - cyclic frequency in hertz; 9 - an initial phase in radians; X(t) - instant value of

process at the time of t; Im - an imaginary part of complex number; j - imaginary unit; X - complex form of record process. Cyclic frequency off (number of cycles in unit of time) (1) is expressed by a formula [8]:

fo = Re

(- jl 2я/dX (t )/it 1 X (t )

(2)

Determine also circular frequency of rao connected with the cyclic frequency of f rao _ 2%fo ratio.

At substitution of expression for X (t ) = Ael(2n /oi+9) from (1) in (2) we will receive identity. We will consider non-stationary process of a look:

X (t ) = Aeaot sin ( 2ft + 0) = Im

Aea'otej( 2nfof+0)] = Im Г AeJ [( 2f - jao)t+0]

(3)

where A (t) = Aea°l - amplitude of sinusoidal process changing in time under the exponential law; a0 -

attenuation coefficient (at a0 > 0 amplitude of A(t) increases in time, a0 < 0 amplitude of A(t) fades, at a0 = 0 process (3) turns into harmonious sinusoidal process with an invariable amplitude). Comparing expression (3) with (1) it is possible to determine complex circular frequency

= 2f - jao = ®o - jao.

(4)

It is possible to give strict definition of complex circular frequency of w according to expression (2), having omitted operation (Re) of capture of a material part:

та =

j

dX (t )/it X (t )

(5)

Substituting in (5) complex form of record of process (3)

X (t ) = Aej(a0" ja0)t+e we will receive:

w = ~j(j[®0 " ja0 ]) = ®0 " jao, at the same time

/0 =(1/2 *) Re w , a0 = - Im w .

(6)

(7)

(8)

Рагозин А.Н., Тележкин В.Ф.

Анализ данных физиологических сигналов на плоскости комплексных частот с использованием процедуры Прони

Thus, the complex circular frequency of to bears information as about the cyclic frequency of non-stationary process off0, and about coefficient of attenuation of a0, at change of amplitude of process under the exponential law. It should be noted that introduction of complex estimates increases their informational content [9]. For harmonious process of to = 2f = ra0, owing to its stationarity the complex frequency of to coincides with the usual frequency of ra0. On the contrary, if process gets change of amplitude in time under the exponential law, then its frequency becomes complex size.

The reflecting non-stationary physiological process it is correct to present a signal through superposition of elementary sinusoidal processes, each of which is characterized by constants in time f0 frequency, the initial phase 0, initial amplitude And yes coefficient of attenuation of a0. In this case the non-stationary physiological signal decays on sinusoidal components with amplitudes changing in time under the exponential law and is presented in the range form on the plane of complex frequencies. At the same time the range on the plane of complex frequencies is generalization of a usual range and can serve for the analysis of the non-stationary processes reflecting including transitional states. As the basic procedure for calculation of such range the Proni procedure is used [1]. It is noted that the Proni method isn't among methods of spectral estimation [1]. However, definition of a concept of complex frequency of to (5) and its physical sense allow the Proni method to carry to a method of the spectral analysis on the plane of complex frequencies.

In the spectral analysis of a physiological signal on the plane of complex frequencies there is an approximation of the sequence from counting of process x (1), x (n) of a linear combination of cosinusoids, fading, increasing or invariable on amplitude:

L 2

X M=X2 Ak exp [ak (n - 1)T]cos[2тfk (n - 1)T + 0k ]

(9)

K=1

where 1 < n < N; L - order of the approximating cosinusoidal model; T - an interval of counting in seconds; Ak and ak - amplitude and coefficient of attenuation (in s-1) k of a cosinusoid; fk and 0k - frequency (in Hz) and an initial phase (in rad) cosinusoid k. Parameters for plotting the spectrum on the complex frequency plane of the analyzing process are: fk, ak, Rk - the power of the k-th varying amplitude cosine wave, the index k varies from 1 to L/2. In Figs. 1 and 2 different variants of spectral analysis in the complex frequency plane of the transition process of the heart rhythm of the person resulting research are considered.

Fig. 1. The first version of the transition process: left - the dynamics of heart rate, right - heart rate spectrum of the transition process in the complex frequency plane

b) a

Fig. 2. Second option of transition process: a) dynamics of a warm rhythm; range of a warm rhythm of transition process on the plane of complex frequencies

In conclusion, we can note the following: 1) the range on the plane of complex frequencies as a frequency image of temporary dynamics of warm rhythm can be a basis at creation of system of classes of transition processes of a warm rhythm; 2) the generalized of complex frequencies parameters allow to estimate quantitatively reactivity of various mechanisms of regulation of warm rhythm in response to functional reflex types [10].

Conclusions

1. The range on the plane of complex frequencies as a frequency image of temporary dynamics of warm rhythm can be a basis at creation of system of classes of transition processes of a warm rhythm.

2. The generalized of complex frequencies parameters allow to estimate quantitatively reactivity of various mechanisms of regulation of warm rhythm in response to functional reflex types [10].

References

1. Marple S.L. Tsifrovoy spectral'nyy analiz i ego prilozhenie [Digital Spectral Analysis and Its Application]. Moscow, World Publ., 1990. 265 p.

2. Ragozin A.N., Kononov D.Yu. [Classic Spectral Analysis, Autoregressive Models, Analysis on the Plane of Complex Frequencies in Assessing the Structure of Heart Rate Oscillations]. Proc. of the Symposium. Slow Oscillatory Processes in the Human Body: Theory, Practical Application in Clinical Medicine and Prevention. Novokuznetsk, 2001, pp. 83-91. (in Russ.)

3. Ritm serdtsa u sportsmenov [A Heart Rhythm at Athletes]. Moscow, Physical Culture and Sport Publ., 1986. 143 p.

4. Kodkin V.L., Aksenov V.V., Usynin A.N., Wagner N.I. [Assessment of a Functional Condition of an Organism of School Students According to Variability of a Rhythm of Heart in the Conditions of Stress Psychoemotional Load]. Tezisy dokladov mezhdunarodnogo simpoziuma "Kompyuternaya elektrografiya na rubezhe stoletiy " [Thesis of Reports of the International Symposium "A Computer Electrography at a Turn of Centuries"]. Russia, Moscow, on April 27-30, 1999, pp. 54-59. (in Russ.)

5. Astakhov A.A., Usynin A.M., Ragozin A.N. [Non-REM Change Display to Users of Blood Circulation in the Sympathicotonia and Parasympathotonia at Children]. Proc. of Scientific Papers of

Рагозин А.Н., Тележкин В.Ф.

Анализ данных физиологических сигналов на плоскости комплексных частот с использованием процедуры Прони

the Symposium. Slow Oscillatory Processes in the Human Body: Theory, Practical Application in Clinical Medicine and Prevention, Novokuznetsk, 1997, pp. 88-95. (in Russ.)

6. Astakhov A.A., Ragozin A.N. Compare Method Heart Variability (HRV) for Preanesthesia Control. Proc. of 9th European Congress of Anesthesiology, Ierusalem, Israel, October 2-7, 1994.

7. Astakhov A.A., Ragozin A.N., Bubnova I.D., Govorov B. M. [Comparison of Absolute Values of Five Blood Circulation Parameters, their Variances, the Data of Spectral Analysis and Complex Frequency in Healthy Young Men at Rest]. Proc. of Scientific Papers of the Symposium. Slow Oscillatory Processes in the Human Body: Theory, Practical Application in Clinical Medicine and Prevention, Novokuznetsk, 1999, pp. 23-33. (in Russ.)

8. Weinstein L.A., Waxman D.Je. Razdelenie chastot v teorii kolebaniy i voln [Division of Frequencies in the Theory of Fluctuation and Waves]. Moscow, Science Publ., 1983. 288 p.

9. Ragozin A.N. [Complex Angular Coordinates of the Radar Center of Rreflection]. Radio electronics, 1991, no. 11, pp. 67-69. (in Russ.)

10. Ragozin A.N., Telezhkin V.F. Processing of Signals as the Sum of the Deter-Mined Function and Realization of Stationary Casual Process. 2nd International Conference on Industrial Engineering, Applications and Manufacturing (ICIEAM), 2016, pp. 154-157. DOI: 10.1109/ICIEAM.2016.7911548

Received 2 June 2017

УДК 612.17(075.8) DOI: 10.14529/^сг180108

АНАЛИЗ ДАННЫХ ФИЗИОЛОГИЧЕСКИХ СИГНАЛОВ НА ПЛОСКОСТИ КОМПЛЕКСНЫХ ЧАСТОТ С ИСПОЛЬЗОВАНИЕМ ПРОЦЕДУРЫ ПРОНИ

А.Н. Рагозин, В.Ф. Тележкин

Южно-Уральский государственный университет, г. Челябинск, Россия

Предложен метод спектрального анализа на плоскости комплексных частот с использованием процедуры Прони, а также графическая форма его представления. В физиологических исследованиях и диагностике функциональных состояний важное место занимает анализ данных физиологических сигналов как колебательных процессов. При этом широко распространенной моделью колебательного процесса является представление исследуемого процесса через суперпозицию (сумму) синусоидальных периодических процессов, каждый из которых характеризуется постоянными во времени амплитудой, частотой и начальной фазой. Набор рассчитанных из исходного сигнала параметров гармоник задает спектр исследуемого процесса. Необходимо отметить, что физиологические сигналы в общем случае представляют собой сложно организованные во времени колебательные процессы. Поэтому в литературе физиологические процессы в основном рассматриваются как случайные процессы. В этом случае в качестве меры распределения по частоте колебательной активности физиологического сигнала используют статистическую оценку спектральной плотности мощности сигнала, которая отражает зависимость распределения в среднем мощности (колебательной активности) сигнала от частоты.

Ключевые слова: временной ряд, метод Прони, комплексные частоты, спектральная плотность, физиологические сигналы, колебательные процессы, случайные процессы, колебательная активность.

Литература

1. Марпл-мл., С.Л. Цифровой спектральный анализ и его приложение / С.Л. Марпл-мл. - М. : Мир,1990. - 265 с.

2. Рагозин, А.Н. Классический спектральный анализ, авторегрессионные модели, анализ на плоскости комплексных частот в оценке структуры колебаний сердечного ритма / А.Н. Раго-

зин, Д.Ю. Кононов // Сборник научных трудов симпозиума. Медленные колебательные процессы в организме человека: теория, практическое применение в клинической медицине и профилактике. - Новокузнецк, 2001. - С. 83-91.

3. Ритм сердца у спортсменов / под. ред. Р.М. Баевского и Р.Е. Матылянской. - М.: Физкультура и спорт, 1986. - 143 с.

4. Оценка функционального состояния организма школьников по данным вариабельности ритма сердца в условиях стрессовой психоэмоциональной нагрузки / В.Л. Кодкин, В.В. Аксенов, А.Н. Усынин, Н.И. Вагнер // Тезисы докладов международного симпозиума «Компьютерная электрография на рубеже столетий», Россия, Москва, 27-30 апреля 1999 г. - М., 1999. - С. 54-59.

5. Астахов, А.А. Медленноволновые изменения показателей кровообращения при симпати-котонии и парасимпатикотонии у детей / А.А. Астахов, А.М. Усынин, А.Н. Рагозин // Сборник научных трудов симпозиума. Медленные колебательные процессы в организме человека: теория, практическое применение в клинической медицине и профилактике. - Новокузнецк, 1997. -С. 88-95.

6. Astahov, A.A. Compare method heart variability (HRV) for preanesthesia control / A.A. Astahov, A.N. Ragozin //Proc. of 9th European Congress of Anesthesiology, Ierusalem, Israel, October 2-7, 1994.

7. Сопоставление абсолютных значений пяти параметров кровообращения, их дисперсии, данных спектрального анализа и комплексных частот у здоровых молодых мужчин в покое / А.А. Астахов, А.Н. Рагозин, И.Д. Бубнова, Б.М. Говоров // Сборник научных трудов симпозиума. Медленные колебательные процессы в организме человека: теория, практическое применение в клинической медицине и профилактике. - Новокузнецк, 1999. - С. 23-33.

8. Вайнштейн, Л.А. Разделение частот в теории колебания и волн / Л.А. Вайнштейн, Д.Е. Ваксман. - М. : Наука, 1983. - 288 с.

9. Рагозин, А.Н. Комплексные угловые координаты радиолокационного центра отражения / А.Н. Рагозин //Радиоэлектроника. - 1991. - № 11. - С. 67-69.

10. Ragozin, A.N. Processing of Signals as the Sum of the Deter-Mined Function and Realization of Stationary Casual Process / A.N. Ragozin, V.F. Telezhkin // 2nd International Conference on Industrial Engineering, Applications and Manufacturing (ICIEAM). - 2016. - P. 154-157. DOI: 10.1109/ICIEAM. 2016.7911548

Рагозин Андрей Николаевич, канд. техн. наук, доцент кафедры инфокоммуникационных технологий, Южно-Уральский государственный университет, г. Челябинск; [email protected].

Тележкин Владимир Фёдорович, д-р техн. наук, профессор кафедры инфокоммуникационных технологий, Южно-Уральский государственный университет, г. Челябинск; telezhkinvf@ susu.ru.

Поступила в редакцию 2 июня 2017 г

ОБРАЗЕЦ ЦИТИРОВАНИЯ

Ragozin A.N. The Analysis of These Physiological Signals on the Plane of Complex Frequencies with Use of the Proni Procedure / A.N. Ragozin, V.F. Telezhkin // Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника». - 2018. - Т. 18, № 1. -С. 69-74. DOI: 10.14529/ctcr180108

FOR CITATION

Ragozin A.N., Telezhkin V.F. The Analysis of These Physiological Signals on the Plane of Complex Frequencies with Use of the Proni Procedure. Bulletin of the South Ural State University. Ser. Computer Technologies, Automatic Control, Radio Electronics, 2018, vol. 18, no. 1, pp. 69-74. DOI: 10.14529/ctcr180108