Интервально-дифференциальные уравнения и метод их решения

УДК 62-50:519.7/8

ИНТЕРВАЛЬНО-ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ И МЕТОД ИХ РЕШЕНИЯ

© В.И. Левин

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

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

1. ВВЕДЕНИЕ

В различных областях науки и техники очень часто встречаются задачи, для решения которых нужно решить одно уравнение или систему уравнений, содержащих производные искомых функций. Такие уравнения называются дифференциальными уравнениями. Теория и разнообразные методы решения таких уравнений разработаны весьма подробно [1]. При этом всегда предполагается, что все фигурирующие в уравнениях функции являются полностью определенными. Однако на практике задачи, сводящиеся к решению дифференциальных уравнений, весьма часто связаны с исследованием неполностью определенных систем. Поэтому в них могут возникать дифференциальные уравнения иного типа. Фигурирующие в них искомые функции являются неполностью определенными. Соответственно этому и производные искомых функций в дифференциальных уравнениях данного типа оказываются неполностью определенными. Настоящая работа является введением в изучение дифференциальных уравнений именно этого типа.

2. МАТЕМАТИЧЕСКИЕ ОСНОВЫ

Будем использовать в качестве математического аппарата алгебру интервальных чисел [2], интервальный анализ [3] и интервально-дифференциальное исчисление [4].

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

а = [а1,а2] = {а | а1 < а < а2}

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

Таким образом, алгебраические операции над интервалами - сложение, вычитание, умножение, деление -вводятся в виде

а + Ь = {а + Ь | а е а, Ь е Ь}, а - Ь = {а - Ь | а е а, Ь е Ь}, к  а = {к  а | а е а},

а  Ь = {а  Ь | а е а,Ь е Ь}, а / Ь = {а/ Ь | а е а,Ь е Ь}.

Из определений (3) вытекают следующие формулы для вычисления результатов алгебраических операций над интервальными числами

у + Ь = [а1, а2] + [Ьх, Ь2] = [а1 + Ь1, а2 + Ь2], у -~ = [аъа2] -[Ь^Ь2] = [«1 -Ь2,а2 -Ь1],

к  а = к [а1,а2] =

[ка1,ка2], к > 0,

1[ка2,ка1], к < 0, ^2] = [тт(аг-  Ъ. X

и'

а/Ь = [аъ а2 ] /[Ь1, Ь2] = [а1, а2] ^[1/Ь2,1/Ь1].

~  Ь = [а1,а2]  [Ь1,Ь2] = [шт( аг  Ь'),шах(аг  Ь')],

г.' ' г,3 '

Объектом изучения интервального анализа являются интервальные функции [3]. Интервальная функция вводится как однозначное отображение множества замкнутых вещественных интервалов {~} вида (1), т. е. ~ = [д, х2 ], на множество замкнутых вещественных интервалов {у} того же вида у = [у, у2 ]. Символически интервальная функция записывается в виде

а о Ь = {а о Ь | а е а,Ь е Ь } .

у = / (^Х где х = [ х2 ] у = [Уl, У2], / (X) = / X), /2( X)],

где х называется интервальной независимон переменной (интервальным аргументом); у - интервальной

зависимой переменной; / - интервальной функцией; / (•) - нижней граничной функцией интервальной

функции /; / (•) - верхней граничной функцией

интервальной функции / .

Базовым понятием интервального анализа является понятие предела интервальной функции, которое вводится таким образом. Независимая интервальная переменная х = X, х ] интервальной функции (5) по определению неограниченно приближается к некоторому интервалу х0 = [х01, х02], если в процессе этого изменения х неограниченно приближается к х , а х - к х02. Символически это записывается так:

Аналогично определяется неограниченное приближение зависимой интервальной переменной у = [ у, у2 ]

функции (5) к интервалу ~0 = [уо1, уо2 ]:

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

~Нт у = уо

или

_1ип / (X) = Уо .

Если интервальная функция (5) непрерывна, т. е. нижняя У1 и верхняя У2 границы зависимой переменной у - непрерывные функции нижней Х1 и верхней Х2 границ независимой переменной X , то предел функции (5) равен значению функции в предельной точке х0 аргумента X , или символически

