Научная статья на тему 'USING EMG PATTERNS FOR HUMAN GAIT CYCLE RECOGNITION'

USING EMG PATTERNS FOR HUMAN GAIT CYCLE RECOGNITION Текст научной статьи по специальности «Медицинские технологии»

CC BY
39
15
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
PATTERN RECOGNITION / LOWER LIMB MUSCLES / SIGNAL FILTERING / FUNCTIONAL DIAGNOSTICS / MOVEMENT ANALYSIS / ELECTROPHYSIOLOGY

Аннотация научной статьи по медицинским технологиям, автор научной работы — Bobe A.S., Konyshev D.V., Vorotnikov S.A.

The article describes the system for human gait cycle recognition based on EMG signal processing. The analysis of electrical signals produced by lower limb muscles during gait is carried out. The correlation between EMG signal characteristics and actual leg movements during gait are researched. The method of signal preprocessing and adaptive segmentation for detection of muscle activation patterns is developed. The most informative muscles for analysis are selected basing on signal to noise ratio. The pattern sequences processing method for gait cycle coding is proposed. The algorithm is implemented in C# and has shown about 90 % recognition rate during testing. The system can be used for exoskeleton control as well as in functional diagnostics and sports medicine applications.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

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

В статье описана система распознавания цикла ходьбы человека на основе обработки сигналов ЭМГ. Проведён анализ особенностей электрических сигналов от различных мышц нижних конечностей в процессе ходьбы. Изучены взаимосвязи характеристик ЭМГ сигнала с реальными фазами движения ног в процессе совершения шага. Разработан метод предобработки и адаптивной бинаризации ЭМГ сигнала для выделения паттернов активации мышц. Выбраны наиболее информативные для анализа мышцы на основе соотношения сигнал-шум. Предложен метод обработки последовательности паттернов для кодирования элементов цикла ходьбы. Реализован алгоритм распознавания фаз ходьбы на языке C#, при тестировании достигнут 90 % уровень верного распознавания элементов цикла ходьбы. Система предназначена для использования при управлении экзоскелетом, но может быть также использована в системах функциональной диагностики и спортивной медицине.

Текст научной работы на тему «USING EMG PATTERNS FOR HUMAN GAIT CYCLE RECOGNITION»

26. Saurbaeva R. A., Yesaydar U. S. Primenenie zarubezhnykh sistem bronirovaniya i ispolzovanie tekhnologiy v gostinichnom biznese [Application of foreign booking systems and use of technologies in hotel business]. Vestnik Almatin-skogo tekhnologicheskogo universiteta [The Bulletin of the Almaty Technological University], 2013, no. 4, pp. 86-89.

27. Skobkin S. S. Industriya gostepriimstva ili turizma? Kto prav? [Industry of hospitality or tourism? Who is right?]. Rossiyskoepredprinimatelstvo [Russian Business], 2012, no. 21, pp. 130-135.

28. Startseva O. M. Razvitie sportivnogo turizma v Astrakhanskoy oblasti [Development of sporting tourism in the Astrakhan region]. Ekologiya Rossii: na puti k innovatsiyam [Ecology of Russia: on the way to innovations]. 2015. № 11. S. 147-150.

29. Trubitsina L.F., Zanozin V.V. Prirodno-rekreatsionnye resursy Limanskogo rayona Astrakhanskoy oblasti i per-spektivy ikh ispolzovaniya dlya razvitiya turizma [Natural and recreational resources of the Limansky district of the Astrakhan region and perspective of their use for tourism development]. Ekologiya Rossii: na puti k innovatsiyam [Ecology of Russia: on the Way to Innovations], 2013, no. 7, pp. 27-30.

30. Tyger L. M., Goryacheva Ye. D., Sultaeva N. L., Kostikova N. V., Rudina L. M. Sovremennye trebovaniya k personalu v sisteme upravleniya chelovecheskimi resursami v industrii gostepriimstva [The modern requirements to staff in human resources management system in the hospitality industry]. Naukovedenie [Science of Science], 2016, vol. 8, no. 2, pp. 80.

