Интервальная производная и интервально-дифференциальное исчисление Текст научной статьи по специальности «Математика»

УДК 62-50:519.7/8

ИНТЕРВАЛЬНАЯ ПРОИЗВОДНАЯ И ИНТЕРВАЛЬНО-ДИФФЕРЕНЦИАЛЬНОЕ ИСЧИСЛЕНИЕ

© В.И. Левин

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

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

1. ВВЕДЕНИЕ

а ° Ь = {а ° Ь | а є а,Ь є Ь} .

(2)

Проектирование и исследование свойств разнообразных систем обычно опирается на тот или иной подходящий математический аппарат. К настоящему времени создано большое число различных видов математического аппарата. Однако при всем их различии почти все они обладают одним общим свойством - применимостью только к полностью определенным (детерминированным) системам. В то же время встречающиеся на практике системы обычно характеризуются той или иной степенью неопределенности (недетерми-нированы). С целью построения и исследования таких систем чаще всего применяют математический аппарат теории вероятностей [1], теории нечетких множеств [2], интервальной математики [3]. В настоящей работе впервые предлагается новый математический аппарат для исследования недетерминированных систем - недетерминистское дифференциальное исчисление. Этот аппарат, в отличие от трех названных выше, нацеленных в основном на исследование статических систем, применим также к изучению динамических систем.

2. ВСПОМОГАТЕЛЬНЫЕ МАТЕМАТИЧЕСКИЕ СВЕДЕНИЯ

Мы будем использовать в качестве вспомогательных сведений прежде всего основные математические сведения из алгебры интервальных чисел [3-4]. В этой алгебре в качестве операндов берутся замкнутые вещественные интервалы, определяемые как множества всех вещественных чисел между нижней и верхней границами интервала, включая сами эти границы:

