CC BY
26
УДК 519.61
ДИФФЕРЕНЦИРОВАНИЕ ИНТЕРВАЛЬНЫХ ЭЛЕМЕНТАРНЫХ ФУНКЦИЙ И ПРИНЯТИЕ РЕШЕНИЙ
В.И. Левин
Пензенский государственный технологический университет, Пенза, Россия vilevin@mail.ru
Аннотация
В статье рассмотрены задачи, связанные с вычислением производных от интервально-определён-ных функций. Эти задачи актуальны при изучении систем с той или иной степенью неопределённости (недетерминированные системы). Рассматриваются простейшие системы, описываемые элементарными интервально-определёнными функциями. Решаются задачи вычисления производных от элементарных интервально-определённых функций. Приведены основные определения, связанные с производными от интервально-определённых функций, и формулы двух типов, которые позволяют вычислять указанные интервальные производные. Формулы первого типа выражают производные в закрытой интервальной форме, которая для вычисления требует использования аппарата интервальной математики. Формулы второго типа выражают производные в открытой интервальной форме, т.е. в форме двух формул, первая из которых выражает нижнюю границу интервала, представляющего производную, а вторая - верхнюю границу интервала. В такой форме вычисление производной от интервально-определённой функции сводится к вычислению двух обычных функций. Получены производные для ряда элементарных интервальных функций: интервальной константы, степенной функции, показательной функции, экспоненциальной функции, логарифмической функции, функции натуральный логарифм, тригонометрических функций, обратных тригонометрических функций. Формулы всех производных представлены в открытой интервальной форме. Показано отличие производных интервальных элементарных функций от производных соответствующих неинтервальных элементарных функций.
Ключевые слова: интервал, интервальная функция, интервальные вычисления, интервальная производная, интервально-дифференциальное исчисление, интервальные элементарные функции.
Цитирование: Левин, В.И. Дифференциальное исчисление интервальных элементарных функций и принятие решений / В.И. Левин // Онтология проектирования. - 2016. - Т.6, №3(21). - С. 340-354. - DOI: 10.18287/2223-9537-2016-6-3-340-354.
Введение
Проектирование и анализ свойств разнообразных систем требует соответствующего, адекватного рассматриваемой задаче, математического аппарата. Встречающиеся на практике системы часто характеризуются той или иной степенью неопределённости (недетерминирова-ны). Для исследования и построения таких систем обычно применяют специализированный математический аппарат - теорию вероятностей, нечёткие множества, интервальную математику [1-3].
В работах [4, 5] автором был предложен новый математический аппарат для проектирования и исследования недетерминированных систем - недетерминистское (интервальное) дифференциальное исчисление. Этот аппарат является аналогом классического дифференциального исчисления Ньютона-Лейбница [6]. Он позволяет переносить основные идеи классического дифференциального исчисления на неполностью определённые функции, задаваемые с точностью до интервалов возможных значений переменных. Однако, несмотря на сходство основных исходных идей двух исчислений, предложенное интервальное дифференциальное исчисление по форме не похоже на классическое дифференциальное исчисление Ньютона-
Лейбница, что является следствием неопределённости интервальных функций, обычно фигурирующих в интервальной математике. Кроме того, предложенное интервальное дифференциальное исчисление, по нашему мнению, более адекватно природным явлениям и процессам, чем классическое дифференциальное исчисление [7].
1 Постановка задачи
Согласно классическому дифференциальному исчислению [6], нахождение производной от любой функции базируется на: 1) представлении этой функции в виде соответствующей суперпозиции так называемых элементарных функций; 2) переходе в полученном представлении от функций к их производным, с использованием для этого основных теорем дифференциального исчисления (производная суммы функций, производная произведения функций и т.д.); 3) подстановке выражений производных элементарных функций.
Процедура, сходная с описанной, может быть полезна и при нахождении производных от интервальных функций, рассматриваемых в интервальном дифференциальном исчислении. При этом необходимо учитывать большое отличие свойств и форм представления обычных и интервальных функций и производных от них. Задачами статьи являются: 1) составление полного набора производных от всех интервальных элементарных функций; 2) выявление различий производных от интервальных элементарных функций и производных от детерминированных элементарных функций.
Интервальное дифференциальное исчисление может быть использовано для исследования поведения неполностью определённых функций, служащих характеристиками недетерминированных систем, при принятии решений на этапе проектирования таких систем. Эта ситуация аналогична использованию классического дифференциального исчисления при исследовании поведения функций, являющихся характеристиками детерминированных систем.
При решении задачи будем использовать в качестве вспомогательных сведений основные математические сведения из алгебры интервальных чисел [3, 4]. В этой алгебре оперируют интервальными числами (операндами), операциями над интервальными числами и интервальными функциями [4, 5, 7].
2 Математический аппарат
Базовое для данной статьи понятие интервальной функции [4, 5, 7] вводится как однозначное отображение множества {Х} замкнутых вещественных интервалов ~ = [х1, х2] на множество {у} замкнутых вещественных интервалов у — [у1,у2] того же вида. Символически интервальная функция записывается в виде
(1) у-/(у) или у-Цух,У2] -[/К^Лта
где у-[х1, Х2] - интервальная независимая переменная, у—У у] - интервальная зависимая переменная, /-/,/2] - интервальная функция, с её нижней /1 и верхней /2 граничными функциями. Второе базовое понятие, используемое далее, - понятие предела интервальной функции.
Независимая интервальная переменная ХХ - [х1, х2] интервальной функции (1) по определению неограниченно приближается к некоторому интервалу (пределу) ЗГ-Х", Х2 ], если в процессе указанного изменения х1 неограниченно приближается к х[ , а х2 к х2 . Символически это записывается в виде
(2) 33 ^33° = (х1 ^ х[ , х2 ^ х2 ).
Аналогично определяется неограниченное приближение зависимой интервальной переменной y =[y1,y2] интервальной функции (1) к интервалу (пределу) y° = [yl,y\]:
(3) y = (yj ^yi, y2 ^у2 ).
Если независимая интервальная переменная у своим неограниченным приближением к интервалу-пределу у° вызывает неограниченное приближение зависимой интервальной переменной у к интервалу-пределу y°, то говорим, что предел интервальной функции (1) при у, стремящемся к у °, равен y°:
(4) lim у = ~ ° или lim f(x) = у °.
Если интервальная функция f (1) непрерывна, т.е. её нижняя f1 и верхняя f2 граничные функции есть непрерывные функции нижней x1 и верхней x2 границ независимой переменной у = [x1, x2], то предел функции f равен её значению в предельной точке у° аргумента у, или символически
(5) lim f(~)=f(~ °).
Основное понятие интервальной производной вводится аналогично понятию обычной производной функции [6]. Рассмотрим произвольную интервальную функцию f вида (1). Будем считать её непрерывной. Зафиксируем значение независимой переменной у = у = [x[, x22 ]. Этому значению, в силу непрерывности функции, соответствует фиксированное значение самой функции У = f(Т). Определим приращение независимой и зависимой переменных функции относительно этих фиксированных значений
(6) Ax = у-у°, Ду = у -у° = f(x) -f(x°) и составим отношение второго приращения к первому
(7) Ду / Дс = (у - y ° )/(У - Г) = (ДУ) - f(F ))/(У - у°).
Предел отношения (7) при неограниченном приближении независимой переменной у к её фиксированному предельному значению у°, если он существует, называется интервальной производной функцией от исходной интервальной функции f (у) (1) в точке у ° и обозначается yL или flc (x). Таким образом,
(8) уу = JL (х) = lim ДУ / Дс, где Д, Ду из (6).
Доказано [4, 5, 7], что для существования у непрерывной интервальной функции у = f(у) в
точке у° интервальной производной необходимо и достаточно, чтобы в этой точке и некоторой её окрестности как независимая ~x , так и зависимая y переменные были существенно интервальными (т.е. не вырождались в точку).
Как и в случае обычной производной, понятие интервальной производной можно обобщить путём повторного выполнения операции взятия производной. При этом из интервальной производной 1-го порядка y~ получается интервальная производная 2-го порядка y~,
а из последней - интервальная производная 3-го порядка уУ и т.д. Согласно введённым определениям интервальной производной любого порядка, все интервальные производные, как и исходная интервальная функция, при любом численном значении аргумента ~x в виде интерва-
ла возможных значений x = [x1, x2 ] также принимают численные значения в виде интервала возможных значений. Поэтому вычисление интервальной функции и интервальной производной от неё любого порядка состоит в вычислении нижних и верхних граничных функций соответствующих интервальных функций. Вычисление интервальной функции f выполняется согласно формуле (1), задающей эту функцию в виде пары «нижняя f и верхняя f2 граничные функции». Вычисление интервальной производной любого n -го порядка yin) = f~n\x) от интервальной функции у = f (f выполняется с помощью формулы [7]
(9) yxn) = fn)x=f^ff^c)],
где fi(~n) и f2( f - соответственно нижняя и верхняя граничные функции интервальной производной f ~n) n -го порядка от исходной интервальной функции y = f (x). Эти граничные функции в формуле (9) выражаются в виде
(10) f!§(x) = -2nAy -yi)/(x2 -xi)n, f^(x) = 2nA(y2 -л)/^ -xi)n,
при этом x1, x2 - нижняя и верхняя границы интервальной независимой переменной х =[x1, x2] в точке взятия производной от функции у = fx, У1, У2 - нижняя и верхняя границы интервальной зависимой переменной у=[y1, y2] этой функции в той же точке. Как следует из выражений (9), (10), интервальная производная любого порядка имеет вид интервала, симметричного относительно нуля. Это позволяет записать выражение интервальной производной любого n -го порядка (9), (10) в более простой форме
(11) yfxn) (У)=fn) (f=f (X), ffn) (у)],
где
(12) ff(X)=2n-1y - У1)/(Я2 - x1)n,
а значения x1, x2, y1, y2 раскрыты в пояснениях к формуле (10).
Формулы (11), (12) позволяют находить интервальные производные любого порядка от любых интервальных функций, в том числе от элементарных интервальных функций.
3 Производные от элементарных интервальных функций
По определению любая элементарная интервальная функция y = f(~x) вида (1) получается из соответствующей элементарной вещественной функции y=f(x) путём преобразования в соответствующие интервальные формы независимой переменной у = [x1, x2], зависимой переменной
У = [У1, y2] и интервальной функции f = [f1, f2]. Определение производных от элементарн^1х интервальн^1х функций начнём с простейшей интервальной функции - интервальной константы. Функция интервальная константа выражается в виде
(13) у = c =[c1,c2], c1 = const, c2 = const, c1 <c2.
Из сравнения выражений (13) и (1) видно, что интервальная константа - это интервальная функция, нижняя и верхняя граничные функции которой имеют следующий вид
(14) У1 = f1(f) = f1(x1, x2) = c1, У2 = f2(X) = f2(x1, x2) = c2.
Подставив (14) в общие формулы интервальных производных (11), (12), получим выражение производной любого n -го порядка от интервальной константы y =c в виде
(15) y(n) -X(n) = [Cl,c2] = [-cXn)(0,c^n)(X)],
где
4И)(~) = 2п-1(с2 -с1)/(х2 -х,)п. Как видно из (15), интервальная производная любого п -го порядка от интервальной константы с (13) не равна нулю (точнее, нулевому интервалу 0=[0,0]), в отличие от классической производной от вещественной константы, равной нулю. Более того, эта производная не является постоянной величиной, а существенно зависит от аргумента ~ = [х,, х2]. Согласно формуле (15), она существует во всех точках х, в которых х1 Фх2, и монотонно убывает при увеличении разности х2 - х1.
Интервальная степенная функция выражается в виде
(16) у = [у,, у2 ] = С" - [х,, Х2Г ,
где х = [х,,х2] - интервальная независимая переменная, у = [у,,у2] - интервальная зависимая переменная. Будем считать, исходя из физических соображений, что в пределах одной (любой!) решаемой задачи переменная величина х = [х,, х2] может быть только неотрицательной или только неположительной, т.е. выполняется условие
(17) х,,х2 > 0 или х,,х2 < 0.
Тогда интервальную степенную функцию (16) можно записать в явном интервальном виде посредством следующей формулы:
х ([х1,х2т], при х,,х2 >0 или х,,х2 <0, т нечётно;
(18) у =[у1,у2] Нг т тп
([х2,х ], при х,,х2 <0, т чётно.
Из сравнения выражения (18) с (1) видно, что интервальная степенная функция есть интервальная функция, нижняя и верхняя граничные функции которой имеют такой вид:
ч ¡хт, при х,,х2 >0 илих,,х2 <0, т нечётно;
у = /1(Х) = х2) Ч т ' -1 2
1х2, при х,,х2 <0, т чётно;
(19) у 2
„ ч ¡хт при х,,х2 >0 илих,,х2 <0, т нечётно;
у2 = /2(Х) = х2) Ч т ' „ 2
(х , при х,,х2 <0, т чётно Подстановка (19) в (11) и (12) даёт выражение производной любого п -го порядка от интервальной степенной функции у = хт (16) в виде
(20) у(п) = (хт )(п) = ([х,, х2]т )(п) = [-/ (п)(х),/(п)(х)],
где
(2п 1(х2Т - хт )/(х2 - х|)п, при х1, х2 > 0 или х1, х2 < 0, т нечётно; (2п-1(хт - хт )/(х2 - х|)п, при х,, х2 < 0, т чётно.
/ (п)(~) =
Как видно из формулы (20), в интервальной производной любого п -го порядка (хт )(п) от
интервальной степенной функции хт (16) при увеличении п показатель степени т не уменьшается, приближаясь к нулю, в отличие от классической производной от вещественной степенной функции. Интервальная производная (хт )(п), согласно формуле (20), существует во всех точках х =[х,,х2], в которых х, Фх2; она монотонно убывает с увеличением разности х2 х,.
Интервальная показательная функция выражается в виде (21) у = [у,,у2] = ах -а[**], а>0,
где 33 - [х1,х2] - интервальная независимая переменная, у -[у1,у2] - интервальная зависимая переменная. Согласно [4, 5, 7] интервальная функция аХ определяется в виде
(22) 0х -{аХ | хех},
здесь аХ - исходная обычная показательная функция, которая монотонно возрастает. Это позволяет записать интервальную показательную функцию (21) в явном интервальном виде с помощью формулы
(23) у = [у!,у2] -[а\аХ2], а>0.
Сравнивая (23) с (1), устанавливаем, что интервальная показательная функция есть интервальная функция, нижняя и верхняя граничные функции которой принимают вид:
(24) у - /1С33) - /1(х1, Х2) - аХ1, у 2 - /2 (Х) - /2(х1, Х2) - аХ2, а > 0. Подставив (24) в (11) и (12), получим выражение производной любого п -го порядка от
интервальной показательной функции у - аХ (21)
(25) у(п) -(аХ)(п) -(а^)^ -[-/(п)(Х),/(п)(Х)],
где
/(п)(33) - 2п-1 (аХ -аХ')/(х2 -х1)п . Из (25) видно, что в интервальной производной любого п -го порядка (аХ )(п) от интервальной
показательной функции аХ (21) при увеличении п перед данной функцией не появляются дополнительные множители 1п а, в отличие от классической производной от вещественной показательной функции, у которой указанные множители появляются. Интервальная производная (а )(п), согласно (25), существует во всех точках х -[х1,х2], в которых х1 Фх2; она монотонно убывает с увеличением разности х2 - х1 .
Интервальная экспоненциальная функция выражается в виде
(26) у -[у1,у2] - еХ = е[Х1,Ч
где 33-[х1, х2] - интервальная независимая переменная, у-[у у2] - интервальная зависимая переменная. Из сравнения (21) и (26) видно, что интервальная экспоненциальная функция (26) есть частный случай интервальной показательной функции (21) при а - е. Таким образом, из выражения (25) производной любого п -го порядка от интервальной показательной функции, положив в нем а - е, получим следующее выражение производной любого п -го порядка от интервальной экспоненциальной функции
(27) у(п) -(еХ)(п) -(е^)^ -[-/(n)(Х),/(n)(Х)], при этом
/(п>(3с) - 2п-1(еХ -еХ1)/(х2 -х1)п .
Из (27) видно, что в интервальной производной (еХ )(п) любого п -го порядка от интервальной экспоненты еХ (26) при увеличении п изменяется множитель перед экспонентой, равный
(28) М- 2п-1 /(х2 -Х])п ,
в отличие от классической производной от вещественной экспоненты, у которой такого множителя нет. Интервальная производная (еХ )(п), согласно формуле (27), существует во всех
точках 33 - [х1, х2], в которых х1 Ф х2. Она монотонно убывает с увеличением разности Интервальная логарифмическая функция выражается в виде (29) У -[у1,у2] - 1ogaХ - 10а,[х1,х2] а > 0,
Х2 Х1
где ~ = [xl,x2] - интервальная независимая переменная, ~ = [y1,y2] - интервальная зависимая переменная. Согласно [4, 5, 7] интервальная функция loga~ определяется в виде
(30) logax -{logax|xe~), a>0,a*1,
при этом log ax - исходная обычная логарифмическая функция, которая монотонно возрастает. Это позволяет записать интервальную логарифмическую функцию (29) в явном интервальном виде посредством формулы
(31) ~ = [У1,У2] - [logaXl,logaX2], a > ° a * 1
Сравнивая (31) с (1), приходим к выводу, что интервальная логарифмическая функция - это интервальная функция, нижняя и верхняя граничные функции которой выглядят как
(32) У1 - /1(~) - f1(x1,x2) - logaX1, У2 - f2(X) - f2(x1,x2) - logaX2, a > 0, a * 1. Подставляя выражения (32) в (11) и (12), получим выражение производной любого n -го порядка от интервальной логарифмической функции ~ - log ax (29) в виде
(33) ~(П) - (l0ga~)(n) - (l0ga[x1,x2])(n) - [-/(n)(~),/(n)(~)],
где
/(n)(~) - 2n-1(l°gax2 -logax1)/(x2 -x1)n, a > 0, a * 1,
или, после потенцирования,
/(n)(~) - 2n-1log a(x2/x1)/(x2 - x1)n, a > 0, a * 1. Как видно из (33), в интервальной производной (log a~)(n) любого n -го порядка от интервальной функции логарифм log a~ исходная логарифмическая функция остаётся логарифмической, в отличие от классической производной от вещественной логарифмической функции log ax, которая равна 1/x ln a, т.е. является рациональной функцией. Интервальная производная (log~)(nl существует во всех точках ~ - [x1, x2], в которых x1 * x2, и монотонно убывает с увеличением разности x2 -x1.
Интервальная натуральная логарифмическая функция выражается в виде
(34) ~ - [y1, У2] - ln~ - ln[x1, x2],
где х - [x1, x2] - интервальная независимая переменная, ~ - [у1, у2] - интервальная зависимая переменная. Из сравнения (29) и (34) видно, что интервальная натуральная логарифмическая функция (34) есть частный случай интервальной логарифмической функции (29) при a - e. Производная любого n -го порядка от интервальной натурально-логарифмической функции представляется в виде
(35) уn -(ln~)(n) -(ln[c1,x2])(n) -[-/(n)(~),/(n)(~)],
где
/(n)(~) - 2n-1 ln(x2 /x,)/(x2 -x1)n . Как видно из (35), интервальная производная любого n -го порядка (lnx)(n) от интервальной натурально-логарифмической функции ln ~ является натурально-логарифмической, в отличие от производной от обычной натурально-логарифмической функции ln x, которая равна 1/ x и
является рациональной функцией. Найденная интервальная производная (lrnf существует во
всех точках ~-[x1, xj], в которых x1 * x2. Она монотонно убывает с увеличением разности
x2 - x1 .
Интервальная тригонометрическая функция синус выражается в следующем виде
(36) ~ = [y1, y2] = sin ~ = sin^Cj^, x2],
где ~=[xj,x2] - интервальная независимая переменная, y = [y1,y2] - интервальная зависимая переменная. Согласно [4, 5, 7] интервальная функция sin ~ определяется в виде
(37) sin x = {sin x| xex},
здесь sin x - исходная обычная тригонометрическая функция синус, которая монотонно возрастает на интервале —ж/2<x<ж/2 и принимает все возможные значения от -1 до 1. Последнее позволяет записать интервальную тригонометрическую функцию синус (37) в явном интервальном виде
(38) ~ =[y1, У2] = sin ~=sin[x1, x2] =[sin x1,sin x2].
Сравнивая (38) с (1), заключаем, что интервальная тригонометрическая функция синус есть интервальная функция, нижняя и верхняя граничные функции которой имеют вид
(39) У1 = /1(~) = f^x2) = sinx1, У2 = /2(~) = /2(x1,x2) = sinx2 . Подставляя (39) в (11) и (12), получаем выражение производной любого n -го порядка от
интервальной тригонометрической функции y = sin x в виде
(40) y(n) = (sin ~)(n) = (sinft, x2])(n) = [—/(n)(X),/(n)(~)],
где
/ (n)(x) = 2n—!(sin x2 — sin xj)/(x2 — xj)n. Из формулы (40) видим, что в интервальной производной любого n -го порядка (sin ~)(n) от интервальной тригонометрической функции sin x исходная тригонометрическая функция синус остаётся синусом, в отличие от производной от обычной тригонометрической функции sin x, которая равна cos x. Интервальная производная (sin x)(n) существует во всех точках ~=[x1, x^], в которых x^x;,. Она монотонно убывает с увеличением разности x2 —x1. Интервальная тригонометрическая функция косинус выражается в виде
(41) y =[y1,y2] = cos ~ = cos [x1,x2],
где x =[x1, x2] - интервальная независимая переменная, y=[>i,y2] - интервальная зависимая переменная. Согласно [4, 5, 7] интервальная функция cos ~ определяется в виде
(42) cos ~ = {cos x|xe~},
здесь cos x - исходная обычная тригонометрическая функция косинус, которая монотонно убывает на интервале 0<x <ж, принимая последовательно все свои возможные значения от 1 до -1. Это позволяет записать интервальную тригонометрическую функцию косинус в виде
(43) ~ =[y1,y2]=cosх = cos [x1,x2]=[cos x2,cosx1].
Сравнивая (43) с (1), получаем, что интервальная тригонометрическая функция косинус является интервальной функцией, нижняя и верхняя граничные функции которой имеют вид
(44) y1 = f1(x) = /1(x1,x2) = cos x2, y2 = f2(x) = /2(x1,x2) = cos x1. Подставляя (44) в (11) и (12), получим выражение производной любого n -го порядка от
интервальной тригонометрической функции y = cos x в виде
(45) y(n) = (cos x)(n) = (cos [x!,x2])(n) = [—/(n)(X),/(n)(x)],
где
/ (n)(3c) = 2n—2(cos xj — cos x2)/(x2 — xj)n. Согласно (45), в интервальной производной любого n -го порядка (cos ~)(n^ исходная тригонометрическая функция косинус остается косинусом, в отличие от производной от вещественной тригонометрической функции cos x, равной — sin x. Интервальная производная
(cos у)(n) существует во всех точках х—xx2], где x1 Фx2. Она монотонно убывает с увеличением разности x2 -x1.
Интервальная тригонометрическая функция тангенс выражается в виде
(46) ~ - [л^] - tg ~ - tg[x1,x2],
где х-[x1,x2] - интервальная независимая переменная, y —[>1,y2] - интервальная зависимая переменная. Согласно [4, 5, 7], интервальная функция tg х определяется в виде
(47) tg x — {tg x | x e ~},
причём tg x - исходная обычная тригонометрическая функция тангенс, которая монотонно возрастает на интервале -л/2<x<л/2, где она последовательно принимает все свои возможные значения от — да до + да. Это позволяет записать интервальную тригонометрическую функцию тангенс (47) в явном интервальном виде
(48) ~ = у y2] — tg ~ — tg [xl, x2] — [tg xi,tg x2].
Сравнив выражение (48) с (1), устанавливаем, что интервальная тригонометрическая функция тангенс есть интервальная функция, нижняя и верхняя граничные функции которой имеют вид
(49) У1 — /1(у) — fl^x2) — tg x1, У2 — /2(у) — /2(x1,x2) — tg x2.
Подставляя (49) в (11) и (12), получим выражение производной любого n -го порядка от интервальной тригонометрической функции у — tg х в виде
(50) у(n) — (tg ~)(n) — (tg[x1, x2 ])(n) — [-/(n)i~),/(n)m,
где
/(n)(~) — 2n—1(tg X2 — tg x1)/(x2 — x1)n . Формула (50) показывает, что в интервальной производной любого n -го порядка (tg x)(n^ от интервальной тригонометрической функции tg ~x исходная тригонометрическая функция тангенс остаётся тангенсом, в отличие от классической производной от вещественной тригонометрической функции tg x, равной 1/cos2 x. Интервальная производная (tg х)(n^ существует во всех точках х—[x1, x2], где x1 Ф x2. Она монотонно убывает с увеличением разности x2 — x1. Интервальная тригонометрическая функция котангенс выражается в виде
(51) ~ — [У1,У2] — ctg~ — ctg[x1,x2],
где х — [ x1, x2] - интервальная независимая переменная, у — [ у1, у2] - интервальная зависимая переменная. Согласно [4, 5, 7], интервальная функция ctg х определяется в виде
(52) ctg х — {ctg x | x e~},
здесь ctg x - исходная обычная тригонометрическая функция котангенс, которая монотонно убывает на интервале 0 < x <л, где она последовательно принимает все возможные значения от + да до — да. Это позволяет записать интервальную тригонометрическую функцию (52) в явном интервальном виде
(53) у = [ л, у2] — ctg ~ — ctg [x1, x2] — [ctg x2, ctg x1 ].
Сравнивая (53) с (1), видим, что интервальная тригонометрическая функция котангенс есть интервальная функция, нижняя и верхняя граничные функции которой имеют вид
(54) у1 — /1(у) — /1(x1, x2) — ctg ^ у2 — /2(У) — А^Ъ x2) — ctg x1.
Подставляя (54) в (11) и (12), получаем выражение производной любого n -го порядка от интервальной тригонометрической функции у — ctg х в виде
(55) у(п) — (ctg ~)(n) — (ctg[x1, x2])(n) — [—/(n)(X),/(n)(X)],
где
f (n>(x) = 2n—1(ctg X! - ctg x2)/(x2 -x). Выражение (55) показывает, что в интервальной производной любого n -го порядка (ctg x)(n) от интервальной тригонометрической функции ctg у исходная тригонометрическая функция котангенс остаётся котангенсом, в отличие от классической производной от обычной тригонометрической функции ctg x, которая равна -1/sin2 x. Интервальная производная (ctg ~)(n) существует во всех точках у=[х, x2], в которых x1 Ф x2. Она монотонно убывает с увеличени-Xo xi
ем разности 2 1.
Интервальная обратная тригонометрическая функция арксинус выражается в виде
(56) у = [y1, y2] = arcsin у=arcsin [x1, x2 ],
где x^[xi,x2] - интервальная независимая переменная, ~ =[yi, y2] - интервальная зависимая переменная. Согласно [4, 5, 7], интервальная функция arcsin x определяется в виде
(57) arcsin у = {arcsin x | xe x},
в котором arcsin x есть исходная обычная обратная тригонометрическая функция арксинус. Последняя монотонно возрастает на интервале —1< x <1, проходя последовательно все возможные значения от —ж/2 до ж 12. Это позволяет записать интервальную обратную тригонометрическую функцию (57) в явном интервальном виде следующим образом
(58) у = [y1, y2] = arcsin у = arcsin [x1, x2]=[arcsin x1, arcsin x2 ].
Сравнивая (58) с (1), устанавливаем, что интервальная обратная тригонометрическая функция арксинус является интервальной функцией, нижняя и верхняя граничные функции которой имеют вид
(59) y1 = f1(x) = f1(x1,x2) = arcsinx1, y2 = f2(x) = f2(x1,x2) = arcsinx2. Подставив (59) в (11) и (12), получим выражение производной любого n -го порядка от
интервальной обратной тригонометрической функции y = arcsin x в виде
(60) y(n) = (arcsin ~)(n) = (arcsin [x1, x2])(n) =[—f (n)(3X),f (n)(3X)],
где
f (n)(x) = 2n—1'(arcsin x2 — arcsin x1)/(x2 — x1)n. Как видно из формулы (60), в интервальной производной любого n -го порядка (arcsin у)(n) от интервальной обратной тригонометрической функции arcsin у исходная обычная обратная тригонометрическая функция арксинус остаётся арксинусом, в отличие от производной от обычной обратной тригонометрической функции arcsin x, равной 1/V1 —x . Производная (arcsinу)(n) существует во всех точках х = [xx,x2], где xx Фx2 и монотонно убывает с увеличением разности x2 — x.
Интервальная обратная тригонометрическая функция арккосинус выражается в виде
(61) у =[y1,y2] = arccos у=arccos[xx,x2],
где у=[х, x2] - интервальная независимая переменная, у =[у, y2] - интервальная зависимая переменная. Согласно [4, 5, 7], интервальная функция arccos у определяется в виде
(62) arccos у = {arccos x|xex},
при этом arccos x - исходная обычная обратная тригонометрическая функция арккосинус, которая монотонно убывает на интервале —1 < x <1, проходя последовательно все возможные значения от ж до 0 . Это позволяет записать функцию арккосинус (62) в интервальном виде
(63) у =\У\,У2]-arccos ~=arccos \х1,х2]=[arccos х2,arccos х1 ].
Сравнив выражение (63) с (1), видим, что интервальная обратная тригонометрическая функция арккосинус является интервальной функцией, нижняя и верхняя граничные функции которой имеют вид
(64) у1 -/1(х) -/1(х1,х2) -arccos х2, у2 -^(~) -^(х1,х2) -arccos х1. Подставив выражения (64) в (11) и (12), получим выражение производной любого п -го
порядка от интервальной обратной тригонометрической функции у - arccos х в виде
(65) у(п) -(агссоБ ~)(п) -(arccos \х1,х2])(п) -\-f (n)(X),f (п)(~)],
где
Уп)(~) - 2n-1(arccos х1 - arccos х2)/(х2 -х1)п. Формула (65) показывает, что в интервальной производной любого п -го порядка (arccos х)(п) от интервальной обратной тригонометрической функции arccos х исходная обычная обратная тригонометрическая функция арккосинус остаётся арккосинусом, в отличие от производной от обычной обратной тригонометрической функции arccos х, которая равна
—х . Интервальная производная (arccos х) существует во всех точках х -\х1,х2], в которых х1 Ф х2. Она монотонно убывает с увеличением разности х2 -х1.
Интервальная обратная тригонометрическая функция арктангенс выражается в виде
(66) у-[у1,у2]-arctg х -arctg \х1,х2],
где х-\х1, х2] - интервальная независимая переменная, у-\>1,у2] - интервальная зависимая переменная. Согласно \4, 5, 7] интервальная функция arctg х определяется в виде
(67) arctg х - {arctg х | х е х},
в котором arctg х - исходная обычная обратная тригонометрическая функция арктангенс. Последняя монотонно возрастает на интервале -да< х <«, проходя последовательно все возможные значения от -я / 2 до я/2, что даёт возможность записать интервальную обратную тригонометрическую функцию арктангенс (67) в виде
(68) ~ =\у1,у2]-arctg~-arctg \х1,х2]-\arctgх^гС^х2].
Сравнив выражение (68) с (1), увидим, что интервальная обратная тригонометрическая функция арктангенс является интервальной функцией, нижняя и верхняя граничные функции которой имеют вид
(69) у -У1(х) -^^х2) - ЭГС^ xl, У2 -У2(х) - f2(x1,х2) - arctg х2 . Подставив (69) в (11) и (12), найдём выражение производной любого п -го порядка от интервальной обратной тригонометрической функции ~ - arctg х в виде
(70) у(п) -(вл^х)(п) -(aгctg\xl,х2])(п) -\-/{п)(~),/{п\х)1
где
f (п)(х) - 2п-11(вг^ х2 -aгctg х1)/(х2 -х1)п. Формула (70) показывает, что в интервальной производной любого п -го порядка ^гС^ х)(п) от интервальной обратной тригонометрической функции arctg х исходная обратная тригонометрическая функция арктангенс остаётся арктангенсом, в отличие от производной от обычной обратной тригонометрической функции arctg х, которая равна 1/(1+х2). Интервальная производная ^г^ х)(п) существует во всех точках х -\х1,х2], в которых х1 Фх2. Она монотонно убывает с увеличением разности х2 -х1.
Интервальная обратная тригонометрическая функция арккотангенс выражается в следующем виде
(71) 3-[у1, у2] - arcctg33 - агсС^^, х2],
где 33 - [х1, х2] - интервальная независимая переменная, 33 - [у1, у2] - интервальная зависимая переменная. Согласно [4, 5, 7] интервальная функция агсС^ х определяется в виде
(72) агс^ х - {агсС^ х | х е х},
при этом агс^ х есть исходная обычная обратная тригонометрическая функция арккотангенс. Последняя монотонно убывает на интервале -да< х и проходит последовательно все возможные значения от л до 0. Это позволяет записать интервальную обратную тригонометрическую функцию арккотангенс (72) в виде
(73) 3=[у1, у2] - агс^33 - агс^ [х1 , х2] - [агс^ х2 , arcctg х1 ].
Сравнив (73) с (1), видим, что интервальная обратная тригонометрическая функция арккотангенс является интервальной функцией, нижняя и верхняя граничные функции которой имеют вид
(74) у - Ж3) - /1 (х1 , х2) - arcctg x2, у2 - /2 (х) - /2 (х1, х2) - агсс^ х1.
Подставив (74) в (11) и (12), получим выражение производной любого п -го порядка от интервальной обратной тригонометрической функции у - агс^ х в виде
(75) у(п) -(агс^х)(п) -Сarcctg [х1,Х2])(п) -[-/{п)(Х),/(п)(Х)],
где
/(п)(х) - 2п-1(агс^ х1 - arcctg х2)/(х2 -х1)п. Формула (75) показывает, что в интервальной производной любого п -го порядка (агс^ 33)(п) от интервальной обратной тригонометрической функции агс^ х исходная обратная тригонометрическая функция арккотангенс остаётся арккотангенсом, в отличие от производной от обычной обратной тригонометрической функции агс^ х, которая равна -1/(1+х2). Интервальная производная (агс^ х)(п) существует во всех точках 33-[х1, х2], в которых х1 Ф х2. Она монотонно убывает с увеличением разности х2 -х1 .
4 Обсуждение
Проблема вычисления производных от интервально-определённых функций существенно отличается от аналогичной проблемы для полностью определённых функций. Это отличие связано с известным отличием процесса вычисления полностью определённой функции от процесса вычисления интервальной функции - во втором случае приходится вычислять параллельно две независимые точно заданные функции - нижнюю и верхнюю граничные функции интервальной функции. В связи с этим отличием вычисление неполностью определённой - интервальной функции оказывается более сложным [8-13]. Эта сложность вычисления распространяется и на производные от интервальных функций, поскольку они также являются интервальными функциями. Поэтому предпринятое в статье детальное изучение производных элементарных функций, сопровождаемое получением простых формул для вычисления этих производных, является важным шагом на пути уменьшения сложности вычисления интервальных функций.
В процессе изучения выявлены различия между интервальными производными от интервальных элементарных функций и классическими производными от детерминированных элементарных функций. Главное различие состоит в том, что производная любого п -го
(п -1,2,...) порядка от любой интервальной функции Р(х) , в частности, элементарной функции, всегда представляет собой функцию класса Р(х)/х" . Поэтому интервальную производную любого п -го порядка с полным правом можно интерпретировать как скорость п -го порядка изменения функции Р(х) относительно её аргумента х . Известно, что классическая производная любого порядка п от любой детерминированной функции Р(х) не обязательно представляет собой функцию класса Р(х)/хп . Поэтому такую производную в общем случае нельзя интерпретировать как скорость п -го порядка изменения функции Р(х) относительно её аргумента. Это означает, в свою очередь, что классическая производная Ньютона-Лейбница, в отличие от интервальной производной, вообще говоря, не может считаться адекватной моделью динамики большинства природных объектов и процессов.
Заключение
В работе получен полный набор производных любых порядков от всех элементарных интервальных функций. Выявлены различия между интервальными производными от интервальных элементарных функций и классическими производными от обычных детерминированных элементарных функций. Показано, что учёт неопределённости исходных функций в форме интервальности их возможных значений делает производные от них адекватными моделями динамики природных процессов.
Список источников
[1] Гнеденко, Б.В. Курс теории вероятностей / Б.В. Гнеденко // - М.: Наука, 2004. - 350 с.
[2] Заде, Л.А. Понятие лингвистической переменной и его применение к принятию приближенных решений / Л.А. Заде // - М.: Мир, 1976. - 160 с.
[3] Алефельд, Г. Введение в интервальные вычисления / Г. Алефельд, Ю. Херцбергер // - М.: Мир, 1987. -356 с.
[4] Левин, В.И. Интервальная производная и начала недетерминистского дифференциального исчисления / В.И. Левин // Онтология проектирования. - 2013. - № 4. - С. 72-84.
[5] Левин, В.И. Интервально-дифференциальное исчисление и некоторые его применения / В.И. Левин // Информационные технологии. - 2014. - № 7. - С. 3-10.
[6] Фихтенгольц, Г.М. Курс дифференциального и интегрального исчисления. Т. 1. / Г.М. Фихтенгольц // - М.: Физматлит, 2001. - 616 с.
[7] Левин, В.И. Дифференциальное исчисление для интервально-определенных функций / В.И. Левин // Эвристические алгоритмы и распределенные вычисления. - 2015. Т. 2. - № 2. - С. 8-25.
[8] Вощинин, А.П. Оптимизация в условиях неопределённости / А.П. Вощинин, Г.Р. Сотиров //. - М: Изд-во МЭИ, 1989. - 224 с.
[9] Ащепков, Л.Т. Универсальные решения интервальных задач оптимизации и управления / Л.Т. Ащепков, Д.В. Давыдов // - М.: Наука, 2006. - 285 с.
[10] Moore R.E. Interval Analysis. - N.Y.: Prentice-Hall, 1966. - 230 p.
[11] Libura M. Integer Programming Problems with Inexact Objective Function // Control and Cybernetics. - 1980. Vol. 9. - № 4. - P. 189-202.
[12] Куржанский, А.Б. Задача идентификации - теория гарантированных оценок / А.Б. Куржанский // Автоматика и телемеханика. - 1991. - № 4. - С. 447-465.
[13] Hyvonen E. Constraint Reasoning Based on Interval Arithmetic: the Tolerance Propagation Approach // Artificial Intelligence. - 1992. Vol. 58. - P. 19.
DIFFERENTIAL CALCULUS OF INTERVAL ELEMENTARY FUNCTIONS AND DECISION-MAKING
V.I. Levin
Penza State Technological University, Penza, Russia vilevin@mail.ru
Abstract
The article deals with problems, related to calculation of derivatives of interval-specified functions. These problems are relevant in study of systems with any level of uncertainty (nondeterministic systems). Specifically we speak about simple systems described by elementary interval-specific functions. Accordingly we solved the problem of calculating derivatives of elementary interval-specified functions. Previously obtained formulas and methods of finding of derivatives of any intervally defined functions are used. Basic definitions, related to derivatives of the interval-specified functions, are given and formulas of two types allowing calculation of interval derivatives are presented. Formulas of the first type express derivatives in the closed interval form, which requires computation using the apparatus of interval mathematics. But formulas of the second type express derivatives in open interval form, i.e. in form of two formulas. Formulas above expresses the lower and the upper limits of the interval representing the derivative. Here the calculation of the derivative of interval-defined function is reduced to computation of two ordinary certain functions. Using the aforementioned mathematical apparatus we find the derivatives of all elementary interval functions: interval constant, interval power function, interval exponential function, interval logarithmic function, interval natural-logarithmic function, interval trigonometric functions (cosine, sine, tangent, cotangent), interval inverse trigonometric functions. Formulas of all the derivatives are shown in form of an open interval. The difference between derivatives of interval elementary functions and the derivatives of normal (i.e. noninterval) elementary functions is discussed.
Keywords: interval value, interval function, interval derivative, interval computations, interval-differential calculus, interval elementary functions.
Citation: Levin VI, Differential calculus of interval elementary functions and decision-making [In Russian]. Ontology of designing. 2016; 6(21): 340-354. DOI: 10.18287/2223-9537-2016-6-3-340-354.
References
[1] Gnedenko BV. Kurs Teorii Veroyatnostey [In Russian]. - Moscow: Nauka, 2004. - 350 p.
[2] Zadeh LA. The Concept of a Linguistic Variable and its Application to Approximate Reasoning // Information Sciences. - 1975. - № 8; 9. - P. 199-249, 301-357; 43-80.
[3] Alefeld G, Herzberger J. Introduction to Interval Computation. - N.Y.: Academic Press, 1983 - 352 p.
[4] Levin VI. The interval derivative and the basis of nondeterministic differential calculus [In Russian]. Ontology of designing. 2013; № 4 (10): 72-84.
[5] Levin VI. Interval-differential calculus and some of its applications [In Russian]. Informacionnye Tehnologii. 2014; № 7: 3-10.
[6] Fihtengolz GM. Differential and integral calculus [In Russian]. V. 1. - M.: Fizmatlit, 2001. - 616 p.
[7] Levin VI. Differential calculus for interval-specific functions [In Russian]. Heuristic algorithms and distributed computing. 2015; 2(2): 8-25.
[8] Voschinin AP, Sotirov GR. Optimization under conditions of uncertainty [In Russian]. - M: MEI Publishers, 1989. -224 p.
[9] Aschepkov LT, Davydov DV. Universal solution of interval governance and optimization problems [In Russian]. -M.: Nauka, 2006. - 285 p.
[10] Moore RE. Interval Analysis. - N.Y.: Prentice-Hall, 1966. - 230 p.
[11] Libura M. Integer Programming Problems with Inexact Objective Function. Control and Cybernetics. - 1980. Vol. 9. - № 4. - P. 189-202.
[12] Kurzhanskiy AB. Identification Problem - Theory of Guaranteed Estimates [In Russian]. Automation and Remote Control. - 1991. Vol. 52. - № 4. - P. 447-465.
[13] Hyvonen E. Constraint Reasoning Based on Interval Arithmetic: the Tolerance Propagation Approach // Artificial Intelligence. - 1992. Vol. 58. - P. 19.
Сведения об авторе
Левин Виталий Ильич окончил Каунасский политехнический ин-т, Открытый ун-т Израиля. Доктор технических наук, профессор, PhD, Full Professor, заведующий кафедрой математики (1975-2000), советник рек-тора по науке (2006-2011) Пензенского гососударственного технологического университетата, профессор Московского университета им. С.Ю. Витте (с 2003 г.). В списке научных трудов сотни работ (в том числе десятки монографий) по логике; математическому моделированию в технике, экономике, социологии, принятию решений; оптимизации; теории надежности; истории науки; проблемам образования. Действительный член МАИ, ЕАИ, МАНЭБ и АСН, заслуженный деятель науки Российской Федерации, лауреат международных премий «Соросовский профессор», Международный эксперт в области социологии конфликта и рейтингования униерситетов. Vitaly I. Levin graduated from Kaunas Politechnical Institute, Open University of Israel, Doctor of Engineering Science, Professor, PhD, Full Professor. Head of Mathematics Department (1975-2000), the scientific counselor of rector (2006-2011) of Penza State Technological University, professor of Moscow University named after S.J. Vitte (since 2003). He is the autor of hundreds of publications (among them dozens of monographs), in logic; mathematical modelling of engineering, economics, sociology; optimization; reliability; history of science; education problems. The member of IIA, EIA, IAELP, ASS. Honoured scientist of Russia, Laureate of International Prizes "Soros Professor", International Reviewer in Sociology of conflicts and Universiry ranking.