31. Shtankova I. S., Bezuglova M. S. Problemy i perspektivy razvitiya rechnogo kruiznogo turizma v Astrakhanskoy oblasti [Problems and perspectives of river cruise tourism development in the Astrakhan region]. Astrakhanskiy vestnik ekologicheskogo obrazovaniya [The Astrakhan Bulletin of Ecological Education], 2014, no. 1 (27), pp. 215-217.

УДК 004.5, 004.93, 612.776

USING EMG PATTERNS FOR HUMAN GAIT CYCLE RECOGNITION1

The article has been received by editorial board 02.07.2016, in the final version — 28.07.2016.

Bobe Anatoly S., Engineer, Neurobotics Ltd., 4922 passage, 2/4 Yuzhnaya promzona, Moscow, Ze-lenograd, 124498, Russian Federation, e-mail: [email protected]

Konyshev Dmitry V., post-graduate student, Bauman Moscow State Technical University, 5/1 Baumanskaya 2nd St., Moscow, 105005, Russian Federation, e-mail: [email protected]

Vorotnikov Sergey A., C.Sc. (Egineering), Associate professor, Bauman Moscow State Technical University, 5/1 Baumanskaya 2nd St., Moscow, 105005, Russian Federation, e-mail: [email protected]

The article describes the system for human gait cycle recognition based on EMG signal processing. The analysis of electrical signals produced by lower limb muscles during gait is carried out. The correlation between EMG signal characteristics and actual leg movements during gait are researched. The method of signal preprocessing and adaptive segmentation for detection of muscle activation patterns is developed. The most informative muscles for analysis are selected basing on signal to noise ratio. The pattern sequences processing method for gait cycle coding is proposed. The algorithm is implemented in C# and has shown about 90 % recognition rate during testing. The system can be used for exoskeleton control as well as in functional diagnostics and sports medicine applications.

Keywords: electromiography, biosignal processing, pattern recognition, lower limb muscles, gait, signal filtering, functional diagnostics, movement analysis, electrophysiology

Graphical annotation (Графическая аннотация)

Движение

___Movement

'Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований, грант №16-07-01080.

ИСПОЛЬЗОВАНИЕ ПАТТЕРНОВ ЭМГ В ЗАДАЧЕ РАСПОЗНАВАНИЯ ЦИКЛА

ХОДЬБЫ ЧЕЛОВЕКА

Бобе Анатолий Сергеевич, инженер, ООО «Нейроботикс», 124498, Российская Федерация, г. Москва, Зеленоград, Южная промзона, проезд 4922, стр. 2, д. 4, e-mail:[email protected]

Конышев Дмитрий Владимирович, аспирант, Московский государственный технический университет имени Н. Э. Баумана, 105005, Российская Федерация, г. Москва, 2-я Бауманская ул., д. 5, стр. 1, e-mail: [email protected]

Воротников Сергей Анатольевич, кандидат технических наук, доцент, Московский государственный технический университет имени Н. Э. Баумана, 105005, г. Москва, 2-я Бауманская ул., д. 5, стр. 1, e-mail: [email protected]

В статье описана система распознавания цикла ходьбы человека на основе обработки сигналов ЭМГ. Проведён анализ особенностей электрических сигналов от различных мышц нижних конечностей в процессе ходьбы. Изучены взаимосвязи характеристик ЭМГ сигнала с реальными фазами движения ног в процессе совершения шага. Разработан метод предобработки и адаптивной бинаризации ЭМГ сигнала для выделения паттернов активации мышц. Выбраны наиболее информативные для анализа мышцы на основе соотношения сигнал-шум. Предложен метод обработки последовательности паттернов для кодирования элементов цикла ходьбы. Реализован алгоритм распознавания фаз ходьбы на языке C#, при тестировании достигнут 90 % уровень верного распознавания элементов цикла ходьбы. Система предназначена для использования при управлении экзоскелетом, но может быть также использована в системах функциональной диагностики и спортивной медицине.

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

Адрес видеоролика (Address of the video)

http://hi-tech.asu.edu.ru/docs/addmaterials/2016-3-1(Gait).mp4

