СХЕМА ВЫЧИСЛЕНИЯ КОВАРИАЦИОННЫХ ФУНКЦИЙ ВЕКТОРОВ СОСТОЯНИЯ НЕСТАЦИОНАРНЫХ ЛИНЕЙНЫХ СТОХАСТИЧЕСКИХ ДИФФЕРЕНЦИАЛЬНЫХ СИСТЕМ С ЗАПАЗДЫВАНИЕМ Текст научной статьи по специальности «Математика»

2019

ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА

Математика. Информатика. Механика

Вып. 2 (45)

УДК 519.2

Схема вычисления ковариационных функций векторов состояния нестационарных линейных стохастических дифференциальных систем с запаздыванием

И. Е. Полосков

Пермский государственный национальный исследовательский университет Россия, 614990, г. Пермь, ул. Букирева, 15 polosk@psu.ru; тел. (342) 239-65-60

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

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

Б01: 10.17072/1993-0550-2019-2-36-45

Введение

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

©Полосков И.Е., 2019

Среди таких примеров отметим [1,2]: 1°. Хлопки зрителей, сидящих на футбольном стадионе. Несмотря на то, что они хлопают в ладоши синхронно, те, кто сидит в противоположных концах стадиона, не слышат синхронного с их собственным хлопания из-за задержки распространения звука с одного конца стадиона до другого.

2°. Течения Эль-Ниньо/Ла-Нинья или Южная осцилляция, где запаздывающая обратная связь представляет собой эффект океанских волн: распространяющиеся на запад волны Россби на океаническом слое

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

3°. Работа оператора системы управления, наблюдающего за состоянием системы и корректирующего параметры функционирования системы на основе поступающих ДЭ.ННЫХ« Поскольку эти корректировки никогда не могут быть выполнены мгновенно, между наблюдением и управляющим действием возникает запаздывание.

4°. Задержки соединений в сетях связанных динамических систем из-за запаздывания потока информации между ними, которые могут привести к появлению множества новых типов поведения: бифуркации с фазовым переворотом, бифуркации Неймарка-Сакера и др. Было доказано, что запаздывания в соединениях улучшают синхронизацию сетей.

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

Обилие приложений стимулирует бурное развитие теории дифференциальных уравнений с запаздыванием (ДУсЗ), которые являются основными математическими моделями систем с ретардацией и представляют собой важнейший подкласс привлекающего в последнее время значительный и все увеличивающийся интерес с различных точек зрения класса функционально-дифференциальных уравнений (ФДУ) и различных частных форм ФДУ, таких как дифференциальные уравнения с отклоняющимся аргументом обыкновенные и в частных производных, в том числе дифференциальные уравнения с запаздыванием (дифференциально-разностные, ДРУ) и нейтральные, а также интегро-дифференциальные уравнения [1,3-16] и др., а сама теория ДУсЗ в настоящее время принадлежит к числу наиболее быстро развивающихся разделов теории ДУ.

В последнее время значительный интерес уделяется стохастическим ФДУ (СФ-ДУ) различных типов [17-28], разработки

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

(вероятностных характеристик, в основном, моментных функций) приближенных численных алгоритмов анализа систем СФ-ДУ.