а = [аі,а2І = {а | а — а — а2} .

(1)

Эти операнды естественно называть интервальными числами. Операции ° над интервальными числами

у = [аг, а2], Ь = [Ьъ Ь2] вводятся как прямые теоретико-множественные обобщения соответствующих операций над вещественными числами а, Ь , т. е.

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

у + Ь = {а + Ь | а є у, Ь є Ь}, у - Ь = {а - Ь | а є а,Ь є Ь}, к ■ а = {к ■ а | а є у}, (3)

у ■ Ь = {а■ Ь | а є а,Ь є Ь}, 07Ь = {а/Ь | а є а,Ь є Ь}.

На основе определений (3) операций над интервальными числами можно вывести следующие формулы для вычисления результатов этих операций [3]:

у + Ь = [а1, а2] + [Ь1, Ь2] = [а1 + Ь1, а2 + Ь2], у-у = [аl,а2] -[Ь1,Ь2] = [а1 -Ь2,а2 -Ь1],

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

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

(4)

1[ка2,ка1], к < 0, у ■ Ь = [а1,а2] ■[Ь1,Ь2] = [шіп( аг- ■ Ь.■ ),тах(аг- ■ Ь.■)],

и' и'

о/у = [al, а2] /[Ь1,Ь2] = [a1, а2] ■ [1/Ь2,1/Ь1].

Также в качестве вспомогательного нам потребуется понятие интервальной функции [5], которая вводится как однозначное отображение множества замкнутых вещественных интервалов {у}, у = [х1,х2] вида (1) на множество замкнутых вещественных интервалов {Xх}, у = [л у2] этого же вида. Символически интервальная функция записывается в виде

У = /(а ) =

(5)

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

Введем теперь понятие предела интервальной функции (5). Рассмотрим независимую переменную х = [хь х2] этой функции. Будем говорить, что переменная х в процессе своего изменения неограниченно приближается к некоторому предельному интервалу Х0 = [д^ц,х02], если в указанном процессе х1 неограниченно приближается к х01, а х2 неограниченно приближается к х02. Символически неограниченное приближение х к х~0 показывается как

(~ ^ Хо) = (х1 ^ Х01, Х2 ^ Х02) = (Ншх1 = Хо1,Ншх2 = Х02). (6)

Совершенно аналогично независимой переменной х зависимая переменная ~у интервальной функции (5) в процессе своего изменения может неограниченно приближается к некоторому предельному интервалу

~0 = [У01У02], т. е.

с- —' Л)) = (y1 —'y01,y2 —'y02) .

(7)

Ay

Ax

.а - а f (а) - f (ao)

-x0

X - Xo

(11)

Возьмем предел отношения (11) при неограниченном приближении независимой переменной ~х к ее фиксированному значению ~0 :

lim Ay

а —>~o Ax

(12)

Предел (12), если он существует, будем называть производной интервальной функцией от исходной интервальной функции (5) в точке уо или, коротко, интервальной производной от функции (5) и обозначать УУ , или (УХ) . Таким образом,

~ ~ ~ Лу

yx0 = f5?0(x) = ~lima та 0 0 x—x0 Ax

(1З)

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

lim у = уо или, по-другому Jim f (x) = уо . (8)

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

alim f (x) = f (xo) .

x —a0

(9)

3. ИНТЕРВАЛЬНАЯ ПРОИЗВОДНАЯ ФУНКЦИЯ

Рассмотрим теперь произвольную интервальную функцию (5). Будем считать ее непрерывной. Зафиксируем некоторое значение х0 = [х01, х02] независимой переменной. Этому значению, в силу непрерывности нашей функции, будет соответствовать некоторое фиксированное значение функции у0 = /(~0) .

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

= x -^, Лу = у -У0 = f(x) -f(x0)

(1G)

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

Условие существования интервальной производной интервальной функции определяется следующей теоремой.

Теорема 1. Для того чтобы в точке ~0 существовала интервальная производная от интервальной функции (5), определяемая формулами (10)-(13), необходимо и достаточно, чтобы в некоторой окрестности этой точки, включая ее саму, все значения независимой переменной х функции были невырожденными интервалами (т. е. интервалами с несовпадающими верхней и нижней границами).

Доказательство. Из выражения (11) интервальной функции, предел которой есть интервальная производная в точке х0 , видно, что эта производная существует только тогда, когда в некоторой окрестности указанной точки, включая ее саму, все возможные значения знаменателя выражения (11) не равны нулю. Но знаменатель выражения (11), согласно формуле (4) разности интервалов, равен интервалу

х _ х0 = [ХЬ х2] _ [x01, х02] = [х1 _ x02, х2 _ х01] .

Правый интервал равен нулю (нулевому интервалу [0,0] ) только при условии

Х1 = ^ х2 = Х01, равносильном следующему

Х1 = Х01 = х2 = х02, или

х = Х0 = [ х, х],

что означает вырожденность интервалов х, Х0 .

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

X — X

X —X

0

0

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

Интервальная производная от интервальной функции (5), определяемая для произвольной точки х 0 формулами (10)-(12) в виде предела, может быть выражена также в конечном виде через значения независимой х и зависимой у = / (х) переменных этой функции в указанной точке. Именно, справедлива следующая теорема.

Теорема 2. Интервальная производная от непрерывной интервальной функции (5), определяемая для произвольной точки / формулами (10)-(13) в виде предела, может быть выражена в конечном виде через значения независимой х и зависимой у = / (х) переменных этой функции в указанной точке следующим образом:

= Д, (х) =

/ (х0) - / (х0) х0 - х0

(14)

Примечание 1. При первом взгляде выражение

(14) может показаться неопределенностью вида 0/0. Но это впечатление неверно, поскольку, согласно теореме 1, у любой существующей в точке х0 интервальной производной интервал х~0 невырожден, и потому, в соответствии с формулой (4) для разности интервалов, разность Х0 - ~0 не равна нулю (нулевому интервалу [0,0] ). Точно также не равна нулю разность

/(х0)-/(х0) .

Доказательство. Согласно определению интервальной производной в точке / (10)-(13) она может быть записана в виде предела

(х) =

/ (х ) - / (~0)

_Нш ~ ~

~^Х0 Х — Х0

(15)

В процессе предельного перехода в правой части равенства (15) ~х неограниченно приближается к ~х0 , а

/ (~х ) , в силу непрерывности интервальной функции / - к /(/). Таким образом, предел в правой части

(15) равен

Нш

/(х) — /(/0) _ /(/0) — /(/0)

(16)

х0 — х0

Подставив полученное значение предела из равенства (16) в выражение (15), будем иметь формулу (14). Что и требовалось доказать.

Итак, выражение интервальной производной функции (14) представляет ее через исходную (первообразную) интервальную функцию (5). При этом, поскольку в (14) точка х0 , в которой определяется производная, произвольна, выражение (14) можно записать в общем виде так:

У~ = Я (х ) =

/ (х ) — / (х )

(17)

где х - произвольное значение независимой переменной непрерывной интервальной функции (5) из ее области определения, при котором интервальная производная этой функции существует. Как видно из формулы (17), интервальная производная выражается непосредственно через исходную (первообразную) интервальную функцию простой алгебраической формой, что, конечно, упрощает процесс вычисления производной. Этот эффект отсутствует при нахождении обычных производных, фигурирующих в традиционном дифференциальном исчислении детерминированных функций [6].

4. ИНТЕРВАЛЬНЫЕ ПРОИЗВОДНЫЕ ВЫСШИХ ПОРЯДКОВ

Производная /~ (х) от интервальной функции

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

ух ~ /х(х) [/х(х)]х

(18)

затем третью производную / ( ) (производную от второй производной / ( ) )

ух _ /х(х) [/Х(х)]Х

(19)

и т. д., вплоть до интервальной производной любого и-го порядка, определяемой как производная от производной п — 1 порядка

у!п) = /х(п)(Х) = [/(^(Х]/

(20)

Условие существования интервальной производной любого и-го порядка определяется нижеследующей теоремой.

Теорема 3. Для того чтобы в точке / существовала интервальная производная и-го порядка (п = 2,3,...) от интервальной функции (5), определяемая формулами (18)-(20), необходимо и достаточно, чтобы в некоторой окрестности этой точки, включая ее саму, все значения независимой переменной функции были невырожденными интервалами.

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

—

0

—

0

0

функция. Поэтому, согласно теореме 1, для существования в точке х0 производной от производной 1 -го порядка функции (5) (т. е. производной 2-го порядка от интервальной функции (5)) необходимо и достаточно выполнения того же условия, что и для существования в точке х 0 производной 1-го порядка от функции (5), т. е. невырожденности всех интервалов, являющихся значениями х функции (5) в некоторой окрестности точки х 0 . Продолжая рассуждения, придем к необходимым и достаточным условиям существования в точке х0 производной п -го порядка от функции (5) в той форме, в которой они даны в теореме 3. Что и требовалось доказать.

Интервальная производная любого п -го порядка от интервальной функции (5), определяемая формулами (18)-(20) итеративно, может быть выражена в конечном виде аналогично выражению (17) интервальной производной 1 -го порядка. Это основано на следующей теореме.

Теорема 4. Интервальная производная п -го порядка (п = 2,3,...) от интервальной функции (5) может быть выражена в конечном виде через интервальную производную (п — 1) -го порядка от этой функции (если она непрерывная) и независимую переменную х в следующем виде:

У~п) - /x(n)(X )=

/~(п-1)( X) - /ІЙ-1)(Х)

