Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

СТАБИЛОМЕТРИЧЕСКИЙ ИНВАРИАНТ УДЕРЖАНИЯ РАВНОВЕСИЯ В НОРМЕ И ПАТОЛОГИИ

Давыдов О.Д. 1 Монтиле А.И. 2 Марчук Ю.В. 3 Монтиле А.А. 4
1 ФГБУ «Уральский научно-исследовательский институт травматологии и ортопедии им. В.Д. Чаклина» Минздрава РФ
2 ФГБОУ ВПО «Уральский государственный лесотехнический университет»
3 ГБУЗ СО ДКБ восстановительного лечения научно-практический центр «Бонум»
4 ФГБУН «Ботанический сад УрО РАН»
В результате клинической апробации нового способа обработки данных стабилометрического тестирования выявлена пространственная конфигурация частоты появления интервалов движения с постоянной скоростью, инвариантная относительно анатомических и функциональных патологий костно-суставного аппарата, мышечной и нервной систем. Параметрические и структурные особенности данной конфигурации обеспечивают дифференциацию нарушений функции равновесия по двигательной (мышечной), управляющей (нервной) и опорной (костно-суставной) подсистемам. Аппроксимация поверхности конфигурации, полученная в результате обработки данных стабилометрического тестирования четырех разновозрастных групп обследуемых без выявленных патологий, может быть использована для замены «нормы» при сравнительной обработке клинических данных конкретных пациентов и в качестве модели идеального удержания равновесия при отсутствии целенаправленного движения для верификации концептуальных формальных моделей организации движения.
компьютерная стабилометрия
нервная система
мышечная система
костно-суставной аппарат
1. Бернштейн Н.А. Физиология движения и активность. – М.: Наука, 1990. – 495 с.
2. Братусь А.С., Новожилов А.С., Платонов А.П. Динамические системы и модели биологии. – М.: ФИЗМАТЛИТ, 2010. – 400 с.
3. Витяев Е.Е. Объяснение Теории Движений Н.А. Бернштейна // VII Всероссийская научно-техническая конференция «Нейроинформатика-2005»: сборник научных трудов. – Часть 1. – М., 2005. – С. 234–240.
4. Гаже П.М., Вебер Б. Постурология. Регуляция и нарушения равновесия тела человека. – СПб.: Изд. Дом СПбМАПР, 2008. – 316 с.
5. Гельфанд И.М., Розенфельд Б.И., Шифрин М.А. Очерки о совместной работе математиков и врачей. – М.: Едиториал УРСС, 2011. – 320 с.
6. Гумбель Э. Статистика экстремальных значений. – М.: Мир., 1965. – 460 с.
7. Монтиле А.А., Шавнин С.А., Монтиле А.И. Кластерный и регрессионный анализ изменений количественных морфометрических признаков // Математическая биология и биоинформатика. – 2007. – Т. 2, № 1. – С. 60–65.
8. Исследование хаотических колебаний в нелинейной системе управления поддержанием вертикальной позы тела человека / М.Г. Розенблюм, Г.И. Фирсов, Р.А. Кууз, В.Н. Суслов // Вестник научно-технического развития: научно-технический журнал. – 2007. – № 3. – С. 32–42.
9. Скворцов Д.В. Клинический анализ движений, стабилометрия. – М.: АОЗТ «Антидор», 2000. – 192 с.
10. Ташкинов А.А., Вильдеман А.В., Бронников В.А. Индивидуальное прогнозирование двигательного развития у больных детским церебральным параличом на основе подходов статистического анализа // Российский журнал биомеханики. – 2010. – Т. 14, № 2 (48). – С. 69–78.
11. Тяжелов А.А. Анализ стабилограмм на основе математической модели тела человека как многозвенной системы / А.А. Тяжелов, Н.Н. Кизилова, В.А. Фищенко [и др.] // Травма. – 2012. – Т. 13, № 4. – С. 17–25.
12. Динамическая стабилизация вертикального положения тела человека / В.И. Усачёв, М.И. Говорун, А.Е. Голованов, М.С. Кузнецов // Известия Южного федерального университета: Технические науки. – 2010. – Т. 110, № 9. – С. 164–169.
13. Юдин Д.Б., Юдин А.Д. Математики измеряют сложность. – М.: Книжный дом «ЛИБРОКОМ», 2014. – 192 с.
14. Пат. 2497451 РФ, МПК А 61 В 5/103. Способ диагностики функциональных нарушений опорно-двигательного аппарата / О.Д. Давыдов, А.И. Монтиле, Ю.В. Марчук, Н.Л. Кузнецова; ФГБУ «Уральский НИИТО им. В.Д. Чаклина» МЗ РФ. – № 2012124615/14; заявл.14.06.2012; опубл.10.11.2013, Бюл. № 31. – 16 с.