К настоящему времени разработан достаточно широкий класс методов решения детерминированных ФДУ [29,30], а также их программных реализаций, которые включают и процедуры из математических пакетов Maple [31], Matlab [32,33] и Mathematica [34]. На основе этих процедур и методов численного решения стохастических дифференциальных уравнений (СДУ) [35-38] интенсивно развиваются приближенные алгоритмы прямого численного интегрирования ('ФДУ различных классов [39-42] в связке с процедурами статистического моделирования (метод Монте-Карло) [43,44] и др.

Важной теоретической и практической задачей является построение схем расчета ковариационных функций векторов состояния систем С ФДУ различных типов. Даже для простых линейных систем ОДУ с одним постоянным запаздыванием использование схем, не связанных с применением статистического моделирования, далеко не прост (см., например, алгоритм из работы [45]).

В статье описана методика и алгоритм исследования во временной области систем линейных стохастических обыкновенНЫХ.

дифференциальных уравнений (СОДУ) с переменными коэффициентами и постоянным запаздыванием, предназначенные для получения ковариационных функций векторов состояния. В предположении, что стохастические возмущения являются гауссовски-ми нестационарными случайными шумами, предложена схема, состоящая из трех этапов. Первый, основанный на классическом методе шагов и расширении пространства состояний (см. [46-48]), позволяет получить

цепочку систем линейных СОДУ без запаздывания; второй, использующий корреляци-оннную теорию, дает возможность построить системы линейных ОДУ для векторных функций математического ожидания и матричных функций ковариации, и, наконец, третий этап дает требуемый результат.

1. Постановка задачи

Рассмотрим следующую систему линейных СОДУ в смысле Стратоновича с запаздыванием (в нестрогой форме):

X(t) = Ao(t) X(t) + co(t) + Go(t) V(t),

to <t < to + t, (1.1)

X (to) = X0, (1.2)

X(t) = A(t) X(t) + Q(t) X(t - t) + + c(t)+ G (t) V (t), to + T<t, (1.3)

где t - время; to и T - начальное и конечное время, to <t ^ T < t > 0 - постоянное запаздывание;

X(t) = (x1(t),x2(t),...,xn(t))J e Rn

- вектор состояния. Начальное условие

x o = (xo,x2o,...,xn )J

- нормальный случайный вектор второго

порядка с плотностью вероятности p°(x),

x e Rn. Его среднее значение равно mX = Xo

DXX = ш[(Xo - mX)(Xo - mX)J] ,

т.е. p o(x) = N (x; mX, DX x). Вход

{V(t) = (Vi(t), V2(t),..., Vm(t))J, t > to}

- m-мерный гауссовскпм случайный белый

o

X

из независимых компонент, так, что

M[V (t)] = 0, rn[v(t) vJ(t')] = Im s(t -1').

В уравнениях (1.1), (1-3) Ao(t) = {aoij(t)}, A(t) = {aij(t)}, Q(t) = {qij(t)}, Go(t) =

{g ij (t)} G (t) = {g ij (t)} - детерминированные вещественные матрицы, co(t) =

{coi(t)} и c(t) = {ci(t)} - детерминированные действительные векторные функции в Мга. Предполагается, что элементы всех этих матриц и векторов являются непрерывными

t

лами Ш[...], Im, ¿(...) и т обозначены оператор математического ожидания, единичная матрица т-го порядка, дельта-функция Дирака и символ транспонирования матрицы соответственно, а точка сверху указывает на производную по времени.

Заметим, что определение вектора состояния на начальном множестве [to, to + т] в форме X (t) = X0 или X (t) = ф(t), где ф({) - заданная детерминированная векторная функция, легко сводится к формулировке задачи, приведенной выше.

Условия существования и единственности для сильных и слабых решений линейных и нелинейных стохастических систем с запаздыванием рассматриваются в

[21, 22, 26, 27]. Поскольку уравнения (1.1),

Xo

ляется нормальным случайным вектором второго порядка, то в соответствии с [26], {X(t), t G [to,T]}, т.е. решение уравнений (1.1), (1.3), является гауссовским непрерывным в среднем квадратичном случайным процессом второго порядка.

Если взглянуть на уравнения (1.1), (1-3) с точки зрения теории случайных процессов, то можно сделать вывод, что вследствие наличия запаздывания случайный векторный процесс X (t), являющийся решением этих уравнений, в общем случае не будет случайным марковским векторным процессом [21]. Поэтому классическая корреляционная теория и известный аналитический аппарат теории марковских процессов, основанный на использовании уравнений Фоккера-Планка-Колмогорова (ФПК-уравнения) для плотностей вероятности p(x, t) и переходных плотностей n(x,t| y, в) векторов состояний линейных или нелинейных динамических систем, не может использоваться для получения статистических характеристик X (t), таких как функции математического ожидания, ковариационные функции и т.д.

С другой стороны, попытки построения других математических моделей явлений, описываемых такими уравнениями, на осно-

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

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

тх(*) = {шхг(*)} = Ш[Х(*)],

а также матричных функций

Ъхх(*) = {Щ(*)} = ж\Х(г) Xт(*)

C XX (Í1,Í2) = {Cij (^)}=Ж [x(ti) X T(Í2) для всех t, ti and t2 in (to, T], где

X(t) = X(t) - mx(t).

Поэтому предлагаемая схема состоит из трех этапов: 1) построение цепочки систем линейных СОДУ без запаздывания; 2) расчет зависимости mx (t) и Dxx (t) от времени; 3) вычисление Cxx(ti,t2)-

2. Уравнения для расширенного вектора состояния и первых моментных функций

Чтобы получить ОДУ для необходимых статистических характеристик вектора X (t) для любого t > to, расширим пространство состояний исследуемой системы путем преобразования немарковской системы в марковскую. Чтобы сделать такое преобразование, воспользуемся следующими обозначениями:

s £ [0,т], tk = to + кт, к = 0,1, 2,..., N + 1,

tN+1 ^ T, sk = s + tk, Afc = (tk, tk+i], Xk(s) = X(s + tk), Xk(0) = Xk_i(T) п.н., Wk(s) = W(s + tk), Wk(0) = Wk_i(T) п.н., Y(s) = Xo, Zo(s) = col (Y(s),Xo(s)), Zi(s) =col (Y(s), Xo(s), Xi(s)) =

= col (Zo(s), Xi(s)), Z2(s) = col (Zi(s),X2(s)), ... Zn(s) = col (Zn_i(s), XN(s)).

Теперь используя пошаговую процедуру, построим цепочку систем линейных СОДУ без запаздывания для расширенных векторных случайных процессов %о (в), % 1 (в), %2(в), ..., (в), принадлежащих семейству

вложенных пространств состоянии

С

R

3n

...

Rn(N+2) ^

0°. Рассмотрим полуинтервал До- Компоненты случайного векторного процесса %о(в), определенного на До, удовлетворяют следующей системе СОДУ:

Г(в) = Оп,

Xо(в) = Ло(во) Xо(в) + Со(во) +

+ Яо(8о) Уо(в), (2.1)

У (0) = Xо, X о(0) = Xо

(здесь и далее точка над символом обозначает производную по в).

Если объединить уравнения (2.1) в виде

Zo(s) = ^+(s) Zo(s)+ co+(s) +

+ G+(s) o(s), Zo(0) = col (Xo, Xo),

(2.2) (2.3)

A+(s)

Or.

n

d

°nxn

^o(so)

co+(s) =

On

co(so)

G+(s) =

со

°nxm

Go(so).

где Ои и - нулевые вектор и матрица

соответствующей размерности, и сравнить структуру систем СОДУ (2.2), (2.3) и (А.1), (А.2), можно заключить, что эти системы будут одинаковы в случае, когда

г = 2 п, д = т, в* = 0, в* = т,

Л(в) = л+(в), с(в) = с+(в), д(в) = ^0+(в).

Кроме того, этот факт дает нам возможность построить явную форму уравнений (А.З), (А.4), определить начальные условия для mz0 (в) и Dz0Zo (в) в форме

mZo (0) DZ0Z0 (0) =

m

x

m

x

oo Dxx Dxx

oo Dxx Dxx

n

соответственно, а затем вычислить mz0 (в) как ш■z(в) и Dz0Zo(в) как Dzz(в), причем рассчитывать mz0 (в) и Dz0z0 (в) как функции в можно как одновременно, так и последовательно, что одновременно даст тх (¿) и ^хх(¿) для ¿о <1 ^ ¿ь

1°. Теперь обратимся к сегментам До и А\. Система СОДУ для вычисления расширенного вектора состояния Zг(в) может быть представлена в следующем виде:

Г(в) = Оп,

Xо (в) = До (во) Xо (в) + со(во) +

+ О о (во) V о(в), X1 (в) = А(вг) Xг(1) + а(вг) Хо(в)+ (2.4)

+ с(вг) + Я(вг) Уг(в), У(0) = Хо, Xо(0) = Xо, X 1(0) = Хо(т).

Действуя как на шаге 0°, мы можем написать краткую форму уравнений (2.4):

^ 1 (в) = А+(в) Zг(в) + с + (в)+

+ (в) V/ +(в), (2.5)

Zг(0) = со1 {Xо(0), Xо(т)), (2.6)

А+(в) =

п

Опхп Опхп Опхп Опхп Ао(во) Опхп

Оп Я(вг) А(вг)

с+(в) =

Оп

Со(во) с(в1)

V+(в)

+ —

О+(в)

п

п

V о(в)

V г(в)

пхт Опхт °о(во) О пхт Опхт О(в1) _

Как и выше, сравниваем структуры уравнений (2.5), (2.6) и (А.1), (А.2) и находим, что эти системы будут идентичны, если

г = 3 п, д = 2 т, в* = 0, в* = т, А(в) = А+(в), с(в) = с+(в), О(в) = д+(в).

Из этого следует, что и в этом случае мы имеем возможность построить подобную шагу 0° форму систем ОДУ на основе уравнений (А.З), (А.4), зафиксировать начальные

условия для mz1 (в) и Dz1z1 (в) в виде:

mz1 (0)

т

х

т

х

(0) =

тХо (т)_

Ъухо (т) Ъухо (т) Ъх0у(т) Ъх0у(т) Ъх0х0 (т)

о

ихх о

ихх

о

ихх о

ихх

и вычислить mz1 (в) как mz(в) и Dz1z1 (в) как Dzz(в), чт° позволит одномоментно получить тх(¿) и ^хх(¿) как тх1 (в) и Х>х1х1 (в) соответственно для ¿1 <Ъ ^ ¿2-

№. На последнем шаге рассматриваем временные полуинтервалы Д^ Д^.., Дм и строим систему СОДУ для вектора состояния Zм (в) в следующей форме:

Г(в) = Оп,

X о (в) = Ао(во) X о(в) + со(во)+ + ОоЫ Vо(в),

XI (в) = А(в1) X 1(1) + а(в1) Хо(в) + + с(в1) + О в) У1(в),

