Научная статья на тему 'К вопросу об устойчивости по начальным данным по части переменных решений линейных импульсных дифференциальных уравнений Ито с последействием'

К вопросу об устойчивости по начальным данным по части переменных решений линейных импульсных дифференциальных уравнений Ито с последействием Текст научной статьи по специальности «Математика»

CC BY
127
32
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
УСТОЙЧИВОСТЬ РЕШЕНИЙ ПО ЧАСТИ ПЕРЕМЕННЫХ / ЛИНЕЙНЫЕ ИМПУЛЬСНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЙ ИТО С ПОСЛЕДЕЙСТВИЕМ / МЕТОД ВСПОМОГАТЕЛЬНЫХ ИЛИ МОДЕЛЬНЫХ УРАВНЕНИЙ / SOLUTIONS' STABILITY ON THE PART OF VARIABLES / LINEAR IMPULSIVE DIFFERENTIAL EQUATIONS WITH AFTEREFFECT / METHOD OF THE AUXILIARY OR MODEL EQUATIONS

Аннотация научной статьи по математике, автор научной работы — Кадиев Рамазан Исмаилович

Исследуются вопросы моментной устойчивости решений по части переменных относительно начальных данных для линейных импульсных систем дифференциальных уравнений Ито с последействием. Исследование проводится методом вспомогательных или модельных уравнений. Получены достаточные условия устойчивости в терминах параметров исследуемых систем.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

THE PROBLEM OF STABILITY ACCORDING TO INITIAL DATA ON THE PART OF VARIABLE SOLUTIONS OF LINEAR IMPULSIVE DIFFERENTIAL EQUATIONS WITH AFTEREFFECT

Problems of moment solutions’ stability on the part of variables regarding initial data for linear impulsive differential equations with the consequence are considered here. The research is conducted by the method of auxiliary or model equations. Sufficient stability conditions in terms of parameters of the researched systems are received.

Текст научной работы на тему «К вопросу об устойчивости по начальным данным по части переменных решений линейных импульсных дифференциальных уравнений Ито с последействием»

УДК 517.929.4+519.21

К ВОПРОСУ ОБ УСТОЙЧИВОСТИ ПО НАЧАЛЬНЫМ ДАННЫМ ПО ЧАСТИ ПЕРЕМЕННЫХ РЕШЕНИЙ ЛИНЕЙНЫХ ИМПУЛЬСНЫХ

ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ИТО С ПОСЛЕДЕЙСТВИЕМ

© Р.И. Кадиев

Ключевые слова: устойчивость решений по части переменных; линейные импульсные дифференциальные уравнений Ито с последействием; метод вспомогательных или модельных уравнений.

Исследуются вопросы моментной устойчивости решений по части переменных относительно начальных данных для линейных импульсных систем дифференциальных уравнений Ито с последействием. Исследование проводится методом вспомогательных или модельных уравнений. Получены достаточные условия устойчивости в терминах параметров исследуемых систем.

Импульсные дифференциальные уравнения Ито с последействием являются хорошей математической моделью для финансовых процессов. Некоторым вопросам устойчивости решений для линейных импульсных систем дифференциальных уравнений Ито с последействием посвящены работы [1—4]. Вопросы устойчивости решений по части переменных относительно начальных данных для линейных импульсных систем дифференциальных уравнений Ито с последействием, по-видимому, ранее не исследовались. Главной целью исследований было развитие метода вспомогательных уравнений применительно к исследованию вопросов устойчивости по начальным данным по части переменным решений линейных импульсных систем дифференциальных уравнений Ито с последействием. Вопросы устойчивости по части переменных для линейных стохастических функционально-дифференциальных уравнений методом вспомогательных уравнений были изучены в работах [5], [6].

Пусть: (Q, F, (F)t^o,P) —стохастический базис; kn —линейное пространство n-мерных Fo — измеримых случайных величин; Bi,i = 2, ...,m — скалярные независимые стандартные винеровские процессы; 1 ^ p < ж ; E —символ математического ожидания; |.| — норма в Rn ; ||.|| —норма l х n -матрицы, согласованная с нормой в Rn ; ß —мера Лебега на [0, ж) .

Рассматривается линейная импульсная система дифференциальных уравнений Ито вида

mi m mi

dx(t) = Aij(t)x(hij(t))dt + Y^ Aij(t)x(hij(t))dBi(t) (t ^ 0), (1)

j=0 i=2 j=0

x(v) = p(v) (v < 0), (2)

x(ßj) = Ajx(ßj — 0), j = 1, 2, 3,... почти наверно (п.н.), (3)

