CC BY
52
Теоретическая механика
УДК 531.31; 519.21
О построении дифференциальных уравнений по заданному интегральному многообразию при наличии случайных возмущений с независимыми приращениями
М. И. Тлеубергенов, Д. Т. Ажымбаев
Лаборатория динамических систем Институт математики ул. Пушкина, д. 125, Алматы, 050010, Казахстан
Строятся уравнения Лагранжа, Гамильтона и Биркгофа по заданным свойствам движения в классе стохастических дифференциальных уравнений типа Ито при наличии случайных возмущающих сил из класса процессов с независимыми приращениями. Полученные результаты иллюстрируются на примере движения искусственного спутника Земли под действием сил тяготения и аэродинамических сил.
Ключевые слова: дифференциальные уравнения, стохастические уравнения, интегральное многообразие, функция, вероятность, вектор.
Введение
Основы теории и общие методы решения обратных задач дифференциальных систем разработаны в [1-8] и др. для детерминированных систем, уравнения которых являются обыкновенными дифференциальными уравнениями (ОДУ). Так, в работе Еругина [1] строится множество ОДУ, которые имеют заданную интегральную кривую. Эта работа впоследствии оказалась основополагающей в становлении и развитии теории обратных задач динамики систем, описываемых ОДУ. В работах [2-8] изложены постановка, классификация обратных задач дифференциальных систем и их решение в классе ОДУ. Следует отметить, что один из общих методов решения обратных задач динамики в классе ОДУ предложен в [7,8].
В работах [9-11] обратные задачи динамики рассматриваются при дополнительном предположении о наличии случайных возмущений из класса винеровских процессов и, в частности, решены:
1) основная обратная задача динамики — построение множества стохастических дифференциальных уравнений второго порядка типа Ито, обладающих заданным интегральным многообразием;
2) задача восстановления уравнений движения — построение множества управляющих параметров, входящих в заданную систему стохастических дифференциальных уравнений второго порядка типа Ито, по заданному интегральному многообразию;
3) задача замыкания уравнений движения — построение множества замыкающих стохастических дифференциальных уравнений второго порядка типа Ито по заданной системе уравнений и заданному интегральному многообразию.
Для разрешения обратных задач широко используется метод квазиобращения, в основе которого лежит Лемма 1 [7, с. 12-13].
Лемма 1. Совокупность всех решений линейной системы
ну = д, Н = ), v = (ук), д = (д^), ^ = 1, т; к = 1 ,п, т < п, (1)
Статья поступила в редакцию 22 августа 2012 г.
где матрица Н имеет ранг, равный т, определяется выражением
V = 8УТ + Vй.
Здесь 8 — произвольная скалярная величина,
(2)
ут = [НС] = [Н ... Нтст+1 ... сп— 1] =
61 . . Сп
Н11 . . Н1п
Нт1 . Н . Нтп
Ст+1,1 . . Ст+1,п
Сп-1,1 . . Сп— 1,п
есть векторное произведение векторов Н^ — (Н^) и произвольных векторов Ср — = (срк), р = т + 1,п — 1; ек — единичные орты пространства Кп, уТ = ), где
0
Н11
Нт1
Ст+1,1
Сп-1,1
1
Н1к
Ншк Ст+1,п
Сп— 1,к
0
Н1п Н
Нтп Ст+1,п
Сп— 1,п
, у" = Н+д,
Н + = Н (НН ) 1, Н —матрица, транспонированная к Н.
Обратные задачи построения уравнений движения в форме Лагранжа, Гамильтона и Биркгофа по заданным свойствам движения в классе ОДУ (при отсутствии случайных возмущающих сил) рассматривалась в [12]. В [13] рассмотрена стохастическая задача Гельмгольца — задача построения по заданному стохастическому уравнению Ито второго порядка эквивалентных стохастических уравнений лагранжевой, гамильтоновой и биркгофиановой структуры. В [14,15] рассматриваемая в данной работе задача исследуется в предположении, что случайные возмущающие силы есть системы независимых винеровских процессов (как частный случай процессов с независимыми приращениями).
В данной работе рассматриваются задачи построения функции Лагранжа, Гамильтона и Биркгофа по заданным свойствам движения в предположении, что случайные возмущающие силы принадлежат более общему классу случайных процессов, а именно, классу процессов с независимыми приращениями.
1. Постановка задачи
По заданному множеству
Л(г): \(х, X,г) = 0, X е кт, х е кп (3)
требуется построить стохастические уравнения лагранжевой, гамильтоновой и биркгофиановой структуры
ё / дЬ \ дЬ ' . • • , -— . -—ч
К [в^) — = ^(Х>Х, (" = 1>п>1 = , (4)
Цк =
дН
дрк,
дН ■
Рк = - + о'ьэ (4,р,£)£:>, (к = 1,п)
(5)
дЩ(г,г) дЩ(г, г)
дхг дzi
Zi -
дВ(г,€) дЩг,г) дх1 №
(г,1 = 1, 2п, р = 1,п + г) (6)
так, чтобы множество Л(1) было интегральным многообразием построенных уравнений. Здесь ..., и {г^\(1,ш), ..., фп+г(Ь,1^)} — системы процессов с независимыми приращениями [4], В = В(г,Ь) — функция Биркгофа, а
, , ГдЕЛх,г) дШг,€)'
№ = (№ц) — тензор Биркгофа с компонентами =
дх1
Поставленная задача в классе обыкновенных дифференциальных уравнений ранее рассматривалась в [12]. В [13] рассмотрена стохастическая задача Гельм-гольца — задача построения по заданному стохастическому уравнению Ито второго порядка эквивалентных стохастических уравнений лагранжевой, гамильтоно-вой и биркгофиановой структуры. В [15,16] поставленная выше задача рассматривается в предположении, что системы ..., (Ь,ш)} и {-ф\(1,Ш), ..., 'фп+г (Ъ,ш)} есть системы независимых винеровских процессов (как частный случай процессов с независимыми приращениями).
В данной работе задачи построения функций Ь, Н и В по заданным свойствам движения рассматриваются в предположении, что случайные возмущающие силы принадлежат классу процессов с независимыми приращениями.
Таким образом, пусть ... и {г^\(1,ш), ..., фп+г— систе-
мы случайных процессов с независимыми приращениями, которые, следуя [16], можно представить в виде суммы процессов:
1) £ = + / с(у)Р°(1,йу), — винеровский процесс; Р0 — пуассоновский процесс; Р0(Ь^у) — число скачков процесса Р0 в интервале [0,Ь], попадающих на множество dy; с(у) — векторная функция, отображающая пространство В?п в пространство значений Кк процесса £(1) при любом ¿;
2) 'ф = ф0 + / с(у)Р0ё7), ф0 — винеровский процесс; Р0 — пуассоновский процесс; Р0(Ь, ё7) — число скачков процесса Р0 в интервале [0,Ь], попадающих на множество ё7; с(у) — векторная функция, отображающая пространство Р2п в пространство значений рп+г процесса ф(Ь) при любом t.
Для решения поставленных задач на первом этапе по заданному множеству методом квазиобращения [7,8] в сочетании с методом Еругина [1] и в силу стохастического дифференцирования сложной функции в случае процессов с независимыми приращениями [16] строится дифференциальное уравнение Ито второго порядка
х = / (х,х + а(х,х ,£)£ (7)
так, чтобы множество Л(1) (3) являлось интегральным многообразием построенного уравнения (7). И, далее, на втором этапе решается стохастическая задача Гельмгольца, а именно, по построенному уравнению Ито (7) строятся эквивалентные ему стохастические уравнения лагранжевой, гамильтоновой, а также бирк-гофиановой структуры.
2. Построение стохастического уравнения лагранжевой структуры (4) по заданным свойствам движения (3)
Предварительно, по правилу стохастического дифференцирования Ито, составляется уравнение возмущенного движения
; дХ дХ • ЭХ . „ „ ЭХ •
Х = + 7Т ж + 51 + ¿2 + ¿з + дъ дх ах ах
(8)
я 1д2Х т д2Х п п т [16]
где ¿1 = -7ТТ7Т : (У(У , а под —тт : п, п = аа± , следуя 1161, понимается вектор,
2 ох2 ох 2
элементами которого служат следы [16] произведений матриц вторых производных соответствующих элементов Х)Л(х,х, 1) вектора Х(х,х, 1) по компонентам х на матрицу п
I т
дЧ
дх2
: п =
...
\ дх2 ) ' ' \ дх2 )
¿2 = [ ^Х(х,х + ас(у), 1) — Х(х,х, ^ + ^ас(у)^ ёу; ¿з = ¡[Х(х,х + ас(у), 1) —
Х(х,х, £)]Р0(1, ёу) и вводятся произвольные функции Еругина [1] А и В, обладающие свойствами А(0,х,х, 1) = = 0, В(0,х,х, 1) = 0, и такие, что имеет место равенство
А = А(х,х,х, г) + в(х,х,х, г) £. (9)
Сравнивая уравнения (8) и (9) приходим к соотношениям
ЭХ . . ЭХ ЭХ .
ттт I = А — — — —х — ¿1 — ¿2 — ¿3,
ОХ 01 ОХ
эх в = В.
ОХ
(10)
Для определения из (10) искомых функций f и а воспользуемся леммой 1, в силу которой вектор-функцию f и матрицу а запишем в виде
= 1
—С дх
+
(&с) (
'А — ^ — ^ — ¿1 — ¿2 — ¿3 01 ОХ
(11)
2
—С дх
+
+ ^дх ) г
(12)
где аг = (а1г, а2г,..., апг) — г-й столбец матрицы а = (<&1>])\V = 1,п, ] = 1,г) ;
Вг == (В1г,В2г,..., Втг)Т — г-й столбец матрицы В = (В^), = 1,т, ] = 1,г), 1 , 2 — произвольные скалярные величины.
Таким образом, из (11), (12) следует, что множество дифференциальных уравнений Ито второго порядка, обладающих заданным интегральным многообразием (3), имеет вид
х = в1
—С дх
(дх\ + /. дх дх. а а а .
+ ( д±) { А — — — —х — ¿1 —¿2 —¿3 ) +
дЬ дх
—С дх
(дхХ п
Паv В1,..., ^
+
—с
дх
Далее раскроем выражение — [ —— ) по правилу стохастического дифферен-
ёъ \ дХу 1
цирования Ито в случае процессов с независимыми приращениями [16]:
ё / д— \ д2— д2— . д2— .. - - -
-Г ^Г- = я- я, + -77^,—Хк + я- я. Хк + Ь\и + Ь2и + Ьзи, (13)
т \дхи) дхудъ дхудхк дхудхк
х 1 д3Ь ~ [(д—(х,Х + ас(у), ^ д—(х,Х, ^} ,
где Ь1и = - . . (7^(7кэ, Ь2и = \-7Г.---7Г.-} ёу,
2 дхудхгдхк ] у дху дху )
д—(х,Х + ас(у), 1) д—(х,Х, 1)
Ьзи =
Р, у).
дху дху
Следовательно, уравнение (4) с учетом (13) принимает следующий вид:
А (д— ^ д—-' ( • ) ¿3 = д2— д2— . д2— .. ё \дХУ ) дху дХУдЬ дХУдхкХк дХУдХкХк
д — '
+ Ь1и + Ь'2„ + Ь3и - — - (х, Х, . (14) Или, учитывая (14) и уравнение (7), имеем
д2— д2— . д2— .. - - - д— ' . . ..
+ —Хк + я- Хк + Ь1и + Ь2и + Ь3и - я--а • (х, х, ^ ^ =
гч ♦ гч. I гч ♦ гч ^ к I «Л • «Л* ^ к I ^ 1У I ^ 2 У I 3 У гч
дхудъ дхУдхк дхУдхк дх
У
= хи - ¡и(х, Х, £) + а^(х, Х, Ь)^. (15)
Из соотношения (15) следуют равенства
д2— ^ д2— д2— . ... д— ~ ~7Г- о хк + Ь1у
дхУдхк
я • Я • = ; Я- Я4- + я- Я-Хк + + Ь2у + - -- = /у,
дхудхк дхидъ дхудхк дху (16)
Тем самым доказана
Теорема 1. Для построения стохастического уравнения лагранжевой структуры (4) по заданному множеству (3) так, чтобы множество (3) было интегральным многообразием построенного уравнения необходимо и достаточно выполнения условий (16).
3. Построение стохастической уравнения гамильтоновой структуры (5) по заданным свойствам движения (3)
Для построения функции Гамильтона предварительно введем новую переменную ук = Хк и перепишем заданное уравнение (5) в виде
( хк = Ук,
I Ук = ! к (х, У, ^ + Окз (х, У, ^
где вектор-функция / = (/1, /2,..., /п)Т и столбцы матрицы а = (а1, а2,..., аг) имеют соответственно вид (11), (12). Затем с помощью замен
*к =
{2 ; °к={I; ^={
0, при ] = 1,п ;
; = \ -
^ п, при ] = п +1,п + 2,...,п + т;
/ \ I 0пхп 0пхт \ / \
^ = (^ко) = 0 _ ; а = (а^)
\ 0ПХП ипхш /
перепишем уравнение (17) в виде
¿к = Ск (г, г) + Ркз (г, ^ 1ру,
(18)
а с помощью некоторой матрицы Г = запишем непрямой вид уравнения (18) в виде
(¿к — Ок — Рк^ ^ = 0. (19)
Далее, стохастическое уравнение гамильтоновой структуры (5) с помощью замены _
к = 1,п,
к = п + 1,п + 2,..., 2п
Vк
Г Чк,
\ Рк-П,
и матриц
V = ) = ^
0ПХП 1 пхп 0
1 пхп 0пхп
) ■
Р = (р кэ) =
I 0пхп 0пхтп \ \ 0пхп апхт )
а также с учетом того, что
( —I ) = №)
перепишем в виде
дН ■
Vк — фки = РкзУз
или, если ввести обратную к (^ки) матрицу
и вектор
(^ки) = (фкн)
%к = =
1
0
I К
пхп 1пхп IпXп X п
—рк, к = 1,п
дк-п, к = п + 1,2п
тогда уравнение (20) преобразуется к эквивалентному уравнению
дН ■
(20)
(21)
Рассмотрим задачу непрямого представления уравнения (18), или, что то же самое, задачу представления уравнения (19) в виде уравнения гамильтоновой структуры (21), т.е. рассмотрим соотношение
дН
(¿к — Ск — = ^ик ¿к — — ^„кРизФз
(
или
дН
С„к¿к - (г;, Ь) - з = ¿к - д--, (22)
где С„к = 7к; (х, г) = ^вк.
Для выполнения тождества (22) требуется выполнение условий
д Н
С„к = ш„к, (г, Ь) = - — , (23)
= ш^кРиз к = 1, 2п, з = 1,п + т) , (24)
7^ = шик ■ (25)
Из соотношений (22) и условий (23)-(25) следует, что ркз = Р^з, а это влечет выполнение равенства а^ = а'к^, (к = 1,п, з = 1,т) . Следовательно справедлива Теорема 2:
Теорема 2. Для непрямого построения стохастического уравнения гамиль-тоновой структуры (5) по заданному множеству (3) так, чтобы множество (3) было интегральным многообразием уравнения (19), необходимо и достаточно выполнение условий (22)-(24).
4. Построение стохастического уравнения биркгофиановой структуры (6) по заданным свойствам движения (3)
Для решения поставленной задачи рассмотрим соотношение
Сик¿к-Ри(г, 3 =
дЕк (г, г) дЕи (г, г)
д г»
д гК
дВ(г, ^ + дЕи (г, ^
д
д
-Тиз1рз ,
(и, к = 1, 2п, з = 1, п + т), которое выполняется тождественно при следующих условиях:
Сик —
дЕК (г, Ь) дЕи (г, Ь)
д
д гк
Ри =
дВ (г, г) + дЕи (г, г)
д г,.
д
р = Т. (26)
Следовательно, имеет место Теорема 3:
Теорема 3. Для построения стохастического уравнения биркгофиановой структуры (6) по заданному множеству (3) так, чтобы множество (3) было интегральным многообразием уравнения (6), необходимо и достаточно выполнение условий (26).
5. Пример
Рассмотрим стохастическую задачу построения функций Лагранжа, Гамильтона и Биркгофа по заданному свойству движения на примере движения искусственного спутника Земли под действием сил тяготения и аэродинамических сил [17].
Пусть свойства движения заданы в виде:
А (г): А = §2 + аг192 + = 0, А е Е1. (27)
Тогда уравнение возмущенного движения (8) примет вид
А = 2№ + 2а1М + 51 + ¿2 + ¿3 = 2М + 2а.1$/ + 31 + 5 + 5з + Ъа^а^, (28)
где
Si = аха2, = j ^2aiac(y)[4'd + ac(y)]^áy, S3 = j {2aiac(y)[M + ac(y)]}p°(t, dy).
Введем функции Н. П. Еругина а = a(X,i,H, t), b = b(X,i,H, t) со свойством a(0, i, i, t) = = b(0, i, i, t) = 0 и такие, что имеет место соотношение
X = a(X,i,í, t) + b(X,i,í, t)l (29)
Из соотношений (27), (28) следует, что множество уравнений (7), в нашем примере имеющее вид
i = f(i,í, t) + a(tí,í, t)Í,
будет обладать интегральным многообразием (26), если f и а будут иметь соответственно вид
. a(X,i,í, t) - 2ii -Si -S2 -S3 b(X,i,í, t)
T = -:-, а = -:-. (3U)
2aiü 2ai$
Уравнение движения искусственного спутника Земли под действием сил тяготения и аэродинамических сил, следуя работе [17], запишем в виде
i = f($,i) + a($,i)i (31)
где i — угол тангажа, а функции f и аимеют вид
f = Qlsin2i -Q[g(i) + r¡í], а = Q5[g(i) + rji]. (32)
Построим лагранжиан по уравнению (31). Предварительно в уравнении (31) следует учесть соотношения (3U), которые обеспечивают интегральность заданного множества (27). Из равенств f = f, а = а вытекает, что четыре параметра Q,5,r¡,l, определяющие динамику движения спутника (31), должны удовлетворять следующим соотношениям
í а(Х, i, i, t) - 2ii -Si -S2 -S3 = 2aii)\Ql sin 2i - Q[g(i) + \ b(X,i,i, t) = 2aiÍQ5 [g(i) + r]i].
Тогда по определению из [18] уравнение (31) допускает непрямое аналитическое представление в терминах стохастического уравнения лагранжевой структуры, если существует функция h такая, что имеет место тождество
Ч S) - ^ , Ы = W -f-а& (33)
Найдем функцию h = h(t) так, чтобы для скалярного уравнения ^ (i, i, t)$ + 12(i, i, t) = 0 выполнялось необходимое и достаточное условие Гельмгольца [18, с.
1ГУ71 912 911 , h911
1U7I существования лагранжиана, которое имеет вид —- = ——+ i^r;.
9i 9t 9i
Этому условию удовлетворяет функция Ь вида Ь = е ь. Подставляя найденное Ь в (33) получим
п™/. д2Ь д2Ь • д2Ь •• - ~ ~ дЬ . • [I? - ! = + —.-$ +--г—$ + 5 + Я2 + 5з - ^т - а'^,
где
5 = -
1д 3ь
а
¿2 =
5з =
2д$3
дЬ($,$ + а с (у), ^ дЬ($
дЬ($,$ + а с (у), ^ дЬ($,$, ^
№ ) <У, $ ) У
Таким образом, искомый лагранжиан строится в виде
Р0^, У).
Ь =
-я(21соъ2$ + в)
0 = 1 (34)
который обеспечивает непрямое аналитическое представление уравнения (31) в виде уравнения лагранжевой структуры
<1 дЬ дЬ _от ■ ■
---г - — = е пща($, $)£.
ИдФ д$ V > ^
(35)
Используя функцию Лагранжа (34) и преобразование Лежандра, определим функ-
дЬ
цию Гамильтона в виде Н = х$ - Ь($,$, 1) . И так как % = —т, то
X = и, следовательно, $ = Тогда каноническое уравнение, соот-
ветствующее стохастическому уравнению лагранжевой структуры (35), примет вид
■ дН
ддн ■ (36)
% = - +&($,х, Ж,
где а = а'($, 1) х1), а функция Гамильтона определяется в виде
Н = 1 еп\2 е-^1 Ь($).
(37)
Рассмотрим теперь на этом примере задачу представления биркгофиана по заданному уравнению (31). По построенному выше каноническому уравнению (36) и
функции Гамильтона (37) с использованием соотношений (26) при С = ^ ( ^
функции Еv(ь = 1, 2), Е = (Е1,Е2) и В определяются в следующем виде Е = {х, (1 + р)$}, В = 2реп™гх2 - рЬ($), где р — произвольная постоянная.
Заключение
Таким образом, по заданному интегральному многообразию методом Еруги-на в сочетании с методом квазиобращения строятся стохастические уравнения в
форме Лагранжа, Гамильтона и Биркгофа в предположении, что случайные возмущающие силы принадлежат классу процессов с независимыми приращениями.
Литература
1. Еругин Н. П. Построение всего множества систем дифференциальных уравнений, имеющих заданную интегральную кривую // ПММ. — 1952. — Т. 16, вып. 6. — С. 659-670.
2. Построение систем программного движения / А. С. Галиуллин, И. А. Муха-метзянов, Р. Г. Мухарлямов, В. Д. Фурасов. — М.: Наука, 1972.
3. Галиуллин А. С. Построение поля сил по заданному семейству траекторий // Дифференциальные уравнения. — 1981. — Т. XVII, вып. 8. — С. 1487-1489.
4. Галиуллин А. С. Об определении силовой функции по заданному интегралу уравнений движения // Дифференциальные уравнения. — 1982. — Вып. 5. — С. 744-748.
5. Галиуллин А. С. Методы решения обратных задач динамики. — М.: Наука, 1986.
6. Галиуллин А. С., Мухаметзянов И. А, Мухарлямов Р. Г. Обзор исследований по аналитическому построению систем программного движения // Вестник Российского университета дружбы народов. Серия прикладная математика и информатика. — 1994. — Вып. 1. — С. 5-21.
7. Мухаметзянов И. А., Мухарлямов Р. Г. Уравнения программных движений. — М.: Изд-во УДН, 1986.
8. Мухарлямов Р. Г. О построении систем дифференциальных уравнений движения механических систем // Дифференциальные уравнения. — 2003. — Т. 39, вып. 3. — С. 343-353.
9. Тлеубергенов М. И. Об обратной задаче динамики при наличии случайных возмущений // Известия МН-АН РК. Серия физико-математическая. Алма-ты. — 1998. — Вып. 3. — С. 55-61.
10. Тлеубергенов М. И. Об обратной задаче восстановления стохастических дифференциальных систем // Дифференциальные уравнения. — 2001. — Т. 37, вып. 5. — С. 714-716.
11. Тлеубергенов М. И. Об обратной стохастической задаче замыкания // Доклады МН-АН РК. Алматы. — 1999. — Вып. 1. — С. 53-60.
12. Туладхар Б. М. Построение уравнений в форме Лагранжа, Гамильтона и Биркгофа по заданным свойствам движения, Автореф. на соиск... к.ф.-м.н. — М.: Изд-во УДН, 1983. — 11 с.
13. Тлеубергенов М. И. Обратные задачи стохастических дифференциальных систем, Автореф. на соиск... д.ф.-м.н. — Алматы: Ин-т математики, 1999. — 33 с.
14. Тлеубергенов М. И., Ажымбаев Д. Т. О построении дифференциального уравнения по заданным свойствам движения при наличии случайных возмущений // Известия НАН РК, Серия физико-математическая. Алматы. — 2007. — Вып. 3. — С. 15-20.
15. Тлеубергенов М. И., Ажымбаев Д. Т. О построении множества стохастических дифференциальных уравнений по заданному интегральному многообразию, не зависящему от скоростей // Украинский математический журнал, Киев. — 2010. — Т. 62, вып. 7. — С. 1002-1009.
16. Пугачев В. С., Синицын И. Н. Стохастические дифференциальные системы. Анализ и фильтрация. — М.: Наука, 1990.
17. Сагиров П. Стохастические методы в динамике спутников // Механика. Период. сб. переводов иностр. статей. М. — 1974. — Вып. 5(147), 6(148). — С. 28-47, 3-38.
18. Santilli R. M. Foundations of Theoretical Mechanics.1. The Inverse Problem in Newtonian Mechanics. — New-York: Springer-Verlag, 1978.
UDC 531.31; 519.21
To the Problem of Closure of Differential Systems with Degenerating Diffusion
M. I. Tleubergenov, D. T. Azhymbaev
Laboratory of Dynamic Systems Institute of Mathematics 125, Pushkin str., Almaty, 050010, Kazakhstan
It is built the Lagrange, Hamilton and Birkhoff equations by properties of movement in the class of stochastic differential Ito equations in the presence of random perturbing forces from the class of processes with independent increments. The received results are illustrated on an example of movement of Earth's satellite under the action of gravitation's and aerodynamics forces.
Key words and phrases: differential equations, stochastic equations, integral manifold, function, probability, vector.