_Нш /(X) = /.

Основным математическим понятием, используемым в данной статье, является понятие интервальной производной функции [3-4]. Оно вводится следующим образом, на базе понятия обычной производной функции [1]. Рассмотрим произвольную интервальную

функцию / в виде (5). Будем считать ее непрерывной. Зафиксируем в ней значение независимой переменной X = X0 = [^1, ]. Этому значению, в силу непрерывности функции, соответствует фиксированное значение функции у = / (X)) . Определим приращения незави-

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

^ = X - Xо, Ау = у - Уо = /

и составим отношение второго приращения к первому

Ау / М = (у - у0)/( x - x0) = = (/(X) - / (^ж X - X,,).

Предел отношения (11) при неограниченном приближении независимой переменной X к ее фиксированному предельному значению X0 , если он существует, называется интервальной производной функцией от исходной интервальной функции /(X) (5) в точке X0

и обозначается или /X (X) :

У-„ = /к (X) = Ау /А

где АX и Ау определяются формулами (10).

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

интервальной функции / (5) существует в любой точке X0, в которой функция / является существенно

интервальной по независимой переменной X .

Как и в случае обычной производной [1], понятие интервальной производной (12) может быть обобщено путем повторного выполнения операции взятия производной. При этом интервальная производная у~ из формулы (12) становится производной 1 -го порядка, производная от у- - производной 2-го порядка у- , производная от у- - производной 3-го порядка у~ и т. д. Вообще, интервальная производная любого п -го порядка у(п) = /~п\X) определяется выражением следующего вида

у^п) = [yXn-1)(X)]' , п = 1,2,3,...,

где обозначение у о означает исходную интервальную

функцию вида у = /(X) , определяемую по формуле (5).

Согласно определениям (12), (13) интервальной производной любого порядка, все интервальные производные, как и исходная интервальная функция (5), при

X

X

о

X

любом численном значении аргумента ~ в виде интервала возможных значений ~ = X, ^ ] также принимают численное значение в виде интервала возможных значений. Поэтому вычисление интервальной функции и интервальной производной от нее любого порядка заключается в вычислении нижних и верхних границ соответствующих интервалов. Вычисление интервальной функции / выполняется по формуле (5), задающей данную функцию в виде пары «нижняя / и верхняя / граничные функции». Вычисление интервальной производной любого п -го порядка уУ") = /(п)(~) выполняется с помощью следующей формулы, выведенной в [3-4]:

от = кп~, у2п~ 1=

2П-1

у = [х1,х2], п = 1,2,3,...

или, по-другому,

/ (п)(х) = [./1(п)(х), /2(п)(х)] =

2п-\т-

(х2 Х1) X = [х, X ], п = 1,2,3,...,

где Уха = //п)(0,У(пу = /^() - нижняя и верхняя граничные функции интервальной производной п -го порядка у~п) = / (и)0 от исходной интервальной функции

~ = [У1,У2] = 7(°) = [/1(у),/2(у)], задаваемой формулой (5).

Как видно из формул (14), (15), интервальная производная любого п -го порядка, в отличие от обычной производной, выражается в явном виде через независимую 5 = [Х1,хг] И зависимую у = [УЬУ2], У! = /¡(°), У2 = /2(°), переменные исходной функции у = /) = [/ ), /2 )]. Эта важная особенность интервальных производных функций принципиально упрощает теорию и методы решения интервально-дифференциальных уравнений.

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

~

виде У , в остальных случаях, когда нас интересуют ее значения во всех точках.

3. ОСНОВНЫЕ ПОНЯТИЯ ОБ ИНТЕРВАЛЬНО-ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЯХ

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

мых функций. Такие уравнения называются дифференциальными. Более точно дифференциальным уравнением называется соотношение, связывающее независимую переменную х , искомую функцию У = /(х) и

ее производные у',у",..., у'я) . Если искомая функция есть функция одной независимой переменной, как в нашем случае, дифференциальное уравнение называется обыкновенным. Если искомая функция является функцией двух или более независимых переменных, дифференциальное уравнение называется уравнением с частными производными. В данной работе рассматриваются только обыкновенные дифференциальные уравнения.

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