где ßj, Aj ,j = 1,2,3,... — действительные числа такие, что 0 = ß0 < ßl < ß2 < ..., lim ßj = ж ; p — случайный процесс, который не зависит от винеровских процессов

Bi, i = 2,..., m , и имеет п.н. ограниченные в существенном траектории; Aj, i = 1,..., m,j = = 0,..., mi — n х n -матрицы, элементы матриц Aij,j = 0,..., mi — прогрессивно измеримые случайные процессы, траектории которых п.н. локально суммируемы, элементы матриц Aij,i = 2,...,m,j = 0,...,mi — прогрессивно измеримые случайные процессы, траектории

которых п.н. локально суммируемы с квадратом; Ъ^ , г = = 0, ...,Шг — изме-

римые по Лебегу функции такие, что Ъ^(Ь) ^ Ь при Ь € [0, то) ^-почти всюду (п.в.), г = 1,..., т, ^ = 0,..., т^.

Ведем следующие обозначения

(Я )(Ь) = / х(Ъ(Ь)), если Ъ(Ь) ^ 0 ^ — п.в., (,) = ( 0, если Ъ(Ь) ^ 0 ^ — п.в.,

(Я*х)(Ь) = \ 0,еслиЪ(Ь) < 0 ^ — п. в., ^(Ь) = \ р(Ъ(Ь)),еслиЪ(Ь) < 0 ^ — п. в..

В дальнейшем будем пользоваться следующей записью уравнения (1)—(3):

¿ж(Ь) = ((^х)(Ь) + /(Ь) (Ь ^ 0), (4)

(те т1 т2

Е (1 — А)5(Ь — ^)х(Ь) + £ (Ь)^.Х)(Ь)), £ ^(ь)^.х)(ь), ...,

3=0 3=0 3=0

тт \ / т1 т2 тт \

Е Атз(Ь)^х(Ь) , /(Ь) = ^ А13(*)),£ А23(Ь)^(Ь), ..., £ Атз(Ь)^(¿Л ,

3=0 ) \3=0 3=0 3=0 )

2(Ь)= С01(Ь, г?2(Ь),..., Вт(Ь)) .

Пусть: — линейное пространство п —мерных прогрессивно измеримых случайных процессов на [0, +то) , траектории которых п.н. непрерывно справа и имеют пределы слева; Чп(2) —линейное пространство п х ш —матриц на [0, +то) , строки которых ш —мерные прогрессивно измеримые случайные процессы локально интегрируемые по 2 ; Ьп — линейное пространство п —мерных случайных процессов на (—то, 0) , которые не зависят от вине-ровских процессов В^, г = 2,..., ш , и имеет п.н. ограниченные в существенном траектории; 7 : [0, +то) ^ Е1 — положительная непрерывная функция; I — некоторое фиксированное натуральное число, удовлетворяющее неравенству 1 ^ I < п; для любого х € Пп введем обозначения у = со1(х!,..., х1) и Ъ = со1(х1+1,..., хп).

Введем следующие обозначения линейных нормированных пространств:

МР = (х : х € Бк, ||х|М ^вир (Е|7(Ь)х(Ь)|р)1/р < то! (МР = Мр);

i р г^ю j

к'п = {а : а € кп, ||а||к™ (Е|а|р)1/р < то} ;

ч; = € £п,|М|ч „га18ИР < то}.

I ^<0 J

Известно, что через любое х(0) € кп проходит единственное решение уравнения (1)—(3) (с точностью до Р —эквивалентности). Обозначим это решение через х(Ь, х(0), . Отметим, что уравнение (1)—(3) называют однородным, если ) = 0 (V < 0) . Через Xf (Ь,х(0)) обозначим решение уравнения (4) такой, что Xf (0,х(0)) = х(0) . Тогда для уравнения (4) также справедлива

Лемма. Для решения Xf (Ь, х(0)) уравнения (4), где х(0) € кп задает начальное условие для этого решения, имеет место представление

Xf (Ь, х(0)) = X(Ь)х(0) + (С/)(Ь) (Ь ^ 0), (6)

где X(Ь)(Х(0) = Е — единичная матрица) — п х п -матрица, столбцами которой являются решения однородного уравнения (4) (фундаментальная матрица), а С : Чп(2) ^ Пп — линейный оператор (оператор Коши) такой, что С/ — 'решение уравнения (4), удовлетворяющее условию (С/)(0) = 0 .

Определение 1. Тривиальное решение однородного уравнения (1)-(3) x(t, 0, 0) = 0 называют:

