Математические и инструментальные методы экономики II TT TT Т Актуальные проблемы экономики и права. 2016. Т. 10, № 4 Mathematical and instrumental methods in economics ¡¿¡^fafa Actual Problems of Economics and Law, 2016, vol. 10, No. 4
........................................................................... ISSN 1993-047Х (Print) / ISSN 2410-0390 (Online) ...........................................................................
УДК 517:330.4 DOI: http://dx.doi.org/10.21202/1993-047X.10.2016.4.77-87
В. В. ТАРАСОВА1 В. Е. ТАРАСОВ1
1 Московский государственный университет имени М. В. Ломоносова, г. Москва, Россия
ДЕТЕРМИНИРОВАННЫЙ ФАКТОРНЫЙ АНАЛИЗ: МЕТОДЫ ИНТЕГРО-ДИФФЕРЕНЦИРОВАНИЯ НЕЦЕЛОГО ПОРЯДКА
Цель: обобщение методов детерминированного факторного экономического анализа, а именно метода дифференциального исчисления и интегрального метода.
Методы: математические методы интегро-дифференцирования нецелого порядка, теория производных и интегралов дробного (нецелого) порядка.
Результаты: сформулированы базовые понятия, и разработаны новые методы, позволяющие учитывать эффекты памяти и нелокальности при количественном описании влияния отдельных факторов на изменение результативного экономического показателя. Предложены два метода интегро-дифференцирования нецелого порядка для детерминированного факторного анализа экономических процессов с памятью и нелокальностью. Показано, что метод интегро-дифферен-цирования нецелого порядка может давать более точные результаты по сравнению со стандартными методами (методом дифференцирования, использующим производные первого порядка, и интегральным методом, использующим интегрирование первого порядка) для широкого класса функций, описывающих результативные экономические показатели. Научная новизна: предложены новые методы детерминированного факторного анализа: метод дифференциального исчисления нецелого порядка и интегральный метод нецелого порядка.
Практическая значимость: основные понятия и формулы статьи могут быть использованы в научной и аналитической деятельности для факторного анализа экономических процессов. Предлагаемый метод интегро-дифференциро-вания нецелого порядка расширяет возможности детерминированного факторного экономического анализа. Новый количественный метод детерминированного факторного анализа может стать началом количественных исследований поведения экономических агентов с памятью, эредитарностью и пространственной нелокальностью. Предлагаемые методы детерминированного факторного анализа могут быть использованы при изучении экономических процессов, подчиняющихся степенным законам, в которых показатели (эндогенные величины) являются степенными функциями факторов (экзогенных величин), включая процессы, описываемые производственной функцией Кобба - Дугласа, поскольку эти методы позволяют точнее описывать суммарное влияние факторов по сравнению со стандартным методом. Предлагаемые методы могут быть использованы при изучении экономических процессов, описываемых уравнениями со степенной нелокальностью в факторном пространстве и в пространстве состояний.
Ключевые слова: математические и экономические методы экономики; факторный анализ; метод дифференциального исчисления; интегральный метод; процессы с памятью; эредитарность; производная нецелого порядка; интегрирование нецелого порядка
Как цитировать статью: Тарасова В. В., Тарасов В. Е. Детерминированный факторный анализ: методы интегро-дифференцирования нецелого порядка // Актуальные проблемы экономики и права. 2016. Т. 10, № 4. С. 77-87. DOI: 10.21202/1993-047Х.10.2016.4.77-87
V_У
Введение ного исчисления и метод интегрального исчисления,
В экономическом анализе активно используются которые базируются на математической теории произ-
методы детерминированного факторного анализа, ко- водных и интегралов целого порядка. В данной работе
торые дают точную характеристику влияния факторов предлагается обобщение метода дифференциального
на изменение результативного показателя. Среди ос- исчисления и интегрального метода факторного ана-
новных методов выделяются метод дифференциаль- лиза. В стандартном подходе к факторному анализу ме-
ISSN 1993-047Х (Print) / ISSN 2410-0390 (Online) ...........................................................................
тодом дифференциального исчисления используются производные целого (первого) порядка. В современной математике, помимо производных и интегралов целых порядков, хорошо известны производные и интегралы нецелых (дробных) порядков [1—5]. Данный математический аппарат позволяет описывать различные типы процессов, характеризующихся памятью и нелокальностью, что активно применяется в настоящее время в естественных науках [6-8]. В данной статье анализируется возможность применения данного математического инструмента для количественного описания влияния отдельных факторов на изменение результативного экономического показателя.
Кратко опишем стандартный метод дифференциального исчисления. Математической основой этого метода описания влияния отдельных факторов на изменение результативного показателя являются операция взятия производной (дифференцирование) и теорема Тейлора. Для функции f (x, y) от двух переменных x и y рассматриваются частные производные первого порядка по x и по у:
W): = M (1)
Метод дифференциального исчисления основан на формуле Тейлора. Например, формула Тейлора для функции от двух переменных z = f (x, у) может быть записана в виде:
fix, у) = /Оо.Уо) + /*(1)Оо-Уо) • А* +
+ /у(1)(*о,Уо)-Ду + К2(х,у),
(2)
где Ух(1)С^О'Уо) и $1}(.хо>Уо) - значения частных производных первого порядка по х и у соответственно в точке (х0, у0); Дх и Ду - факторные приращения соответствующих переменных (Дх := х — х0 , Ду := у — у0); Я2(,х,у') - остаточный член. Например, в асимптотической форме (форме Пеано) остаточный член имеет вид Я2(*.у): = о + (Ду)2)
является бесконечно малой величиной (бесконечно малой функцией в окрестности точки (х0 , у0)).
В результате, пренебрегая остаточным членом, полное приращение функции йг ■•— f(x1,y1) — f(x0,y0) записывают через приращение факторов Ах и Ау в виде:
Аг * (хо, Уо) • Ах + /у(1) (х0, у0) • Ду. (3)
Влияние фактора х на обобщающий показатель вычисляется по формуле:
Агх = ^(хо^-Ах, (4)
а влияние фактора у - с использованием формулы:
Агу=$1\ха,уа)-Ау, (5)
где х0 и у0 - базисные (плановые) значения факторов х и у, оказывающих влияние на результативный показатель; х иу1 - фактические значения этих факторов Дх: = х-1~ х0, и Ду: = у1—у0 - абсолютные изменения (отклонения) факторов х и у соответственно.
В методе дифференциального исчисления предполагается, что общее приращение функций Дг (результирующего показателя) разлагается на слагаемые йгх и Лгу. Значение каждого из этих слагаемых вычисляется как произведение соответствующей частной производной на приращение переменной (фактора), по которому берется эта производная. В этом методе неразложимый остаток (остаточный член R2 (х, у)), которым пренебрегают, интерпретируется как логическая ошибка метода дифференциального исчисления. Отбрасывание неразложимого остатка является одним из недостатков этого метода, поскольку для экономических расчетов часто требуется точный баланс изменения результативного показателя и алгебраической суммы влияния всех факторов.
В стандартном подходе к факторному анализу в методе дифференциального исчисления используются производные целого (первого) порядка. В современной математике, помимо производных целых порядков, хорошо известны производные нецелых (дробных, произвольных) порядков [1-5]. Производные нецелого порядка обладают набором нестандартных свойств. Например, нарушение стандартного правила Лейбница (правила дифференцирования произведения) является характеристическим свойством всех типов производных нецелого порядка. Стандартное выражение для дифференцирования сложной функции также нарушается. Производные и интегралы нецелых порядков активно применяются для описания различных физических процессов, характеризующихся памятью и нелокальностью [6-8]. Производные и интегралы нецелого порядка недавно были использованы для описания различных финансовых и экономических процессов в статьях [9-18] и наших работах [19-24].
В настоящее время известны различные виды производных произвольного (дробного, нецелого) порядка, предложенные Риманом, Лиувиллем, Сониным, Лет-никовым, Маршо, Риссом, Вейлем, Адамаром, Капуто
........................................................................... ISSN 1993-047Х (Print) / ISSN 2410-0390 (Online) ...........................................................................
[1, 3-5]. Для построения метода дробного интегро-дифференцирования произвольного (нецелого) порядка наиболее удобными являются производные Капуто [3, с. 90-99]. Основной отличительной особенностью этих производных является то, что их действие на постоянную функцию дает ноль. Это приведет нас к нулевому влиянию соответствующего фактора на обобщающий показатель в эредитарном случае. Приведем определения производной Капуто, данное в книге [3, с. 92].
Определение. Левосторонняя и правосторонняя производные Капуто порядка а > 0 определяются формулами:
Киа+гях) ■ г(п _ а) ^ (х _ ?)а_п+1, V
(па ЛОО := Г ^^ , (х<Ь), (7)
Г(п - а) I а-х)а~п+1 v 7
где Г (а) - гамма-функция, а < х < Ь, а f (п) (х) - производная целого порядка п = [а] + 1 от функции f (х) по переменной х. Предполагается, что функция f (х) имеет производные вплоть до (п-1) порядка, которые являются абсолютно непрерывными функциями на интервале [а, Ь].
Из формул (6) и (7) видно, что производная Капуто, несмотря на использование термина «производная», фактически является интегро-дифференциальным оператором со степенным ядром для нецелых значений а. Отметим, что для целых положительных значений а = п производная Капуто совпадает [3, с. 92] с обычной производной порядка п:
Ф1Мх)=^\х),фиШ) = (-1)"/Сп)(*). (8) В силу свойства (8) производные Капуто являются обобщением стандартных производных положительного целого порядка. При этом эти стандартные производные являются частными случаями производных дробного (произвольного) порядка.
Успешное применение в естественных науках производных Капуто нецелых порядков для описания процессов и систем с нелокальностью и памятью [6-8] позволяет ожидать, что данный математический аппарат позволит расширить возможности анализа финансовых и экономических процессов. В недавних работах [19-24] производные Капуто были применены для описания экономического поведения потребителей [19, 22, 24], вычисления предельной полезности для экономических процессов с памятью [20], описания
ценовой эластичности спроса с памятью [21], анализа внебиржевого кассового оборота валютного рынка [23]. В данной работе исследуются возможности дифференциального (и интегрального) исчисления дробного (нецелого) порядка по реализации детерминированного факторного анализа в экономике.
В данной статье предлагаются методы детерминированного факторного анализа, которые могут быть использованы при анализе экономических показателей (эндогенных величин), являющихся степенными функциями факторов (экзогенных величин) и описываемых уравнениями со степенной нелокальностью. Степенные закономерности и функции широко представлены в экономическом анализе, включая хорошо известные функции Кобба - Дугласа [25] и множество других законов, описанных в обзорах [26, 27]. Нелокальные свойства могут описываться интегральными и интегро-дифференциальными уравнениями. Для описания нелокальности степенного типа применяются уравнения с интегралами и производными нецелого порядка в силу содержания в этих операторах степенных ядер. Отметим, что применение интегрирования нецелого порядка к степенной функции Кобба - Дугласа дает функцию Кобба - Дугласа с другими показателями степени. Для экономических показателей (эндогенных величин), резко меняющихся в пространстве факторов, следует учитывать влияние показателей в удаленных точках этого пространства на экономическую среду и состояние экономики, в данной точке факторного пространства. Нелокальные экономические процессы в пространстве состояний (state space) рассматриваются, например, в работах [15, 17, 18].
Результаты исследования
1. Метод дифференциального исчисления произвольного (нецелого) порядка
Поскольку метод дифференциального исчисления базируется на формуле Тейлора, то для обобщения этого метода за счет использования теории производных дробного (нецелого) порядка нам понадобится обобщение формулы Тейлора для производной Ка-путо. Рассмотрим сначала обобщение ряда Тейлора, которое было предложено в статье [28].
Функцию f (х) для х > a можно разложить [28, с. 289] в обобщенный ряд Тейлора порядка 0 < а < 1 с левосторонними производными Капуто:
ISSN 1993-047Х (Print) / ISSN 2410-0390 (Online) ...........................................................................
+ ЯЫа(х,а+1 (9)
где Кца(х,а +) - остаточный член, который можно представить в виде:
RNaix, а +) =
T(Na +1)
ix-a)Na, (10)
где a < < x, и r(z) - гамма-функция. Например, для N = 2 имеем:
fix) = fia) + Щ^у ■ (х - а)а + R2a(x, а +). (11)
Для x < a формулу (11) нельзя применить, если параметр а не является целым. Данную проблему можно решить, если использовать правостороннюю производную Капуто порядка а. Используя свойства правосторонней производной [3, с. 95], функцию f (x) для x < a, если она определена в этой области, можно представить в виде в обобщенного ряда Тейлора порядка 0 < а < 1 с правосторонними производными Капуто:
/со=ш KVllT ■ о - +а -),(12)
где RNa(x, а—) - остаточный член, который можно представить в виде:
а = ((№ + 1)")"(а " ж)№' (13)
и x < < b. Для N = 2 имеем:
fix) = fia) + ■ О - х)а + R2aix, а -). (14)
Однако остается другая проблема для формул (9-14), обусловленная совпадением начальных и конечных значений в производных iD /) (а) (£"-/) (а). Данный недостаток приводит к ограничениям на применимость формул (9-14) в расчетах влияния факторов на изменение результативного показателя. Для решения этой проблемы предлагается обобщение формулы Тейлора (11) на случай, когда производная Капуто вычисляется для произвольной точки х0 > а. Рассмотрение ряда Тейлора (9) для функции и для самих производных Капуто в точке х = х0 позволяет получить требуемую формулу:
(££+/)(* о)
fix) = fix0)+-
Г( сс + 1)
разности R2aix, а—) и R2aix0, а—), и мы использовали обозначение:
Аах •■= ix - а)а - ix0 - а)а (16)
Аналогичные формулы можно получить и для ряда Тейлора правосторонней производной Капуто. Для упрощения можно будет рассматривать только случай левосторонней производной с начальной точкой а = 0, поскольку многие факторы могут описываться положительными действительными числами. В этом случае формулу (15) можно записать в виде:
(Доа+/)Оо)
fix) = fix о)+-
Г(а + 1)
Aax + R2aix,x0,0+), (17)
где х0 > 0, х > 0 и Аах := ха-хд .
Для функции от двух переменных г = /(х, у) формула Тейлора имеет вид:
fix,У)= ixo,Уo)+—т.. , .,4--Аах +
Г(а + 1)
(Рр+/)(*о.Уо)
■ Г(Р+1) ЛрУ+Д2а,2Р(^У,0+), (18) где К2а,2р(х>У> 0 +) - остаточный член.
Для а = в = 1 формула Тейлора (18) дает стандартную формулу (2).
В результате, пренебрегая остаточным членом Я2а,2^х,у), полное приращение функции Дг := fix1,y1) - - fix0,y0) запишем через приращение факторов х и у в виде:
Л (До"+Ж*о,Уо) . , К+/) (хР'Уо) (19) Лг Ж Г(а + 1) • ^+ Г(р + 1) ' ^ (19)
где х0, у0 - базисные (плановые) значения факторов х и у, оказывающих влияние на результативный показатель; х у1 - фактические значения этих факторов; Аах = х? — Хд и Дру = у" — Уд - обобщенные (степенные) абсолютные изменения (отклонения) факторов х и у соответственно.
В результате влияние фактора х на результативный показатель будет вычисляться по формуле:
Л (PqVXWO) л
Azr„ - ————--Аах,
• Дах + R2aix, х0, а +), (15)
Г(а + 1)
а влияние фактора у - с использованием формулы:
(20)
где х0 > а, х > а, остаточный член R2aix, х0, а +) равен
(DÜ+f)ix0,y0)
гср + 1) "Ару-
(21)
ISSN 1993-047Х (Print) / ISSN 2410-0390 (Online) ...........................................................................
Для а = в = 1 формулы (20) и (21), влияние факторов х и у, принимают стандартный вид (4) и (5) соответственно. Известно, что производная Капуто порядка а > 0 (и в > 0) постоянной функции равна нулю. В результате для функции f (х, у) = const мы получим нулевое влияние обоих факторов Azxa = Azyp = 0. Кроме этого, для функции, не зависящей от одного из факторов, влияние этих факторов будет равно нулю.
2. Сравнение со стандартным методом дифференциального исчисления
Отметим, что формулы Тейлора (2) без остаточного члена не могут дать точных результатов для степенных функций с нецелыми показателями. Формулы Тейлора (16) с производными нецелого порядка даже без остаточного члена являются более точным инструментом для аппроксимации степенных функций. Чтобы проиллюстрировать это утверждение, рассмотрим нелинейную степенную функцию fix) = саха + с0 , где а - нецелое положительное число. Если мы воспользуемся стандартной формулой Тейлора для разложения этой функции в окрестности точки х0 > 0, то получим бесконечный степенной ряд. Для этой функции выражение обобщенной формулы Тейлора (17) с производными Капуто можно получить, используя следующие формулы левосторонней производной Капуто для степенной функции [3, с. 95], имеющие вид:
„ ГГВ +1) „
{п-Ка<п,$>п-1), (22)
D%+(x-a)k = 0,(fc = 0,l.....n-1). (23)
В частности, имеем D%+1 = 0 и D%+(x — а)а = = Г(а + 1). Используя эти формулы, получим: Фо+П(.хо) = саГ(а+ 1). При этом производные Капуто более высокого порядка будут равны нулю, т. е. ((Do+)kf)(x) = 0 для k = 1, 2, ... Подставляя (0"+/)Оо) = саГ(ос +1) и /(х0) = сах$ + с0 в формулу Тейлора (17), получим:
/О) = caxg + с0 + са • Аах = саха + с0. (24)
В результате ряд Тейлора (24) дал саму функцию fix) — саха + с0. Следовательно, обобщенная формула Тейлора такой степенной функции дает точный результат (саму функцию), а не приближенное выражение.
В результате можно сделать вывод, что применение в факторном анализе метода дифференцирования
дробного (нецелого) порядка может давать более точные результаты, чем стандартные методы.
Примером степенной функции двух переменных является производственная функция Кобба - Дугласа P = P (L, C). Эта функция описывает зависимость объема производства P (production) от значений затрат труда L {labor) и капитальных затрат C (capital). Функция Кобба - Дугласа имеет вид: P(L, С) ■= А • La • С", где A - совокупная производительность факторов, показатели a и b интерпретируются как эластичности по капиталу и по труду соответственно (эластичности капитала и рабочей силы). Значения константы 0 < a < 1 и 0 < b < 1 определяются имеющимися технологиями. Функция Кобба - Дугласа была впервые предложена в статье [25, с. 151] для промышленности США за период 1899-1922 гг. в виде P (L, С) ■■= 1,01 • L3/4 • С1/4, т. е. A = 1,01, a = 0,75 и b = 0,25. Используя уравнение (20) для функции Кобба - Дугласа P(L, С) ~ А • La • Cb, мы получаем:
(DZ+P)(K,C) :=
(D<¡+P)(L,C) :=
А-Па + Ï) ,„_„ „ь
Г(а - а + 1) А-Г(Ь + 1)
¿а-а . сЬ> (25)
La ■ С*""3. (26)
Г(Ь - р + 1)
Если порядок производных Капуто брать равным показателям функции Кобба - Дугласа (а = а, в = Ь), то формулы (25) и (26) дают выражения ф£+Р)(Ь, С) := А ■ Г(а + 1) • Ср и (Оо+Р) {I, С):= А-Ьа- Г((3 + 1). Если порядок производных Капуто брать равным единице, то получаем стандартные выражения.
В результате влияния факторов Ь и С на результативный показатель Р, вычисленные по формулам (20) и (21), имеют вид:
ЛРга =
СDS+P)iL0,C0)
Г(а + 1) А ■ Г(а + 1) Г(а — а + 1) • Г (а + 1)
AaL =
1Га-С0й-Ла1, (27)
ДРс, р =
_ (о0р+р)а0,со) г(Р+1)
АрС =
А-Г(Ь+1)
r(b-p+i)-r(p+i)
L%-Cb0 Р'Л«С, (28)
где AaL = и ДрС = С? - сЦ.
........................................................................... ISSN 1993-047Х (Print) / ISSN 2410-0390 (Online) ...........................................................................
Если рассматривать порядки производных а = а и в = Ь, то тогда формулы (27) и (28) можно записать в виде: АРиа = А-СЦ- АаЬ и АР^ =А • Ь% ■ АЬС.
В качестве численного примера рассмотрим произ-
3 1
водственную функцию в виде Р(Ь, С) := 1,01 • Ь* ■ С* , т. е. А = 1,01, а = 0,75 и Ь = 0,25. Для простоты базисные (плановые) и фактические значения факторов выберем в виде четвертных степеней: Ь0 = 6,14, Ьг = 6,24 и С0 = 0,24, С-у = 0,54. В этом случае формулы (27) и (28) приводят к следующим значениям:
^,<х=а = А-С$- Аа1 = 2,2920940, (29) АРСЗ=Ь = А-Щ-АЬС = 68,775243. (30)
Суммарное влияние факторов Ь и С равно:
АРЬа + АРС# = 71,067337. (31)
При этом прирост АР = Р(/,1( Сг) — Р(!0, С0) результативного показателя Р, задаваемого производственной функцией Р{Ъ, С) •= А- Ьа • Сь, равен:
АР = А-С?-А-Ь^-СЦ = 74,505478. (32)
Суммарное влияние факторов, вычисленное методом дифференцирования дробного порядка, отличается от реального изменения результативного показателя АР приблизительно на 4,6 %. Дополнительный прирост результативного показателя от взаимодействия факторов определяется выражением:
6а>(3: = АР- {АРКа + АРС#) = 3,438141. (33)
Сравним теперь полученные результаты со стандартным методом дифференциального исчисления. Если порядок производных Капуто в формулах (27) и (28) приравнять единице (а = в = 1), то получим стандартные выражения стандартного метода дифференциального исчисления, задаваемые формулами (4) и (5). Для рассматриваемой задачи с производственной
3 1
функцией Кобба - Дугласа Р(Ь, С) := 1,01 • Ь* • С* влияние факторов Ь и С принимают значения:
АРи1 = А-а- 1а0-х ■ С0Й ■ АЬ = 2,310983484, (34)
АРС д =А-Ъ-Ь%- СЦ-1 • АС = 436,2929478, (35)
а соответствующее суммарное влияние факторов Ь и С равно:
АРЬ1 + АРСЛ = 438,6039313. (36)
Суммарное влияние факторов, вычисленное стандартным методом дифференциального исчисления, отличается от реального изменения результативного
показателя АР приблизительно в 4,9 раза, т. е. на 390 %, тогда как в предлагаемом методе отличие составляет 4,6 %. Стандартное выражение дополнительного изменения результативного показателя от взаимодействия факторов определяется выражением: 5 = 51Д: = АР- (ЛРМ + АРС1) = -364,0984533. (37)
Из сравнения значений и выражений видно, что ошибка стандартного метода больше ошибки нового метода более чем в 100 раз ОЗц/б^р! = 105,8998026).
3. Интегральный метод произвольного (нецелого) порядка
Математической основой стандартного интегрального метода являются операция интегрирования первого порядка, фундаментальная (основная) теорема математического анализа и классическая формула Ньютона - Лейбница:
£Г™Шх=т-Па). (38)
Интегральный метод является одним из наиболее общих методов факторного анализа, позволяющих разложить общий прирост результативного показателя по факторным приращениям. Приведем пример стандартных формул, описывающих связь между приращением функции и приращением факторных признаков. Для упрощения рассмотрим функцию двух переменных г = f (х, у). Формула интегрального метода, позволяющая вычислять влияние фактора х на результативный показатель, имеет вид:
^Х = Гг^Чх.у^х = Г/^ЪУЖ (39)
•'Хо •'Ха
Влияние фактора у вычисляется по формуле:
Агу= Г^1\х,у)ах= Г/^лЖ (40)
•'Уо •'Уо
Интегральный метод позволяет получать точные оценки факторных влияний и не предполагает разделение факторов на количественные и качественные.
В данном разделе исследуется возможность применения интегрального исчисления дробного (нецелого) порядка в факторном экономическом анализе. Как было доказано в работе [29], обратной операцией для производной Капуто является интегрирование Рима-на - Лиувилля порядка. Приведем определение этого интегрирования, данное в книге [3, с. 69].
Определение. Левосторонний и правосторонний интегралы Римана - Лиувилля порядка а > 0 определяются формулами:
ISSN 1993-047Х (Print) / ISSN 2410-0390 (Online) ...........................................................................
(«ъЛМ-^да. (х > а) (41)
сълм-^г^р*(х < Ь). (42)
где Г (а) - гамма-функция, а < х < Ь, а функция/(£) предполагается измеримой на интервале (а, Ь) и удовлетворяет условию < оо.
Интегрирование по Риману - Лиувиллю (39) и (40) является обобщением стандартной операции я-кратного интегрирования. Отметим, что интегралы Римана -Лиувилля (41) и (42) для порядка, равного единице а = 1, равны стандартному интегралу первого порядка:
(/¿+ям= Гтжаихх)= Саш (43)
•'а •'х
Таким образом, операция интегрирования, используемая в базовых формулах интегрального метода факторного анализа, является частным случаем операции интегрирования дробного порядка.
Поскольку интегральный метод базируется на формуле Ньютона - Лейбница и фундаментальной теореме математического анализа, то для обобщения этого метода нам необходимо иметь обобщение этой теоремы и формулы Ньютона - Лейбница на случай операторов дробного (нецелого) порядка. Фундаментальная (основная) теорема теории инте-гро-дифференцирования дробного (нецелого) порядка была сформулирована в статье [29] и монографиях [8, с. 247-248; 7, с. 263-267]. Некоторые дополнительные аспекты этой теоремы обсуждаются в статье [30]. Приведем обобщение формулы Ньютона - Лейбница на случай интегралов и производных нецелого порядка. Для левосторонних операторов выполняется следующая обобщенная формула Ньютона - Лейбница:
(/«+^+/)(Ь) = /(Ь) -/(а) - <0* (44)
где я-1 < а < я, а для правосторонних операторов формула имеет вид:
(М_/)(а) = /(а) -/(Ь) -Ш^Г^ - а)*.(45) Используя выражения:
(10ь-/)(а) = - С^Шх.^аЖМЮ =
'а
получаем, что формулы (42) и (43) для а = 1 дают стандартную формулу Ньютона - Лейбница (38).
Для применения в факторном анализе метода интегрирования нецелого порядка нам понадобятся формулы левостороннего и правостороннего интегралов Римана - Лиувилля порядка а > 0 для степенной функции [3, с. 71], имеющие вид: ГГВ +11
1*+(х - а)Р = г(р ^а + 1) (* " а)Р+а- (Р > 0), (47) ГГВ +11
lUb - *)р = т + а+1) (Ь - *)р+а. (Р > 0). (48)
Использование метода интегро-дифференцирова-ния нецелого порядка позволяет получить более точные результаты для влияния факторов по сравнению с методом дифференциального исчисления дробного порядка. Это обусловлено тем, что дополнительный прирост результативного показателя, возникающий из-за взаимодействия факторов, распределяется между ними в равной пропорции. В качестве примера приведем формулы для степенной функции от двух переменных, аналогичной рассмотренной выше функции Кобба - Дугласа. Для функции fix,у) = А-ха • ур, где а > 0 и в > 0, получаем выражения:
в А
Azx,a = А ■ yl ■ Лах + - ■ Лах ■ Ару =
= 2-Л«х-(У1+Уо)> (49)
AzyP =А-х£-Ару + --Аах- Ару =
= А-х£- (yf + у0Р). (50)
В этих формулах фактически учтено, что 8а,р: = А-Аах- Лру.
Для производственной функции С):=
з 1
:= 1,01 ■ L* • С* плановых и фактических значений L0 = 6,14, = 6,24 и С0 — 0,24, Сх — 0,54 факторов формулы (49) и (50) приводят к следующим значениям:
ДРь,а=а =A-CZ-AaL + ~ AaL ■ ЛрС =
= 4,0111645, (51)
J
Ja
= I fx(1)ix)dx,
(46)
АРс#=ъ = А • Ь% • АЪС + - ■ АаЬ • АрС =
= 70,4943135. (52)
В результате А Риа + АРС>Ъ = АР = 74,5054780.
Для а = в = 1 эти формулы принимают стандартный вид формул для метода дифференциального исчисле-
........................................................................... ISSN 1993-047Х (Print) / ISSN 2410-0390 (Online) ...........................................................................
ния с производными целого порядка для мультипликативной модели вида f (х, у) = А • х • у.
Помимо (47) и (48) существует набор формул интегрирования нецелого порядка, приведенных в табл. 9,1-9,2 в монографии [1, с. 140-142]. Для применения предлагаемого обобщения интегрального метода требуется использование этих табличных формул, которые позволяют разработать конечные рабочие формулы для наиболее распространенных видов факторных зависимостей и сделать этот метод более доступным.
Выводы
Предлагаемые методы интегро-дифференцирования нецелого порядка расширяют возможности детерминированного факторного экономического анализа. Метод дифференциального исчисления нецелого порядка может давать более точные результаты по сравнению со стандартным методом (методом дифференциального исчисления, использующим производные целого порядка), для широкого класса функций, включающего функции степенного типа. Кроме этого, методы инте-гро-дифференцирования нецелого порядка позволяют учитывать эффекты памяти и нелокальности в экономических процессах [19, 20, 23, 24, 32].
Предлагаемые дифференциальные и интегральные методы детерминированного факторного анализа могут быть использованы при анализе экономических или финансовых процессов, в которых показатели (эндогенные величины) являются степенными функциями
факторов (экзогенных величин) или описываемых уравнениями со степенной нелокальностью. Эти методы могут быть использованы для изучения процессов, описываемых степенными законами и функциями [26, 27], включая производственные функции Кобба - Дугласа [25], поскольку позволяют точнее описывать суммарное влияние факторов по сравнению со стандартным методом. Нелокальные экономические процессы часто связаны с эффектами запаздывания, которые описываются интегральными и интегро-дифференциальными уравнениями [31, с. 23-29]. Для описания нелокальности степенного типа применяются уравнения с интегралами и производными нецелого порядка, в силу содержания в этих операторах степенных ядер. Экономические процессы, являющиеся нелокальными в пространстве состояний (state space), описываются в статьях [15, 17, 18].
При применении этих методов в факторном анализе следует учитывать, что производные нецелого порядка обладают рядом нестандартных свойств, включающих нарушение стандартного правила дифференцирования произведения [1, с. 214-219] и нарушение стандартного правила дифференцирования сложной функции [4, с. 91-92]. С целью практического применения методов интегро-дифференцирования нецелого порядка для факторного экономического анализа необходимо разработать конечные рабочие формулы на основе табличных интегралов [1, с. 140-142] для различных видов факторных зависимостей, чтобы эти методы стали доступными для широкого круга аналитиков.
Список литературы
1. Самко С. Г., Килбас А. А., Марычев О. И. Интегралы и производные дробного порядка и некоторые приложения. Минск: Наука и Техника, 1987. 688 с.
2. Diethelm К. The Analysis of Fractional Differential Equations: An Application-Oriented Exposition Using Differential Operators of Caputo Type. Berlin, Heidelberg: Springer-Verlag, 2010. 247 p.
3. Kilbas A. A., Srivastava H. M., Trujillo J. J. Theory and Applications of Fractional Differential Equations. Amsterdam: Elsevier, 2006. 540 p.
4. Podlubny I. Fractional Differential Equations. San Diego: Academic Press, 1998. 340 p.
5. Samko S. G., Kilbas A. A., Marichev O. I. Fractional Integrals and Derivatives Theory and Applications. New York: Gordon and Breach, 1993. 1006 с.
6. Учайкин В. В. Метод дробных производных. Ульяновск: Артишок, 2008. 512 с.
7. Тарасов В. Е. Модели теоретической физики с интегро-дифференцированием дробного порядка. М.: Институт компьютерных исследований, 2011. 568 с.
8. Tarasov V. E. Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media. New York: Springer, 2011. 505 p.
9. Cartea A., Del-Castillo-Negrete D. Fractional diffusion models of option prices in markets with jumps // Physica A. 2007. Vol. 374. № 2. Рр. 749-763.
........................................................................... ISSN 1993-047Х (Print) / ISSN 2410-0390 (Online) ...........................................................................
10. Gorenflo R., Mainardi F., Scalas E., Raberto M. Fractional calculus and continuous-time finance III: the diffusion limit // In: M. Kohlmann, S. Tang, (Eds.) Mathematical Finance. Trends in Mathematics. Basel: Birkhauser, 2001. Рр. 171-180.
11. Kerss A., Leonenko N., Sikorskii A. Fractional Skellam processes with applications to finance // Fractional Calculus and Applied Analysis. 2014. Vol. 17. № 2. Рр. 532-551.
12. Laskin N. Fractional market dynamics // Physica A. 2000. Vol. 287. № 3. Рр. 482-492.
13. Mainardi F., Raberto M., Gorenflo R., Scalas E. Fractional calculus and continuous-time finance II: The waiting-time distribution // Physica A. 2000. Vol. 287. № 3-4. Рр. 468-481.
14. Scalas E., Gorenflo R., Mainardi F. Fractional calculus and continuous-time finance // Physica A. 2000. Vol. 284. № 1-4. Рр. 376-384.
15. Skovranek T., Podlubny I., Petras I. Modeling of the national economies in state-space: A fractional calculus approach // Economic Modelling. 2012. Vol. 29. № 4. Рр. 1322-1327.
16. Tenreiro Machado J., Duarte F. B., Duarte G. M. Fractional dynamics in financial indeces // International Journal of Bifurcation and Chaos. 2012. Vol. 22. № 10. Article ID 1250249. 12 p.
17. Tenreiro Machado J. A., Mata M. E. Pseudo phase plane and fractional calculus modeling of western global economic downturn // Communications in Nonlinear Science and Numerical Simulation. 2015. Vol. 22. № 1-3. Рр. 396-406.
18. Tenreiro Machado J. A., Mata M. E., Lopes A. M. Fractional state space analysis of economic systems // Entropy. 2015. Vol. 17. № 8. Рр. 5402-5421.
19. Тарасова В. В., Тарасов В. Е. Критерии эредитарности экономического процесса и эффект памяти // Молодой ученый. 2016. № 14 (118). С. 396-399.
20. Тарасова В. В., Тарасов В. Е. Предельная полезность для экономических процессов с памятью // Альманах современной науки и образования. 2016. № 7 (109). C. 108-113.
21. Тарасова В. В., Тарасов В. Е. Ценовая эластичность спроса с памятью // Экономика, тециология и право. 2016. № 4-1. С. 98-106.
22. Тарасова В. В., Тарасов В. Е. Предельные величины нецелого порядка в экономическом анализе // Азимут Научных Исследований: Экономика и Управление. 2016. Том 5. № 3 (16). С. 197-201.
23. Тарасова В. В., Тарасов В. Е. Эластичность внебиржевого кассового оборота валютного рынка РФ // Актуальные проблемы гуманитарных и естественных наук. 2016. № 07-1 (90). С. 207-215.
24. Tarasova V. V., Tarasov V. E. Elasticity for economic processes with memory: Fractional differential calculus approach // Fractional Differential Calculus. 2016. Vol. 6. № 2. Рр. 219-232.
25. Cobb C. W., Douglas Р. H. A theory of production // American Economic Review. 1928. Vol. 18 (Supplement). Рр. 139-165.
26. Gabaix X. Power laws in economics and finance // Annual Review of Economics. 2009. Vol. 1. № 1. Pp. 255-293.
27. Gabaix X. Power laws in economics: An introduction // Journal of Economic Perspectives. 2016. Vol. 30. №2 1. Рр. 185-206.
28. Odibat Z. M., Shawagfeh N. T. Generalized Taylor's formula // Applied Mathematics and Computation. 2007. Vol. 186. № 1. Рр. 286-293.
29. Tarasov V. E. Fractional vector calculus and fractional Maxwell's equations // Annals of Physics. 2008. Vol. 323. № 11. Рр. 2756-2778.
30. Grigoletto E. C., De Oliveira E. C. Fractional versions of the fundamental theorem of calculus // Applied Mathematics. 2013. Vol. 4. Рр. 23-33.
31. Allen R. G. D. Mathematical Economics. Second edition. London: Macmillan, 1960. 812 p.
32. Тарасова В. В., Тарасов В. Е. Экономический показатель, обобщающий среднюю и предельную величины // Экономика и предпринимательство. 2016. № 11-1 (76-1). С. 817-823.
Дата поступления 06.09.2016 Дата принятия в печать 21.10.2016 Дата онлайн-размещения 29.12.2016
© Тарасова В. В., Тарасов В. Е., 2016
Информация об авторах
Тарасова Валентина Васильевна, магистрант, Высшая школа бизнеса, Московский государственный университет им. М. В. Ломоносова
Адрес: 119991, Москва, ГСП-1, Ленинские горы, д. 1, стр. 52А, ВШБ МГУ, тел.: + 7 (495) 939-59-89
E-mail: v.v.tarasova@mail.ru
ORCID: http://orcid.org/0000-0003-4410-8470
Researcher ID: http://www.researcherid.com/rid/J-4141-2016
ISSN 1993-047Х (Print) / ISSN 2410-0390 (Online)
Контактное лицо:
Тарасов Василий Евгеньевич, доктор физико-математических наук, ведущий научный сотрудник, Научно-исследовательский институт ядерной физики им. Д. В. Скобельцына, Московский государственный университет им. М. В. Ломоносова Адрес: 119991, Москва, Ленинские горы, д. 1, стр. 2, НИИЯФ МГУ, тел.: + 7 (495) 939-59-89 E-mail: v.e.tarasov@bk.ru, tarasov@theory.sinp.msu.ru ORCID: http://orcid.org/0000-0002-4718-6274 Researcher ID: http://www.researcherid.com/rid/D-6851-2012
DETERMINISTIC FACTOR ANALYSIS: METHODS OF INTEGRO-DIFFERENTIATION OF NON-INTEGRAL ORDER
Objective: to summarize the methods of deterministic factor economic analysis, namely the differential calculus and the integral method. Methods: mathematical methods for integro-differentiation of non-integral order, the theory of derivatives and integrals of fractional (non-integral) order.
Results: the basic concepts are formulated and the new methods are developed that take into account the memory and non-locality effects in the quantitative description of the influence of individual factors on the change in the effective economic indicator. Two methods are proposed for integro-differentiation of non-integral order for the deterministic factor analysis of economic processes with memory and non-locality. It is shown that the method of integro-differentiation of non-integral order can give more accurate results compared with standard methods (method of differentiation using the first order derivatives and the integral method using the integration of the first order) for a wide class of functions describing effective economic indicators.
Scientific novelty: the new methods of deterministic factor analysis are proposed: the method of differential calculus of non-integral order and the integral method of non-integral order.
Practical significance: the basic concepts and formulas of the article can be used in scientific and analytical activity for factor analysis of economic processes. The proposed method for integro-differentiation of non-integral order extends the capabilities of the determined factorial economic analysis. The new quantitative method of deterministic factor analysis may become the beginning of quantitative studies of economic agents behavior with memory, hereditarity and spatial non-locality. The proposed methods of deterministic factor analysis can be used in the study of economic processes which follow the exponential law, in which the indicators (endogenous variables) are power functions of the factors (exogenous variables), including the processes described by the Cobb - Douglas production function, since these methods allow to more accurately describe the total influence of the factors in comparison with the standard method. The proposed methods can be used in the study of economic processes described by equations with a power-law non-locality in factor space and in state space.
Keywords: Economic and mathematical methods of Economics; Factor analysis; Method of differential calculus; Integral method; Processes with memory; Hereditarity; Derivative of non-integral order; Integration of non-integral order
1. Samko, S. G., Kilbas, A. A., Marychev, O. I. Integrals and derivatives of non-integral order and some applications, Minsk: Nauka i Tekh-nika, 1987, 688 p. (in Russ.).
2. Diethelm, K. The Analysis of Fractional Differential Equations: An Application-OrientedExposition Using Differential Operators ofCaputo Type, Berlin, Heidelberg: Springer-Verlag, 2010, 247 p.
3. Kilbas, A. A., Srivastava, H. M., Trujillo, J. J. Theory andApplications of Fractional Differential Equations, Amsterdam: Elsevier, 2006, 540 p.
4. Podlubny, I. Fractional Differential Equations, San Diego: Academic Press, 1998, 340 p.
5. Samko, S. G., Kilbas, A. A., Marichev, O. I. Fractional Integrals and Derivatives Theory andApplications, New York: Gordon and Breach, 1993, 1006 p.
6. Uchaikin, V. V. Technique of derivatives of non-integral order, Ulyanovsk: Artishok, 2008, 512 p. (in Russ.).
7. Tarasov, V. E. Models of theoretical physics with integro-differentiation of non-integral order, Moscow: Institut komp'yuternykh issledo-vanii, 2011, 568 p. (in Russ.).
8. Tarasov, V. E. Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields andMedia, New York: Springer, 2011, 505 p.
9. Cartea, A., Del-Castillo-Negrete, D. Fractional diffusion models of option prices in markets with jumps, PhysicaA, 2007, vol. 374, No. 2, pp. 749-763.
10. Gorenflo, R., Mainardi, F., Scalas, E., Raberto, M. Fractional calculus and continuous-time finance III: the diffusion limit. In: M. Kohlmann, S. Tang, (Eds.) Mathematical Finance. Trends in Mathematics, Basel: Birkhauser, 2001, pp. 171-180.
11. Kerss, A., Leonenko, N., Sikorskii, A. Fractional Skellam processes with applications to finance, Fractional Calculus andApplied Analysis, 2014, vol. 17, No. 2, pp. 532-551.
V. V. TARASOVA1 V. E. TARASOV1
1 Moscow State University named after M.V. Lomonosov, Moscow, Russia
References
........................................................................... ISSN 1993-047Х (Print) / ISSN 2410-0390 (Online) ...........................................................................
12. Laskin, N. Fractional market dynamics, PhysicaA, 2000, vol. 287, No. 3, pp. 482-492.
13. Mainardi, F., Raberto, M., Gorenflo, R., Scalas, E. Fractional calculus and continuous-time finance II: The waiting-time distribution, Physica A, 2000, vol. 287, No. 3-4, pp. 468-481.
14. Scalas, E., Gorenflo, R., Mainardi, F. Fractional calculus and continuous-time finance, Physica A, 2000, vol. 284, No. 1-4, pp. 376-384.
15. Skovranek, T., Podlubny, I., Petras, I. Modeling of the national economies in state-space: A fractional calculus approach, Economic Modelling, 2012, vol. 29, No. 4, pp. 1322-1327.
16. Tenreiro Machado, J., Duarte, F. B., Duarte, G. M. Fractional dynamics in financial indeces, International Journal of Bifurcation and Chaos, 2012, vol. 22, No. 10, Article ID 1250249, 12 p.
17. Tenreiro Machado, J. A., Mata, M. E. Pseudo phase plane and fractional calculus modeling of western global economic downturn, Communications in Nonlinear Science and Numerical Simulation, 2015, vol. 22, No. 1-3, pp. 396-406.
18. Tenreiro Machado, J. A., Mata, M. E., Lopes, A. M. Fractional state space analysis of economic systems, Entropy, 2015, vol. 17, No. 8, pp. 5402-5421.
19. Tarasova, V. V., Tarasov, V. E. Criteria of hereditiarity of economic process and memory effect, Molodoi uchenyi, 2016, No. 14 (118), pp. 396-399 (in Russ.).
20. Tarasova, V. V., Tarasov, V. E. Marginal utility for economic processes with memory, Al'manakh sovremennoi nauki i obrazovaniya, 2016, No. 7 (109), pp. 108-113 (in Russ.).
21. Tarasova, V. V., Tarasov, V. E. Price elasticity of demand with memory, Ekonomika, cotsiologiyaipravo, 2016, No. 4-1, pp. 98-106 (in Russ.).
22. Tarasova, V. V., Tarasov, V. E. Marginal values of of non-integral order in economic analysis, Ekonomika i Upravlenie, 2016, vol. 5, No. 3 (16), pp. 197-201 (in Russ.).
23. Tarasova, V. V., Tarasov, V. E. Elasticity of over-the-counter cash turnover of the Russian currency market, Aktual'nye problemy gumani-tarnykh i estestvennykh nauk, 2016, No. 07-1 (90), pp. 207-215 (in Russ.).
24. Tarasova, V. V., Tarasov, V. E. Elasticity for economic processes with memory: Fractional differential calculus approach, Fractional Differential Calculus, 2016, vol. 6, No. 2, pp. 219-232.
25. Cobb, C. W., Douglas, P. H. A theory of production, American Economic Review, 1928, vol. 18 (Supplement), pp. 139-165.
26. Gabaix, X. Power laws in economics and finance, Annual Review of Economics, 2009, vol. 1, No. 1, pp. 255-293.
27. Gabaix, X. Power laws in economics: An introduction, Journal of Economic Perspectives, 2016, vol. 30, No. 1, pp. 185-206.
28. Odibat, Z. M., Shawagfeh, N. T. Generalized Taylor's formula, AppliedMathematics and Computation, 2007, vol. 186, No. 1, pp. 286-293.
29. Tarasov, V. E. Fractional vector calculus and fractional Maxwell's equations, Annals of Physics, 2008, vol. 323, No. 11, pp. 2756-2778.
30. Grigoletto, E. C., De Oliveira, E. C. Fractional versions of the fundamental theorem of calculus, AppliedMathematics, 2013, vol. 4, pp. 23-33.
31. Allen, R. G. D. Mathematical Economics. 2nd edition, London: Macmillan, 1960, 812 p.
32. Tarasova, V. V., Tarasov, V. E. Economic indicator summarizing the average and the marginal values, Ekonomika i predprinimatel'stvo, 2016, No. 11-1 (76-1), pp. 817-823 (in Russ.).
Valentina V. Tarasova, Master student, Higher School of Business, Moscow State University named after M.V. Lomonosov
Address: 1 Leninskiye gory, building 52A, 119991 Moscow, tel.: + 7 (495) 939-59-89
E-mail: v.v.tarasova@mail.ru
ORCID: http://orcid.org/0000-0003-4410-8470
Researcher ID: http://www.researcherid.com/rid/J-4141-2016
Contact:
Vasiliy E. Tarasov, Doctor of Physics and Mathematics, Leading Researcher, Scientific-Research Institute of Nuclear Physics named after D.V. Skobel'tsyn, Moscow State University named after M.V Lomonosov
Address: 1 Leninskiye gory, building 2, 119991 Moscow, tel.: + 7 (495) 939-59-89 E-mail: v.e.tarasov@bk.ru, tarasov@theory.sinp.msu.ru ORCID: http://orcid.org/0000-0002-4718-6274 Researcher ID: http://www.researcherid.com/rid/D-6851-2012
For citation: Tarasova V. V., Tarasov V. E. Deterministic factor analysis: methods of integro-differentiation of non-integral order, Actual Problems of Economics and Law, 2016, vol. 10, No. 4, pp. 77-87 (in Russ.). DOI: 10.21202/1993-047X.10.2016.4.77-87
Received 06.09.2016 Accepted 21.10.2016 Available online 29.12.2016
© Tarasova V. V., Tarasov V. E., 20i6
Information about the authors
V
y