Одним из современных и перспективных методов исследования функционального состояния двигательной системы человека является стабилометрия [4, 9]. Актуальность рассматриваемой тематики, на фундаментальном уровне относящейся к разработке математических моделей физиологии движения [1, 4], обусловлена в первую очередь потребностями клинической практики, связанными с потенциальной возможностью использовать результаты стабилометрических исследований для решения задач диагностики, дифференциальной диагностики и оценки эффективности лечебных мероприятий. Проблема заключается в невозможности решения вышеперечисленных задач на базе единой методики обработки данных стабилометрических тестов [4], создание которой тормозится крайне узким применением стабилометрии в клинической практике лечебных учреждений любого профиля, несмотря на высокий по международным стандартам уровень разработанной аппаратуры, математического и программного обеспечения и невысокую по сравнению с другими современными методами диагностики стоимость оборудования.

Основной причиной сложившейся ситуации, на наш взгляд, является не столько отсутствие у врачей специальных знаний в соответствующей области, в частности математических методов обработки данных, сколько давно отмеченный на методологическом уровне разрыв между «теоретической» медициной, предполагающей существенный уровень статистического обобщения результатов наблюдений, и клинической практикой, с необходимостью, ориентированной как раз на учет индивидуальных особенностей организма пациента, специфики форм и заболеваний [5]. Например, при травмах в разной степени поражаются структурные и функциональные элементы костно-суставной, мышечной и нервной систем, при этом диагностика травм, определение стратегии и тактики лечения предполагает различение патологических процессов и определение степени влияния каждого из повреждений на общее состояние человека с учетом анамнеза.

Непосредственно получаемым результатом отдельного стабилометрического теста является упорядоченный по номерам измерений, осуществляемым через фиксированный промежуток времени, набор пар чисел: фронтальная и сагиттальная координаты центра давления стоп на стабилометрическую платформу. Стандартным способом отображения результатов являются: фронтальная стабилограмма (изменение фронтальной координаты во времени), сагиттальная стабилограмма (изменение сагиттальной координаты во времени) и статокинезиограмма (траектория перемещения центра давления по плоскости опоры). На основе регистрируемых данных вычисляются стандартные, векторные и частотные показатели. В качестве критерия для оценки функции равновесия используются либо отдельные показатели, либо один производный критерий, значение которого вычисляется алгоритмически или по некоторой содержательно обосновываемой формуле [9, 10, 12].

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

Ситуация усугубляется отсутствием объясненных на уровне анатомии и/или физиологии в норме и патологии (т.е. клинически интерпретируемых) связей между изменениями значений подавляющего большинства вышеперечисленных наборов показателей и наблюдаемой врачом структурной (лучевые методы диагностики) и функциональной (функциональные методы диагностики) динамикой состояния двигательной системы больного. Что обусловлено, в первую очередь, исторически сложившимся «стохастическим» подходом к анализу и интерпретации результатов наблюдений, инициированным естественной комбинаторной сложностью изучаемой системы [13], при любом способе моделирования наблюдаемых процессов. Для иллюстрации на рис. 1 приведены примеры статокинезиограмм, фиксирующих перемещение центра давления по платформе при прохождении «простейшего» 20-секундного стабилометрического теста для обследуемого без выявленных патологий (рис. 1, а) и пациентов с различными патологиями (рис. 1 б, в, г, д), демонстрирующие «очевидную» хаотичность движений.

Новая система показателей

С целью преодоления вышеописанных недостатков был разработан новый набор показателей и способ их получения [14], с 2012 года проходящий клиническую апробацию в УНИИТО им. В.Д. Чаклина. Теоретической основой предложенной системы показателей служит модель движения, учитывающая в явном виде фактически выявленную на всех уровнях организации движения дискретность процессов [1, 4]. Второй особенностью модели, по сравнению с единственной известной в различных модификациях кинематической моделью Бернштейна – Гельфанда, является первичность моделирования взаимодействующих мышечной и нервной подсистем. Костно-суставная система рассматривается как «источник» динамических ограничений.