^ (X, у, у', у",...,у(п)) = 0,

где F() - некоторая функция от переменных в скобках. В частных случаях в уравнение (16) могут и не входить переменные X, у и отдельные производные от

функции у порядка ниже, чем п, но это не изменит порядка указанного уравнения, который равен п. Например, уравнения 2x + 3у - 4у' = 0 и 3у - 4у' = 0 имеют порядок, равный 1, а уравнения IX + 3у - 4 у' + у" = 0 и 3у + у' = 0 - порядок, равный 2.

Любая функция у = / (X), которая при подстановке в уравнение (16) обращает его в тождество, называется решением этого уравнения. Например, функция у = вх является решением уравнения у - 2у+у" = 0 , т. к. она при подстановке в это уравнение обращает его в тождество.

Будем теперь рассматривать задачи, решение которых в конечном итоге сводится к решению уравнений (одного или нескольких), содержащих интервальные производные искомых интервальных функций, которые были введены в § 2. Такие уравнения будем называть интервально-дифференциальными. Более точно интервально-дифференциальным уравнением будем называть соотношение, связывающее независимую интервальную переменную у = [^ X2 ], интервальную

искомую функцию у = /(у) вида (5) и ее интерваль-

о, о,, а(п) Т-П

ные производные у ,у ,...,у . Если искомая интервальная функция является функцией одной независимой интервальной переменной, как в рассматриваемом случае, интервально-дифференциальное уравнение назовем обыкновенным. Если искомая интервальная функция является функцией двух или более независимых интервальных переменных, интервально-диффе-ренциальное уравнение назовем уравнением с частными интервальными производными. В этой работе будем рассматривать только обыкновенные интервально-дифференциальные уравнения.

Как и в случае обычных (детерминированных) дифференциальных уравнений (16), порядок старшей производной, которая входит в любое интервально-

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

X (X, у, у', у'',...,у(п)) = X

где Е(•) - интервальная функция от переменных в скобках; х = [Xl, X2] - независимая интервальная переменная; у = [ У1, У2] - зависимая интервальная переменная, находящаяся в функциональной зависимости у = /(X) от независимой переменной х; у"у',...,у(п) -интервальные производные порядка 1,2,...,п от интервальной функции у = /(X); а = [а, ] - числовой интервал.

Аналогично случаю детерминированных дифференциальных уравнений (16), в частных случаях в ин-тервально-дифференциальное уравнение (17) могут и не входить интервальные переменные X, у и отдельные интервальные производные от функции у порядка ниже, чем п, но это не изменит порядка данного уравнения. Любая интервальная функция у = / (X), которая при подстановке в уравнение (17) обращает его в тождество, называется решением уравнения.

Наша задача заключается в нахождении систематического метода и алгоритма решения интервально-дифференциального уравнения (17).

4. ЭЛЕМЕНТЫ ТЕОРИИ ИНТЕРВАЛЬНО-ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

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

Е, связывающая переменные в интервально-дифференциальном уравнении (17), является алгебраической. В этих условиях оказывается справедливой следующая основная теорема.

Теорема 1. Любое интервально-дифференциальное уравнение вида (17) эквивалентно некоторой системе из двух детерминированных алгебраических уравнений вида

Е1( xl, x2, Уь У2) = а1 Е2( xl, x2, Уь У2) = а2

в которой X = [X), X2] - интервальная независимая переменная искомой интервальной функции у = /(X); у = [у, у2 ] - интервальная зависимая переменная указанной выше интервальной функции; Е1 и Е2 - некоторые детерминированные функции, представляющие собой нижнюю и верхнюю границы

интервальной связывающей функции Е уравнения (17) (т. е. Е = ЕЕ2]).

Доказательство. Согласно нашему предположению, интервальная функция Е от переменных X,у,у',у1",...у(п) в уравнении (17) является алгебраической. Поэтому она имеет вид суперпозиции элементарных операций над указанными интервальными переменными. В свою очередь, часть указанных интервальных переменных, а именно, производные у",у",...,у(п) по формуле (15) представляют собой суперпозиции элементарных алгебраических операций над интервальными переменными X = [X], X2], у = [У1, У2], где У1 = /l(x), У2 = /2^). Таким образом, интервальная