The development of the systems, based on the biofeedback control, has been a subject of interest of many researchers during the last several years. In particular, significant progress has been made in elec-tromiogram (EMG) signals recognition for such applications as the limb prostheses operation, functional diagnostics and sports medicine. Most of the recognition systems utilize machine learning approach (neural networks [8], Support Vector Machines [5]) to recognize the certain movements and simple gestures of a limb. This approach has been widely used and proved its efficiency, but it is actually not suitable for recognizing such complex, extended in time and subject-specific movements as strides during human walking. We have found no evidence that EMG informationfrom lower limbs has ever been successfully used for movement recognition and/or biofeedback control, such as exoskeleton [3, 6, 12]. At the same time, these applications have been implemented for the upper limb movement.

The aim of this work was to develop a robust and efficient algorithm for human stridemovement recognition. Another issue was to minimize the time delay of the recognition system, so that it could be able to operate as close to real-time regime as possible.

General characteristic of article theme. Constraint for time delay of the recognition system is essential for any biofeedback control applications. Assuch type of movement as stride is rather extended in time; there is no possibility to learn the whole history of the corresponding EMG signal in advance before processing. That's why we abandoned the idea to use machine learning methods to solve the problem and concentrated on empirical pattern recognition approach.Our method for resolve the problem can be summarized as follows.

• Filtered signals from each of the considered muscles are thresholdedinto activation and deactivation phases.

• A stride movement is described by activation pattern [5], consisting of several stages, which correspond to different combinations of activated and deactivated muscles.

• Threshold values for each muscle activation, as well as generally stable stages and the order of their occurrence for the strides, are adaptively calculated for each subject during training phase.

• During testing phase the system recognizes a stride as a set of activation patterns, occurring in the defined order.

This approach can be used for the muscle activity research for different types of movement, as well as for training a classifier for automatic recognition of these movements or their sequences, such as in [2, 10].

EMG signal processing tools. All our experiments were carried out using the Kardi3/9 wireless 9-channel amplifier with sampling rate of 1000 Hz. A subject performed series of 5 to 9 separate strides with a pause of 2-3 seconds between them. The obtained signals were passed through bandpass filter with low and high frequencies of 30 and 200 Hz respectively. Then the signal root mean square (RMS) was calculated. For gait analysis we do not need the high frequency information. So we applied a lowpass Butterworth filter with cutoff of just 1 Hz to the resulting signal. This way smooth curves, characterizing the general muscle activation in scope of several seconds, were achieved (Fig. 1).

At the first stage of our experiments the goal was to find the most informative (in the relation to the signals recognition) muscles. A wide number of muscles are involved in gait movements, but very few of them actually provide stable and reliable signals for stride phases recognition. Some of the measured signals are presented on Fig. 2. It is obvious, that although different muscles produce different activation intensities and pattern forms during strides, the actual activation of each muscle almost always happens at the very beginning of the stride and deactivation happens at its end.Exception are the muscles, which have non-zero base level - as they are activated during standing, such as Gluteus Medius. Under «activation» we mean the level of signal6 which significantly distincts from noise for this muscle (over 3c deviation).

Figure 1 - Muscle activation signal sample

We have chosen Soleus and Tibialis Anterior muscles as the most informative for our task [4, 9]. The patterns generated by these muscles have proved to be rather stable during all our experiments. The former provides a single-peak activation pattern during the leg pushing off the ground.The latter remains active during the whole leg movement providingone local minimum around the push off stage and,for certain subjects,another one before the start of the support stage (when a leg is approaching the ground). The correspondence between the actual stride stages and the muscle patterns, which is shown on Fig. 3, was determined using accelerometer, attached to one of the legs.As we passed to recognition in our later experiments we utilizedonly the signals from these two muscles, from both legs.

In our experiments average voltage during activation was around 40^Vfor Tibialis Anterior and 100^V for Soleus. With average fluctuations during rest state under 7 ^V and 20 ^V, both muscles provided signal to noise ratio of about 5, which is sufficient for the task.

As presented on the figures above, in ideal case it is possible to precisely extract the time periods of six stages of a stride. In practice, however, signals may contain some other local peaks or artifacts which make it difficult to judge about the gait cycle utilizing the signal form. A much more consistent, though more rough approach, is to apply adaptive signal thresholding and extract the binary activation patterns of each muscle as described below.