п = 2,3,... (21)

Доказательство теоремы получается с помощью формулы (17), если подставить в нее в качестве функции /(х) производную (п — 1) -го порядка /~п—1)(х) и

учесть еще, что производная 1-го порядка от /~(п—1)(Х) согласно определению (20), есть производная п -го порядка /Хп\х).

С помощью теорем 2 и 4 можно последовательно получить конечные выражения для интервальной производной 2-го порядка:

Ух - /у(х) =

/X(х)- Я(х)

/(X ) - /(~ ) /(~ ) - /(х )

для интервальной производной 3-го порядка:

(22)

УХ - /"(X) -

/Х( X) - Жу) _

х-х ______ __ (23)

/( )-/( ) /( )-/( ) /( )- /( ) /( )- /( )

X — X_______X — X

-

x — x_______X-X

-

и т. д. Отметим, что все эти выражения аналогичны конечному выражению (17) для интервальной производной 1-го порядка.

Примечание 2. Выражения (22), (23),... интервальных производных высших порядков кажутся с первого

взгляда суперпозициями неопределенностей вида 0/0 , однако, по причине, указанной в примечании 1, не являются таковыми. Заметим также, что эти выражения не могут быть преобразованы к более простому виду с помощью эквивалентных преобразований, поскольку в алгебре интервальных чисел [3] не выполняются эквивалентности типа