функция Е в уравнении (17) в целом может быть представлена в виде суперпозиции элементарных алгебраических операций только над интервальными переменными X = [X!, X2] и у = [У1, У2].

Учитывая, что этими элементарными операциями являются сложение, вычитание, умножение и деление, определяемые выражениями (3), и используя формулы (4) для выполнения этих операций, мы получим выражение левой части уравнения (17) в виде интервала, нижняя и верхняя границы которого представляют собой некоторые суперпозиции границ интервальных переменных х = [Xl,X2], у = [У1,У2]. Уравнение (17), таким образом, перепишется в явной интервальной форме

[F1

где Е1 и Е2 - некоторые детерминированные функции. Но два интервала равны, только если равны их одноименные границы [2]. Поэтому из (19) получаем эквивалентное ему условие (18). Но условие (19), как следует из приведенного доказательства, эквивалентно уравнению (17), следовательно, условие (18) тоже эквивалентно уравнению (17), что и требовалось доказать.

Как следует из теоремы 1, для решения произвольного интервально-дифференциального уравнения (17) следует действовать согласно алгоритму.

Шаг 1. Заменить в уравнении (17) все вхождения

интервальных производных у",у" ,...,у(п) их выражениями (15) через искомую интервальную функцию х = [У1 = /¡(х), У2 = /2(х)]. В результате интервально-дифференциальное уравнение (17) перейдет в интер-вально-алгебраическое уравнение с неизвестной переменной у .

Шаг 2. Использовать выражения (4) для выполнения элементарных алгебраических операций над интервалами, выразить левую часть получившегося уравнения (17) в виде интервала, нижняя и верхняя границы которого имеют вид суперпозиций (функций) нижних и верхних границ интервальных переменных

у = 1УиУ2] и X = ^x2] .

Шаг 3. Используя результаты шага 2, необходимо представить уравнение (17) в явной интервальной форме (19).

Шаг 4. Приравняв в (19) одноименные границы левой и правой частей, перейти к эквивалентной заданному интервально-дифференциальному уравнению (17) системе двух детерминированных алгебраических уравнений (18).

Шаг 5. Решить систему алгебраических уравнений

(18). Найденное решение этой системы в виде у = /(X) ,

где у = ^= у У2] (у = /¡(x), У2 = /2(у)), будет также искомым решением интервально-дифференциаль-ного уравнения (17).

5. ЗАКЛЮЧЕНИЕ

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

но-дифференциальное уравнение (17) легко сводится к эквивалентной системе двух алгебраических уравнений (18), решение которых дает решение исходного уравнения (17).

ЛИТЕРАТУРА

1. Фихтенгольц Г.М. Курс дифференциального и интегрального исчисления: в 3 т. М.: Наука, 2005.

2. Алефельд Г., Херцбергер Ю. Введение в интервальные вычисления. М.: Мир, 1987. 360 с.

3. Левин В.И. Интервальная производная и начала недетерминистского дифференциального исчисления // Онтология проектирования. 2013. № 4 (10). С. 72-85.

4. Левин В.И. Интервально-дифференциальное исчисление и некоторые его применения // Информационные технологии. 2014. № 7. С. 3-10.

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

Levin V.I. THE INTERVAL-DIFFERENTIAL EQUATIONS AND METHOD OF ITS SOLUTION

The generalization of ordinary differential equations on the interval case is considered. The general concept of interval-differential equation, its order, as well as the solution of interval differential equations are given. An approach to the solution of interval differential equations is developed. It is proved that any interval-differential equation can be reduced to a system of two ordinary algebraic equations. An algorithm for solving equations is given.

Key words: interval-differential equation; interval derivative; system; uncertainty.

Левин Виталий Ильич, Пензенская государственная технологическая академия, г. Пенза, Российская Федерация, доктор технических наук, профессор, советник ректора по науке, заслуженный деятель науки Российской Федерации, e-mail: levin@pgta.ru

Levin Vitaly Ilyich, Penza State Technological Academy, Penza, Russian Federation, Doctor of Technics, Professor, Science Advisor of Rector, Honored Worker of Science of Russian Federation, e-mail: levin@pgta.ru