- •Введение
- •Методы моделирования и анализа детерминированных систем на основе анализа конечных изменений
- •Понятие и классификация факторного анализа функции многих переменных
- •Общая характеристика анализа конечных изменений
- •Модели лагранжева анализа конечных изменений
- •Методы лагранжева анализа конечных изменений
- •Графоструктурное моделирование лагранжева анализа конечных изменений
- •Прикладные аспекты лагранжева анализа конечных изменений
- •Основные понятия и определения
- •Различительные исчисления
- •Построение рабочих моделей анализа конечных изменений функций
- •Лагранжев анализ конечных изменений
- •Лагранжев анализ конечных изменений функций специальной структуры
- •Цепной лагранжев анализ конечных изменений и его модификации
- •Цепной анализ конечных изменений. Связи с интервальным анализом и экономическим факторным анализом
- •Цепной экономический факторный анализ
- •Зависимость от матрицы факторов
- •Зависимость от индивидуальных наборов параметров разной длины
- •Решение прикладных задач методами лагранжева анализа конечных изменений
- •Лагранжев анализ конечных изменений в оценке финансового состояния организации
- •Лагранжев анализ конечных изменений для структурной оптимизации объекта управления дорожным движением
- •Использование новых методов цепного лагранжева анализа конечных изменений для решения прикладных задач
- •Заключение
При этом коэффициенты факторного влияния равны:
1 |
= |
1 |
; |
|
(4.22) |
||
|
|
||||||
3 + 3 + 1 |
|||||||
|
|
|
|
|
|
||
2 |
= |
3 + 3 |
; |
|
(4.23) |
||
3 + 3 + 1 |
|||||||
|
|
|
|
|
|
||
3 = |
2 |
+ 2 − 1 − 1 |
. |
(4.24) |
|||
|
|||||||
|
|
( 3 + 3 + 1)2 |
|
||||
Таким образом, в соответствии с обозначениями для коэффициента
получаем следующие искомые выражения:
|
|
|
|
= |
|
|
|
1 |
|
|
; |
|
|
|
|
(4.25) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
+ + 1 |
|
|
|
|
|
|||||||
|
|
|
|
= |
|
+ |
; |
|
|
|
|
|
(4.26) |
||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
+ + 1 |
|
|
|
|
|
|||||||||
|
|
3 |
= |
+ − (1 − )( + ) |
, |
|
(4.27) |
||||||||||||
|
|
|
|
|
( + + 1)2 |
|
|
|
|
|
|||||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= √ |
|
|
|
|
|
|
− |
|
|
|
|
|
и = . |
(4.28) |
|||||
|
|
|
|
+ 1) |
+ 1) |
|
|||||||||||||
|
|
( |
+ 1)( |
+ |
|
( |
|
|
|
|
|
||||||||
Разработанные модели определения факторного влияния при детерминированной структуре моделей с помощью лагранжева анализа конечных изменений для некоторых финансовых показателей могут быть применены для исследования вклада аргументов функции в ее изменение, что определяет большое прикладное значение полученных результатов. Полученные общие формулы могут быть использованы не только для экономической, но и для производственной сфер исследования.
4.2Лагранжев анализ конечных изменений для структурной оптимизации объекта управления дорожным движением
Всвязи с ростом автомобилизации населения, а, следовательно,
иего транспортной мобильности, проблемы, связанные с управлением
53
транспортными потоками, становятся все более актуальными. Одним из видимых проявлений увеличения количества личного транспорта являются транспортные заторы на пересечении городских магистралей.
Для решения проблемы транспортных заторов на регулируемых перекрестках обычно в качестве показателя качества их функционирования используется величина средней транспортной задержки регулирования. Транспортная задержка с точки зрения математики — это функция одного или нескольких аргументов. В одном из подходов к оптимизации минимизация этой функции по соответствующим аргументам дает оптимальное значение длительности цикла и её оптимальное распределение по фазам регулирования.
Интересной является задача исследования поведения этой функции при изменении её аргументов [8; 10; 23]. Ответ на такой вопрос может быть получен при помощи лагранжева анализа конечных изменений функции агрегированной транспортной задержки.
Рассмотрим класс рациональных функций специального вида:
|
|
|
|
|
|
|
|
|
|
|
( 1 |
, ..., ) = |
+ |
∑ |
|
∑ |
2 |
, ≥ 0, |
0, , ≥ 0. |
(4.29) |
|
|
||||||||||
|
0 |
+ |
|
|
|
|
||||
|
|
|
|
=1 |
|
=1 |
|
|
|
|
|
|
|
|
∑ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
=1 |
|
|
|
|
|
Начальным значениям аргументов 0, = 1, ..., отвечает начальное значение функции ( 01, ..., 0 ). Аргументы претерпевают изменения
, которые порождают новые значения аргументов 0 + и функции
( 01 + 1, ..., 0 + ).
Частные производные указанной функции:
|
|
|
|
|
|
|
|
|
|
|
= |
|
+ |
|
+ |
∑ |
|
2∑ |
2 |
|
|
|
|
. |
||||||||
∂ ( 1, ..., ) |
+ 2 |
0 |
|
+ |
|
|||||
|
|
|
=1 |
|
=1 |
|
|
|||
∂ |
|
=1 |
|
|
|
( =1 ) |
|
|
|
|
|
|
∑ |
|
|
|
∑ |
|
|
|
|
Непосредственное приращение функции может быть найдено как
( 10 + 1, ..., 0 + ) − ( 1, ..., ), |
(4.30) |
54
а с другой стороны лагранжево приращение
|
∂ |
|
||
∑ |
|
|
|
|
∂ ( 10 + · 1, ..., 0 + · ) · . |
(4.31) |
|||
=1 |
||||
|
|
|
||
Приравнивая выражения (4.30) и (4.31) и решая полученное уравнение относительно , получаем:
− =1 + |
√ |
|
|
|
|
− =1 − |
√ |
|
|
|
|
|
=1 |
=1( 0 |
+ ) |
|
=1 |
=1( 0 |
+ ) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1 = |
∑ |
|
∑ |
∑ |
|
, 2 = |
∑ |
|
∑ |
∑ |
|
. |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
∑ |
|
|
|
|
|
∑ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
=1 |
|
|
|
|
=1 |
|
|
|
|
Всилу предположения о положительности значений аргументов
иих приращений только значение 1 удовлетворяет условию теоремы Лагранжа о принадлежности параметра промежуточной точки промежутку (0, 1).
Процесс последовательного использования оптимизации и анализа системы, основанных на использовании функции транспортной задержки в качестве критерия оптимизации и численной оптимизации и анализа с помощью теоремы Лагранжа о промежуточной точке как математического аппарата, представлен в виде структурной схемы на рис. 4.
Рис. 4. Оптимизация и анализ процесса функционирования регулируемого перекрестка
При этом на этапе численной оптимизации находится оптимальный набор времени цикла регулирования и его распределения по фазам. Затем на этапе лагранжева анализа конечных изменений проводится (на основе статистической информации, собранной на данном объекте регулирования) предсказательный анализ изменения интенсивности в крити-
55
ческих направлениях фаз с целью выявления фазы, оказывающий наибольшее воздействие на изменение задержки на перекрестке. Особенность анализа состоит в том, что найденная фаза не обязана обладать наибольшим среди всех рассматриваемых в системе отношением интенсивности к пропускной способности.
В качестве примера рассмотрим двухфазный регулируемый перекресток (рис. 5). Фаза 1 включает в себя направления движения 1,2,5,6; критическим направлением в этой фазе является направление 2. Фаза 2 объединяет направления 3,4,7,8; критическим направлением в этой фазе является направление 4.
Рис. 5. Пример двухфазного регулируемого перекрёстка
Алгоритм получения средней задержки для регулируемого пере-
∑
( ) ·
крёстка приводит к модели вида = |
=1 |
|
. Анализ этой функции, |
|
|
||
|
|
|
∑
=1
где в качестве каждого ( ) выступает выражение
( ) = ( − )2 +
( − )
+ 4 |
|
+ |
|
|
|
|
1 |
|
+ 14400 2 2 |
, |
|||||
|
|
( |
|
) |
√( |
|
|
− ) |
2 |
|
|
(4.32) |
|||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
56
является вычислительно сложной задачей. В выражении (4.32) — время цикла светофорного регулирования (с), — длительность фаз светофорного регулирования (с), — коэффициент, определяющий тип прибытия транспортных средств, — интенсивность прибытия транспортных средств на направлении (прив. авт. / ч), — поток насыщения направления (прив. авт. / ч).
Таблица 1
Результаты анализа функции задержки регулирования
1 |
2 |
1 |
2 |
1 + · 1 |
2 + · 2 |
, c |
1, % |
2, % |
|
|
15 |
5 |
1007,465 |
72,488 |
2,485 |
57 |
43 |
|
70 |
25 |
1007,431 |
82,385 |
197,180 |
48 |
52 |
|
|
|
|||||||
|
75 |
5 |
1036,824 |
72,455 |
329,072 |
58 |
42 |
|
|
|
|||||||
1000 |
|
25 |
1036,663 |
82,221 |
444,307 |
49 |
51 |
|
|
|
|||||||
|
15 |
5 |
1007,467 |
122,489 |
2,313 |
29 |
71 |
|
|
|
|||||||
|
120 |
25 |
1007,434 |
132,390 |
287,127 |
26 |
74 |
|
|
|
|||||||
|
75 |
5 |
1036,853 |
122,457 |
159,473 |
30 |
70 |
|
|
|
|||||||
|
|
25 |
1036,698 |
132,233 |
389,049 |
27 |
73 |
|
|
|
|
||||||
|
|
15 |
5 |
1507,476 |
72,492 |
2,271 |
83 |
17 |
|
70 |
25 |
1507,453 |
82,421 |
100,154 |
71 |
29 |
|
|
|
|||||||
|
75 |
5 |
1537,034 |
72,469 |
398,044 |
85 |
15 |
|
|
|
|||||||
1500 |
|
25 |
1536,921 |
82,307 |
437,816 |
72 |
28 |
|
|
|
|||||||
|
15 |
5 |
1507,477 |
122,492 |
2,494 |
44 |
56 |
|
|
|
|||||||
|
120 |
25 |
1507,454 |
132,424 |
218,672 |
40 |
60 |
|
|
|
|||||||
|
75 |
5 |
1537,048 |
122,470 |
267,177 |
45 |
55 |
|
|
|
|||||||
|
|
25 |
1536,938 |
132,313 |
434,851 |
41 |
59 |
|
|
|
|
Для упрощения анализа применим к функции (4.32) разложение в ряд Маклорена (до трёх членов). В результате такого приема получаем модель, структурно схожую с (4.29). Здесь:
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
2 |
|
|
|
|
|
0 = 0, |
|
= |
|
|
|
(1 |
− |
|
) |
|
, |
||||||
|
|
|
2 |
|
|
|||||||||||||||
= |
1 |
|
( − )2 |
+ |
1 |
|
( |
|
− |
|
( − 7200) |
, = 1, ..., . |
||||||||
|
|
|
|
|
||||||||||||||||
2 |
|
4 |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
) |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
2 |
|
|
|
|
Для некоторых начальных значений интенсивностей и их приращений проведён полный факторный эксперимент (значение потока насы-
57