(а + Ь )с = а с + Ьс;

а + Ь

а Ь

(24)

5. ВЫЧИСЛЕНИЕ ИНТЕРВАЛЬНЫХ ПРОИЗВОДНЫХ

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

Теорема 5. Интервальная производная от непрерывной интервальной функции (5), выражаемая для произвольной точки формулой (17) в виде суперпозиции операций над интервальными переменными, может быть также представлена явно в виде интервала, нижняя и верхняя границы которого выражены через нижнюю и верхнюю границы интервальных независимой ~ = [х1,Х2] и зависимой / = [_ц,У2] переменных этой функции в точке следующим образом:

УЇ - УК-О =

У2 - У1 У2 - У1 х2 - хг’ х2 - Х1

(25)

Доказательство. Будем исходить из выражения интервальной производной от интервальной функции (5) формулой (17). Представим указанную функцию в явном виде интервала

/(л0 = /2(x)] - [УьУ2] - У ■

Здесь

У1 = /l(x), У2 = /2(х)

(26)

(27)

есть соответственно нижняя и верхняя границы интервальной зависимой переменной функции (5). Аналогично представим интервальную независимую переменную :

x = [xl,x2] .

(28)

X - X

X - X

X - X

X - X

X - X

После проведения необходимых подстановок выражение (17) интервальной производной от интервальной функции (5) примет вид

уу - Я(У) =

У2 - У1 У2 - У1

X2 - X: X2 - X!

(36)

ух - /х(а) =

У -У [Уl,У2] -[УъУ2]

X - X ^1, X2] - [Xl, X2]

(29)

Разности интервалов в числителе и знаменателе выражения (29) по формуле (4) представим в виде интервального числа

[Уи У2] - [Уи У2] = [У1 - У2, У2 - Уl], [x1, x2] - [x1, х2] = [Х1 - ^ x2 - x1].

(30)

Подставив выражения (30) в (29), получим представление интервальной производной в виде частного двух интервалов

УХ - /У(x):

[У1 - У2, У2 - У1] [х1 - x2, x2 - х1]

(31)

В свою очередь, выражение (31) можно, согласно формуле (4), представить в виде следующего произведения двух интервалов:

Ух - Я (х) =

= [У1 - У2, У2 - У1] ■ [1/(х2 - х1),1/(х1 - х2)].

(32)

Перемножая интервалы в правой части (32) по соответствующему правилу (4), получим

Ух - Жх) =

У1 -У2 У1 - У2 У2 - У1 У2 - У1

х2 - Х1 Х1 - х2 х2 - Х1 Х1 -х2

г У1 - У2 У1 - У2 У2 - У1 У2 Л - У1

х2 - Х1 ’Х1 - х2 х2 - Х1 ’Х1 - х2

(33)

что и требовалось доказать.

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

Теорема 6. Интервальная производная любого и-го порядка от непрерывной интервальной функции

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

= [УіП-1), У(п-1)] от интервальной функции (5) следующим образом:

у£° - /(п)(х) =

У(п-1) - У(п-1) У(п-1) - У(п-1)

Х2 - Х1

Х2 - Х1

(37)

Доказательство теоремы 6 аналогично доказательству теоремы 5, с той разницей, что в качестве исходного выражения искомой интервальной производной п -го порядка от интервальной функции (5) берется формула (21), тогда как в случае теоремы 5 исходным выражением искомой интервальной производной 1 -го порядка от функции (5) служит такая же по форме формула (17).

Используя формулу (37), нетрудно получить простые явные выражения типа (25) для интервальных производных 2-го, 3-го и т. д. порядка от интервальной функции (5). Действительно, для интервальной производной 2-го порядка из (37) имеем:

Одинаково подчеркнутые члены в выражении (33) равны. Оставив из каждых двух равных членов по одному, представим выражение (33) в более простом виде:

УУ - /х(Х) =

У2 - У1 У2 - У1

Х2 - X Х0 - X-

2-

х1

(38)

Ух - /х(а) =

У2 - У1 У2 - У1 Х9 - Х1 Х1 - X'

х2

2-

Х1 Х1 ■

х2

.(34)

