СОЦИАЛЬНАЯ ФИЗИКА

К разделам сайта:          Домой        Примеры       Теория индексов      Банки      Физика

        матричная квантовая экономика    

 



 
 

 

 

ОЦЕНКИ И СРАВНЕНИЕ МАССИВОВ В АНАЛИЗЕ ДАННЫХ.

В данной работе рассмотрен механизм построения количественной и качественной оценки больших многоэшелонированных структур, при описании которых применяются многомерные массивы, или многотомные книги. Такие массивы можно получить, например, бесконечной дезагрегацией (с фиксацией уровней) целеустремлённой системы (градуированный пучёк [1, 2]), применяемой при описании и оценке действия сложной иерархически структурированной организации [4], нейронных сетей [5], описываемых в категориях [6], методы решение которых с применением римановой метрики даются в работе [7], а с применением метрики симплектических пространств - в работе [8].

Терминологию и методы решения задач будем заимствовать из MATLAB [9]. 

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

Прямые методы анализа опираются на линейные модели, которые дают возможность проводить анализ, не отступая от смысловой интерпретации задачи. Однако решение на основе линейных моделей в большинстве случаев не даёт необходимой точности результатов. Возникает необходимость в построении нелинейных моделей, первым приближением которых, например, является применение моделей математической статистики. Если линейные задачи решаются в пространстве суммируемых функций (L1), то решение вторых задач строится в пространстве функций, суммируемых с квадратом (L2).  Моделирование задач в пространство L2 даёт возможность не только использовать хорошо разработанные математические методы (например, такие, как методы гильбертова пространства), но и позволяет воспользоваться геометрическими методами интерпретации полученных решений.

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

Как правило, наблюдения поведения объекта отражается в числовой информации, в которой состояния отражены в системе факторов его описания. При множестве наблюдений получаем прямоугольную числовую таблицу – двумерный массив действительных натуральных чисел, матрицу. Предполагая, что процесс наблюдения носит регулярный характер (к примеру, бухгалтерские отчёты деятельности некоторой фирмы), отчёты накапливаются в папку, страницы которой нумеруются. Получаем книгу, которая уже составляет трёхмерный массив. Постепенно полка заполняется папками – полку с папками можно полагать четырёхмерным массивом. Со временем такими книгами заполняются все полки данного стеллажа. Формируется массив пятой размерности и т.д. Приходим к понятию многомерных числовых массивов, в которых накапливается громадная информация. Естественно, что справиться с обработкой подобной информации можно, если её перевести на машинные носители. Для обработки такой информации, разработаны соответствующие математические методы.

Пусть двумерный числовой массив  имеет m строк и n столбцов. Его будем записывать в виде матрицы Am´n, где m = |M|, n = |N|. Такие матрицы-сводки могут зависеть от временного фактора  t = T. Для трёхмерного массива введём обозначения

A =A(T, M, N).

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

a = A(2, 3, 4).

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

Матрица – как двумерный числовой массив. Для построения оценок массивов рассмотрим широко распространённую  табличную форму различных производственных отчётных материалов, сформулировав её в виде некоторой затратной модели (см. табл. 1).

 

Таблица 1.   Характеристики производственного процесса.

 

Затраты

на производство

продукции

 

 

Расходы ресурсов

 

Стоимость

 

         1               2           …          j           …           n

 

Производства

 

1

2

i

m

 

       A11                 A12         …        A1 j         …          A1n

      A21                 A22         …        A2j          …          A2n

     . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

A21                 A22         …        A2j          …          A2n

. . . . . . . . . . . . . . . . . . . . . . . . . . .. . . . . . . . . . . . . . .

    A21                 A22           …        A2j           …          A2n

 

q1

q2

qi

qm

 

 

Цена

 

p1                 p2           …        pj           …          pn

          Q

   P

       

 

Это модель некоторой производственной системы, конечным результатом деятельности которой является стоимость продукции.

Введём матрицу-столбец цен ресурсов

p = (pj | j Î N),

и матрицу-столбец объектов

q = (qi| i Î M).

Здесь введены множества N = (1, 2, …, n) и M = (1, 2, …, m).

Если ввести матрицу затрат

A = ((Aij)),

где Aij – затраты j-го ингредиента на выпуск изделий i-го вида, то данную модель производственного процесса можно записать в матричной форме

q = Ap.                                                                                   (1)     

Вектор q выступает в роли вектора оценок строк матрицы A в пространстве Xm. Его можно нормировать в пространстве l1(Xm) суммируемых последовательностей, либо в пространстве l2(Xm) последовательностей, суммируемых с квадратом.  Но, вектор p также можно рассматривать в качестве вектора оценок столбцов этой же матрицы. Имеем

                             q*q = p*(A*A)p,                                                                              (2)