Signal Thresholding Method. In this paper we propose a novel algorithm for signal thresholding, which was designed to automatically distinguishthe signal features observed in our experiments.

The general aim of signal thresholding in our system is to localize the most significant and stable EMG signal phases. At the same time these signals are subject-specific as the walking movements are generally not the same for different persons [7]. Thus the threshold values for the muscles are calculated given a training sample, consisting of several strides made by subject.

The criteria for the threshold value is the following: the number of activations/deactivations should be the same for each step and should be stable enough to variations of a threshold value.

Tibialis Anterior

Vastus mediaiis

Semitendinosus

Biceps femoris

TO

n

o M

S

O

TO g

o d ft

o C3 Q*

Ci' TO

0 C3

1

5'

TO

Gluteus medius

Gracilis

Vastus lateralis

Soleus

(SC) £ '9102 'saiSojouipax il§ïH PUB IOJ^uoo ¡IVNHllOr NVIdSVO

tz

: Hi1

Figure 3 - Correspondence between the stride stages and the muscle activation

Given a signal s(t) within a period of time 7and number of steps in a sample n a threshold value xt corresponds to a maximum value offx, Ax) at x = xt:

x+Ax

f (x, Ax) = | R(x)g(x) dx ^ max ,

x-Ax

T

d (sign(s(t ) - x))

dt,

0 dt fl, if g(x) = k * n, k = 2,4

[ 0 otherwise

where Ax is a window size. The Ax value is chosen with respect to the desired stability of a threshold chosen.

In other words, if we display the number of crossings of a signal and a threshold level line as a function of threshold value g(x), the optimal value would correspond to a middle of the longest plateau of this function (Fig. 4). Such approach is not of much use in terms of muscle activity detection, but gives a good representation of some basic patterns of a particular movement.

During our experiments threshold levels were calculated independently for each muscle, using a sample signal containing from 7 to 10 strides, made by test subject.

Signal segmentation. After signal thresholding a set of activation patterns are obtained for each step as shown on Figure 5. We can introduce a stage as a particular combination of activated and deactivated muscles. A stage can be therefore presented in a form of a binary code. Note that the same stage may occur several times during single movement.

We have developed an algorithm, which performs automatic stage extraction from a set of activation patterns.This algorithm learning the order of patterns occurrence, their duration and their repetitiveness from stride to stride. A stage is excluded from consideration in following cases.

• If its average duration is below a specified limit (we set it as a 3% of average overall stride duration).

• If its position in stride patternsis inconsistent.

• If it is not detectable inmore than 20 % strides within the training signal.

The stages extracted during training phase are associated with particular physical movements of the legs, according to the scheme, presented on Figure 3. Depending on the test subjects from 4 to 6 consistent stages were usually detected, successfully covering the major leg movements.

As for the test mode, our system was configured to detect the stride movement as a set of stages in predefined order. During the neutral position (0000) the first stage is awaited. When its occurrence is detected, the second stage is awaited and so on. If a wrong stage occurs at some point of the stride the system ignores it. If a neutral stage signal 0000 occurs, the system always resets to the neutral position (an interrupted stride is declared).

IO №

x 10

o >

CO

NN >

53

Ch

O

C »

53 >

r

o o 3

r+ >1 O i—»

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

V 3 0.

ffi

ll

-3 n o tf 3 o

lit

Figure 4 - Signal thresholding (Tibialis Anterior). A sample signal containing 9 steps can be thresholded to detect either 1 activation/deactivation per stride (18 crossings overall) or 2 (36 crossings). The second variant gives more phases for each stride but is less stable

Is) O

CO

CO 01

Figure 5 - Activation Patterns from two pairs of muscles (R-right; L-left; S-Soleus; Ta-Tibialis Anterior). A set of 5 consistent stages can be extracted judging on these 3 strides. Our algorithm will discard a 0111 stage between 2 and 3 as it is too short, as well as the final 0100 stage as it is almost not detectable within the 3rd stride

We tested our system on 5 subjects, each performing from 7 to 10 strides for training and about 50 strides for testing. During the experiments we achieved the recognition rate about 90 %. Actually the beginning of the movement was correctly detected for 98 % of the strides.The rest of errors were due to false interruption or mis-detection of one of the stages.The system time delay did not exceed 0,3 seconds, which we consider satisfactory for our task (it is small in comparison with the durationof theentire stride movement).