Два члена в круглых скобках формулы (34) различаются лишь знаком, причем (с учетом того, что Х2 > Х1, У2 > У1) левый член положителен, а правый -отрицателен. Так, из (34) получается простейшее выражение интервальной производной от интервальной функции (5):

Подставляя в соотношении (38) значения нижней и верхней границ интервальной производной 1 -го порядка из (25) в виде

У2 -У1 У2 -У1

У1 =-----------, У2 =--------------------

(39)

с учетом необходимых эквивалентных преобразований, получим:

Ух- /у (х) =

У2 - У1 У2 - У1

Х1'

х2 х2 "

х1

(35)

С учетом того, что Х1 - х2 =-(х2 - Х1), выражение (35) окончательно переписывается в виде

ух- Я (х) =

2(У2 - У1) 2(У2 - У1) (х2 - х1)2 (х2 - х1)2

(40)

Аналогично для интервальной производной 3-го порядка из (37) имеем:

2 Л1

2

Таблица 1

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

Функция Обозначение функции Явное выражение функции

Исходная интервальная функция У = к(хХ х = [хі5 х2]У = [УьУ2] [уі>у2] = /[хі5Х2] или [Уі5 У2] = [/і (хі5 Х2)5к2(хі5 Х2)]

Интервальная производная функция 1-го порядка ух = /Х(х) или У(і) = /і)(Х) Г У2 - Уі У2 - Уі _ Х2 - Хі ’ Х2 - Хі _

Интервальная производная функция 2-го порядка ух - 1х (х) или уХ2) - /Х2)(х) 2(У2 - Уі) 2(У2 - Уі) 2 5 2 _ (х2 -Хі) (х2 -Хі) _

Интервальная производная функция 3-го порядка уі - /?(Х) или уі) - /3)(Х) 4(У2 - Уі) 4(У2 - Уі) 3 5 3 _ (Х2 -Хі) (Х2 -Хі) _

Интервальная производная функция п -го порядка УХп) - 7~п)(Х) Г 2п-і( У2 - Уі) 2п-і( У2 - Уі)" _ (х2 -Хі)П (х2 -Хі)П _

ух - к(х) =

У 2 - Уі У2 - У1

Х2 - X Х2 - Хі

(41)

Как и в предыдущем случае, подставляя в формулу (41) значения границ интервальной производной 2-го порядка из (40) в виде

Уі =

2(У2 - Уі) (х2 - х1)2

У2 =

2(У2 - Уі) (х2 - хі)2

(42)

будем, после необходимых эквивалентных преобразований, иметь:

ух -IX(Х)=

4(У2 - Уі) 4(У2 - Уі) 3 5 3

(х2 -хі) (х2 -хі)

(43)

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

Теорема 7. Интервальная производная любого п -го порядка от непрерывной интервальной функции

у = / (X) вида (5), выражаемая для произвольной точки X формулой (21) в виде суперпозиции операций над интервальными переменными, может быть представлена в указанной точке также в явном виде интервала с явно выраженными нижней и верхней границами вида (25), (40), (43)

уХп) - /~п)( х )=

2п-і( У2 - Уі) 2п-і( У2 - Уі) (х2 - хі)П (х2 - хі)П

-Ь

(44)

Здесь х1, х2 - нижняя и верхняя границы интервальной независимой переменной х = [х1, х2] в точке

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

Сводка полученных явных выражений интервальных производных функций различных порядков приведена в табл. 1.

6. ВОЗМОЖНЫЕ ПРИМЕНЕНИЯ

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

ных значений данных параметров; обычные производные в исходной детерминированной системе заменяются соответствующими интервальными производными, в результате получаем математическую модель собственно изучаемой системы, все параметры которой определяются с точностью до интервалов возможных-значений; в этой модели изучаемые характеристики системы принимают интервальные значения и могут быть представлены в виде суперпозиции интервальных операций (3) над интервалами - значениями параметров звеньев системы; 3) используя соотношения (4), выражающие результаты элементарных преобразований интервалов, представляем интервальнозначные характеристики недетерминированной системы в явной форме интервала ¥ = к /2], нижняя /і и верхняя /2 границы которого являются обычными (детерминированными) функциями от нижних и верхних границ интервалов - параметров звеньев системы. После этого математическое моделирование изучаемой недетерминированной системы по любой ее интервальной характеристике ¥ = [/і5 /2] сводится к анализу двух обычных детерминированных функций /і5 /2 .