— p -устойчивым относительно первых l компонент, если для любого е > 0 найдется такое 5(е) > 0, что при любых ж(0) € kn, f € L^ и ||ж(0)||кп + ||f||ln < ¿(е) будет выполнено неравенство (E|y(t, ж(0), f)|p)1/p ^ е для любого t ^ 0 ;

— асимптотически p -устойчивым относительно первых l компонент, если оно p -устойчиво, и, кроме того, для любых ж(0) € кП, f € L^ и ||ж((0)||кп + ||f||ln < ¿(е) будет

lim (E|y(t, ж(0), f)|p)1/p = 0 ;

— экспоненциально p -устойчивым относительно первых l компонент, если существуют некоторые положительные числа с, ß такие, что для любых ж(0) € fcp , f € Ln справедливо неравенство (E|y(t,ж(0), f)|p)1/p < c(||x(0)||fc™ + ||f||L™)exp{-£t} .

Определение 2. Уравнение (1)-(3) назовем Mp y -устойчивым, если для любых ж(0) € kpn , f € L^ для решения уравнения (1)-(3) x(.,x(0),f) имеем y(.,x(0),f) € Mp и выполнено неравенство

||у(.,ж(0), f)|MY < c(l|x(0)||fci? + ||f Nl-) (7)

для некоторого положительного числа с .

Очевидно, что:

— из Mpу -устойчивости уравнения (1)-(3) следует p -устойчивость тривиального решения однородного уравнения (1)-(3) относительно первых l компонент;

— из Mp у -устойчивости уравнения (1)-(3) (где Y(t) ^ 5 > 0 (t ^ 0) и lim 7(t) = следует асимптотическая p -устойчивость тривиального решения однородного уравнения (1)-(3) относительно первых l компонент;

— из Mp у -устойчивости уравнения (1)-(3) (где 7(t) = exp{^t} , ß — некоторое положительное число) следует экспоненциальная p -устойчивость тривиального решения однородного уравнения (1)-(3) относительно первых l компонент.

Так как ж = col(y, h) и Dn = Dl х Dn-1 , то уравнение (4) эквивалентно системе вида

dy(t) = [(V1y)(t) + (V2h)(t) + fy(t)]dZ(t) (t ^ 0), dh(t) = [(V3y)(t) + (V4h)(t) + fh(t)]dZ(t) (t ^ 0),

hmi^m ^ (8)

где VI : Б1 ^ ) , : Бп-1 ^ Ь1 (Я) , : Б1 ^ Ьп-1(^) , : Бп-1 ^ Ьп-1 (Я) суть линейные операторы, определяемые оператором V, /у € Ьг(^) , /н € Ьп-1 (^) , / = = со1(/у,/Л).

В силу того, что через любое ж(0) € кп проходит единственное решение уравнения (4), каждое из уравнений системы (8) в отдельности будеть иметь единственное решение при любых фиксированных у(0) € к1, Н € Пп-1 и Н(0) € кп-1, у € Б1 соответственно. Тогда, в силу леммы, второе уравнение системы (8) эквивалентно уравнению

Н(*) = Н(¿)Н(0) + (С/ + ^3у))(^ (* ^ 0),

где Н —фундаментальная матрица, а С1 — оператор Коши для второго уравнения системы (8). Следовательно, из первого уравнения системы (8) получим

¿у(*) = [№)(*) + ^(НН(0)))(^ + ^С/)(*) + /у(*) (* ^ 0), (9)

где = VI + V2ClVз .

Отсюда следует, что уравнение (1)-(3) Мр у -устойчиво тогда и только тогда, когда при любых ж(0) € кП , р € ЬП решение уравнения (9) у(.,ж(0),р) принадлежит пространству Мр и для него выполнено неравенство (7).

Для установления принадлежности решения уравнения (9) у(.,х(0),р) пространству Мр при любых х(0) € крп , р € ^ и выполнимости для него неравенства (7) воспользуемся Ш —преобразованием, т.е. эквивалентным преобразованием уравнения (9). Для описания Ш —преобразования уравнения (9) рассмотрим модельное уравнение, асимптотические свойства решений которого известны. Пусть модельное уравнение имеет вид

