Файл: Шепелев, И. Г. Математические методы планирования и управления в строительстве конспект лекций.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 29.10.2024
Просмотров: 48
Скачиваний: 0
С помощью датчика случайных чисел вырабатывается слу чайная величина РСП (0,1). Допустим, выработалась h = = 0,750, сравниваем ее с начальными вероятностями матрицы. Так как gi > 0,556, но меньше 0,756, то выбирается состояние 2, Пусть следующая случайная величина | 2 = 0,223. Сравни ваем с вероятностями второй строки той же матрицы, аналогич
но выбирается состояние 2 .
Третье случайное |
число |
|з |
(принимаем |
g3 = 0,538) |
будет |
|||||||||
сравниваться с вероятностями второй |
строки матрицы П, так |
|||||||||||||
как |
выбор |
второго столбца |
в предыдущей |
матрице |
означает |
|||||||||
переход на вторую строку последующей. |
|
выбирается |
||||||||||||
Так как |
больше 0; |
0,166; |
0,332, |
но меньше 1, |
||||||||||
состояние 4 и т. д. для всех пяти испытаний серии. |
результаты |
|||||||||||||
Всего |
проигрывается |
10 |
серий из 5 испытаний, |
|||||||||||
приведены в табл. 2 0 . |
что наибольшее число раз (5) |
повторил |
||||||||||||
Из таблицы |
видно, |
|||||||||||||
ся следующий набор состояний: |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
2 —2 —4—4—6 . |
|
|
|
|
|||
Эти состояния и принимаются в качестве результата |
моде |
|||||||||||||
лирования. |
|
|
|
|
|
|
|
|
|
|
|
|
||
Обращаясь к табл. 17, получаем следующие значения поста |
||||||||||||||
вок: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( П |
|
|
( 2 ) |
( 3 ) |
|
( 4 ) |
( 5 ) |
|
|
|
||
|
|
q, = |
20%, qt = 0, q, = |
40%, qi = 0, q, = 40%. |
|
|||||||||
Результаты |
можно |
рассматривать |
как прогноз |
выполнения |
||||||||||
недельного заказа. |
|
|
|
|
|
qi =? 42,47 тыс. |
шт., |
мож |
||||||
Зная |
величину недельного заказа |
|||||||||||||
но определить абсолютную величину поставок: |
|
|
|
|||||||||||
( 4 ) |
8,49 тыс. шт.; |
( 2 ) |
|
( 3 ) |
|
|
|
( 4 ) |
(5 ) |
|
||||
qt = |
qt = |
0;qt = |
16,98 тыс. шт. qt = 0; |
q, — 16.98 |
||||||||||
тыс. шт. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Формула |
(8.5.7) |
для расчета Уп может быть преобразована |
||||||||||||
следующим образом |
|
|
|
|
|
|
|
|
|
|
||||
— |
|
1 |
/ |
|
|
( 1 ) |
|
( 2 ) |
( 3 ) |
( 4 ) |
( 5 ) |
|
|
\ |
Уп = |
2 . 5 |
\9 Уп—1 + 9 qn + |
7 qn + 5 qn-f- 3qn -j- qn — 4 Sn + |
ynj. |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
(8.5.10) |
|
Подставляем в |
(8.5.10) |
вычисленные значения поставок и задан |
||||||||||||
ное значение потребности (табл, |
17), вычисляем величину |
|
У( = ( 9 - 0 + 9 3 , 4 9 + 7 . 0 + 5Л6,98+3.0+16,98—4.0+42,47) =22,08
тыс. шт.) |
(8.5.11) |
121
Значение функции затрат h (42,47; 42,47) = 333,12 определя ется подстановкой вычисленного значения У1 в формулу (8.5.8).
Это значение функции заносится в первую строку табл. рис. 25 (столбец У) (Уо), так как для первого шага
L1 (Y0) = min[/1 (q1 -y1) + 0] = /1 (q1; у,). |
(8.5.12) |
Аналогично вычисляются остальные значения функции при других начальных уровнях запасов (табл. рис. 25).
ч |
16 |
|
17,6 |
27,6 |
37,6 |
47,6 |
W |
|
V |
|
|||||||
\ |
|
|
|
|
|
|
|
|
0 |
349,61 327,46 |
327,42 |
352,21 385,62 |
660,54 |
||||
+ |
|
|
|
333,12 |
333,12 |
27,6 |
||
|
333,12 333,12 |
333,12 |
|
|||||
ю |
327А |
355,80 |
355,25 |
4 2 6 ,9 3 '///// |
660,52 |
|||
+ |
|
|
|
333,12 |
|
7,6 |
||
|
333,12 333,12 |
333,12 '///А |
|
|||||
20 |
354,91 |
434,14 |
433,07 У /^ |
У /// |
686,03 |
|||
333,12 |
333,12 |
333,12 |
|
7,6 |
||||
|
У Л ' / ' / а У /. |
|
||||||
30 |
432,41 |
|
|
|
|
|
7,6 |
765,53 |
+ |
|
3 |
3 |
1 2 / / / |
|
|||
|
зззЦ 1 3 |
/ / / / / / Л |
|
|||||
40 |
563,50 7 |
/7 |
? |
У / / / / / У |
/ / v |
916,62 |
||
4- |
|
|
|
|
|
76 |
||
|
333,12 |
У |
|
/ ш |
/ Ш |
/ , |
|
|
|
Рис. 26. |
Таблицы потерь за вторую неделю |
Структура табл. рис. 26 несколько отличается от предыдущей. Каждому начальному уровню запасов соответствует не одно, а пять значений поставок, так как уровень запасов' на конец неде ли уже не определяется однозначно.
Числа в таблице представляют собой сумму потерь рассмат риваемой недели, оптимальных для последующей.
Для каждого начального уровня запасов у\ значение функ
ции L2 (//i) представляет собой минимальную из всех сумм в клетках рассматриваемой строки. Так, например, если у\ — О,
q2 = 7,6 тыс. шт., то г/г = 0 при S2 = 7,6 тыс. шт. и у (7; 6 ; 0 ).
Вторым слагаемым в этой клетке является значение функции ЬДО), выбранное из табл. рис. 25, так как остаток на конец дан ной недели (у2 = 0 ) равен остатку на начало следующей {у2 ~ у /
Поэтому в табл. рис. 26 выбирается значение функции Li(y0), соответствующее уо = 0 . Аналогично заполняются и другие
122
клетки первой строки табл. рис. 26. Минимальные потери при данном уровне начального запаса равны 660,54, а соответствую щая им поставка равна 27,6 тыс. шт.
Это значение потерь и соответствующее ему значение поставки приведены в крайних правых столбцах таблицы рис. 26. Затем заполняется вторая строка и т. д.
Заштрихованные клетки в таблице соответствуют запрещен ным вариантам, т. е. значения поставок и остатков не отвечают ограничению, заданному формулой (8.5.2)..
Рис. 27. Таблицы потерь за третью неделю
/
Затем строятся таблицы (рис. 27—34).
После построения таблиц можно приступать к выбору опти мального плана. Выбор осуществляется по таблицам, начиная с рис. 34. Они содержат только одну строку, соответствующую на чальному уровню г/g = г/ср = 35 тыс. шт.
При таком начальном уровне запаса минимальным потерям для десятинедельного планового периода соответствует поставка qm = 0,9 тыс. шт., и остаток на конец недели г/ю = 20 тыс. шт. (формула 8.5.5).
Остаток на начало следующей недели г/о = 2 0 тыс. шт. Ему со ответствует третья строка табл. рис. 33. Минимальным потерям
за оставшиеся девять недель соответствует поставка |
qo = 7,6 |
тыс. шт. Затем переходим к табл. рис. 33 и т. д. |
который |
В результате получаем оптимальный план поставок, |
|
приведен в табл. 2 1 . |
|
123
Рис. 28.~Таблицы потерь за четвертую неделю
|
|
V s |
4,4 |
14,4 |
24,4 |
Ъ |
W |
) |
|
|
|
40 (955,81536/5 |
|
|
4,4 |
2Щ95 |
|
||
|
|
Рис. |
29. Таблицы |
потерь |
за пятую |
|
|||
|
|
|
|
неделю |
|
|
|
|
|
|
|
|
|
|
|
|
|
Т а б л и ц а 21 |
|
% |
Че |
W |
|
|
|
|
|
План |
План |
|
Месяц |
Недели |
поставок |
потребности |
|||||
10 |
56,4 |
2964J5 |
|
|
|
1 |
|
0,9 |
15,9 |
|
|
|
|
|
IV |
2 |
|
32,7 |
52,7 |
20 |
46,4 |
2976,2 |
|
|
3 |
|
56,82 |
46,8 |
|
|
|
|
4 |
|
41,4 |
41,4 |
|||
|
|
|
|
|
|
5 |
|
56,4 |
26,6 |
30 |
36,4 |
2998,8 |
|
Итого |
|
|
188,22 |
183,22 |
|
|
|
|
6 |
|
4,4 |
4,4 |
|||
|
|
|
|
|
|
|
|||
40 |
26,4 |
зовщ |
|
|
|
7 |
|
9,2 |
9,2 |
|
|
V |
8 |
|
|||||
|
|
|
|
|
|
9 |
|
7,6 |
7,6 |
Рис. 30. Таблицы потерь |
|
|
|
10 |
|
12,47 |
— |
||
|
Итого |
|
|
33,67 |
21,2 |
||||
за шестую |
педелю |
|
|
|
124
Рис. 31. Таблицы потерь за седьмую неделю
Рис. -32. Таблицы потерь за восьмую неделю
Рис. 33. Таблицы потерь за девятую неделю
0,9 |
10,9 |
20,9 |
30,9 |
Яю |
B w (V |
|
|
|
|
||
328,56 |
327,60 |
328,36 |
336,67 |
0,9 |
4284,88 |
35 |
4- |
3976,09 |
3992,36 |
||
3056,3 |
3974,0 |
|
|
Рис. 34. Таблицы потерь за десятую неделю
Задача решается в скользящем режиме, т. е. оптимальный план ежемесячно корректируется, что позволяет более оператив но планировать поставки, постоянно уточняя их в ходе решения задачи.
§ 8.6. Моделирование и прогнозирование процесса поставок строительных материалов
Процесс поставки строительных материалов является сложным стохастическим процессом. Это вызывается большим количеством поставщиков и потребителей, сложностью производства и достав ки материалов, влиянием климатических, организационных и других факторов. Поэтому даже в, условиях централизованного планирования и фондового распределения материалов часто на блюдаются отклонения поставок от плановых как по объемам, так и по срокам. Для прогнозирования поставок можно: а) ап проксимировать процесс потавок факторными корреляционными
126
моделями; б) считать процесс поставок случайным процессом, определив параметры этого процесса.
Рассмотрим случай аппроксимации процесса поставок фак торными моделями. В качестве зависимых переменных необходи мо рассматривать некоторые переменные, являющиеся разностью между фактическими и плановыми '(договорными) сроками или объемами поставок строительных материалов.
Анализ фактических данных показывает, что плановая равно мерность поставок по большинству строительных материалов не подтверждается практикой.
Устойчивые многолетние данные показывают, что существует зависимость между отклонениями в сроках и объемах поставок, с одной стороны, и рядом перечисленных выше факторов, с дру гой. Однако большинство из этих факторов очень трудно учесть, так как многие не имеют количественного выражения. Анализ показывает, что ряд климатических и организационных факторов статистически связан с календарным временем года. Сезонность особенно существенно сказывается на работе транспорта, так как периодически он отвлекается на перевозки сельскохозяйственных продуктов,-
В гл. II приведен пример аппроксимации процесса поставок цемента тригонометрической функцией времени. Формула (2.7.6) является выражением зависимости отклонений поставок цемента от времени года. Как видно из оценки уравнения (2.7,6), где
г| =0,97 и е = 15,8, в этом случае получена достаточно надежная для прогнозирования модель. Часто при учете функции задержки поставок в выражении (8.4.5) можно ограничиваться функцией только одного переменного, времени. Если функция получена ме тодом корреляционного анализа, то все неучтенные факторы вы ражаются через случайную компоненту этой функции I. Характер этой компоненты и метод получения доверительных интервалов к уравнению рассмотрены в §3.5.
При учете в модели помимо фактора времени также и других факторов следует-ожидать повышения точности модели, умень шения дисперсии случайной компоненты / и уменьшения довери тельного интервала к-модели. Тогда можно строить более надеж ные модели, пригодные для целей планирования и управления реальными процессами в строительстве и его материально-техни ческом обеспечении. Так, например, для одного из трес тов была построена модель отклонения фактических поставок це мента от плановых. В качестве фактора в этой модели рассматривался фактор времени. Корреляционное отношение для этой формулы было равно 0,4. Для повышения тесноты связи в
127
модель ввели еще один фактор — расстояние между поставщиком и потребителем цемента. При этом была получена формула:
Y = (— 9,41 -}- 6,42 sin x t + 4,26 cos x t -(- 4,87 sin 2 x t —
— 1,99 cos 2 x x— 1,45 sin 3^! — 6,69 cos З л , — l,74sin4.v1 +
+ 0,01 cos 4*,) ( — 2,30+ 15,95 x 2), |
(8.6.1) |
где У — отклонения фактических поставок цемента от |
плано |
вых, %; |
|
Xi — календарное время года, в днях от начала года; Хо — средневзвешенное расстояние поставок, км.
Средневзвешенное расстояние поставок определяется из фор мулы
1=П
М m i 4i |
|
x 2 = -ig ----- . |
(8 .6 .2 ) |
‘fq.
i=1
где Ш] — расстояние между трестом и i-ым заводом-поставщиком,
км (i = 1 , 2 , 3,..., п);
q, — объем поставок с i-ro завода за рассматриваемый пе риод.
При аппроксимации поставок двухфакторной моделью точность модели повысилась, корреляционное отношение повысилось с 0,4 до 0,97. Такая модель может реально применяться при планиро вании поставок, при этом можно с высокой вероятностью прогно зировать поставки материалов по сравнению с фондами или заяв ками в разные периоды года, с учетом отдаленности поставщиков, что даст возможность более рационально планировать заказы и строительное производство.
Однако факторные, в том числе и динамические, модели типа (8.6.1) не являются единственными. Процесс поставок некоторых материалов можно отнести к стационарным процессам с задан ными законами распределения. Тогда поставки могут характери зоваться параметрами распределения вероятностей этих законов. Например, нормальным законом, биноминальным и т. д.
Среди различных классов случайных процессов большой инте рес представляют марковские процессы, названные так в честь русского математика А. А. Маркова, впервые исследовавшего их свойства.
Выше (см. §8.5) рассмотрен метод оптимизации поставок строительных материалов с определением среднего запаса или при помощи цепей Маркова. Здесь рассмотрим методику аппрок симации поставок Марковским процессом.
128