Базовым понятием новой системы показателей является интервал неизменного движения – двухкомпонентная величина (delta_ti, Vi), характеризующаяся значениями двух переменных: длительность интервала (delta_ti) и постоянная для интервала скорость (Vi). Рассматриваются интервалы неизменного движения при перемещении во фронтальной и сагиттальной плоскостях, а также для статокинезиограммы. В последнем случае определяются интервалы с постоянными линейными и угловыми скоростями. Помимо вышеперечисленных «элементарных» интервалов вводятся «обобщенные» интервалы неизменного движения, характеристической функцией принадлежности к множеству которых является одновременное постоянство скорости изменения фронтальной и сагиттальной координат или линейной и угловой скоростей. Стабилограммы интерпретируются как последовательности элементарных интервалов. Траектория перемещения проекции центра давления по стабилометрической платформе (статокинезиограмма) рассматривается как последовательность обобщенных интервалов неизменного движения. Границам интервалов сопоставляются мгновенные изменения соответствующих скоростей.

pic_19.tif аpic_22.tif бpic_20.tif в
pic_23.tif гpic_24.tif дpic_21.tif е

Рис. 1. Статокинезиограммы для обследуемых с различными диагнозами: а – физиологическая норма; б – сенситивная атаксия; в – ДЦП, спастическая диплегия; г – двусторонний остеоартроз тазобедренных суставов; д – посттравматический остеоартроз таранно-пяточного сустава; е – идиопатический сколиоз II степени

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

В процессе клинической апробации на заболеваниях как с доминированием патологий отдельных систем, так и сочетанных, связанной с отработкой методик диагностики, контроля за ходом лечения и оценки его результатов, был обнаружен фундаментальный, с точки зрения физиологии движения и системного анализа в первую очередь теории управления сложными системами, факт – наличие однотипной формы поверхности, отображающей частоту появления интервалов неизменного движения с определенными длительностями и скоростями при удержании равновесия, если отсутствует целенаправленное движение. Различным анатомическим и функциональным патологиям соответствуют специфичные параметрические и визуальные изменения формы, более точно – специфические отклонения от «нормальной» формы, присущей обследуемым без выявленных патологий, которые уменьшаются при успешном лечении. Для иллюстрации на рис. 2 приведен пример первичной обработки результатов прохождения стабилометрического теста пациентом с идиопатическим сколиозом I ст.; для сравнения выводятся результаты тестирования здорового обследуемого.

pic_25.tif

Рис. 2. Частота появления интервалов движения длительностью delta_ti (ось Х) с постоянной скоростью Vi (ось Y) во фронтальной плоскости у обследуемого без выявленных патологий (зеленый цвет) и обследуемого с идиопатическим сколиозом I ст. (красный цвет)

Материалы и методы исследования

В обследовании приняло участие 324 человека в возрасте от 6 до 75 лет. Из них 140 человек (контрольная группа) без выявленной патологии (в том числе: 15 – в возрасте 6–8 лет; 42 – в возрасте 18–22 года; 31 – в возрасте 30–45 лет; 28 – в возрасте 50–60), 26 больных с посттравматическим остеоартрозом голеностопных суставов, 36 больных с деформирующим остеоартрозом коленных суставов, 54 больных с деформирующим остеоартрозом тазобедренных суставов, 16 больных ДЦП, 22 больных идиопатическим сколиозом I–II степени, 40 больных продольным плоскостопием I–II степени, 12 пациентов с дисметаболическими полинейропатиями нижних конечностей.

Исследования проводились с использованием компьютерного стабилоанализатора «Стабилан-01» (ЗАО ОКБ «Ритм», г. Таганрог). Обследование включало проведение стандартного стабилометрического теста продолжительностью 20 секунд. Для повышения достоверности получаемых данных при каждом обследовании тест повторялся три раза с перерывами в 3–5 минут.

Для больных, проходящих лечение в УНИИТО им. В.Д. Чаклина, тестирование проводилось до операции, на 12 сутки после проведения операции, а также через 6 и 12 месяцев.

Для фронтальной и сагиттальной стабилограмм, полученных при обследовании, были определены значения новых показателей и сформированы необходимые для дальнейшего анализа выборки результатов. По каждой из стабилограмм для фронтальной и сагиттальной координат были выделены интервалы движения центра давления с постоянной скоростью и определены длительность каждого интервала (delta_ti) и скорость на нем (Vi), а также количество интервалов с различными длительностями и скоростями. Результаты были экспортированы в программу обработки биометрических данных «Радикал» [7], с помощью которой были построены диаграммы рассеяния (delta_t, V), и в «Statistica v.6.0», с помощью которых были получены гистограммы, характеризующие распределения интервалов (delta_ti, Vi), длительностей (delta_ti) и скоростей (Vi) и определены значения числовых параметров.