Conclusion. We have developed a system for single stride recognition, based on analysis of a filtered EMG signal of lower limb muscles. We have localized the most informative muscles and determined the relation between their EMG signal and corresponding actual leg movements during stride. Alsoauthors have developed an algorithm for adaptive signal thresholding - to generate the pattern of muscle activation/deactivation. The patterns, obtained during training, were automatically analyzed to determine their most consistent stages - the combinations of activated/deactivated muscles. A complex movement was interpreted as an ordered set of stages. This approach allowed us to efficiently recognize with relatively small time delay such a long and complex movement as a stride. The configuration, described in this paper, was used for the stride detection.However the system can also be trained and configured for any other complex movement recognition, using a desired number of EMG channels, without any major changes. This makes our approach general and flexible for different tasks.

We consider our system as a promising step in the field of biofeedback-based control system development. It could also be useful in medical applications such as functional diagnostics and rehabilitation.

In our future work we plan to use more EMG channels for obtaining more robust and accurate detection of different stride stages. Another issue is to develop a system, which could rapidly adapt for recognition of any complex movement of lower or upper limb by automatically choosing the most informative EMG channels and produce the general «map» of the particular movement and the muscles, involved in it.

Список литературы

1. Benedetti M. G. Muscle Activation Patterns During Level Walking and Stair Ambulation / M. G. Benedetti, V. Agostini, M. Knaflitz, P. Bonato // Applications of EMG in Clinical and Sports Medicine. - Available at: http://www.intechopen.com/books/applications-of-emg-in-clinical-and-sports-medicine/muscle-activationpatterns-during-level-walking-and-stair-ambulation (accessed 20.06.2016).

2. Brumshteyn Y. System analysis of approaches and functional capabilities of human body movement research methods / Y. Brumshteyn, J. Aksenova, M. Ivanova // Caspian Journal: Control and High Technologies. - 2014. - № 3. -P. 80-98.

3. Gopura R. A. R. C. Application of Surface Electromyographic Signals to Control Exoskeleton Robots / R. A. R. C. Gopura, K. Kiguchi // Applications of EMG in Clinical and Sports Medicine. - Available at: http://www.intechopen.com/ books/applications-of-emg-in-clinicaland-sports-medicine/application-of-surface-electromyographic-signals-to-control-exoskeleton-robots (accessed 20.06.2016).

4. He H. A Study on Lower-Limb Muscle Activities during Daily Lower-Limb Motions / H. He, K. Kiguchi, E. Hori-kawa // International Journal of Bioelectromagnetism. - 2007. - Vol. 9, № 2. - P. 79-84.

5. Hernandez R. SVM-Based Volitional Artificial Leg Control Via Ubiquitous Small and Low Power Architectures / R. Hernandez // University of Rhode Island, Open Access Dissertations, Paper 206. - 2014. - Available at: http://digitalcommons.uri.edu/oa_diss/206 (accessed 20.06.2016).

6. Medved V. Kinesiological Electromyography / V. Medved, M. Cifrek // Biomechanics in Applications. - 2011. -P. 349-366.

7. Schmitz A. Differences in lower-extremity muscular activation during walking between healthy older and young adults / A. Schmitz, A. Silder, B. Heiderscheit, J. Mahoney, D. G. Thelen // Journal of Electromyography and Kinesiology. -2009. - № 19. - P. 1085-1091.

8. Shahmoradi M. H. Combined Classification and Regression for Simultaneous and Proportional EMG Control of Wrist Forces / M. H. Shahmoradi, M. A. Akhaee, M. S. Mirian // 23rd European Signal Processing Conference (EUSIPCO). -2015. - P. 2471-2475.

9. Srivastava J. An Electromyography analysis of lower limb muscles for differentlocomotion activities / J. Srivastava, A. Srivastava, R. P. Tewari // InternationalJournal of Engineering Research & Technology (IJERT). - 2012. - Vol. 1, issue 3. -P. 1-3.