... ... ... ... ...

Xм(в) = А(вм) Xм(¿) + я.(вм) Хм-1(в) + + с(вм) + О (вм) Ум (в),

У (0) = Xо, X о(0) = Xо,

X1 (0) = Xо(т), ... Xм(0) = Xм-1(т).

Дальнейшие действия аналогичны тем, которые выполнены в рамках шагов 0°, 1°, ..., (N-1)°:

1) запись уравнений (2.7) в сжатой форме:

^м (в) = Ам (в) Zм (в) + с+(в) +

+ ом (в) ум (в), (2.8)

Zм(0)=со1 ^м-1(0), Xм-1(т)), (2.9)

где используемые матрицы и векторы рассчитываются так:

ам (в) =

Ам~ 1(в)

о

см+(в) =

(м +1)пхп

[Оп мп Я(вм)] А(вм)

™+-1(в)

™ м (в)

сМ-1(в) с(вм)

™ м (в)

дN (в) =

-1(в) °(М+1)пу °пхИт д(вМ)

2) сравнение уравнений (2.8), (2.9) с уравнениями (А.1), (А.2) и вывод о том, что эти системы будут подобны, если

г = (Ж + 2) п, д = (Ж + 1) т, в* = 0, в* = т, Л(в) = Л+ (в), с(в) = с+(в), д (в) = д + (в);