Пример. Рассмотрим простейшее применение предложенного математического аппарата в экономике. Пусть у = /(х) - некоторая детерминированная функция, показывающая зависимость одного экономического показателя у от другого х [7]. Тогда эластичностью ЕХ (у) этой функции называется предел относительного приращения функции У к относительному приращению аргумента (независимой переменной) х при Ах ^ 0 , т. е.

Ех (У) = 1іт I Ау /— 1 = Х 1іт Х • у'х-

Ах^-ОІ у I х ) у Ах^-0 Ах у

(45)

Здесь у'х - производная от функции у = /(х) по переменной х . Эластичность функции показывает приближенно, на сколько процентов изменится значение функции у = /(х) при изменении значения аргумента х на 1 %. Например, если функция у = /(х) показывает зависимость спроса у на товар от его цены х , то эластичность этой функции Ех (у) есть коэффициент, приближенно показывающий, на сколько процентов изменится спрос на товар при изменении его цены на 1 %. Аналогично, если функция у = /(х) показывает зависимость издержек производства у от объема выпускаемой продукции х , то эластичность этой функции Ех(у) показывает приближенно, на сколько процентов изменятся издержки производства при изменении объема выпускаемой продукции на 1 %. Наконец, если изучаемая нами функция у = / (х) есть зависимость себестоимости единицы продукции у от стоимости всей выпускаемой продукции х , то эластичность данной функции Ех (у) показывает приближенно, на сколько процентов изменится себестоимость единицы продукции при изменении стоимости всей продукции на 1 % и т. д.

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

(интервальную) функцию ~ = / (X) вида (5), где ~ = [хі5 х2] - независимая интервальная переменная (интервальный аргумент), а у = [у^ у2] - зависимая

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

1. Исходная математическая модель идеальной (детерминированной) системы уже есть, она описывается функцией эластичности системы (45) с точными значениями переменных х и у .

2. Заменяем в математической модели, установленной на шаге 1, точные значения параметров х и у соответствующими интервалами их возможных значений X = [хі5 х2] и у = [уі5 у2], а детерминированную производную 1-го порядка у'Х - соответствующей интервальной производной Ух , явное выражение которой в виде интервала дается формулой (25). В результате получаем математическую модель искомой недетерминированной экономической системы - ее интервальную функцию эластичности - в виде следующей суперпозиции интервальных операций над интервалами -параметрами системы:

Ех (у )=Х • Ух = ^

У [Уі5 У2]

У2 - Уі У2 - Уі

Х2 - X Х2 - Хі

(46)

3. После вычисления выражения (46) по формулам (4) элементарных преобразований интервалов с учетом положительности значений хьХ25Уі5У2 в экономических системах получим окончательно следующую математическую модель искомой системы

Ех (У) =

х2( У2 - Уі) х2( У2 - Уі) Уі(х2 -Хі)5 Уі(х2 -Хі)

(47)

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

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

ЛИТЕРАТУРА

1. Гнеденко Б.В. Курс теории вероятностей. М.: Наука, 2004.

2. Заде Л.А. Понятие лингвистической переменной и его применение к принятию приближенных решений. М.: Мир, 1976.

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

4. Левин В.И. Интервальные методы оптимизации систем в условиях неопределенности. Пенза: Изд-во Пензен. технол. ин-та, 1999.

5. Левин В.И. Оптимизация в условиях интервальной неопределенности. Метод детерминизации // Автоматика и вычислительная техника. 2012. № 4.

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

7. Кремер Н.Ш. и др. Высшая математика для экономистов. М.: ЮНИТИ, 2001.

Levin V.I. INTERVAL DERIVATIVE AND INTERVAL DIFFERENTIAL CALCULUS

The possibility of generalization of the classical differential calculus to functions which variables have interval uncertainty is shown. The concept of the derivative of function is introduced. The formulas representing explicitly interval derivatives of any order are given.

Key words: interval value; interval function; interval derivative; interval computation; nondeterministic differential calculus.

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

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

Levin Vitaliy Ilich, Penza State Technological Academy, Penza, Russian Federation, Doctor of Technics, Professor, Honored Scientist of RF, е-mail: [email protected]