Результаты исследования и их обсуждение

Полученные в результате обработки данных пассивных стабилометрических тестов для обследуемых без выявленных патологий («здоровые», «норма») из четырех разновозрастных групп, трехмерные гистограммы распределения интервалов движения с постоянной скоростью (delta_ti, Vi) имеют повторяющуюся при всех измерениях форму, одинаковую как для фронтальной (X), так и для сагиттальной (Y) координат (рис. 3, а, б).

pic_26.tif аpic_27.tifб
pic_28.wmfв pic_29.wmf г pic_30.wmf д

Рис. 3. Частота появления интервалов движения с постоянной скоростью при трехкратном проведении стабилометрического теста обследуемым без выявленных патологий для сагиттальной (а) и фронтальной (б) плоскостей, гистограммы распределения скоростей (в), длительностей интервалов (г) и диаграмма рассеяния интервалов (д)

Основными визуальными признаками экспериментально выявленной формы являются колоколообразный вид сечений распределения по delta_ti = const и гистограммы скоростей в целом (рис. 3, в), монотонное убывание значений в сечениях по Vi = const, и гистограммы длительностей интервалов в целом (рис. 3, г); симметрия фигуры относительно секущей V = 0. Одинаковы по виду и симметричны относительно оси delta_t диаграммы рассеяния интервалов (delta_ti, Vi), являющиеся проекцией гистограммы на плоскость (delta_ti; Vi), см. рис. 3, д.

Множество точек диаграммы рассеяния ограничено сверху и снизу монотонно убывающими гиперболическими функциями. Точный вид распределения, наиболее подходящего для аппроксимации гистограмм скоростей для фиксированных длительностей delta_ti – первое предельное распределение экстремальных значений [6].

Для аппроксимации зависимости, выявленной для здоровых обследуемых, в настоящее время используется функция вида

μ(V, Δt) = k1∙exp(–k2∙(V∙delta_t)2)/delta_t,

где k1 и k2 определяются индивидуально.

Пример поверхностей гистограммы, полученной при тестировании одного обследуемого, и аппроксимирующей функции приведены на рис. 4.

pic_31.tifа pic_32.tifб

Рис. 4. Поверхности: а – гистограммы частоты появления интервалов движения с постоянной скоростью у здорового обследуемого; б – аппроксимирующей функции

pic_33.tif аpic_34.tifбpic_35.tifв
pic_36.tif гpic_37.tif дpic_38.tif е

Рис. 5. Частота появления интервалов движения с постоянной скоростью во фронтальной плоскости для обследуемых с различными диагнозами: а – сенситивная атаксия; б – ДЦП, спастическая диплегия; в – ДЦП, спастический гемипарез; г – деформирующий остеоартроз коленных и тазобедренных суставов; д – посттравматический остеоартроз таранно-пяточного сустава; е – идиопатический сколиоз II степени

Предлагается использовать найденную и параметризованную для обследуемых без выявленных патологий функциональную зависимость в качестве феноменологической модели идеального поддержания равновесия при отсутствии целенаправленного движения. Данная модель может быть использована как в качестве эталона для сравнения с результатами тестирования пациентов с нарушенными движениями, что снимает неопределенность понятия «нормы», так и для верификации концептуальных моделей [1, 2, 3, 4, 8, 11].

Визуально общий вид формы огибающей поверхности сохраняется и для пациентов со всеми рассмотренными диагнозами при различной тяжести проявления функциональной патологии (рис. 5), вплоть до 2-х предельных случаев с точки зрения удержания равновесия (рис. 5, а, б). Пациенты (выраженная сенситивная атаксия и гемиплегическая форма ДЦП) с трудом могли стоять в течение 20 секунд. Фактически для них ни разу не удалось при одном обследовании повторить тест, хотя стандартное обследование подразумевает трехкратное прохождение данного теста с интервалом от 3 до 30 минут.

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

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

pic_39.tif аpic_40.tifб
pic_41.tif вpic_42.tif г

Рис. 6. Изменение характеристик процесса удержания равновесия: а – частичное восстановление симметрии диапазонов скоростей относительно оси delta_t после 4 операций эндопротезирования коленных и тазобедренных суставов; б – частичное восстановление количества длительностей интервалов у пациента З. с гемиплегической формой ДЦП после периартериальной криосимпатодеструкции; в и г – сечения Vi = 0,02 до и после периартериальной криосимпатодеструкции пациента М. с ДЦП (спастическая диплегия). Красный цвет – состояние до лечения, зеленый – после