3) построение явной формы уравнений

(А.З), (А.4);

4) установка начальных условий для mzN (в) и DzNzN (в) в виде:

mzN (0)

mZN-1 (0) mXN-1 (т)

Т>

ZN ZN

(0) =

DzN-1ZN-1 (0) V

Т>

[21]

ZN ZN

[12]

ZN ZN

(0)

(0) DXN-1 (Т )

[12] (0) = ZN ZN '

DYXN-1 (т ) DYXN-1 (т ) DXoXN-1 (т )

VXN-2XN-1 (т ).

(0)= -1Y(т) VXN-1Y(т)

DXN-lXо (т) ••• DXN-lXN-2 (т

5) расчет mzN (в) как mz (в) я DzNzN (в) как Dzz(в), чт° позволит одномоментно получить mх (Ъ) и &хх (Ъ) как mхN (в) и DхNхN (в) соответственно для Ъм-1 < Ъ ^

3. Уравнения для ковариационных функций

Расчет матрицы ковариационных функций Схх(£,£') представляет собой более сложную задачу. Дело в том, что наличие запаздывания приводит к необходимости использования нестандартных инструментов для вычисления значений Схх(£,£') в верхних треугольниках, примыкающих к прямой Ъ = Ъ' при переходе от области к области по горизонтали, и в нижних треугольниках, примыкающих к прямой Ъ = Ъ' при перехо-