где A* - эрмитово сопряжённая матрица матрицы A.

Из (1) находим

p = A+q.

Здесь A+ - псевдообратная матрица матрицы A. Умножая обе части данного равенства на произведение A*A и учитывая свойство псевдообратной матрицы [10, стр. 35]

A*AA+ = A*,

приходим к соотношению

(A*A)p = A*q,

подставляя которое в (2), будем иметь

q*q = p*A*q.                                                                                     (3)

Если полагать, что векторы p и q ортонормированы, то из равенства (3) получаем соотношения

p*A*q = q*q = p*p = 1,                                                                            (4)

из которого следует аналогичная (1) обратная зависимость

   p = A*q                                                                                       (5)

и которые вместе с нормой данных векторов в соответствующих пространствах

||q|| = (q*q)1/2,     ||p|| = (p*p)1/2

дают возможность вычислить норму оператора A

||A|| = (lim||p||->1(p*A*Ap/p*p))1/2.

Из соотношения

l = lim||q||->1(q*AA*q/q*q)

находим, что

l = ||A||2.

Отсюда следует, что векторы оценок столбцов p и строк q прямоугольной квадратной матрицы A с неотрицательными членами являются собственными векторами операторов A*A и AA* - симметричные квадратные матрицы с неотрицательными членами, отвечающие одному и тому же наибольшему собственному значению этих матриц, числу l, т.е. имеют место равенства

(A*A)p = lp,     (AA*)q = lq.                                                                    (6)

Бинарность (5) оценок одного и того же оператора  A вытекает из его материальной сущности - экономической (в данном случае, см. табл. 1), которая наделяет данный оператор двумя (и только двумя!) фундаментальными одинаково существенными и взаимно дополнительными структурными свойствами материи в любой её форме самоорганизации – количества и качества, и, таким образом, выражает симметрию этих свойств в виде равенств (1) и (5) [11, стр. 44].

Механизм построения оценок массива, заданного в виде прямоугольной матрицы, и сравнение двух массивов продемонстрируем на действительных матрицах A1 и A2, взятых из публичной печати,  и представленные таблицами 1 и 3 (см. Приложение 1).

Прежде всего, в данных массивах определим совокупность факторов  описания в их относительно независимом виде. Однако следует отметить, что поскольку система замкнута, то все факторы связаны друг с другом и изменение одного влияет на изменение другого. Но, в данном представлении в фактор "активы" включены "кредиты", а в фактор "депозиты" включены "вклады населения". Поэтому, пользуясь языком MATLAB, из столбцов исходной матрицы  Ai (i = 1, 2) сформируем новый массив A:

A = [Ai (:,1)  Ai (:, 3)  Ai (:, 2) – Ai (:, 3)  Ai (:, 5)  Ai (:, 4) – Ai (:, 5)  Ai (:, 6)].

Здесь факторы выступают в последовательности: "Капитал банка", "Кредиты",  "Другие активы", "Вклады населения", "Другие депозиты", "Прибыль". Получаем новую таблицу 1 , с другим описанием факторов и с новым массивом их значений в системе. Аналогично получаем табл.4 (см. приложения 1) с новым массивом данных.

Количественные оценки матрицы. Любой прямоугольный числовой массив A характеризуется двумя собственными  последовательностями индикаторов – последовательностью индикаторов оценок строк v и последовательностью индикаторов оценок  столбцов u, которые связаны с абсолютными величинами p и q зависимостями

u = q / ||A||,   v = p / ||A||,

а поскольку уравнения (6) однородные, то v и u будут собственными векторами операторов A*A и AA*, соответственно, отвечающие их наибольшим собственным значениям. После решения уравнений (6) искомые абсолютные оценки находим по формулам

p = ||A||v,    q = ||A||u.                                                                             (7)

  Таблица 1. Итоги деятельности банков Ростовской области на 01.07.04 (тыс. руб.)

Наименование

 банка    

Факторы

    Капитал   Кредиты   Др.активы     Вклады   Др.депозиты   Прибыль                                                         

 1. Волгодонский  гкб

 2. Дон-Тексбанк

 3. Донактивбанк

 4. Донбанк

 5. Донинвест

 6. Донкомбанк

 7. Донской народный банк

 8. Донхлеббанк

 9. Земельный КБ

10. Капиталбанк

11. ИБРР-банк 

12. Кредит Экспресс

13. МеТраКомБанк

14. Морозовсккомагробанк

15. Рост.  универсальный банк

16. РостПромстройбанк

17. Сельмашбанк

18. Стелла-Банк

19. Таганрогбанк

20. Центр-инвест

21. Южный регион

22. Южный регион. банк 