Таким образом, появляется возможность детализации имеющих место двигательных расстройств, связанной с их декомпозицией на функциональную и/или анатомическую патологию костно-суставной, мышечной и нервной подсистем. Анатомическим и обуславливаемым ими функциональным нарушениям костно-суставного аппарата соответствует асимметрия относительно нулевых значений середин диапазонов изменения скоростей для интервалов определенных длительностей (рис. 6, а) или характерные нарушения формы трехмерных гистограмм в целом (рис. 5, г, д).

Функциональным нарушениям в нервной системе соответствуют изменения формы гистограмм, связанные с длительностями интервалов, причем уменьшение количества длительностей интервалов свидетельствует о функциональной недостаточности нервной системы в целом (рис. 5, в; 6, б), а отсутствие интервалов определенной длительности, либо нарушение колоколообразной формы гистограмм скоростей для них свидетельствует об отсутствии рефлекторного кольца с соответствующим периодом или их параметрической патологии (рис. 5, а, б; 6, в, г). Функциональным нарушениям в мышечном аппарате соответствуют изменения диапазонов разброса значений скоростей как для всех интервалов, так и для отдельных длительностей интервалов, причем уменьшение диапазонов для всех интервалов относительно нормы свидетельствует о мышечном дефиците.

Заключение

Форма поверхности частоты появления интервалов движения с постоянной скоростью, за которой в процессе исследований закрепилось название – «нормалоид», в настоящее время параметрически описывается набором из 93 показателей, 11 показателей для каждого из 8 сечений по оси delta_t (delta_t = 0,02; 0,04; …; 0,16 с) и 5 показателей фиксирующих количества интервалов различных типов для стабилограммы (фронтальной либо сагиттальной) в целом. Используемый набор показателей обеспечил решение всех задач, поставленных в процессе клинической апробации (например: диагностика начальных стадий идиопатического сколиоза; дифференциация и оценка тяжести функциональных нарушений при ДЦП; контроль лечебного процесса при эндопротезировании тазобедренных и коленных суставов; оценка результатов консервативного и хирургического лечения при посттравматическом остеоартрозе таранно-пяточного сустава), различными методами параметрической обработки данных.

Опыт клинической апробации позволяет сделать вывод, связанный с особенностями использования профильными специалистами различных способов предъявления результатов компьютерной обработки, а именно – визуальное сравнение гистограмм, поверхностей и диаграмм рассеяния, не только наглядно и понятно, благодаря легкости интерпретации (скорости и их изменения – мышцы; локальные дефекты и асимметрии сечений формы – костно-суставный аппарат; количество длительностей интервалов и колоколообразная форма – нервная система), но и позволяет учитывать мелкие структурные особенности поверхности, характеризующие индивидуальную специфику пациента и заболевания, которые не могут быть отображены параметрически в рамках известных на сегодня математических методов обработки данных. Появилась возможность получить статическое, в виде одного трехмерного изображения, представление динамики процесса удержания равновесия, а в случае обработки активных стабилометрических тестов, подразумевающих целенаправленное движение, и процесса движения в целом. Проводя аналогию с рентгенографией, можно сказать, что форма нормалоида играет роль «скелета» процесса удержания равновесия.

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

Рецензенты:

Герасимов А.А., д.м.н., профессор, зав. кафедрой медицины катастроф, ГБОУ ВПО «Уральская государственная медицинская академия» Минздрава РФ, г. Екатеринбург;

Борзунов И.В., д.м.н., зам. декана лечебно-профилактического факультета, ГБОУ ВПО «Уральская государственная медицинская академия» Минздрава РФ, г. Екатеринбург;

Пучиньян Д.М., д.м.н., профессор, зам. директора по науке, ФГБУ «Саратовский научно-исследовательский институт травматологии и ортопедии» Министерства здравоохранения и социального развития Российской Федерации, г. Саратов.

Работа поступила в редакцию 15.07.2014. 


Библиографическая ссылка

Давыдов О.Д., Монтиле А.И., Марчук Ю.В., Монтиле А.А. СТАБИЛОМЕТРИЧЕСКИЙ ИНВАРИАНТ УДЕРЖАНИЯ РАВНОВЕСИЯ В НОРМЕ И ПАТОЛОГИИ // Фундаментальные исследования. – 2014. – № 7-5. – С. 931-939;
URL: https://fundamental-research.ru/ru/article/view?id=34803 (дата обращения: 28.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674