CC BY
15
УДК519.711.3 : 519.722
О ДИФФЕРЕНЦИАЛЬНО-НЕАВТОНОМНОМ ПРЕДСТАВЛЕНИИ ИНТЕГРАТИВНОЙ АКТИВНОСТИ НЕЙРОПОПУЛЯЦИИ БИЛИНЕЙНОЙ МОДЕЛЬЮ ВТОРОГО ПОРЯДКА С ЗАПАЗДЫВАНИЕМ
© 2021 A.B. Данеев1, A.B. Лакеев2, В.А. Русанов2, П.А. Плеснёв1
1 Иркутский государственный университет путей сообщения, Иркутск, Россия 2 Институт динамики систем и теории управления имени В.М. Матросова СО РАН, Иркутск, Россия
Статья поступила в редакцию 13.02.2021
Для нейроморфных процессов, заданных поведением локальной нейропопуляции (например, процессами, индуцированными интерфейс-платформой "мозг-машина" типа Neuralink), исследуется разрешимость задачи существования дифференциальной реализации этих процессов в классе билинейных нестационарных обыкновенных дифференциальных уравнений второго порядка (с запаздыванием) в сепарабельном гильбертовом пространстве. Данная постановка относится к типу обратных задач для аддитивной комбинации нестационарных линейных и билинейных операторов эволюционных уравнений в бесконечномерном гильбертовом пространстве. Метаязыком развиваемой теории служат конструкции тензорных произведений гильбертовых пространств, структуры решеток с ортодополнением, функциональный аппарат нелинейного оператора Релея-Ритца и принцип максимума энтропии. При этом показано, что свойство сублинейности этого оператора, позволяет получить условия для существования таких дифференциальных реализаций; попутно обосновываются метрические условия непрерывности проективизации данного оператора с вычислением фундаментальной группы его компактного образа.
Ключевые слова: обратные задачи нелинейной нейроморфной динамики, тензорный анализ, билинейная неавтономная дифференциальная реализация второго порядка с запаздыванием, энтропийный оператор Релея-Ритца. DOI: 10.37313/1990-5378-2021-23-2-115-126
Работа выполнена при финансовой поддержке РФФИ (проект № 19-01-00301).
ВВЕДЕНИЕ
Многомерные нейродинамические модели эволюции состояния синапсов, изменения скорости нейроимпульсов и т.п. в локальных популяциях нейронов [1] могут иметь форму систем интегро-дифференциальных уравнений [2]. Для данного дифференциального представления нейроморфных процессов актуальна задача разрешимости таких бесконечных систем. При этом используются различные техники современного функционального анализа, включая теоремы сравнения, монотонные итерационные методы, метод усечения, а также топологические методы фиксированной точки. Другой нейродинамический подход описывает нейропроцессы в терминах нейросетевой протомодели "черного ящика" с дальнейшей реализацией [3, с. 21] этих процессов в классе бесконечномерных дифференциальных моделей [4] с запаздыванием и без. В целом можно заключить, что нелинейные неавтономные дифференциальные уравнения в бесконечномерных гильбертовых пространствах весьма продуктивны в теоретико-прикладных задачах математической нейробиологии. При этом факторы нестационарности, полилинейности и запаздывания таких дифференциальных моделей, в том числе высших порядков [5], отражают интегративный процесс перестройки физико-химических взаимосвязей в синергетических комплексах синапсов/рецепторов [6] исследуемых (моделируемых) нейропопуляций.
В контексте второго подхода отметим, что основоположник качественной теории дифференциальной реализации (КТДР) Р. Калман, констатируя, что в теории динамических систем задача реализации играет центральную роль [3, с. 267], сформулировал следующий подход [3, с. 286]: рассматривать задачу реализации как попытку угадать уравнения движения динамической системы по поведению ее входных и выходных сигналов или как построение физической модели, объясняющей экспериментальные данные типа "вход-выход" (протомодель "черного ящика"). Настоящий период развития КТДР в бесконечномерной постановке связан с созданием специального математического метаязыка - энтропийной теории расширений M2- операторов [4]. Данная теория существенным образом пере-
Данеев Алексей Васильевич, доктор технических наук, профессор, профессор кафедры «Информационные системы и защита информации». Е-mail: [email protected] Лакеев Анатолий Валентинович, доктор физико-математических наук, ведущий научный сотрудник. Е-mail: [email protected]
Русанов Вячеслав Анатольевич, доктор физико-математических наук, доцент, старший научный сотрудник. Е-mail: [email protected] Плеснёв Павел Александрович, аспирант. Е-mail: [email protected]
строила и укрепила теоретико-системные основания КТДР и обеспечила гармоничную связь геометрических идей M2-продолжимости [4] с методами теории апостериорного математического моделирования дифференциальных уравнений высших порядков [5] в бесконечномерных пространствах, причем делая упор на приложения [2, 6], а не на достижение максимальной общности изложения КТДР-результатов.
Последнее обстоятельство не могло не наложить отпечаток на содержание данной работы. А именно, так как во многих задачах математической нейрофизиологии модель реализации дифференциального представления нейроморфных процессов типа "ответная реакция синапсов на раздражающее воздействие рецепторов" требует учитывать нелинейную взаимосвязь как от самой ответной реакции и скорости ее протекания, так и от раздражающего воздействия, то ниже основное внимание сосредоточено на исследовании модели дифференциальной реализации, зависящей от пяти нестационарных билинейных структур. Первая из них зависит от реакции (как вектор-функции времени), второй билинейный оператор зависит от этой реакции и скорости ее протекания, третий билинейный оператор зависит только от скорости протекания данной реакции и два других учитывают эти переменные в связи с влиянием на них фактора запаздывания нейроимпульсов синапсов и вектор-функции нейросигналов возбуждаемых рецепторов.
Кроме того, КТДР-методы, рассматриваемые в духе бесконечномерной постановки обратных задач математической физики [4], сложнее, интереснее, глубиннее в приложениях и очень важны для понимания основных репрезентативных свойств самих дифференциальных моделей. Их конструкции могут служить отправными точками современного развития общей (аксиоматической) теории динамических систем (в духе [3]), попутно создавая этим конструкциям репутацию весьма полезного математического инструмента в прецизионном апостериорном нелинейном моделировании сложных бесконечномерных динамических моделей.
1. ТЕРМИНОЛОГИЧЕСКИЙ СЛОВАРЬ И ПОСТАНОВКА ЗАДАЧИ
Далее (X,|| ■ ||X ), (Y,|| • ||Y ) - вещественные сепарабельные гильбертовы пространства (т.е. нормы удовлетворяют "условию параллелограмма" [7, с. 162]); при этом ниже используем [7, с. 176] линейную изометрию (сохраняющую норму) E: Y ^ X пространств Y и X. Как обычно L(B\B") - банахово пространство (с операторной нормой) всех линейных непрерывных операторов для двух банаховых пространств B' и B", Z(X2,X) - пространство всех непрерывных билинейных отображений из декартового квадрата XXX в пространство X, ниже активно используем линейную изометрию [7, с. 650] пространств Z(X2,X) и L(X,L(X,X)).
Обозначим через T := [i0, ij отрезок числовой прямой R с мерой Лебега Ц, через - О -алгебру всех Ц-измеримых подмножеств из T, запись S с Q (S,Qер^) означает p,(S\Q) = 0.
mod ц
Сверх того примем, что ACl{T,X) - множество всех функций ф: T ^X, первая производная которых является абсолютно-непрерывной функцией (относительно меры Ц) на интервале T.
Если ниже (B, || ■ ||) - банахово пространство, то через Lp (T, B), p е [1,будем обозначать банахово пространство классов Ц-эквивалентности интегрируемых по Бохнеру [8] отображений f: T ^ B с нормой (j"T || f (т) ||p yi{dx))иp <«>, соответственно через L^ (T, B) - банахово пространство из этих классов с нормой ess supT || f ||< ^. В данном контексте условимся, что
L2 := L2(T, L(X, X)) X L2 (T, L(X, X)) x L2 (T, L(Y, X)) X X L2(T, Z( X2, X)) X L2(T, Z( X2, X)) X L2(T, Z( X2, X)) X L2(T, Z( X2, X)) X L2(T, Z( X2, X)), L := L(X, X) X L(X, X) X L(Y, X) X Z(X2, X) X Z(X2, X) X Z(X2, X) X Z(X2, X) X Z(X2, X).
Далее, считаем, что на временном интервале T фиксировано a posteriori поведение исследуемой (моделируемой) локальной нейропопуляции в виде не ограниченного по мощности нелинейного
пучка N наблюдаемых нейродинамических процессов в пространстве состояний1 X, представленных парами (x, ы) вектор-функций типа "реакция, воздействие", т.е. нейроморфной бихевиористи-ческой системой в виде N -"черного ящика", или формально:
N с {(x, ы): x е AC\T, X), ы е L2(T,Y)}, Card N < exp Х0,
где (в терминологии [3]) (x, ы) - пара "траектория (состояние синапсов), управление (воздействие на
рецепторы)", Х0 - алеф нуль, ехр К0 - континуум; термин "нелинейный пучок" означает, что для
траекторий пучка a priori не предполагается наличие принципа суперпозиции, - когда зависимость выходных х-переменных от входных н-воздействий суть линейная [3].
Кроме того, примем, что дифференциальная реализация N -системы должна иметь второй порядок2, т.е. содержать член AH, где оператор-функция A: T ^L(X,X) обладает свойствами:3
A е L„ (T, L(X, X)),
^{t е T: A(t) = 0 е L(X,X)} = 0; при этом нарушение условия эквивалентности нормальной системе, а именно:
^{t е T: Ker A(t) Ф {0} с X} Ф 0, необременительно (допустимо); в данном контексте см. ниже замечание 1.
Рассмотрим задачу: для фиксированной пары (N, A) определить необходимые и достаточные условия, выраженные в терминах нелинейного пучка динамических процессов N и оператор-функции A, существования упорядоченного набора из восьми оператор-функций
(Aj,4,,B, Dj,D2,D3,D4,D5)е L2,
для которого осуществима билинейная дифференциальная реализация (БДР) второго порядка с запаздыванием X = const > 0, имеющая аналитическое представление вида:4
A x + Ax x + Aq x — = Bu + Dj(x, x) + D 2( x, x) + D3( x, x) + (1)
+ D4(E(u),y) + D5(E(u),y), V(x,u)е N,
f x(t - x), если tn + x < t < t;
t ^ y(t):= \ v h 0 _ (1*)
[0 е X, если t0 < t < 10+x.
Условимся, что если моделируемые операторы БДР-системы (1)-(1*) предполагается искать в классе стационарных, то будем их строить в классе непрерывных, и при этом писать
AЛ,B, Dl5D2,D3,D4,D5)е L*.
В связи с предложенной математической постановкой разрешимости задачи существования нейроморфной динамической модели отметим, что каждая область математики, как правило, содержит свои ведущие проблемы, которые настолько трудны, что их полное решение даже и не ожидается, но которые стимулируют постоянный поток работ и служат главными вехами на пути прогресса в этой области. В рамках КТДР-исследований такой проблемой является проблема классификации непрерывных бихевиористических систем (в том числе нейродинамических протомоделей типа "черный ящик"), рассматриваемых так, как если бы они точно (прецизионно) совпадали с решениями идеализированных дифференциальных моделей, в том числе высших порядков с фактором запаздывания. В наиболее сильной форме она предполагает модельную классификацию таких си-
1 Выбор сепарабельного гильбертова пространства в качестве пространства состояний уравнений динамики исследуемой нейропопуляции обеспечивает интегративный характер нейроморфного моделирования с большим числом переменных [23]. Так, для счетного числа нейронов их дифференциальные уравнения второго порядка в гильбертовых пространствахХ,- можно свести к одному уравнению (второго порядка) в расширенном гильбертовом пространствеX = ® X, г = 1, 2, ..., где © -прямая сумма пространств [17, с. 54]; см. также ниже пример 5.
2 Использование в нейродинамике дифференциальных моделей второго порядка позволяет моделировать в качестве эндогенного фактора (координат вектора состояния) изменение скорости исследуемых нейроимпульсов.
3 В случае отсутствия (в контексте математического моделирования нейродинамики) априорной информации об оператор-функцииЛ можно считать, что Л представляет оператор гомотетии с коэффициентом 1.
4 Равенство в дифференциальном уравнении (1) рассматривается как тождество в пространстве X).
стем с точностью до соответствующего класса дифференциальных реализаций, в частности, класса БДР-моделей (1)-(Г), в обосновании чего ниже служат КТДР-теоремы 1-2 и их следствия. Данная конструкция может стать отправной точкой в развитии геометрической теории нестационарных полилинейных векторных полей [9], параллельно развивая новый математический метаязык для прецизионного дифференциального моделирования нелинейных бесконечномерных динамических систем.
2. ХАРАКТЕРИСТИЧЕСКИЙ ПРИЗНАК БДР-МОДЕЛИ: ЭКВИВАЛЕНТНЫЕ ФОРМУЛИРОВКИ
Опишем теперь аналитическую схему решения вопроса о разрешимости (или неразрешимости) БДР-задачи (1)-(Г). Итак, пусть Z := X ® X - гильбертово тензорное произведение [10, с. 54] гильбертовых пространств X и X с кросс-нормой || , определяемой скалярным произведением [8] в
X . Сверх того, введем новые обозначения:
U := X х X xY х Z х Z х Z х Z х Z,
L2 :=L2 (Г,AtX^a(T,L(X,X))xL2(T,L(Y + X))х х L2(T, L(Z, X)) х L2(T, L(Z, X)) х L2(T, L(Z, X)) х L2(T, L(Z, X)) х L2 (T, L(Z, X)); ясно, что функциональное пространство L2 (с топологией произведения) линейно гомеоморфно банахову пространству L2(T,L(U,X)) .
Обозначим через П - универсальное билинейное отображение я:XxX^X®X; на языке категорий морфизм П определяет тензорное произведение как универсальный отталкивающий объект [10, с. 40]. Универсальность билинейного отображения П состоит также в том, что
я: X х X ^ X ® X,
(х15х2) ^ я(х15х2) = Xj ®х2, ||xj ®х2||Z =||xJJ|х2|X ; данные соотношения важны для аналитического представления нелинейного функционального оператора Релея-Ритца в части конкретизации (в его конструкции) нормы Ци .
Далее считаем, что декартов квадрат X2 = X х X наделен нормой (||-||x +||j|x)12 • В данной постановке я е Х(X2,Z) и, с учетом теоремы 2 [7, с. 245], для любого билинейного отображения D е Х(X2,X) всегда найдется линейный непрерывный оператор D е L(Z, X) такой, что D = D о п, при этом для любой пары (x,u) е N будут выполняться (согласно (1*)) включения: n(x,x), п(x,X), п(Х,X)е L^(T,Z), n(E(u),y), n(E(u),y))е L2(T,Z). Данные построения подытоживает следующее утверждение. Лемма 1. Для любого набора (Al,Л^,B, D15D2,D3,D4,D5)е L2 и отображения
F: L2(T, X) х L2(T, X) х L2(T,Y) х х L2(T, X2) х L2(T, X2) х L2(T, X2) х L2(T, X2) х L2(T, X2) ^ L^T, X),
(У1, y2, Уз, y4> y5> Уб> Ут y8) ^ F (У1» y2> Уз, y4> y5> Уб> Ут y8):=
= A,y, + A y2 + By3 + Dy + D2 y5 + D3 y6 + D4 y7 + B5 y8
существуют единственный кортеж оператор-функций
(a, d2, D3, a, ds, а, а, а) е l2
и, соответственно, единственное линейное отображение
M :L2(T,U) ^ Lj(T,X), имеющее аналитическое представление вида
(Z1 >Z 2' Z3> Z4 ' Z5 ' Z6 ' ZT Z8 ) ^ M (Z1'Z 2' Z3> Z4 ' Z5 ' Z6 ' ZT Z8 ) •= = D1Z1 + A Z2 + D3 Z3 + D4 Z4 + D5 Z5 + D6 Z6 + D7 Z7 + D8 Z8 >
такое, что выполняется следующее функциональное равенство:
(У1, У2> Уз» Ул> ys> Уб» y7. y8) ^ F (У1' У 2' Уз. Ул> ys> Уб» y7> y8) = =M (У1, У2» УзЖ УЛЯ( У5)»я( Уб)»я( У7)»я( У8))>
которое, в свою очередь, индуцирует для оператор-функций из конструкций отображений F и M следующие операторные соотношения:
A1 = Д, A0 = D2, B = D3, Dj = D4 о D2 = D5 оn,
D3 = D6 о к, D4 = D7 О D5 = D8 О л .
Везде далее (в контексте БДР-задачи (1)-(1*)) примем
VN := Span{(x,x,u,n(x,x), n(x,x), n(x,x), n(E(u),y),n(E(u),y))е L2(T,U): (x,u)е N}.
Следующая лемма обобщает бихевиористическое условие (7) [11].
Лемма 2. Пусть
S := {t е T: (g(t), w(t), v(t), q(t), s(t), h(t),i?CtX"(t)) = 0 е U}, Q := {t е T: g(t) = 0 е X}, где (g,w,v,q,s,h,u,u)е VN. Тогда имеет место S с Q.
mod ц
Далее, пусть (L + (T,R),<L) - положительный конус [12] классов Ц-эквивалентности всех вещественных неотрицательных Ц-измеримых на интервале T функций с квазиупорядочением <L, при котором <L в том и только в том случае, если ^'(t) < ^"(0 Ц-почти всюду в T. При этом для заданного подмножества WсL+ (Г,R через supL W обозначим наименьшую верхнюю грань подмножества W, если эта грань существует в конусе L + (T, R) в структуре квазиупорядочения <L, в частности, легко установить, что имеет место соотношение
В данной постановке рассмотрим решетку с ортодополнением [17, с. 339] вида:
R(W):= е L + (T,R):§ <L supL W}.
Тогда (R(W),<L) - решетка с наименьшим %0е L + (T,R) и наибольшим supL Wе L + (T,R) элементами; здесь и далее X0 - "нуль-функция" конуса L + (T, R). В контексте теоремы 17 [7, с. 68] и следствия 1 [7, с. 69] несложно извлечь более общее утверждение; ниже infL - наибольшая нижняя <L -грань.
Лемма 3. Решетка R(W) - полная, т.е.
infL V, supLV е R(W) VV с R (W).
Пусть ¥: VN ^ L + (T, R) - энтропийный оператор Релея-Ритца [4]:
t » ¥(<p)(t):=
A(t) g (t) /||Ф<7 | „, если 9(t ) Ф 0 eU ;
0 e R, если ç(t ) = 0 eU ; где ф := ( g, w, v, q, s, h, û, u) eVN. Ясно, что имеет место равенство
||( g (t ), w(t ), v(t ), q(t ), s(t), h(t ), ïîtt ), û(t ))|| U :=
=4 g (t t +1 |w(t t -HKtf-+И1; +1И1; +| h(t )||; +|;(о||; y»,
при этом функции
( g, w, v) м- g, w, v) :=|| g, w, v,0,...,0 II U, (q,s,h,û,û) ^ ^(q,s,h,û,û) :=|| 0,0,0,q,s,h,û,û ||U, назовем, соответственно, линейной и билинейной характеристиками оператора Релея-Ритца. В силу леммы 2 на интервале времени T выполняется
supp ¥(ф) = supp A g (mod ц);
здесь в определении supp-конструкции носителя функции следуем [7, с. 137] (т.е. носитель определяется с точностью до множества меры нуль).
Энтропийный оператор ¥ удовлетворяет весьма простым (но важным) соотношениям
X 0 -l ¥(<р) = ¥(/ ф), / е R*:= R\{0}, Ф е Vn ; ниже будем различать в обозначениях образ точки ¥(ф) и образ множества ¥[{ф}].
Теория оператора Релея-Ритца нуждается в точном функционально-геометрическом языке, что заставляет уделять этому языку особое внимание. Поэтому прежде чем идти дальше, будет полезно ввести дополнительную терминологию. А именно, функциональный оператор ¥ индуцирует отображение P¥ : PN ^ L + (T, R), которое, по сложившейся в теории представлений традиции [10, с. 239], назовем проективизацией оператора Релея-Ритца:
P¥(y):=¥[y], у е Pn (у с Vn), где PN - вещественное проективное пространство, ассоциированное с линейным многообразием VN (с топологией, индуцированной из пространства L2(T,U)); т.е. PN есть множество орбит мультипликативной группы R , действующей на VN \ {0}. В данной геометрической трактовке ключевым моментом являются топологические свойства пространства PN, dimPN , разумеется, в первую очередь (в контексте предложения 1), его компактность, в частности, если имеет место dimVN = 3, то компактное 2-многообразие PN устроено как лист Мёбиуса, к которому по его границе приклеен круг [13, с. 162]. Попутно отметим, что на PN можно ввести геометрическую структуру CW-комплекса [13, с. 140], что, в свою очередь, упрощает рассмотрение вопроса о геометрической реализации многообразия PN - теорема 9.7 [13, с. 149].
Теорема 1. Каждое из следующих трех первых условий влечет за собой два других:
(i) БДР-задача (1)-(1*) разрешима относительно (Л15Л0,B, D15D2,D3,D4,D5) е L2;
(ii) 39е L2(T,R) :¥(ф) -L 0, V9е Vn ;
(iii) 3supLP¥[Pn]: supLP¥[Pn]е L2(T,R).
При этом, для выполнения (A1, Aq, B, Dj , D2, D3, D4, D5) е L необходимо чтобы
(iv) R(P¥[Pn ]) с L^ (T,R) .
Замечание 1. Теорему 1 можно рассматривать как начальный этап в изучении проблемы5, когда пучку управляемых траекторных кривых N от неявного дифференциального уравнения высшего порядка требуется сопоставить явную нестационарную билинейную дифференциальную систему второго порядка с тем же пучком управляемых траекторных кривых N.
Доказательство. Будем пользоваться идеями работы [4]. Придерживаясь определения 1 [4], введем в рассмотрение конструкцию M2 -оператора M : L2(T,U) ^ ЬДГ,X) вида
3 (Dj,D2,D3,D4,Ds,D6,D7,Ds)е L2:M(g,w,v, q, s,h, u, u):= = Dj g+D2 w+D3v+D4 q+Ds s+D6 h+DjU+Dsii V(g,w,v,q,s,h,u,u)е L2(T,U).
Остальные детали доказательства с небольшими уточнениями (с учетом приведенных выше лемм 1-3 для решетки R(P¥[PN ])) содержит аналитическая схема M2 -продолжимости в форме следствия 2 [4] и теоремы 3 [4] (принцип максимума энтропии).
При этом необходимое условие (iv) для
5 В частности, данная постановка уместна, когда разрешимую для пары (N, А1) БДР-задачу необходимо редуцировать к
разрешимой БДР-задаче для пары А2) в положении, когда еТ :Кег а41({) Ф 0е X} Ф 0, А2 - оператор гомотетии с коэффициентом 1; характер сопутствующих вычислений проиллюстрирован в примерах 1, 2.
(A, A, B, Dj, D2, D3, D4, D5) е L*
устанавливается модификацией доказательства теоремы 3 [11]. ■
Замечание 2. Необходимо отметить, что даже в случае 1 < Card N <Х0 имеет место положение Card PN = exp ; но можно показать, что существует (теорема 17 [7, с. 68]) такое счетное множество G с PN, что, если в пространстве L + (T, R) лежит supL P¥[PN ], то вещественную функцию := supL P¥[PN ] осуществляет следующая sup -конструкция:
t ^ £(t) = sup {P¥(y)(t) е R: у е G}.
Замечание 3. Доказательство теоремы 1 легко модифицировать так, чтобы сформулировать аналог теоремы 3 [12], выражающий в терминах углового расстояния в гильбертовом пространстве условия существования билинейной системы (1), реализующей пучки Nx, N2, каждый из которых обладает своей БДР-моделью, в том числе, в математической постановке [11], когда моделируемые операторы дифференциальной системы (1) суть стационарные, т.е.
(Д, A,B, Dj,D2,D3,D4,D5)е L*,
в частности, с оптимальной нормой [14].
В конкретных рассуждениях важен также следующий частный случай:
Следствие 1. Если dim VN < Х0, ¥[VN ] с L2(T,R) и найдется p е [1, <»), при котором
¥(Ф1 + ф2) -L p¥(ф1) + р¥(ф2), (ф15ф2)е Vn х Vn , то БДР-задача (1)-(1*) разрешима.
Заметим, что при p = 1 данное свойство в (контексте квазиупорядочения —L) сродни свойству "сублинейности" [12] функциональных операторов.
3. НЕПРЕРЫВНОСТЬ ОПЕРАТОРА РЕЛЕЯ-РИТЦА В АНАЛИЗЕ РАЗРЕШИМОСТИ БДР-ЗАДАЧИ
В случае компактности проективного многообразия PN (т.е. dim PN < Х0) естественно попытаться связать это свойство с задачей построения решетки R(P¥[PN ]) в контексте условий непрерывности проективизации оператора Релея-Ритца (см. также [15]). Ниже при выборе в предложении 1 метрической структуры в конусе L + (T,R) , прибегли к теоремам 15, 16 [7, сс. 65, 67]; в данной постановке L + (T, R) - полное сепарабельное метрическое пространство.
Предложение 1. Пусть dimPN < Х0 и конус L+ (T,R) наделен топологией, индуцированной сходимостью номере Д, или, что эквивалентно, инвариантной6метрикой
Ртf/2):= \I /i(x) -/2(х)|(1 +1 /г(т)-/2(х)|)"V(rfx), fx, / е L+ (T,R).
T
Тогда оператор P¥ : PN ^ L + (T,R) будет непрерывным, если пучок N таков, что
V9е Vn \{0}: supp^lU = T (mod ц), (2)
в частности, если
Vy е PN : suppP¥(у) = T (mod ц). (3)
Отметим, что предложение 1 является развитием теоремы 3 [16], что подтверждает её методологическую важность в дифференциальном моделировании сложных динамических систем [3-5]. Одним из приложений этого результата является следующее утверждение.
Следствие 2. Если при выполнении (2) или (3) оператор P¥ взаимно-однозначный, то P¥ - гомеоморфизм, а фундаментальная группа метрического пространства (P¥[PN],pT) изоморфна аддитивной группе целых чисел Z при dim Span N = 2 и группе вычетов Z 2 при dim Span N > 3, при-
6 Инвариантность означает рт (у + qí g + д) = рт (у; g), у; g; д е ь + (Т, Я) ■ Вариант неинвариантной неполной метрики, обеспечивающей непрерывность оператора Релея-Ритца при нарушении условия (2) или (3), рассмотрен в [15].
чем, пространство (P¥[PN],pT) ориентируемо, если размерность линейной оболочки Span N четная, и неориентируемо, если эта размерность нечетная.
Принимая во внимание, что непрерывная вещественная функция на компактном пространстве достигает своих наибольшего и наименьшего значений, приходим к заключению, что при выполнении (2) или (3) и учетом теоремы 5 [7, с. 28], для случая 1 < dimPN < К0 и при наличии supL P^[PN]
найдутся такие точки у', у"е PN, что будут верны оценки:
Pt(P^(Y'),X0) = sup(Pt(P¥(y),x0):ye Pn} < pr(supLP^P],%0) < ц(Т),
Pt (P^(y"),SupL Pn ]) = inf {Pt (P^(y), SUPl PF[Pn ]) : y e Pn } >0 .
Заметим, что включение P^(y') e L2(T,R) не гарантирует вложение R(P¥[Pn]) с L2(T,R) (см. пример 1 из [16]). При этом отметим, что dim PN = 0 приводит к положению
supL P4"[P„ ] = P4«[Pn ] = ||A + Mx + ll€ + iMív +
+HxI4x +1x||x +|Е<иШ +|вдМ(4)
В контексте теоремы 1 и предложения 1 можно уточнить условия существования решетки R(P^[PN]). В качестве отправной точки введем вспомогательную конструкцию: для натурального
n обозначим через Wn конечное n_1 -плотное подмножество в метрическом пространстве
(P^[PN],pT); подмножество Wn найдется в силу предложения 1. Ниже Limp {£,n} означает предел
последовательности {<^n} с L+ (T,R) в топологии, индуцированной метрикой pT .
Теорема 2. Пусть W . W = с P^P ] и fn := ^ v ... v^ Л = Ci v... v С,,
1 < i < n. Тогда конус L + (T, R) содержит решетку R(P^[PN ]), если и только если
Pt (fn, fm) ^ 0 (n, m , причём БДР-разрешимость приобретает вид: БДР-реализация (1)-(1*) существует тогда и только тогда, когда Lim р {fn} e L2 (T, R), что равносильно R(P^[PN ]) с L2 (T, R).
С учетом теоремы Радона-Никодима [17, с. 38] получаем
Следствие 3. Дифференциальная реализация (1)-(1*) существует, если и только если функционал S ^ J s (supL P^[PN](х)) ^R суть абсолютно р-непрерывная мера на .
В завершении приведем примеры, снимающие возможное представление, что выше особое ударение мы всюду делали исключительно на идейном аспекте каждого понятия, тем самым невольно пренебрегая его рассмотрением с вычислительной точки зрения; всюду ниже считаем, что моделирование осуществляется с нулевым запаздыванием х = 0 (т.е. y(-) = x(-)).
Пример 1. Пусть T = [0,10], Y := X, A - оператор гомотетии [12] с единичным коэффициентом, A = 0e L(X,X), Ц = D3 = D4 = D5 = 0e Z(X2,X), ee X, || e ||X = 1,
t ^ x(t) = (tsint)e, t ^u(t)=0e L2(T,X).
2 2 2 1/2
Тогда функция f := sup^P¥(Pn) = || x||x (||xц^^ +1| xц^у || X ц^у)" (см. рис. 1) не принадлежит пространству L2(T, R) и, следовательно, согласно теореме 1 и формуле (4), реализация (1) для неуправляемого процесса N = {(x, u)j не существует.
Пример 2. Изменим постановку примера 1 тем, что
t ^ u{t) = (t sin21 + 2-112 sin2t + cost )e.
Тогда (см. рис. 2) f := supL P^(PN) = || x\\X (||x||x +|| x||x || X|x +||u||y)"1/2 eL2(T, R) и значит реализация (1) для управляемого процесса N = {(x,u)} существует; нетрудно установить, что
x + x = 2u - 2D 2 (x, x), где D2 =< •,• >X e, < v > - скалярное произведение в X.
Отметим, что для более сложных вариантов задания пары (г ^ х{г) , г ^ и )) символьные вычисления функции /2( ) (аналогичной рис. 1, 2) можно проводить средствами компьютерной алгебры в математической физике [18]. В данном контексте анализ разрешимости БДР-задачи примеров 1, 2 можно модифицировать для качественного анализа редукции точных многомерных решений диффузии со степенными нелинейностями к задаче Коши для счетной системы обыкновенных дифференциальных уравнений с полилинейной структурой.
Следующий пример характеризует факт наличия свойства нестационарности БДР-модели, как эндогенного фактора дифференциальной реализации пучка N (см. пункт (¡у) теоремы 1).
Рис. 1. f 2(t) = (2 cos t -1 sin t)2 х х ((t sin t)2 + (t sin t)2 (sin t +1 cos t)2)-1
Рис. 2. f 2(t) = (2 cos t -1 sint)2 х(( t sint)2 + +(t sin/)2 (sint+1 cost)2 + (t sin2/+2-1/2 sin2i+cost)2)-
Пример 3. Пусть T = [0,1] и коэффициент гомотетии оператора A равен 5, при этом
N = {(x, u)}, t ^ x(t) = t1,6e, t ^ u(t) = iT (t)e. Тогда (как несложно установить) функция f := supL P¥(PN) согласно формулы (4) удовлетворяет
неравенствам
0,61-0'4 - supL P¥(PN) - 4,81~
где второе неравенство гарантирует (в силу пункта (ш) теоремы 1) существование некоторой БДР-модели, реализующей процесс-пучок N, при этом первое неравенство, согласно пункта (¡у) теоремы 1, показывает, что данная модель не может быть стационарной.
Приведем пример, иллюстрирующий положение, когда свойство линейности дифференциальной модели реализации недостаточно для прецизионного построения уравнений динамики пучка N, что важно при моделировании синергетического фактора нейропопуляции.
Пример 4. Пусть Т = [1, 2], А - оператор гомотетии с коэффициентом 1, при этом
N = {(*!, и1),(х2, и2),(х3, и3)}, г ^ хх{г) = (г2 + 2)в, t^и(t) = 2-1хТ(t)в, t ^ х2(г) = гв, г^и2(0=%Т(г)в,
г ^ х3 (г) = (2 - 4Л)гв, г ^ и3 (г) = %Т (0е.
Вначале покажем, что динамический пучок N не может иметь линейную дифференциальную реализацию. С этой целью достаточно установить, что линейная характеристика , отвечающая оператору Релея-Ритца, индуцирует функцию г ^ Ль^'с нулем порядка 2. Характеристическим условием данного положения будет следующая система уравнений (относительно параметров а, Р):
|| х (г) + ах2 (г) + рх3 (г) ||Х = ((г2 + 2) + ш + р(2 - 4л/2> )2 = 0,
у(а,р)=а+р(2-4л/2) >
у(а,р)=а+р(2-4,/2)
II х (г) + ах 2 (г) + рх3 (г) ||Х = (2г + а% т (г) + Р(2 - 4л/2)х т (г ))2 = О, || их (г) + аи 2(г) + ри3 (г) ||Х = (2-1 г + аг + рг )2 = О, откуда, несложно вычислить, что а = -1, Р = 0,5, при этом точкой нуля будет г = у/2 .
Теперь покажем, что БДР-модель (1) для пучка N существует. Для этого достаточно установить, что сумма линейной и билинейной характеристик (параметризированных по а, Р) оператора Релея-Ритца для пучка N ограничена снизу некоторой функцией вида:
г ^ Г1Т (г), г е (0, На этом пути легко видеть, что при любых г е [1,2], а,Р е Я будет
%0,а,р) + Лв(г,а,Р) > >|| хх(г) + ах2(г) + рхз(г) \ \2Х + || их{г)хх(г) + аи2(г)х2(г) + ри3(г)хъ{г) \\2Х =
= (2г + а + Р(2 - 4л/2))2 + (г2 + аг + Р(2 - 4л/2)г )2 = = [(2г + у(а,Р))2 + г 2(г + у(а,р))2]
> [(2г + у(а,Р))2 + (г + у(а,р))2],
= [(г + 1{г,у))2 + 12(г>гН+г(а>Р) >
>[(г + 9(г))2 + 92(г)]|е(0=_г/2 = г2/2> 0,5хт(г).
Можно проверить, что БДР-модель (1) для пучка N имеет аналитическое представление:
х + х = В5(и, х),
где В5(у) = -<•,• >х е.
Пример 5. Конструкция из сноски 1, позволяет построить общую дифференциальную реализацию (с вектором состояний X © X © X) для динамических пучков из примеров 2-4.
Непосредственно из примера 4, как его методологическое обобщение, вытекают два важных теоретико-модельных вывода, с аналитической проверки которых по существу должен начинаться любой качественный анализ задачи БДР-моделирования (1)-(Г); при этом позволяющих сократить на базе чисто инженерных идей (основанных на опыте и здравом смысле) разрыв между теорией и практикой. Первый вывод состоит в том, что конечный динамический пучок N может иметь на интервале Т интерполяционное представление в классе полиноминальных сплайн-функций [19]. Второй вывод выражает следующее утверждение:
Предложение 2. Пусть динамический пучок N удовлетворяет условию:
V,д,к,щи)е VN : || А^ ||х Оу/п^^^у) + X^)-1 * Ь2(Т,Я), (5)
где Хз - характеристическая функция множества Sg = Т\suppg.
Тогда дифференциальная реализация данного пучка не может иметь аналитическое представление (1)-(1*), в котором все билинейные операторы суть нулевые.
Условие (5) можно ослабить (конструктивно усилить), сведя его к поиску средствами компьютерной алгебры нулей функции + , параметризированной коэффициентами интерполяционного представления пучка N в классе полиноминальных сплайн-функций; см. также полезные предложения 2.6.1 и 2.6.4 [19, сс. 67, 69] о расположении нулей полиномов.
ЗАКЛЮЧЕНИЕ
Предложенные выше результаты подтверждают, что перманентно развиваемая КТДР может представлять универсальный язык для междисциплинарного обмена научными результатами в рамках энтропийного анализа. Действительно, она достаточно обща для того, чтобы не вносить своих собственных ограничений в математическое описание процессов апостериорного моделирования нелинейных бесконечномерных динамических систем, и в тоже время, в силу своей строгости, она устраняет возможность весьма распространённых разночтений в толковании такого ключевого понятия теории идентификации систем, как структурная идентифицируемость. При этом не позволяя этому
понятию окрасится в субъективные тона, а напротив, придав ему полностью математически формализованный смысл [20], обеспечивая устойчивый интерес со стороны тех, кто ставит себе в КТДР не самые узкие цели. В целом можно убедиться, что понимание природы полилинейных моделей, как формы представления уравнений динамики сложных систем (в том числе нейроморфных), способствует прояснению интуитивных представлений об их динамике и мотивировке всего обсуждения. К тому же, не откажем себе в удовольствии заметить, что "билинейный аппарат" КТДР, получивший продвижение выше, доставляет, по меньшей мере, эстетическое удовлетворение; аналитик, преследующий чисто эстетические цели, содействует созиданию метаязыка, более приспособленного к тому, чтобы удовлетворить математические запросы нейромоделирования [21-23].
Для того чтобы пойти дальше, не сложно указать то теоретико-системное направление, которое может составить алгебраическую основу (без чрезмерного нагнетания средств алгебраической геометрии) следующего этапа развития качественной теории дифференциальной реализации высших порядков [5], в том числе, с неполным измерением вектора состояний [24]. А именно, переход в моделируемой нейроморфной динамике от билинейной структуры нелинейных связей к полилинейным7 связям. Методологически этот переход состоит в использования геометрического языка тензорных структур пространств Фока [17] и проективных представлений [10] в контексте исследования метрический свойств энтропийного операторов Релея-Ритца [15] средствами компьютерной алгебры в математической физике [18].
СПИСОК ЛИТЕРАТУРЫ
1. Поляков Г.И. О принципах нейронной организации мозга. М.: МГУ, 1965.
2. Brzychczy S., Poznanski R. Mathematical Neuroscience. Academic Press. 2013.
3. Калман P., Фалб П., Арбиб M. Очерки по математической теории систем. М.: Мир, 1971.
4. Rusanov V.A., Banshchikov A.V., Daneev A.V., Lakeyev A.V. Maximum Entropy Principle in the Differential Second-
Order Realization of a Nonstationary Bilinear System // Advances in Differential Equations and Control Processes. 2019. Vol. 20. No. 2. P. 223-248.
5. Rusanov V.A., Daneev A.V., Lakeyev A.V., Sizykh V.N. Higher-Order Differential Realization of Polylinear-Controlled
Dynamic Processes in a Hilbert Space // Advances in Differential Equations and Control Processes. 2018. Vol. 19. No. 3. P. 263-274.
6. Савельев A.B. Источники вариаций динамических свойств нервной системы на синаптическом уровне в нейрокомпьютинге // Искусственный интеллект. НАН Украины. 2006. № 4. С. 323-338.
7. Канторович Л.В. Акилов Г.П. Функциональный анализ. М.: Наука, 1977.
8. Иосида К. Функциональный анализ. М.: Мир, 1967.
9. Rusanov V.A., Daneev A.V., Lakeyev A.V., Linke Yu.E. To Existence of a Nonstationary Quasi-Linear Vector Field Re-
alizing the Expansion of a Control Trajectory Bundle in Hilbert Space // WSEAS Transactions on Systems. 2020. Vol. 19. P. 115-120.
10. Кирилов A.A. Элементы теории представлений. M.: Наука, 1978.
11. Rusanov V.A., Daneev A.V., Lakeev A.V., Linke Yu.E. On the Differential Realization Theory of Nonlinear Dynamic Processes in Hilbert Space // Far East Journal of Mathematical Sciences. 2015. Vol. 97. No. 4. P. 495-532.
12. Эдварде P. Функциональный анализ: Теория и приложения. М.: Мир, 1969.
13. Прасолов В.В. Элементы комбинаторной и дифференциальной топологии. М.: МЦНМО, 2014.
14. Русанов В.А., Данеев А.В., Линке Ю.Э. К оптимизации процесса юстировки модели дифференциальной реализации многомерной системы второго порядка // Дифференциальные уравнения. 2019. T. 55. № 10. С. 14321438.
15. Rusanov V.A., Daneev A.V., Lakeyev A.V., Linke Yu.E. On the Theory Differential Realization: Criterions for the Continuity of the Nonlinear Rayleigh-Ritz Operator // International Journal of Functional Analysis, Operator Theory and Applications. 2020. Vol. 12. No. 1. P. 1-22.
7 Полилинейные дифференциальные уравнения занимают в теории динамических систем особое место, поскольку согласно одной из основных идей нелинейного многомерного анализа является положение, что всякая достаточно гладкая (относительно производных Фреше) функция многих переменных в окрестности каждой своей точки хорошо аппроксимируется полилинейной функцией - формула Тейлора [25, с. 491] для отображений в бесконечномерных банаховых пространствах. Возникающая, таким образом, операция полилинеаризации приводит к полилинейным дифференциальным моделям в качестве эффективного прецизионного приближения при исследовании свойств нелинейных дифференциальных уравнений вблизи каких-либо их решений. Именно этот факт дает возможность проводить анализ поведения любых (репрезентативных) нелинейных нейроморфных дифференциальных моделей высших порядков (с запаздыванием и без) высокой размерности на базе полилинейных структур динамических уравнений исследуемой нейропопуляции. При этом в перспективе претендуя, в процессе обработки данных от интерфейс-платформы "мозг-машина" типа ЫеигаНпк [23], на выявление специальных (типических) динамических свойств этих моделей, отвечающих за когнитивные функции моделируемой нейропопуляции.
16. Русанов В.А., Данеев А.В., Линке Ю.Э. К геометрическим основам дифференциальной реализации динамических процессов в гильбертовом пространстве // Кибернетика и системный анализ. 2017. Т. 53. № 4. С. 71-83.
17. Рид М., Саймон Б. Методы современной математической физики: Функциональный анализ. М.: Мир, 1977.
18. Банщиков А.В., Иртегов В.Д., Титоренко Т.Н. Программный комплекс для моделирования в символьном виде механических систем и электрических цепей // Свидетельство о государственной регистрации программы для ЭВМ. №2016618253 от 25.07.2016. Федеральная служба по интеллектуальной собственности (РОСПАТЕНТ).
19. Лоран П.-Ж. Аппроксимация и оптимизация. М.: Мир, 1975.
20. Rusanov V.A., Daneev A.V., Lakeev A.V., Linke Yu.E., Vetrov A.A. System-Theoretical Foundation for Identification of Dynamic Systems // Far East Journal of Mathematical Sciences. 2019. Vol. 116. No. 1. P. 25-68.
21. Савельев A.B. Философия методологии нейромоделирования: смысл и перспективы // Философия науки. 2003. № 1. С. 46-59.
22. Савельев А.В. На пути к общей теории нейросетей. К вопросу о сложности // Нейрокомпьютеры: разработка, применение. 2006. № 4. С. 4-14.
23. Valle G. An Integrated Brain-Marine Interface Platform with Thousands of Channels // J. Med. Internet Res. 2019. Vol. 21. No. 10: e16194. (doi: 10.2196/16194).
24. Русанов B.A., Шарпинский Д.Ю. К теории структурной идентификации нелинейных многомерных систем // Прикладная математика и механика. 2010. Т. 74. Вып. 1. С. 119-132.
25. Колмогоров А.Н., Фомин С.В. Элементы теории функций и функционального анализа. М.: Наука, 1976.
ON DIFFERENTIAL-NON-AUTONOMOUS REPRESENTATION INTEGRATIVE ACTIVITY OF NEUROPOPULATION BILINEAR SECOND-ORDER MODEL WITH A DELAY
© 2021 A.V. Daneev1, A.V. Lakeev2, V.A. Rusanov2, P.A. Plesnev1
1 Irkutsk State Transport University, Irkutsk, Russia 2 Matrosov V.M. Institute for System Dynamics and Control Theory of SB RAS, Irkutsk, Russia
For neuromorphic processes spedfied by the behavior of a local neuropopulation (for example, processes induced by a brain-machine interface platform of the Neuralink type), we study the solvability of the problem of the existence of a differential realization of these processes in the class of bilinear nonstationary ordinary differential equations of the second order (with delay) in separable Hilbert space. This formulation belongs to the type of inverse problems for an additive combination of nonstationary linear and bilinear operators of evolution equations in an infinite-dimensional Hilbert space. The metalanguage of the theory being developed is the constructions of tensor products of Hilbert spaces, lattice structures with orthocompletion, the functional apparatus of the nonlinear Rayleigh-Ritz operator, and the principle of maximum entropy. It is shown that the property of sublinearity of this operator, allows you to obtain conditions for the existence of such differential realizations; along the way, metric conditions for the continuity of the projectivization of this operator are substantiated with the calculation of the fundamental group of its compact image. This work was financially supported by the Russian Foundation for Basic Research (project no. 19-01-00301). Keywords: inverse problems of nonlinear neuromorphic dynamics, tensor analysis, bilinear non-autonomous differential realization of the second order with delay, the Rayleigh-Ritz entropy operator. DOI: 10.37313/1990-5378-2021-23-2-115-126
Aleksey Daneev, Doctor of Technical Sciences, Professor, Professor of the Department of Information Systems and Information Security. Е-mail: [email protected] Anatoly Lakeev, Doctor of Physical and Mathematical Sciences, Leading Researcher. E-mail: [email protected]
Vyacheslav Rusanov, Doctor of Physical and Mathematical Sciences, Associate Professor, Senior Researcher. E-mail: [email protected] Pavel Plesnev, Postgraduate Student. Е-mail: [email protected]