23. Южный Торговый Банк  

     8340     12520    21867     7290       1000           634

   35274   188500  183456 198064     14700           913

  105662   186120  175495    28365   10000         2473

    55125   208985  366410  166129   20500         7036

  237856  683559   473929  256491   72498         3111

    64042  199017     96677    72223   18620         6829  

   90780   611069   307656  638906   74324       13839

 103735   709324   360219  511942  100342        7772

     9015     31043     13738    14127      5557          531

   26415     52053      32691   35907           0             76

   14179     30345      25504   17379      3277           677

 161368     11849    151079           1            0           154

 282247 1391678  1405994  401210  551851      46256

   15959     27820      45381        161            0          444

   26888     71500      35225    33941     24900         455

 200677   982897  1040477   703540   160488    13387

   35001     49375    458397     19860       6200      3653

   84558   441313    130529   163146     52793      4570

   51661   107172     52176      15401             0      2825

 685032 4691216 2769207  2414785 1288418    74988

 109262   235970   328530    235625     27411      3471

 156950   304918  240781     16551    52500        1425

 330653   987970  650148   495952  228222      30475

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Известно, что если существует хотя бы одна последовательность, например, оценок столбцов матрицы, то можно построить, опираясь на множество выпуклых функций, множество эквивалентных последовательностей оценок столбцов данного массива. Их эк��ивалентность можно обнаружить проекцией этих шка�� на одну из них. В таблицах 2 и 3 представлены оценки столбцов и строк  массива, приведённого в табл.1

Таблица 2. Оценки объектов.

 

Объекты

(Банки)

Оценка в пространствеl1

Оценка в пространствеl2

Оценка

(u)

Ранг

Оценка

(u)

Ранг

 

 

 1. Волгодонский  гкб

 2. Дон-Тексбанк

 3. Донактивбанк

 4. Донбанк

 5. Донинвест

 6. Донкомбанк

 7. Донской народный банк

 8. Донхлеббанк

 9. Земельный КБ

10. Капиталбанк

11. ИБРР-банк 

12. Кредит Экспресс

13. МеТраКомБанк

14. Морозовсккомагробанк

15. Рост.  универсальный банк

16. РостПромстройбанк

17. Сельмашбанк

18. Стелла-Банк

19. Таганрогбанк

20. Центр-инвест

21. Южный регион

22. Южный региональный банк

23. Южный Торговый Банк   

 

0.001409

0.018372

0.014682

0.023986

0.0518

0.013785

0.051447

0.054783

0.0021999

0.0042123

0.0025894

0.0060085

0.1186

0.0026327

0.0053865

0.091889

0.016014

0.027951

0.0069057

0.35881

0.026125

 0.022349

0.078062

 

23

12

14

10

7

15

6

5

22

19

21

18

2

20

17

3

13

8

16

1

9

11

4

 

  0.0034237 

  0.044642  

  0.035676  

  0.058283  

   0.12587   

  0.033497  

   0.12501   

   0.13312   

 0.0053455 

   0.010235 

  0.0062919

     0.0146   

    0.28819  

  0.0063972

   0.013089 

    0.22328  

   0.038913 

   0.067916 

    0.01678  

    0.87186  

   0.063481 

   0.054306 

    0.18968  

 

23

12

14

10

6

15

7

5

22

19

21

17

2

20

18

3

13

8

16

1

9

11

4

 

Таблица 3. Оценки факторов.

 

 

Факторы

 

 

 

 

Оценка в пространстве l1

 

Оценка в пространстве l2

Оценка

(v)

 

Ранг

Оценка

(v)

 

Ранг

 

 

1.     Капитал

2.     Кредиты

3.     Др. активы

4.     Вклады населения

5.     Др. депозиты

6.     Прибыль

 

0.063848

0.38197

0.24893

 

0.19646

0.10209

0.0067018

 

5

1

2

 

3

4

6

  

0.063848

0.38197

0.24893

 

0.19646

0.10209

0.0067018

 

5

1

2

 

3

4

6

 

в пространствах суммируемых последовательностей и суммируемых последовательностей с квадратом. Из табл.3 и 4 видим, что хотя численные значения на данных шкалах различны, их проекции на ранговую шкалу совпадают.

 

Таблица 4.  Оценка факторов в пространстве l2.

 

 

Факторы

 

Системы

A1

A2

 

Оценка

Вероят-

ность

 

 

Ранг

 

Оценка

Вероят-

ность

 

Ранг

1. Капитал

2. Кредиты

3. Др.

    активы

4. Вклады населения

5. Др. депозиты

6. Прибыль

 0.12497

 0.74766

 

 0.48724

 

  0.38454

 0.19984

0.013118

0.015618

0.559

 

0.2374

   

0.14787

0.039936

0.0001720

 

 

 6

0.15316 

0.78759 

 

0.41157 

 