йу(Ь) = [(0У)(Ь) + (Ь) (Ь ^ 0), (10)

где Q : Б1 ^ Ьг(2) — линейный оператор; 2 определен ранее; д € Ь1 (2) . Предполагается, что через любое у(0) € к1 проходит единственное (с точностью до Р—эквивалентности) решение уравнения (9). Тогда, в силу леммы, для решения у этого уравнения имеет место представление у(Ь) = и(Ь)у(0) + (Шд)(Ь), Ь ^ 0 , где и — фундаментальная матрица, Ш — оператор Коши для уравнения (10).

Уравнение (9) при помощью модельного уравнения (10) перепишем в виде

у(Ь) = и(Ь)у(0) + (Ш(V, — Q)y)(t) + (ШУ2(ЯЪ(0)))(Ь) + (Ш(^С/Л + /у))(Ь) (Ь ^ 0).

Обозначив Ш(У5 — Q) = © , получим

((/ — ©)у)(Ь) = и(Ь)у(0) + (ШУ2(НЪ(0)))(Ь) + (Ш(^С/ + /у))(Ь) (Ь ^ 0).

Теорема. Пусть иу(0) + Ш^(ЯЪ(0)) + Ш(^С/ + /у) € М^ для любых х(0) € € кп , р € ьп и | |и у(0) + Ш^(ЯЪ(0)) + + /У)|М < с(||х(0)||к? + |М|^) для

некоторого положительного числа с, а оператор © действует в пространстве Мр . Тогда, если оператор (I — ©) : Мр ^ Мр непрерывно обратим, то уравнение (1)-(3) Мр у -устойчиво.

Пользуясь этой теоремой получены достаточные условия устойчивости решений различных классов уравнений вида (1)—(3) в терминах параметров этих уравнений.

ЛИТЕРАТУРА

1. Кадиев Р.И. Устойчивость решений нелинейных функционально-дифференциальных уравнений с импульсными воздействиями по линейному приближению // Дифференц. уравнения. 2013. Т. 49. № 8. С. 963-970.

2. Kadiev R., Ponosov A. Stability of impulsive stochastic differential linear functional equations with linear delays // J. of Abstract Differential Equations and Applications. 2012. V. 2. № 2. P. 7-25.

3. Кадиев Р.И., Поносов А.В. Устойчивость решений линейных импульсных систем дифференциальных уравнений Ито с ограниченными запаздываниями // Дифференц. уравнения. 2010. Т. 46. № 4. C. 486-498.

4. Кадиев Р.И., Поносов А.В. Устойчивость решений линейных импульсных систем дифференциальных уравнений Ито с последействием // Дифференц. уравнения. 2007. Т. 43. № 7. С. 879-885.

5. Кадиев Р.И. Достаточные условия устойчивости по части переменных линейных стохастических систем с последействием // Изв. вузов. Математика. 2000. № 6. С. 75-79.

6. Кадиев Р.И. Допустимость пар пространств по части переменных для линейных стохастических функционально-дифференциальных уравнений // Изв. вузов. Математика. 1994. № 4. С. 1-10.

Поступила в редакцию 1 июня 2015 г.

Kadiev R.I. THE PROBLEM OF STABILITY ACCORDING TO INITIAL DATA ON THE PART OF VARIABLE SOLUTIONS OF LINEAR IMPULSIVE DIFFERENTIAL EQUATIONS WITH AFTEREFFECT

Problems of moment solutions' stability on the part of variables regarding initial data for linear impulsive differential equations with the consequence are considered here. The research is conducted by

the method of auxiliary or model equations. Sufficient stability conditions in terms of parameters of the researched systems are received.

Key words: solutions' stability on the part of variables; linear impulsive differential equations with aftereffect; method of the auxiliary or model equations.

Кадиев Рамазан Исмаилович, Дагестанский научный центр, Дагестанский государственный университет, г. Махачкала, Российская Федерация, доктор физико-математических наук, профессор кафедры прикладной математики, e-mail: kadiev_r@mail.ru

Kadiev Ramazan Ismailovich, Dagestan Research Center, Dagestan State University, Makhachkala, the Russian Federation, Doctor of Physics and Mathematics, Professor of the Applied Mathematics Department, e-mail: kadiev_r@mail.ru

УДК 517.977.5

ОБ УПРАВЛЕНИИ С ПОВОДЫРЕМ В ЗАДАЧЕ ОПТИМАЛЬНОГО ВЫВЕДЕНИЯ РАКЕТЫ-НОСИТЕЛЯ

© И.Н. Кандоба, И.В. Козьмин, Е.К. Костоусова, В.И. Починский

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

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

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

Рассмотрим задачу оптимального управления нелинейной динамической системой, описывающей движение ракеты-носителя (РН) от точки старта до момента выхода РН на заданную околоземную эллиптическую орбиту. Здесь, в отличие от [4, 5], математическая модель управляемого движения РН включает уравнения поступательного движения его центра масс и уравнения вращательного движения носителя как твердого тела. Двигательная установка РН состоит из одного жестко закрепленного основного двигателя, создающего тягу вдоль оси симметрии РН, и четырех подвижных рулевых двигателей. Движение

i Надоели баннеры? Вы всегда можете отключить рекламу.