де от области к области по вертикали (см. рис.1). Получение же значений Схх(¿,¿0 в остальных областях хорошо подстраивается под схему, представленную в предыдущем разделе.

Ниже опишем необходимую пошаговую процедуру, но сначала укажем необходимые соотношения для вычисления Схх(£,£') в диагональных областях (рис.1). Стандартные ОДУ для вычисления этой матрицы в нижнем и верхнем треугольниках при движении по горизонтали и вертикали соответ-

ственно на основе применения корреляционной теории [49,50] можно представить так:

д C X0X0 (t,t)

dt

to ^ t' <t ^ ti

д C XiXi (t,t')

= Ao(t) CX0X0(t,t'), (3.1)

dt

= A(t) CXiXi (t,t),

(3.2)

ti < t'<t < ti+i, i = 1, 2,...,N; d C X0X0 (t,t')

dt'

t0 ^ t < t' ^ t1; d C XiXi (t,t')

= CX0X0(t,t') AO(t'), (3.3)

dt

= CXiXi(t,t') AT(t'), (3.4)

где Ro(t,t'), R(t, t') - матричные функции Грина, удовлетворяющие уравнениям

dRo(t,t')

dt

dR(t, t')

Ao(t) Ro(t,t'), (3.11)

= A(t) R(t,t'), (3.12)

R(t',t')

дЬ

Ъо(1/,1/ )= I, ЩЪ^ )= I; Но М) = По^в) до(в),

тт М) = щ^в) д(в).

С учетом этого в остальных областях Схх(Ъ, Ъ )

так:

dCXiXj (t,t') dt

= A(t) C (t,t')+

и < КЪ < ^+1, г = 1,

К этим уравнениям необходимо добавить начальные условия в следующей форме:

CхiXi ^М' )= Ъхх (Ъ), (3.5) С хх (Ь,Ь)= Ъхх (Ъ), (3.6)

г = 0,1,

Теперь обратимся к формулам для вычисления матрицы Схх (Ъ^') в нижнем и верхнем треугольниках, но при движении по вертикали и горизонтали соответственно на основе применения корреляционной теории [49], которые можно записать так:

Схохо (Ъ Ъ) = Щ(Ъ Ъо) Ъхохо (Ъ) ЩЦ Ъо) + + / Я-оМ) ЩЩТ(Ъ', в) йв, Ъ > Ъ, (3.7)

■По

СXX ^) = Щ^Ь) Ъх1 х (Ъ) П1^^ + [ ТТ(Ъ,в) П1^'^) йв, Ь>Ь', (3.8)

■По

Схохо (Ъ Ъ) = ЩЪ, Ъо) Ъхохо (Ъ) Я1 (Ъ Ъо) + + Щ(Ъ,в) тго(Ъ,в) йв, КЪ', (3.9)

■По

С XX ^) = R(Ь,Ьо) Ъх1 х (Ъ) Щ^^оК + / щ (Ъ,в) тт Т(Ъ', в) йв, кЪ, (3.10)

■По

г = 1,2,..,Ж,

+ G (t) Cx- Xj (t,t'), t>t', 0 < j <i < N; dCXiXj (t,t')

(3.13)

dt = CXix (t,t') AT(t')+

+ CXiXj-i(t,t) GT(t'), t<t', 0 < i < j < N, (3.14)

ti ^ t ^ ti+1, tj ^ t' ^ tj+1.

Заключение

В работе была представлена многошаговая трехэтапная схема расчета моментных функций векторов состояния систем линейных стохастических обыкновенных дифференциальных уравнений с переменными коэффициентами и одним постоянным запаздыванием. Основное предназначение схемы - получение матриц ковариационных функций для этих векторов. Основной проблемой для практической реализации метода является расчет ковариационных функций по формулам (3.7)—(3.10) в небольшом числе подобластей, где невозможно применять стандартную форму уравнений для этих функций. В рамках вычисления присутствующих в формулах (3.7)—(3.10) интегралов можно рекомендовать аппроксимацию матричной функции Грина итерационным методом [51].

Несложно установить, что изложенный алгоритм может быть перенесен на область анализа подобных систем, но с кратными постоянными запаздываниями.

Приложение А. Общие ОДУ для первых моментных функций

Построим общую систему ОДУ, управляющую во времени поведением векторной функции математических ожиданий

mz(в) = {тг(в)} = Ш(в)],

матрицы функций ковариации

Dzz(в) = {Щ(в)} = М [¿(в) г\в)

и матрицы ковариационных функций

Czz(в, в') = {Сгз(в, в')} = Ж [¿(в) г\в')

вектора состояния ^(в) = {Д¿(в)} € Кг, для которого ¿(в) = Z(в) — mz(в), линейной стохастической дифференциальной системы с аддитивными шумами в форме Стратоно-вича:

¿(в) = А(в) 2(в) + с(в) + О(в) V(в),

в* < в < в*, (А.1)

Z (в*)= £ *, (А.2)

где основная часть обозначений аналогична введенным в первых разделах статьи, а точка сверху обозначает производную по в. Пусть начальный вектор^^ = (Д *1,Д ■■■, Д является случайным и имеет математическое ожидание m*z = ® [Ж *] и матрицу ко-

вариаций D*zz = ® ¿Т .Вход {V (в) = {V ¿(в)}, в* ^ в ^ в*}- д-мерный векторный гауссовский белый шум с независимыми компонентами. В уравнениях (А.1) А(в) = {а^(в)}, О(в) = {д ..(в)} и с(в) = {сг(в)} -детерминированные матрицы и вектор соответственно. Предполагается, что элементы этих матриц и вектора являются непрерыв-

в

Применяя корреляционную теорию (см., например, [49, 50]), мы можем получить необходимые системы ОДУ для mz (в) И Dzz (в) в следующей форме:

т z (в) = А(в) mz (в)+ с(в), (А.з)

Рzz(в) = А(в) Dzz(в) + [А(в) Dzz(в)]т + + О(в) От(в), (А.4)

в* < в' <€ в*

с начальными условиями

mz (s*) = m*z,

Dzz (s*) = D*zz соответственно.

(A.5) (A.6)

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

1. Эльсгольц Л.Э., Норкин С. Б. Введение в теорию дифференциальных уравнений с отклоняющимся аргументом. М.: Наука, 1971. 296 с.

2. Lakshmanan М., Senthilkumar D.V. Dynamics of nonlinear time-delay systems. Berlin, Heidelberg: Springer, 2010. XVII, 313 p.

3. Азбелев H.B., Максимов В.П., Рахматул-лина Л.Ф. Элементы современной теории функционально-дифференциальных уравнений. Методы и приложения. М.: Институт компьютерных исследований, 2002. 384 с.

4. Мышкис А.Д. Общая теория дифференциальных уравнений с запаздывающим аргументом // Успехи математических наук. 1949. Т. 4, вып. 5 (33). С. 99-141.

5. Рубапик В. П. Колебания квазилинейных систем с запаздыванием. М.: Наука, 1969. 288 с.

6. Baker С. Т.Н. Retarded differential equations // Journal of Computational and Applied Mathematics. 2000. Vol. 125, № 1-2. P. 309-335.

7. Breda D., Maset S., Vermiglio R. Stability of linear delay differential equations: A numerical approach with MATLAB. New York: Springer, 2015. XI, 158 p.

8. Erneux T. Applied delay differential equations. New York: Springer, 2009. XII, 204 p.

9. Fridman E. Introduction to time-delay systems: Analysis and control. Basel: Birkhauser, 2014. XVIII, 362 p.

10. Guglielmi N. Delay differential equations // Encyclopedia of Applied and Computational Mathematics / B.Engquist (ed.). Berlin, Heidelberg: Springer, 2015. P. 334-338.

11. Hale J.K., Lunel S.M.V. Introduction to functional differential equations. New York: Springer, 1993. X, 447 p.

12. Kim A. V., Ivanov A.V. Systems with delays: Analysis, control, and computations. Hoboken: John Wiley, 2015. 184 p.

13. Kolmanovskii V., Myshkis A. Introduction to the theory and applications of functional differential equations. Dordrecht: Springer, 1999. XVI, 648 p.

14. Lakshmikantham V., Rao M.R. Theory of in-tegro-differential equations. Reading: Gordon and Breach, 1995. 384 p.

15. Smith H. An introduction to delay differential equations with applications to the life sciences. New York: Springer, 2011. XI, 172 p.

16. Wu J. Theory and applications of partial functional differential equations. New York: Springer, 1996. X, 432 p.

17. Лидский Э.А. Об устойчивости движений системы со случайными запаздываниями // Дифференциальные уравнения. 1965. Т. 1, № 1. С. 96-101.

18. Митрополъский Ю.А., АньН.Д. Случайные колебания в квазилинейных системах стохастических дифференциальных уравнений с запаздыванием // Украинский математический журнал. 1986. Т. 38, № 2. С. 181-183.

19. Потапов В.Д. Устойчивость движения стохастической вязкоупругой системы // Прикладная математика и механика. 1993. Т. 57, вып. 3. С. 137-145.

20. Рубапик В. П. Колебания сложных квазилинейных систем с запаздыванием. Ми.: Изд-во "Университетское", 1985. 143 с.

21. Царьков Е.Ф. Системы стохастических дифференциальных уравнений с запаздыванием // Изв. АН Латв. ССР. Сер. физ. и техн. наук. 1966. № 2. С. 57-64.

22. Царьков Е.Ф. Случайные возмущения дифференциально-функциональных уравнений. Рига: Зинатне, 1989. 421 с.

23. Чайковский М.В., Янович Л.А. О численном нахождении корреляционных функций решения систем линейных интегро-дифференциальных уравнений со случайно возмущенной правой частью // Дифференциальные уравнения. 1987. № 2. С. 328-338.

24. Kushner H.J. Numerical methods for controlled stochastic delay systems. Boston: Birkhau-ser, 2008. XIX, 281 p.

25. Lafuerza L.F., Toral R. Stochastic description of delayed systems // Philosophical Transactions of the Royal Society. Ser. A. Mathematical Physical and Engineering Sciences. 2013. Vol. 371, № 1999. DOI: 10.1098/rsta.2012.0458. P. 1-16.

26. Мао X. Stochastic differential equations and applications. 2nd ed. Cambridge, UK: Wood-head Publishing, 2011. XVIII, 422 p.

27. Mohammed S.E.A. Stochastic functional differential equations. Boston, London: Pitman Publishing, 1984. IX, 245 p.

28. Tsokos C.P., Padgett W.J. Random integral equations with applications to stochastic systems // Lecture Notes in Mathematics,

Vol.233. Berlin: Springer-Verlag, 1971. VIII, 176 p.

29. Bellen A., Zennaro M. Numerical methods for delay differential equations. Oxford: Oxford University Press, 2003. XIV, 395 p.

30. Zhang W. Numerical analysis of delay differential and integro-differential equations: A thesis. St. John's, NL, Canada: Memorial University of Newfoundland, 1998. VIII, 138 p.

31. Heck . I. Introduction to Maple. 3rd ed. New-York: Springer, 2003. XVI, 828 p.

32. Sham,pine L.F., Gladwell I., Thompson S. Solving ODEs with Matlab. Cambridge: University Press, 2003. 272 p.

33. Sham,pine L.F., Thompson S. Delay differential equations with constant lags // CODEE Journal. 2012. Vol.9, Art. 10. 5p.

34. Mangano S. Mathematica cookbook. Cambridge: O'Reilly, 2010. XXIV, 800 p.

35. Кузнецов Д.Ф. Стохастические дифференциальные уравнения: теория и практика численного решения. 4-е изд. СПб.: Изд-во Политехи, ун-та, 2010. XXX, 786 с.

36. Kloeden Р.Е., Platen Е., Schurz Н. Numerical solution of SDE through computer experiments. Berlin: Springer, 2003. XIV, 292 p.

37. Milstein G.N., Tretyakov M.V. Stochastic numerics for mathematical physics. Berlin: Springer, 2004. XIX, 594 p.

38. Platen E., Bruti-Libera,ti N. Numerical solution of stochastic differential equations with jumps in finance. Berlin: Springer, 2010. XXVIII, 856 p.

39. Buckwar E. Introduction to the numerical analysis of stochastic delay differential equations // Journal of Computational and Applied Mathematics. 2000. Vol.125, № 1-2. P. 297307.

40. Buckwar E., Kuske R., Mohammed S.-E., Sha-rdlow T. Weak convergence of the Euler scheme for stochastic differential delay equations // London Mathematical Society Journal of Computation and Mathematics. 2008. Vol.11. P. 60-99.

41. Kiichler U., Platen E. Strong discrete time approximation of stochastic differential equations with time delay // Mathematics and Computers in Simulation. 2000. Vol. 54, № 13. P. 189-205.

42. Мао X., Sabanis S. Numerical solutions of stochastic differential delay equations under local Lipschitz condition // Journal of Computational and Applied Mathematics. 2003. Vol. 151. P. 215-227.

43. Кельт,он, В., Лоу А. Имитационное моделирование. Классика CS. 3-е изд. СПб.: Питер; К.: BHV, 2004. 847 с.

44. Allen Е. Modeling with Ito stochastic differential equations. Dordrecht: Springer, 2007. XII, 230 p.

45. Разевиг В.Д. Корреляционный анализ систем с запаздыванием // Автоматика и телемеханика. 1973. № 9. С. 42-48.

46. Полосков И.Е. Расширение фазового пространства в задачах анализа дифференциально-разностных систем со случайным входом // Автоматика и телемеханика. 2002. № 9. С. 58-73.

47. Полосков И.Е. Численно-аналитические схемы анализа динамических систем с последействием // Вестник Пермского ун-та. Математика. Механика. Информатика. 2011. Вып. 2 (6). С. 51-58.

48. Poloskov I.E Numerical and analytical methods of study of stochastic systems with delay // Journal of Mathematical Sciences. 2018. Vol. 230, № 5. P. 746-750.

49. Астапов Ю.М., Медведев B.C. Статистическая теория систем автоматического регулирования и управления. М.: Наука, 1982. 304 с.

50. Полосков И.Е. Стохастический анализ динамических систем [Электронный ресурс]: монография. Пермь: Изд-во Перм. ун-та, 2016. 772 с.

51. Полосков И.Е. Об одном методе приближенного анализа линейных стохастических интегро-дифференциальных систем // Дифференциальные уравнения. 2005. Т. 41, № 9. С.1276-1279.

The scheme for calculating the eovarianee functions of the state vectors for non-stationary linear stochastic differential systems with delay

I. E. Poloskov

Perm State University; 15, Bukireva st., Perm, 614990, Russia polosk@psu.ru; (342) 239 65 60

The paper describes the method and algorithm of analysis in the time domain of linear systems of stochastic ordinary differential equations (SODE) with variable coefficients and constant delays, designed to obtain the eovarianee functions of state vectors, in the case of assumption that stochastic perturbations are Gaussian non-stationary random noises. To solve the problem, we propose a multi-step scheme. Each of steps consists of three stages. The first based on the classical method of steps, expands the state space which ultimately allows us to obtain a chain of linear SODE systems without delay; on the second step, using the correlation theory, systems of linear ODE are constructed for vector functions of expectation and matrix eovarianee functions with the necessary initial conditions; at the third stage, we apply the correlation theory too and the results of the second stage and achieve the required result of the current step.

Keywords: linear equation; stochastic ordinaty differential equation; delay; moment functions; eovarianee function.