0.3583  

0.23905 

0.036385

0.023458

0.6203

 

0.16939

     

0.12838

0.057145

0.001324

5

1

 

2

 

3

4

6

 

X1 = Y1= 7.0302e+006                   Y2 = X2 = 1.3798e+007

Таблица 5.  Оценка объектов системы в пространстве l2.

 

Объекты

(банки)

Системы

A1

A2

Оценка

Ранг

Оценка

Ранг

 1. Волгодонский  гкб

 2. Дон-Тексбанк

 3. Донактивбанк

 4. Донбанк

 5. Донинвест

 6. Донкомбанк

 7. Донской народный банк

 8. Донхлеббанк

 9. Земельный КБ

10. Капиталбанк

11. ИБРР-банк 

12. Кредит Экспресс

13. МеТраКомБанк

14. Морозовсккомагробанк

15. Рост.  универсальный банк

16. РостПромстройбанк

17. Сельмашбанк

18. Стелла-Банк

19. Таганрогбанк

20. Центр-инвест

21. Южный регион

22. Южный региональный банк 

23. Южный Торговый Банк   

  0.0034237  

  0.044642   

  0.035676   

  0.058283   

   0.12587    

  0.033497   

   0.12501    

   0.13312    

 0.0053455  

   0.010235  

  0.0062919 

     0.0146    

    0.28819   

  0.0063972 

   0.013089  

    0.22328   

   0.038913  

   0.067916  

    0.01678   

    0.87186   

   0.063481  

   0.054306  

    0.18968

23

12

14

10

6

15

7

5

22

19

21

17

2

20

18

3

13

8

16

1

9

11

4

0.0024567   

0.027049     

0.032271     

0.044261     

0.089214       

0.038868     

0.11553        

0.097855       

0.0022528   

0.0067426  

0.030307  

0.013799     

0.13064      

0.0018411 

0.013864   

0.12298      

0.019298   

0.043535    

0.0094427   

0.95646      

0.039744    

0.029841   

0.10335

21

15

12

8

7

11

4

6

22

20

13

18

2

23

17

3

16

9

19

1

10

14

5

 

Эту же систему рассмотрим в другой момент времени, состояние в котором представлено массивом данных A2 (табл. 3 приложения 1).

Переходя к новым факторам, получаем массив B2 (табл.4 прилож.1).  Оценки факторов массивов B1 и B2 в пространстве l2 приведены в таблице 4, где кроме оценок факторов помещены вероятности обнаружения этих факторов в общей деятельности системы.

Из табл. 4 находим, что хотя ранговые оценки факторов остались без изменения, в их относительных оценках произошёл сдвиг.   Существенные изменения произошли в активах системы и основном капитале. Заметные изменения произошли в прибыли. Но, можно отметить, что основные финансовые потоки системы проходят через активы (80%) и, хотя в активах произошли существенные изменения, вероятность  их присутствия в системе практически не изменилась, снизилась лишь на 0,84%:

(B1):  y2 + y3 = 0.7964,   (B2):   y2 + y3 = 0.78969.

Обратившись к табл. 5 заключаем, что в системе объектов произошли ранговые подвижки. Своё присутствие в системе укрепил банк "Центр-инвест", который более чем на 90% определяет состояние финансовой системы Ростовских банков. Чтобы оценить изменения в данных, нужно иметь, по крайней мере, два массива данных, из которых и можно выбрать измерительный эталон. В данном случае измерительным эталоном может служить как массив B1, так и массив B2, а агрегированные оценки структурных сдвигов можно строить с помощью индексов типа Ласпейреса

     Ly = y2*y1 / y1*y1,         Lx = x2*x1 / x1*x1,                                                    (7)   

или индексов типа Пааше

Yy = y2*y / y2*y1,          Yx = x2*x2 / x2*x1.                                                      (8) 

Учитывая равенства (3), получаем

Ly = y2*y1,              Lx = x2*x1.                                                                  (9)

Но,

s(x1) = s(x2) = s(y1) = s(y2) = 1

и, следовательно,

px = Lx2,             py = Ly2

суть вероятности  обнаружения сдвигов в оценках факторов и объектов за данный период наблюдения.

Из выражений находим

Yy = 1 / Ly ,            Yx = 1 / Lx.

Принимая массив B1 за эталон, и используя оценки (9), будем иметь

Ly = 0.99455,      Lx = 0.97223,

т.е.

px = 0.94523,        py = 0.98913.

Видим, что при переходе системы из одного состояния в другое произошли незначительные (1,1%) качественные изменения. Эти изменения не отразились на шкале рангов (табл. 4). В то же время, произошли значительные количественные изменения, вероятность обнаружения которых равна 5,5%. В табл. 5 это отразилось в подвижки рангов.

 Чтобы отследить количественные изменения, разобьём диапазон вариации оценки объектов системы на количественной шкале, например, на четыре интервала. При этом объекты по количественной вероятностной оценке разобьём на четыре группы, как это показано в табл.6 для системы в состоянии А1 (определённым в табл.2).

