I. БИОМЕДИЦИНСКИЕ НАУКИ И СИНЕРГЕТИКА
DOI: 10.12737/5514
ИДЕНТИФИКАЦИЯ ПАРАМЕТРОВ ПОРЯДКА В ПСИХОФИЗИОЛОГИИ
М.А. ФИЛАТОВ, Д.Ю. ФИЛАТОВА, Д А. СИДОРКИНА, С М. НЕХАЙЧИК
ГБОУ ВПО «Сургутский государственный университет ХМАО-Югры», пр. Ленина, д.1,
г. Сургут, Россия, 628412
Аннотация. При использовании набора диагностических признаков x;, которые представляют общий вектор состояния системы x=x(t), возникает задача выбора из общего числа m таких признаков в наиболее значимых числом l<<m. На сегодня общее решение таких задач отсутствуют и в работе представлен метод ее решения для оценки различий в психофизиологических состояниях двух больших групп испытуемых (находящихся в условиях изменения факторов экосреды) с использованием нейроэмуляторов. При решении задачи бинарной классификации на основе анализа сенсомоторных параметров двух обследуемых групп предложена процедура, обеспечивающая реальное ранжирование x; и выделение из них наиболее значимых. Делаются аналогии в такой работе нейроэмулятора с работой нейросетей гиппокампа при эвристической деятельности мозга.
Ключевые слова: биологические динамические системы, фазовое пространство состояний, теория хаоса - самоорганизации.
ORDER PARAMETER IDENTIFICATION IN PSYCHOPHYSIOLOGY
M.A. FILATOV, D.YU. FILATOVA, D.A. SIDORKINA, S.M. NEKHAICHIK
Surgut State University, Lenina, 1, Surgut, Russia, 628412
Abstract. While using a set of diagnostic characters x;, which form a total system state vector x=x(t), there emerges problem of selections such characters with more significant number l<<m from the total number m. Today such problems are still not solved, therefore, the paper presents method of its solving for assessment of distinctions in psychophysiological state of two big groups of subjects (under changing ecological factors) using neuroemulators. While binary classification problem solving based on analysis of sensomotor parameters of two groups the procedure, providing real ranking x; and choosing more significant l, is produced. The work of neuroemulator is similar to work of neural networks of hippocampus while heuristic brain activity.
Key words: biological dynamic systems, phase space of states, theory of chaos and self-organization.
Общеизвестно, что различные внешние воздействия, в том числе и действие экологических факторов внешней среды, способны существенно изменят психофизиологический статус испытуемого. Более того, в последние годы мы говорим о состоянии экологического стресса у жителей северных регионов РФ, которые находятся в
особых условиях климата, продолжительности светового дня и в гипокинезии. Попытки изменения негативного психофизиологического состояния жителей Севера РФ, путем их перемещения в лечебно-оздоровительные зоны РФ, приводят к определенным изменениям в параметрах экологического стресса. Однако качественная и количественная
оценка эффективности подобных оздоровительных мероприятий наталкивается на определенные трудности. Главная из них - невозможность применения статистических методов для такой оценки.
Очень часто эта ситуация демонстрируется незначительными статистическими различиями в параметрах вектора состояния психофизиологических функций человека (ВСПФЧ) х=х(^=(хь Х2,...,Хщ)Т до и после проведения лечебно-оздоровительных мероприятий или при смене экофакторов среды (например, при переезде с севера на юг и обратно). Для объективного установления таких различий нами сейчас разработаны два новых подхода: на базе теории хаоса-самоорганизации (ТХС) (путем анализа параметров квазиаттракторов) и на базе нейроэмуляторов. Последняя и рассматривается в настоящем сообщении.
Отметим, что идентификация параметров порядка (наиболее важных диагностических признаков из всех т, которые изучаются) - основа системного синтеза, а значит и одна из главных проблем естествознания вообще. Один из ответов на этот фундаментальный вопрос естествознания мы предлагаем в рамках решения задачи бинарной классификации с помощью нейроэмулятора (ней-ро-ЭВМ или НЭВМ), как некоторой модели принятия правильного решения мозгом человека, его нейронными сетями. При этом для нейроэмулятора создавались определенные условия в виде многократного представления на вход НЭВМ двух одинаковых обучающих выборок, которые нейросеть устойчиво разделяла на два разных объекта, т.е. выполняла процедуры бинарной классификации. Однако, сама нейросеть каждый раз (в исходном состоянии, до начала деятельности) получала хаотичный набор весов признаков х; (перед настройкой) для конкретного вектора состояний х=х(1) исходной (исследуемой) биологической системы. В общем случае в качестве последней брали группы больных и здоровых людей, или группы больных одинаковым заболеванием, но для которых применялось два вида лечения, или бралась одна экспериментальная группа до воздействия и после воздействия. Всего было обследовано более 20-ти тысяч человек в виде наблюде-
ний и экспериментов в рамках задачи бинарной классификации, из них более 8-ми тысяч - по параметрам психофизиологических функций. Результаты во всех случаях были получены однотипные в отношении общих закономерностей. Сами эти закономерности в работе НЭВМ при идентификации параметров порядка мы представляем в нашем сообщении.
Применение НЭВМ для ранжирования значимости диагностических признаков, т.е. для понимания размерности фазового пространства с т до 1, где 1<<т, используется в физиологии, психологии, медицине уже на протяжении последних 3035 лет. Однако, однократное решение задачи бинарной классификации представляет только одну возможную реализацию весов признаков (компонент х, вектора состояния системы). Разовая диагностика может быть ошибочной (и обычно она ошибочна!). Поэтому мы предлагаем особую процедуру такого ранжирования, которая подобна работе мозга человека [7,8].
Общеизвестно, что любая когнитивная деятельность человека сопровождается многократными реверберациями в нейронных сетях гиппокампа. Возникает закономерный вопрос: зачем нужны мозгу эти реверберации и каково их значение при наших попытках градации умного человека (с хорошей памятью и наличием признаков эвристической деятельности) от глупого человека (у него эти характеристики отсутствуют). При этом, в различных многочисленных исследованиях (более 20-ти тысяч испытуемых) мы доказали, что одинаковых состояний в эфферентных и афферентных системах человека не существует. Тремор, кардиоинтервалы, теппинг, миограммы и электроэнцефалограммы постоянно изменяют свои параметры. Тогда очевидно, что и входные потоки информации (биоэлектрические сигналы, поступающие в нейро-сети мозга) не будут одинаковы, они будут постоянно изменяться, оставаясь в пределах некоторых квазиаттракторов.
Фактически, нейросети мозга постоянно находятся в хаотическом режиме своих динамических характеристик (возбудимость, параметры биопотенциалов и т.д.), но
при этом из этого хаоса рождается порядок. Мы мыслим, принимаем решения, создаем цивилизацию. Однако главное во всем этом - это ответ на вопрос: как мы выбираем параметры порядка, как мы решаем что главное, а что - второстепенное, как из огромного числа х; (параметров внешней среды) выбрать их необходимое число 1? Для ответа на этот вопрос рассмотрим как работает аналог мозга - НЭВМ в режиме малого числа итераций /, когда I=1,2,...,Р и Р<<100.
1. Нейросеть при малом числе повторов работает с низкой эффективностью. Если использовать генератор чисел в рамках равномерного распределения на интервале (0,1) и для каждого конкретного значения весов wi признаков х;,, которые образуют т-мерное фазовое пространство состояний (г=1,2,...,т), задавать конкретные значения весов wio хаотическим образом (т.е. 0^0<1 из генератора равномерного распределения в интервале (0,1)), то после решения задачи бинарной классификации, на выходе нейроэмулятора, мы будем получать наборы Wij весовых значений каждого диагностического признака хг, где г - номер координаты, а,- - номер итераций (-=1,2,...,р). Эти наборы после каждой ,-й итерации будут различными и они образуют матрицу Ж= №£=1,..,т!=1,..,р. Если таких р будет невелико (р<100), то мы будем получать конкретные выборки этих весовых коэффициентов хг, (для каждого отдельного р, где I -номер серии итераций из числа р). Однако попытки выделения из них наиболее значимых признаков (т.е. ранжирование х) обречены на провал. Веса диагностических признаков Wij будут хаотически изменяться для разных серий итераций 1=1,2,.,.,г и мы не сможем выделить главные из них. При малых р<100 мы для каждого значения I будем получать свой набор ранжирования wij и свои значения средних весов <^п> для каждой выборки I. Задача идентификации параметров порядка не будет решена.
Состояние психофизиологических функций зависит от большого комплекса воздействий природного, техногенного и социально-психологического характера. Особенно в последнее десятилетие мы наблюдаем массовое возрастание информа-
ционного «загрязнения» и прессинга на сенсорные системы (интернет, СМИ) человека, которое по большей своей части содержит не ликвидную и не значимую информацию. Такое «загрязнение» вызывает сложные ответные реакции психофизиологических функций (изменения в поведении, памяти, пороге восприятия информации) человека. У современной молодежи наблюдается зависимость от потребности постоянного прослушивания музыки, просмотра телевизора, игры в телефоне и на планшете (выполнение домашнего задания «под телевизор», что непосредственно снижает усвоение и закрепление новой информации). Физиологические или патологические смещения функциональных систем организма (ФСО) особенно четко проявляются у растущего организма под воздействием факторов техногенного, экологического характера. В пубертатном периоде эмоциональная компонента играет ведущую роль в работе сенсорных и психофизиологических систем восприятия информации. Известно, что «информационными» структурами мозга является комплекс гиппокамп-миндалина-гипоталамус, которые выполняют существенную роль в организации восприятия, памяти и переработки информации, соответственно смещение движения вектора сенсорных и психофизиологических функций в т-мерном фазовом пространстве в условиях широтных перемещений, указывает на косвенное изменения в работе этих функций [2,3]. Таким образом, решающую роль в деятельности любой адаптирующейся системы играет способность человека корректировать свои реакции согласно изменениям воздействующего стимула на основе уже имеющегося опыта, а также накопления и хранения вновь поступающей информации. Именно поэтому приспособление (усвоение новой информации) возможно лишь при условии сохранности и четкой работы психофизиологических механизмов.
В рамках применения нейро-ЭВМ были поставлены задачи по идентификации различий в работе психофизиологических функций в разных экологических условиях групп детей (мальчики и девочки) до отъезда в санаторий «Юный Нефтяник» и после
приезда в г. Сургут. В табл. 1-2 представлены результаты 30-ти 100 и 1000-кратного обучения нейросети для определения ранговой значимости признаков (Р4-Р7).
интервал вариации весовых коэффициентов сужается. Следует отметить, что при 50-ти кратном, 100-кратном обучении и более, в результате обучения нейросеть будет выдавать тоже разные результаты ранговой Таблица 1
Усредненные значения отдельных координат весов признаков Wi вектора состояния системы (параметры Wj - показатели психофизиологических функций детей (Р4-Р7), отъезд из Сургута на юг) полученные в процессе идентификации параметров порядка нейроэмулятором после р > 1000 итераций (настроек ЭВМ) в режиме бинарной классификации
Нейросети c p < 400=4x100
Расчеты итераций по выборкам (р > 1000) Средние значения весов признаков для координат вектора состояния системы х1 по наибольшим и наименьшим весам
ZK=<w1> CHIS=<w2> OP=<w3> SIMB=<w4>
р=4000 1=(1,...,100) 1 0,6009 0,5264 0,2524
р=1000 1=(0.....1000) 1 0,6122 0,5196 0,2526
р=1000 1=(1000,...,2000) 1 0,6007 0,5245 0,2461
р=1000 1=(2000,...,3000) 1 0,5930 0,5329 0,2527
р=1000 |=(3000,...,4000) 1 0,5976 0,5288 0,2579
Интервалы изменений <Д^1> 0 0,0192 0,0133 0,0118
Таблица 2
Усредненные значения отдельных координат весов признаков wj вектора состояния системы (параметры wj - показатели психофизиологических функций детей (Р4-Р7), приезде в Сургут из санатория) полученные в процессе идентификации параметров порядка нейроэмулятором после р > 1000 итераций (настроек ЭВМ) в режиме бинарной классификации
Нейросети c p < 400=4x100
Расчеты итераций по выборкам (р > 1000) Средние значения весов признаков для координат вектора состояния системы х1 по наибольшим и наименьшим весам
ZK=<wi> CHIS=<w2> OP=<ws> SIMB=<w4>
р=4000 1=(1.....100) 1 0,4384 0,3343 0,2891
р=1000 )=(0,.,1000) 1 0,4474 0,3296 0,2889
р=1000 )=(1000,.,2000) 1 0,4343 0,3339 0,2919
р=1000 |=(2000,...,3000) 1 0,4366 0,3374 0,2828
р=1000 )=(3000,.,4000) 1 0,4351 0,3365 0,2927
Интервалы изменений <Д^1> 0 0,0131 0,0078 0,0099
Как видно из табл. 1-2 с увеличением количества итераций обучений нейросети
значимости, т.е. вклад каждого признака в процедуру идентификации важных диагностических признаков при каждой итерации (настройке нейроэмулятора) будет различным. Однако, при увеличении числа повторов итераций p, т.е. повторов решения задачи бинарной классификации в рамках исходного задания хаотичного набора весов признаков хю, и при переходе к 102<p<103, картина начинает существенно изменяться.Веса Wj признаков Xj продолжают демонстрировать хаотические вариации (их законы распределения не совпадают!), но эти хаотические изменения выходных значений весов xj (после разделения двух групп) проявляются в пределах некоторого квазиаттрактора для всех компонент вектора х. Иными словами возникает некоторая устойчивая первая цифра, а последующие две цифры (для p>100) могут изменяться от минимума до максимума, образуя интервал изменения tswj=wj max-wj mjn для каждого набора итераций l, где l=1,2,...,N. Характер их изменения подобен изменению весов признаков при p<100. Это значит, что отдельные наборы (выборки весов признаков w\j, где j - номер порядка итерации, а l - номер набора (выборки) итераций) демонстрируют ненормальное распределение (хотя нормальные законы распределения тоже возможны). Возникающие непараметрические распределения при этом постоянно изменяются (при переходе отp1 (l=1) к другому p2 (l=2). Для каждой выборки pl, где l - номер серии повторов настройки нейроэмулятора для задачи
бинарной классификации, мы будем иметь свои функции распределения /(р¡), то есть при N повторах (¡=1,...,Ы), мы можем получить N различных функций распределения.
Ниже мы представляем объединенные результаты применения новых биофизических методов в рамках ТХС и нейро-ЭВМ в сравнительном аспекте для психофизиологических функций в 4-х мерном фазовом пространстве состояний (ФПС). Как показано выше, на основе полученных экспериментальных данных было доказано, что 30-ти кратная итерация нейрокомпьютерного анализа не позволяет идентифицировать параметры порядка (наиболее важные диагностические признаки), а новые биофизические методы в рамках ТХС (метод исключения отдельных признаков) обеспечивают идентификацию параметров порядка, предложен метод одиночного исключения для эффективного использования нейроэмулятора.
Рис. 1. Средние значения весов диагностических признаков нейронной сети (30-ти кратное обучение) параметров психофизиологических (Р4-Р7) функций учащихся в условиях широтных перемещений (столбцы соответствуют порядку этапов тестирования
1.2 п
и
0,711
0,6 0,64 0,6
0,52
0.43
0,58
0,45 I— 0,33
0,55
0,25 0,29 0,28
□ 0,25 0,28
- И.ГГ и.
Р4
Р5
Р6
Р7
Рис. 2. Средние значения весов диагностических признаков нейронной сети (4000 обучений нейросети) параметров психофизиологических (Р4-Р7) функций учащихся в условиях широтных перемещений (столбцы соответствуют порядку этапов тестирования
Наиболее значимый весовой коэффициент идентифицированный в нейро-ЭВМ до отъезда в санаторий является Р4, однако на четвертом этапе при 30-ти кратном обучении происходит смещение значимости и параметром порядка является Р5 (при р=1000 обучении этого не происходит - рис. 1 и 2. Оба теста направлены на концентрацию внимания, но в пятом тесте добавляется новая задача в распознавании четных и нечетных чисел, соответственно для решения поставленной задачи задействованы память, мышление. Как видно из рис. 1 и 2, инверсия весовых коэффициентов при увеличении количества итераций обучения нейросети происходит для параметра Р4 во втором и четвертом измерении, для параметра Р5 в 2-4 измерении; для параметра Р6 в 1-м измерении, для параметра Р7 во всех измерениях веса коэффициентов значимости принципиально не изменились. Определение параметром порядка в нейро-ЭВМ - показателя Р5 говорит о том, что после приезда в г. Сургут у учащихся повысился контроль над произвольным вниманием и решение поставленной задачи сопровождается высокой устойчивостью внимания и внутренней мотивацией.
Рис. 3. Результаты 30-ти итераций (крупный пунктир), 1000 итераций (мелкий пунктир) усредненные значения)) нейронной сети и результаты системного синтеза (идентификация параметров порядка, сплошная линия) параметров психофизиологических функций учащихся в условиях широтных перемещений (до отъезда в санаторий). Примечание: Р4 - ЛП реакции на появление зеленого/красного квадрата в произвольном месте экрана, (с); Р5 - ЛП распознавания чётных и нечётных чисел, (с); Р6 - ЛП распознавания символа, (с); Р7 - сосредоточенность внимания, (с)
ол -
о -I-,-,-,-
Р-1 PÏ Р6 Р7
Рис. 4. Результаты 30-ти обучений (усредненные значения) нейронной сети (сплошная линия) и результаты системного синтеза (идентификация параметров порядка) параметров психофизиологических функций учащихся в условиях широтных перемещений (приезд в г. Сургут)
Сравнительный анализ весовых коэффициентов и результатов идентификации параметров порядка с помощью новых методов на базе ТХС в условиях широтных перемещений представлены на рис. 1 и 2.
На рис. 3 и 4 представлены усредненные весовые коэффициенты полученные в результате 30-ти кратного обучения и 1000-кратного. В сравнении с результатами применения методов ТХС в этом случае мы видим диаметрально противоположные результаты. Однако в случае исключения признаков параметр P6 изначально был идентифицирован как параметр определяющий ВСС.
Таким образом, полученная детальная обработка числовых массивов новыми биофизическими методами позволила получить точные количественные данные о динамике изменения психофизиологических функций человека в условиях широтных перемещений. Подобный подход целесообразно производить перед направлением в оздоровительный комплекс, чтобы иметь конкретные рекомендации для каждого ребенка с целью проведения специфических и строго дозированных индивидуально физиотерапевтических, кинезиологических, психологических (разгрузочные тренинги) и других лечебных мероприятий в санатории.
Представлены результаты 30-ти кратного обучения нейросети для определения ранговой значимости признаков (P1-P7).
Для иллюстрации представим основные результаты 4-го блока исследований настоящей работы, в частности, пример ре-
шения одинаковой задачи бинарной классификации (например, различие параметров сенсомоторных и психофизиологических функций мальчиков и девочек в 1-м измерении), которая повторялась 30 раз, и 30 раз мы получали разные наборы весовых коэффициентов и разные наборы значимости xj.
Из приведенных выше примеров следует, что 30-ти кратная настройка нейросети не даёт ответ на вопрос об идентификации параметров порядка из всех Xj. В случае с нейроэмулятором для p<100 мы будем получать наборы выборок для каждого xj, но функции распределения (для каждого набора p1=p2=...=pN) будут получаться разными. При переходе к p>100 возникает устойчивость в значениях весов признаков xj в первой значащей цифре (после запятой) и далее, с ростом p (p>103, p>104), мы будем иметь хаотические вариации двух последних значащих цифр. Конкретный пример для p=100 мы представляем в табл. 1, где производится сравнение пяти выборок по p=100 итераций в каждой. При переходе к p>1000 мы получаем сходную устойчивость уже по первым двум значащим цифрам после запятой (табл. 2) для всех четырёх значимых цифр после запятой. Очень часто этими хаотическими изменениями в биологии пренебрегают, если мы оперируем с большими цифрами.
2. Нейро-ЭВМ идентифицирует параметры порядка при большом числе итераций. В качестве примера реализации таких процедур мы предлагаем результаты ранжирования весовых значений <wj> признаков xj для решения задачи разделения (бинарной классификации) группы испытуемых из 19-ти человек, которые проходили обследование по параметру сердечнососудистой системы (ССС) до работы в ночную смену и после такой работы (сходные данные были получены и для психофизиологических параметров человека). При многократных повторах процедуры бинарной классификации (разделения двух обучающих выборок, т.е. состояний xj для ССС перед ночной сменой и после смены) ней-росеть каждый раз четко обучается в этом разделении по 15-ти диагностическим признакам xj (î=l,2,3,...,15) указанных двух
групп (до и после смены). При этом каждый раз выдает свой уникальный (и неповторимый) набор весов признаков Xj, то есть Wj, где j - номер итерации (j=1,...,p), а i - номер признака (компартмента xi для ВСС). В табл. 1 мы представляем результаты расчёта четырёх выборок по p=100 итераций в каждой и их общую выборку (p=400) для наглядности динамики изменения <wil> при p=100.
Результаты одной серии из 50-ти повторов представлены на рис.1, для m=11 (были выбраны наиболее важные параметры из биологических признаков). Очевидно, что веса признаков (все xi) изменяются при каждой настройке и эти вариации (распределения каждого x;) НЭВМ не попадают под закон нормального распределения. Особым образом себя ведут и средние значения весов для разных наборов выборок, т.е. с увеличением l и p для всех признаков (на рис. 1 их число m=11) мы получаем разные законы распределения и для разных выборок из x; мы тоже имеем существенные различия. Подобную картину описала группа ученых Стэнфорда при анализе выполнения движения при якобы произвольном (целевом) движении (reaching) [1]. Любое движение человека повторить невозможно. Более того постуральный тремор - это тоже неповторимое (невоспроизводимое) движение, которое произвольное по сути (есть цель!), но оно выполняется непроизвольно, т. е. хаотически. Этот хаос начинается с афферентов, с хаоса в нейросетях мозга и заканчивается хаосом в эфферентных сигналах и в эффек-торном органе (в мышцах для тремора).
Более того, непараметрические распределения для каждой серии (из p=50) тоже получаются разными. Мы имеем просто некоторые наборы выборок, которые генерируются как бы из разных генеральных совокупностей. Подобную картину мы наблюдаем при анализе треморограмм (на интервале Т=5 сек каждый интервал т=1 сек, где Т]=т2=тз=Т4=тз=Т/5 генерирует свою выборку) и здесь статистические функции распределения не совпадают, они как бы получаются из разных генеральных совокупностей [4-6]. Аналогичную закономерность имеем и для теппинграмм, кардио-
грамм, электроэнцефалограмм, миограмм и даже из анализа биохимических параметров больных. При длительном мониторинге параметров крови имеются. В общем случае мы сейчас в рамках ТХС говорим о хаотической динамике всех параметров го-меостаза любого жителя планеты Земля. В целом, мы можем говорить, что при малых повторах неопределенность значимости хг велика, а для НЭВМ определить параметры порядка дляр< 100 невозможно!
3. Переход хаоса в порядок реализуется за счёт ревербераций в нейросетях мозга. Из предыдущего параграфа следует, что разовая настройка нейросети не даёт ответ на вопрос об идентификации параметров порядка из всех хг. При увеличении числа повторов итераций р, т. е. повторов решения задачи бинарной классификации в рамках исходного задания хаотичного набора весов признаков хг0, и при переходе к
2 3
10 <р<10 , картина начинает существенно изменяться. Веса wij признаков хг продолжают демонстрировать хаотические вариации (их законы распределения не совпадают!), но эти хаотические изменения выходных значений весов х; (после разделения двух групп) проявляются в пределах некоторого квазиаттрактора для всех компонент вектора х. Иными словами возникает некоторая устойчивая первая цифра, а последующие две цифры (для р>100) могут изменяться от минимума до максимума, образуя интервал изменения /^wi=wi maX-wi тгп для каждого набора итераций ¡, где ¡=1,2,...,N. Характер их изменения подобен изменению весов признаков при р<100. Это значит, что отдельные наборы (выборки весов признаков Wlij, где - - номер порядка итерации, а I
- номер набора (выборки) итераций) демонстрируют ненормальное распределение (хотя нормальные законы распределения тоже возможны). Возникающие непараметрические распределения при этом постоянно изменяются (при переходе от р] (¡=1) к другому р2 (¡=2). Для каждой выборки р1, где ¡ - номер серии повторов настройки нейроэмулятора для задачи бинарной классификации, мы будем иметь свои функции распределения /(р1), то есть при N повторах (¡=1,...,Щ, мы можем получить N различ-
ных функций распределения. Подобное мы имеем и для тремора (теппинга, кардиоин-тервалов и т.д.), если считать, что функция регуляции позы работает подобно нейро-эмулятору в режиме неопределенности (постоянного изменения весов синаптических связей между нейронами). Постоянная и хаотичная самоорганизация нейросетей мозга порождает хаотичную динамику теп-пинга, тремора, кардиоритма и других параметров гомеостаза. В этом смысле гомео-стаз отдельного организма подобен гомео-стазу экосистем (непрерывный хаос, но в пределах квазиаттрактора).
Действительно, именно такая картина наблюдается и при регистрации тремора. Если мы будем пытаться удержать палец с пластинкой на некотором расстоянии от датчика (рис. 2), то это будет (по факту наличия цели в таком движении) произвольное движение. Однако, это движение (удержание пальца в пространстве) будет реализо-вываться непроизвольно (т.е. хаотически). Каждая секунда тремора будет давать свою функцию распределения (могут возникать даже параметрические распределения!). Однако, повторов распределения мы произвольно получить не можем в этом простом опыте. Тремор будет генерировать хаотический набор функций распределения для каждого отдельного интервала измерения т, что подобно генерации нейроэмулятором весов признаков для набора итераций р1 (I -номер выборки пор итераций в каждой).
При этом, регистрация тремора будет демонстрировать существование некоторого квазиаттрактора в ФПС х=х(0=(х1,х2,х3)т, где х1 - реальная координата пальца по вертикали, х2=dx1/dt - скорость перемещения и xз=dx2/dt - ускорение. Движение пальца будет происходить с точностью до метра или даже сантиметра без изменений, а вот миллиметры и более мелкие единицы измерения будут непрерывно и хаотически изменяться. При этом для подобных сложных систем под хаотической динамикой их поведения мы понимаем не положительные константы Ляпунова (расхождение двух фазовых траекторий) и не стремление автокорреляционной функции к нулю, а фундаментальное определение хаоса. Оно звучит так: задание началь-
ных параметров системы в виде начального значения вектора состояния системы х^0) не определяет дальнейшую траекторию и конечное состояние х(При этом х^0) не определяет дальнейшую траекторию её развития и её конечное состояние (х(к) для наших особых систем тем более, так как х(0 нельзя дважды повторить, любое состояние биосистемы изменяется непрерывно. Аналогия с тремором здесь весьма сильна, так как начальное положение пальца в пространстве мы не можем задать повторно и определенно (это хаотичное движение в квазиаттракторе) и поэтому любое х0 будет выбираться для рецепторов хаотически из некоторого квазиаттрактора. А это означает, что и веса синап-тических связей меняются непрерывно. Это легко наблюдать по нейрограммам афферен-тов и эффекторов (они весьма различны и непрерывно изменяются).
В случае с нейроэмулятором для р< 100 мы будем получать наборы выборок для каждого х1, но функции распределения (для каждого набора р1=р2=...=ры) будут получаться разными. При переходе к р>100 возникает устойчивость в значениях весов признаков х1 в первой значащей цифре (после запятой) и далее, с ростом р (р>103, р> 104), мы будем иметь хаотические вариации двух последних значащих цифр. Конкретный пример для р=100 мы представляем в табл. 1, где производится сравнение пяти выборок по р=100 итераций в каждой. При переходе к р>1000 мы получаем сходную устойчивость уже по первым двум значащим цифрам после запятой (табл. 2) для всех четырёх значимых цифр после запятой. Очень часто этими хаотическими изменениями в биологии пренебрегают, если мы оперируем с большими цифрами. Например для тремора, если мы регистрируем положение пальца с точностью до 0,5 см, хаотической динамикой можно пренебречь. Тогда можно утверждать, что человек держит произвольно свою конечность в данной точке пространства (но эта «точка» определяется с точностью до амплитуды тремора - 0,5 см). Аналогичные высказывания можно делать и в отношении кардиоитервалов, теппинга, электромио-грамм и т. д. Для обозначения границ всех
таких фазовых координат (а это может быть х1 - координата пальца в пространстве, х2=ёх1/ё1 и х3=ёх2/ё1;) мы вводим понятие квазиаттрактора, объемом Уе. Квазиаттрактор представляет некоторый объем фазового пространства, внутри которого непрерывно и хаотически движется ВСС х(1;).
Характерно, что в подобных нейросе-тевых системах не возникают флуктуации, которые закономерно возникают в стохастических системах. Например, при бросании монетки 1000 раз мы можем получить отклонение уже в первой значащей цифре, то есть вместо 0,5... можно получить 0,6... или 0,4... и при этом легко рассчитать вероятность такого отклонения. При этом, первая значащая цифра после запятой (в частоте события р (А)) будет всё-таки варьировать. В стохастике мы имеем флуктуации, которые затухают с увеличением К, но редкие (и большие!) флуктуации всё-таки возможны (хоть и с малой вероятностью). Для систем с хаосом и самоорганизацией (роль хаоса играет генератор равномерного распределения для х;о, а роль самоорганизации выполняет нейроэмулятор и нейронные сети мозга, которые контролируют тремор, теппинг, кар-диоинтервалы, электромиограммы и гомео-стаз человека в целом), картина получается другой. Все наблюдаемые распределения получаются или нормальными, или ненормальными и при малых выборках (числе итераций), т.е. N<100, мы будем иметь непрерывно изменяющиеся функции распределения, которые для каждой выборки итерации имеют свои значения. При бросании монетки для N<100 мы всегда (как и для других стохастических процессов) будем иметь определённые (и неизменные!) законы распределения и даже большие флуктуации их не могут существенно изменить.
Анализируя результаты, полученные на модели, которая демонстрирует переход хаоса в порядок можно сделать ряд важных выводов, которые следует учитывать всем, кто использует НЭВМ для диагностических целей в медицине, психологии и биологии.
Выводы:
1. Если число повторов итераций в работе нейросети будет невелико (например, у человека число ревербераций в гип-покампе будет невелико), то когнитивная
деятельность реализуется с низкой эффективностью, параметры порядка будут неправильно выбраны нейросетью. Это будет нейроэмуляторная модель некомпетентного человека, который не умеет находить параметры порядка. Все такие попытки (при малых N) обречены на ошибки и полученные ранее результаты следует пересмотреть, т.к. мы имеем каждый раз одну особую реализацию из многих возможных.
2. При увеличении числа итераций (N>100) мы повышаем точность идентификации параметров порядка и результаты когнитивной деятельности при многократных реверберациях нейросетей мозга могут быть весьма высокими. Такие нейросети способны к тонким оценкам.
3. Очевидно, что гениальные люди отличаются от обычных людей именно из-за числа повторов решения одной и той же задачи (анализа данных, процессов, явлений). Гений должен многократно возвращаться к решению возникающей задачи даже в режиме исходно хаотического состояния нейросетей мозга (воспроизводить задачу как новую каждый раз). Эвристическая деятельность мозга (как НЭВМ при идентификации параметров порядка) невозможна без многократных повторов решения задачи (при минимальной неопределенности в состоянии сети и при одинаковых обучающих выборках). В нашей модели отсутствует память, которая может резко снизить число необходимых итераций.
Литература
1. Churchland M.M., Cunningham J.P., Kaufman M.T. and others. Neural population dynamics during reaching // Nature.- 2012.-V.487.- P.51-56.
2. Eskov V.M. Modeling of the hierarchical respiratory neuron networks // Neurocomputing.- 1996.- V.11.- P. 203-226.
3. Eskov V.M., Eskov V.V., Filato-va O.E. Characteristic features of measurements and modeling for biosystems in phase spaces of states // Measurement Techniques (Medical and Biological Measurements).-2011.- V.53, №12.- P. 1404-1410.
4. Eskov V.M., Gavrilenko T.V., Kozlo-va V.V., Filatov M.A. Measurement of the dy-
namic parameters of microchaos in the behavior of living biosystems // Measurement Techniques.- 2012.- V.55, №9.- P.1096-1100.
5. Eskov V.M., Eskov V.V., Filato-va O.E., Filatov M.A. Two types of systems and three types of paradigms in systems philosophy and system science // Journal of Biomedical Science and Engineering.- 2012.- V.5, №10.- P. 602-607.
6. Haken H. Principles of brain functioning: a synergetic approach to brain activity, behavior and cognition (Springer series in synergetics). Springer, 1995.- 349 p.
7. Prigogine I. The Die Is Not Cast // Futures. Bulletin of the Word Futures Studies Federation.- 2000.- V.25, №4.- P.17-19.
References
1. Churchland MM, Cunningham JP, Kaufman MT and others. Neural population dynamics during reaching. Nature. 2012;487:51-6.
2. Eskov VM. Modeling of the hierarchical respiratory neuron networks. Neuro-computing.1996;11:203-26.
3. Eskov VM, Eskov VV, Filatova OE. Characteristic features of measurements and modeling for biosystems in phase spaces of states. Measurement Techniques (Medical and Biological Measurements).2011;53(12):1404-10.
4. Eskov VM, Gavrilenko TV, Kozlo-va VV, Filatov MA. Measurement of the dynamic parameters of microchaos in the behavior of living biosystems. Measurement Techniques. 2012;55(9):1096-100.
5. Eskov VM, Eskov VV, Filatova OE, Filatov MA. Two types of systems and three types of paradigms in systems philosophy and system science. Journal of Biomedical Science and Engineering. 2012;5(10):602-7.
6. Haken H. Principles of brain functioning: a synergetic approach to brain activity, behavior and cognition (Springer series in synergetics). Springer, 1995.
7. Prigogine I. The Die Is Not Cast. Futures. Bulletin of the Word Futures Studies Federation.2000;25(4):17-9.