10. Starchenko I. Mimic gesture recognition based on electromyography signal processing / I. Starchenko, R. Budko // Caspian Journal: Control and High Technologies. - 2016. - № 1. - P. 39-50.

11. Todd Farrell M. A method to determine the optimal features forcontrol of a powered lower-limb prostheses / Todd Farrell M., H. Herr // 2011 Annual International Conference of the IEEE Engineering in Medicine and Biology Society. - P. 6041-6046.

12. Vorotnikov S. Biometric system for arm prosthesis control / S. Vorotnikov, V. Strunin, N. Vybornov // Caspian Journal: Control and High Technologies. - 2013. - № 3. - P. 147-162.

13. Winter D. A. EMG profiles during normal human walking: stride-to-stride andinter-subject variability / D. A. Winter, H. J. Yack // Electroencephalography and clinical Neurophysiology. - 1987. - № 67. - P. 402-441.

References

1. Benedetti M. G., Agostini V., Knaflitz M., Bonato P. Muscle Activation Patterns During Level Walking and Stair Ambulation. Applications of EMG in Clinical and Sports Medicine. Available at: http://www.intechopen.com/books/ applications-of-emg-in-clinical-and-sports-medicine/muscle-activationpatterns-during-level-walking-and-stair-ambulation. (accessed 20.06.2016).

2. Brumshteyn Y., Aksenova J., Ivanova M. System analysis of approaches and functional capabilities of human body movement research methods. Caspian Journal: Control and High Technologies, 2014, no. 3, pp. 80-98.

3. Gopura R. A. R. C., Kiguchi K. Application of Surface Electromyographic Signals to Control Exoskeleton Robots. Applications of EMG in Clinical and Sports Medicine. Available at: http://www.intechopen.com/books/applications-of-emg-in-clinicaland-sports-medicine/application-of-surface-electromyographic-signals-to-control-exoskeleton-robots (accessed 20.06.2016).

4. He H., Kiguchi K., Horikawa E. A Study on Lower-Limb Muscle Activities during Daily Lower-Limb Motions. International Journal of Bioelectromagnetism, 2007, vol. 9, no. 2, pp. 79-84.

5. Hernandez R. SVM-Based Volitional Artificial Leg Control Via Ubiquitous Small and Low Power Architectures. University of Rhode Island, Open Access Dissertations, Paper 206, 2014. Available at: http://digitalcommons.uri.edu/ oa_diss/206 (accessed 20.06.2016).

6. Medved V., Cifrek M. Kinesiological Electromyography. Biomechanics in Applications, 2011, pp. 349-366.

7. Schmitz A., Silder A., Heiderscheit B., Mahoney J., Thelen D. G. Differences in lower-extremity muscular activation during walking between healthy older and young adults. Journal of Electromyography and Kinesiology, 2009, no. 19, pp. 1085-1091.

8. Shahmoradi M. H., Akhaee M. A., Mirian M. S. Combined Classification and Regression For Simultaneous and Proportional EMG Control of Wrist Forces. 23rd European Signal Processing Conference (EUSIPCO), 2015, pp. 2471-2475.

9. Srivastava J., Srivastava A., Tewari R. P. An Electromyography analysis of lower limb muscles for differentlocomotion activities. InternationalJournal of Engineering Research & Technology (IJERT), 2012, vol. 1, issue 3, pp. 1-3.

10. Starchenko I., Budko R. Mimic gesture recognition based on electromyography signal processing. Caspian Journal: Control and High Technologies, 2016, no. 1. , pp. 39-50.

11. Todd Farrell M., Herr H. A method to determine the optimal features forcontrol of a powered lower-limb prostheses. 2011 Annual International Conference of the IEEE Engineering in Medicine and Biology Society, pp. 6041-6046.

12. Vorotnikov S., Strunin V., Vybornov N. Biometric system for arm prosthesis control. Caspian Journal: Control and High Technologies, 2013, no. 3, pp. 147-162.

13. Winter D. A., Yack H. J. EMG profiles during normal human walking: stride-to-stride andinter-subject variability. Electroencephalography and clinical Neurophysiology, 1987, no. 67, pp. 402-441.

i Надоели баннеры? Вы всегда можете отключить рекламу.