В табл. 6 система представлена двухуровневой иерархической структурой. На первом уровне она дана разбиением на группы по количественному признаку (мощности): малые, средние и крупные. В четвёртую группу выделяется объект, по мощности действия несоизмеримый с остальными, который определяет по внутренним оценкам системы 76% действия всей системы в целом.

Таблица 6. Группировка объектов системы в состоянии А1.

№ группы

п./п.

Объект

Оценка объекта

в группе

(u)

Вероятностная оценка

группы (p)

п./п.

Объекты группы

u1

p1

 

 

 

 

1.

1.

Волгодонский ГКБ

0,013146

 

 

 

 

0,0008916

2.

Земельный КБ

0,032048

3.

Капиталбанк

0,117500

4.

ИБРР-банк

0,044400

5.

Кредит-Экспресс

0,239070

6.

Морозовсккомагробанк

0,045900

7.

Рост.универсальный банк

0,192140

8.

Таганрогбанк

0,315800

 

 

 

 

2.

п./п.

Объекты группы

u2

p2

1.

Дон-Тексбанк

0,095398

 

 

 

 

0,02089

2.

Донактивбанк

0,060925

3.

Донбанк

0,162610

4.

Донкомбанк

0,053710

5.

Сельмашбанк

0,072483

6.

Стелла-Банк

0,220800

7.

Южный регион

0,192900

8.

Южный региональный банк

0,141170

 

 

 

3.

п./п.

Объекты группы

u3

p3

1.

Донинвест

0,726480

 

 

 

0,21807

2.

Донской народный банк

0,716610

3.

Донхлеббанк

0,081256

4.

МеТраКомБанк

0,380840

5.

РостПромстройбанк

0,228610

6.

Южный торговый банк

0,164980

 

4.

п./п.

Объекты группы

u4

p4

1.

Центр-инвест

1

0,76014

 

Ниже, на рис.1, даётся структурная схема количественного расслоения системы в метрическом пространстве l2, из которой следует, что

u = piui

и

||u||2 = pi||ui||2,

где pi – вероятностная оценка присутствия соответствующей группы в системе.

Рис.1. Расслоение системы по количественным признакам в пространстве l2.

Оценки качества. Множество объектов (при этом, рассматривая множество состояний одного объекта как различные объекты) в n-мерном факторном пространстве будем рассматривать как аффинное пространство E. Если в этом пространстве фиксировать некоторый объект a в качестве начала координат, то оно превратится в n-мерное векторное пространство R, элемент x которого определяется как разность объектов x и a  аффинного пространства E

x = x  a,                                                                                     (10)

что определяет внешний закон композиции

x = a + x.

Наделяя окрестность фиксированной точки евклидовой структурой определением на элементах x, y векторного пространства R положительно-определённой билинейной симметрической формой m(x, y), можно ввести расстояние между точками x, y аффинного пространства

r(x, y) = ||y  x|| = ||x|| = m1/2(x, x),                                                              (11)

связав, метрическим соотношением (10) аффинное и векторное пространства и поставив всякому объекту, достаточно близкому к фиксированному объекту аффинного пространства, во взаимно однозначное соответствие элемент векторного пространства.

Аффинный репер {a, x1, x2, …,xn} в аффинном пространстве E порождает базис {x1, x2, …, xn} присоединённого векторного пространства Rтакой, что координаты (a1, a2, …, an) точки x множества E в данном репере, где введены обозначения xi = xi  a, i  =1, 2, …, n, одновременно являются координатами вектора x = x  a в данном базисе

x = aixi.                                                                                   (12)

Если обозначить Ea векторное пространство, образованное элементами аффинного пространства E из окрестности произвольной точки a множества E, то отображение fa: E -> R порождает естественный изоморфизм между векторным пространством R и векторизацией области Ea аффинного пространства E в точке a.

Установленный изоморфизм позволяет каждому объекту x аффинного пространства области Ea поставить во взаимно-однозначное соответствие его образ в векторном пространстве – вектор x, а паре объектов x, y из Ea – ковариацию их образов m(x, y). Например, путём сопоставления траекторий эволюции этих объектов из состояния a по угловому расхождению касательных к этим траекториям в их общей точке a. Таким образом, каждому состоянию x объекта подпространства Ea пространства E ставится в соответствие вектор

x = qY                                                                                          (13)

 из R, как количественная величина q = ||x|| определённого качества

Y  = x / ||x||

в ориентированном ��епером пространстве Ea. При этом dim R = dim Ea = n. 

         Очевидно,

D(Y) = m(Y, Y) = 1                                                                           (14)

 и при произвольном Y  уравнение (14) определяет в пространстве R единичную сферу S с центром в точке a.

Построим  последовательность нормированных базисных векторов W = (Yj: j = 1, 2, ..., n). Если начала всех векторов системы W совместить с точкой a, то концы всех векторов Yj  будут лежать на сфере S. Качественное сходство и расхождение любых двух векторов, исходящих из точки a будет определяться кратчайшей  дугой сферы, соединяющей концы их масштабных элементов (единичных векторов), в том числе и базисных векторов множества W.

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

z = Yja j.

Конец данного вектора также лежит на сфере S.

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

Пусть dij – длина хорды, стягивающей две точки Yi, Yj сферы S, а aij – длина дуги в радианах.  По теореме косинусов

dij2 = D(Yi) + D(Yj) - 2s(Yi)s(Yj) cos aij.

Поскольку

D(Yk) = s 2(Yk) = 1,   k = i, j,

то имеем

cos aij = 1 – ½ dij2.

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

s(bij) = ||bij|| = |Yi ˄Yj| = s(Yi)s(Yj) sin aij = sin aij.

Из основного тригонометрического тождества, вводя соответственно меры сходства и различия объектов сферы

m(Yi, Yj) = cos aij,    n(Yi, Yj) =  sin aij,

приходим к метрическому тождеству

m 2(Yi, Yj) + n 2(Yi, Yj) = 1

на объектах сферы, а умножая данное тождество на произведение дисперсий

D(x, y) = D(x)D(y),

где

x = qY ,     y = pF,    Y, F Î W,

произвольные элементы метрического пространства R, приходим к основному метрическому тождеству (которое представим тождеством Пифагора)

s 2(x, y) = m 2(x, y) + n 2(x, y),

в котором внутреннее (скалярное) произведение векторов

m (x, y) = y*x

характеризует меру их сходства, а внешнее произведение

n 2 (x, y) = |x ˄ y| 2 = b*b = D(b),   b = x ˄ y,

меру их расхождения.

Рассмотрим n-мерное аффинное пространство E. Произвольно выберем в нём некоторую точку a и в окрестности Ea этой точки зафиксируем точки x, xo, y, yo. Пусть в этой окрестности  из точки xo существует непрерывный путь в точку x, а из точки yo  в точку y. Тогда можно соответствующие переходы записать в виде

x = xo + x,   y = yo + y,

где величины x и y будут элементами некоторого другого присоединённого к Ea  пространства R. В пространстве Ea всегда можно ввести топологию, которая обращает присоединённое пространство в векторное пространство. Эта топология индуцируется любой евклидовой метрикой

r(xo, x) = s(x) >= 0,

где s(x) понимается как статистическая характеристика среднеквадратического отклонения состояния x от состояния xo, а выражение

D(x) = s 2(x)

представляет дисперсию этого отклонения, которая вводится как внутреннее в R произведение

D(x) = x*x,

определяемое на RxR с помощью положительно определённой симметрической билинейной формы и которую символически будем записывать в виде

m(x, y) = y*x.

Таким образом, в основном метрическом тождестве реализуется единство таких диалектических антиномий как сходство и различие, качество и количество, форма и содержание и т.п.

Рассмотрим пространство B бивекторов. Учитывая свойства

x˄x = 0,      y˄x = - x˄y

можно определить его размерность

dim B = n(n/2 - 1)

и, если элементами векторного пространства присоединённого к евклидову пространству Ea  являются векторы, которые геометрически можно интерпретировать как направленные отрезки прямой, то геометрической интерпретацией элемента пространства B будет ориентированная площадь параллелограмма,  натянутая на бинарное соответствие (x, y) векторов пространства R.  Такое пространство называется симплектическим.

Как и в евклидовом пространстве в симплектическом пространстве мера порождается функцией

D(b) = b*b = (± s(b)) 2,

где s(b) – площадь соответствующего параллелограмма, а знак определяется ориентацией этой площади. Эта мера задаётся выбором евклидовой структуры в окрестности фиксированной точки a на аффинном пространстве E, которое при надлежащем выборе точек a превращается в риманово многообразие. Действительно, при выборе функции m(x, y) имеем

 D(b) = D(y) D(x) - m 2(x, y) = D(x˄y) = y*Kx

и можно ввести метрику m(a, b) в пространстве BхB на бивекторах a, b с помощью расширения основного метрического тождества

D(b)D(a) = m 2(a, b) + D(a˄b),

или, в общем случае,

D(a˄b) = b*Ka,                                                                                  (15)

где ядро представления является оператором со следующими конструктивными особенностями

K = K(a, b) = M(b, a) – M(a, b) = M* - M.

Здесь M – некоммутативный билинейный оператор,

M(a, b) = ab*,    M* = ba* = M(b, a),

который в квантовой физике [12, 13] определяется как произведение символов и несёт селективные свойства. Очевидно, если b = la, оператор Mстановится эрмитово-самосопряжённым и обращает оператор K в нуль, и наоборот, если K равен нулю, то M* = M.

Если элементам a поставить в соответствие обобщённое векторное пространство R, соответствующее аффинному пространству A, то обращение в нуль выражения (15) приводит к гомотетии [14, стр. 57]

b = a + la.  

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

Разложим направляющий вектор оси конуса по базисным векторам векторного пространства присоединённого к пространству состояния объектов системы. В барицентрических координатах qj, j = 1, 2,  ..., n,  это разложение принимает вид

z = q jYj.

Тогда отклонение качества базисного вектора Yj  от среднего качества системы z можно оценить функцией качественного сходства

m(z, Yi) = Yi*z = rijq j = ri ,

где для коэффициента вариации базисных векторов введено обозначение

rij = m(Yj, Yi).

Учитывая ограничение

D(z) = q*Rq = 1,

где R = ((rij)) – матрица коэффициентов ковариации, приходим к следующей задаче квадратического программирования

max{ri = m(z, Yi): D(z) =1, z >= 0}.

Предположим, что r = max (ri: i = 1, 2, ..., n) решение данной задачи. Тогда подстановка

q = ry

сводит предыдущую задачу максимизации к задаче минимизации

min{F = y*Ry: Ry = 1, y ³ 0},

после решения которой находим

r = 1/(F)1/2,    q = ry,  z = X*q,

где введено обозначение для матрицы нормированных данных

X = (Y1;Y2; …, Ym)T.

Для решения задачи квадратического программирования воспользуемся функцией quadprog MATLAB.

При постановке задачи квадратического программирования в общем виде

min{F = f *y + ½ y*Hy: Ay <= b, Aegy = beg, lb <= y <= ub}

решение получаем, если воспользуемся программой

[y, F] =quadprog(H, f, A, b, Aeg, beg, lb, ub).

В нашем случае

H =2R, f = 0, F = - R, b = - 1, Aeg = 0, lb = 0  

и

[y, F] = quadprog(2*R, [], - R, - b, [], [], lb)                                                               (16)

Находим вектор z 

z = (0.6839   0.1829   0.5414   0.3330   0.0810   0.2972)                                                (17)

и его корреляцию с векторами системы Y:

     1.  0.9668

     2.  0.8563

     3.  0.9268

     4.  0.9315

     5.  0.9431

     6.  0.8819

     7.  0.7961

     8.  0.8482

     9.  0.9211

    10.   0.9037

    11.   0.9729

    12.   0.7961

    13.   0.8214

    14.   0.9330

    15.   0.8688

    16.   0.9313

    17.   0.7961

    18.    0.8785

    19.    0.8869

    20.    0.8371

    21.    0.9844

    22.    0.8818

    23.    0.9170 

Видим, что на конической поверхности шарового сектора лежат объекты 7, 12 и 17. Эти объекты наиболее удалены по качественным признакам от среднего состояния системы. Наименьшее качественное отклонение от среднего имеют объекты 1, 11 и 21.

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

1.     Группа гармоничного обслуживания: Gr_1={Центр-инвест,  Рост.  универсальный банк}. Это банки, попадающие в подгруппы 4 и 1, соответственно. Их графики имеют вид

2.                Группа по обслуживанию кредитов и вкладов: Gr_2 = {Дон-Тексбанк, Донской народный банк, Донхлеббанк, Земельный КБ, Капиталбанк, Стелла-Банк,  Таганрогбанк}. Графики действия этих банков имеют вид

3.     Группа обслуживания активов: Gr_3={Волгодонский  ГКБ, Донактивбанк, Донбанк, Кредит Экспресс, Морозовсккомагробанк, РостПромстройбанк, Сельмашбанк}. Их графики изображены ниже.

 

4.     Группа обслуживания активов и вкладов: Gr_4={Донинвест,  Донкомбанк, ИБРР-банк, Южный регион, Южный Торговый Банк}.  Их графики изображены на следующем рисунке

 

И последняя группа в качественном расслоении – это группа

5.     Группа обслуживания активов и депозитов юрлиц: Gr_5= {МеТраКомБанк,  Южный региональный банк}. Графики этих банков изображены на следующем рисунке

 

Естественно, что каждый объект имеет индивидуальные качественные особенности и занимает своё место в рассматриваемой финансовой системе.

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

   Сравним два массива A1 и A2, заданных данными в табл. 2 и табл. 4 приложения 1. Массив A1 был исследован выше. В качестве эталона был взят, удовлетворяющий условиям (6), вектор p оценки столбцов. Этот вектор и выберем в качестве эталона при сравнении данного множества массивов. С массивом  A2 сделаем те же операции, что и при анализе массива  A1. Элементы каждой её строки делим на соответствующие элементы эталона. Получаем новую матрицу X. А нормируя строки последней, получаем матрицу Xu, с помощью которой находим корреляционную матрицу объектов множества A2

R = Xu'*Xu.

Используя теперь решение задачи квадратического программирования, находим, что объекты 2, 14, 15, 20 наиболее удалены от оси

z2 =(0.6839    0.1829    0.5414    0.3330    0.0810    0.2972)               (18)

нового конуса соцелей. Объектами же качественно близкими к среднему системы становятся объекты 5, 6, 18 и 21.

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

Таблица 4. Вероятности присутствия факторов в системе.

 

Система

Факторы

Капитал

Кредиты

Др. Акт.

Вклады

Др. деп.

Прибыль

A1

A2

    0.4677    0.0334      0.2931     0.1109      0.0066      0.0883 

    0.4131    0.1058      0.0864     0.1915      0.0065      0.1966 

 

Качественное расхождение систем A1 и A2 определяется коэффициентом корреляции векторов (16) и (17) zz2 = 0.9423, что в угловом измерении составляет приблизительно 20o.

Рассмотрим корреляцию объектов систем с соответствующими направляющими осевых векторов (табл. 5).

Из табл.5 видим, что объекты качественно сблизились (конус качественного расхождения сузился на 5,5%) и ось конуса сместилась на 3,8о, что свидетельствует о том, что произошедшие качественные сдвиги в структуре объектов незначительны. Эти сдвиги можно наблюдать на соответствующих графиках (см. прилож. 3).

Таблица 5. Коэффициенты корреляции субъектов систем с их направляющими векторами осей симметрии.

 

Объекты (банки)

Коэффициенты корреляции систем

A1

A2

 1. Волгодонский  гкб

 2. Дон-Тексбанк

 3. Донактивбанк

 4. Донбанк

 5. Донинвест

 6. Донкомбанк

 7. Донской народный банк

 8. Донхлеббанк

 9. Земельный КБ

10. Капиталбанк

11. ИБРР-банк 

12. Кредит Экспресс

13. МеТраКомБанк

14. Морозовсккомагробанк

15. Рост.  универсальный банк

16. РостПромстройбанк

17. Сельмашбанк

18. Стелла-Банк

19. Таганрогбанк

20. Центр-инвест

21. Южный регион

22. Южный региональный банк 

23. Южный Торговый Банк   

0.9668

0.8563

0.9268

0.9315

0.9431

0.8819

0.7961

0.8482

0.9211

0.9037

0.9729

0.7961

0.8214

0.9330

0.8688

0.9313

0.7961

0.8785

0.8869

 0.8371

 0.9844

 0.8818

 0.9170

 

0.9357

0.8400

0.9339

0.9189

0.9683

0.9703

0.9151

0.9058

0.9164

0.8800

0.8558

0.8414

0.8865

0.8400

0.8400

0.9392

0.9074

0.9787

0.9413

 0.8400

 0.9702

 0.8509

 0.9139

 

Литература.

1.     Бредон Г. Теория пучков /М., Наука, 1988.

2.     Бондал А.И. Представления ассоативных алгебр и когорентных пучков /Изв. АН СССР. Сер. матем., 1989,         том 53, выпуск 1, стр. 25-44.

3.     Месарович М., Мако Д., Такахара И. Теория иерархических многоуровневых систем /М., Мир, 1973.

4.     Кондрашов И.Е., Королёв C. MATLAB как система программирования научно-технических расчётов /М., Мир, 2002.

5. Хайкин С. Нейронные сети //М., С.-П., 2008.

6. Новиков Б.В. Лекции по теории категорий //Луганск, 2004.

7. Нейман И. Математические основы квантовой механики //М., Наука, 1964.

8. Дирак П.И.А. Лекции по квантовой механике//Ижевск, Ижевская республиканская типография, 1998.

9. Кондрашов В.Е., Королёв С., MATLAB, М., Мир, 2002.

10.     Гантмахер Ф.Р. Теория матриц //М., Наука, 1988.

11.     Идлис Г.М. Единство естествознания по Бору и единообразные взаимосвязанные периодические системы физики, химии, биологии и психологии /Исследования по истории физики и механики, 1990 //М., Наука, 1990.

12.     Грин Х. Матричная квантовая механика /М., 2009.

13.     Кемпфер Ф. Основные положения квантовой механики /М., 2007.

14.     Берже М. Геометрия, том первый М., Мир, 1984.

 

 

 механике //