Добавил:
ivanov666
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Применение логистико-математических моделей в лесной отрасли. Монография
.pdf
∞
g
С учетом (2.57) имеем:
∞
11
=− − − =
Hdw
∫
−∞
πσ πσ
22
∑∑
⎧⎫ ⎧⎫
ww ww
() ()
⎪⎪ ⎪⎪
exp log exp
⎨⎬ ⎨⎬
⎪⎪ ⎪⎪
⎩⎭ ⎩⎭
22
∞
11
=− − +
log exp
πσ πσ
22
∞
e
log
+− =πσ+
2
πσ
∞∞
et e t
log log
+−=πσ+−+−=
ππ
22 22
2
tdt te dt
∫∫
−∞ −∞
∫
∑
−∞
⎧⎫ ⎧⎫
exp log 2 exp
⎨⎬ ⎨⎬
⎩⎭ ⎩⎭
⎧⎫
⎪⎪
exp log 2
⎨⎬
⎪⎪
⎩⎭
22
22
∫
∑∑
−∞
ww ww
−−
()()
22
σσ
22
−−
22
σσ
∑∑
⎧⎫
ww
()
⎪⎪
⎨⎬
⎪⎪
⎩⎭
22
22
∑∑
∑
2
dw
dw dt
(
−
2
σ
∑
xw
−
σ
∑
=σ
−
2
dw
∑
t
=
∑
2
t
2
−∞
1
e
=πσ+ π=πσ+ =πσ=
log 2 2 log 2 log log 2
log 1
∑∑
π
22
⎛⎞
⎜⎟
log 2
=πσ
⎜⎟
⎝⎠
nm
ex
∑∑
ji
11
==
⎛⎞
2
⎜⎟
jij
⎝⎠
2
eez
1
2
2
⎛⎞
⎜⎟
⎜⎟
⎝⎠
n
∑
=
j
1
22
jj
2
(2.58)
В соответствии с этим критерий приобретает вид:
RX e x x x
=πσ + −
()
2
1
⎛⎞
⎜⎟
lo
⎜⎟
⎝⎠
nm mn
⎛⎞
2.
2
∑∑ ∑∑
ji ij
⎜⎟
jij ijij
11 11
== ==
⎝⎠
2
()
0
()
2
(2.59)
Третий этап. На этом этапе выбирается начальный опорный
план задачи и реализуется его улучшение. В качестве начального
плана задачи может быть использован модальный набор, обнуляющий второе слагаемое критерия. Дальнейшая итерационная процедура улучшения плана представляет собой редукцию к двумерному случаю процедуры улучшения плана [60].
81

2.5.2. Транспортная задача со случайными
=
⎧
⎪
⎨
⎪
x
R
стоимостями перевозок
Рассмотренная выше транспортная задача, в которой случайным
был спрос на товар, значительно усложняется, если случайными являются стоимости транспортировок. Рассмотрим эту задачу. Примем,
что параметры целевой функции задачи есть случайные величины с
известными плотностями распределения. На стоимости транспортировок товара реально влияет большое число разнообразных факторов.
При этом обычно считают, что
их совместное влияние в соответствии
с центральной предельной теоремой теории вероятностей приводит к
нормальности распределения случайных стоимостей перевозок [60].
В соответствии с этим введем набор плотностей:
σ
2
2
ij
, 1, 2, ...,im
, 1, 2, ...,jn=, (2.60)
ϕ= −
c
()
ij ij
2
πσ
j
случайных значений
1
⎪⎪
exp
⎨⎬
⎪⎪
⎩⎭
c
.
ij
⎧⎫
−
cm
()
ij ij
2
Тогда плотность распределения случайной величины суммарной
mn
стоимости транспортировки
fR fLX
== −
() ( )
()
RLX cx
==
()
1
∑∑
11
ij
==
()
exp
2
πσ
Σ
⎩⎭
имеет вид
ij ij
−
Rm
Σ
2
2
σ
Σ
2
⎫
⎪
, (2.61)
⎬
⎪
где
mn
mmx
=
∑∑
Σ
==
`1
ij
ij ij
σ= σ
,
mn
222
∑∑
Σ
==
`1
ij
. (2.62)
ij ij
В такой постановке задача состоит в отыскании плана транспор-
тировок
Xx=
, доставляющего экстремальное значение какой-либо
()
ij
естественным образом выбранной функции с учетом плотности (2.61) и
удовлетворяющего ограничениям (2.2)–(2.4). Эта задача рассматривалась в работе [67], где критерий был определен следующим образом.
Задается некоторое пороговое значение суммарной стоимости транспортировки
, превышение которого расценивается как признак неэф-
n
фективности соответствующего плана транспортировки. Тогда в качестве критерия предлагается выбрать вероятность того, что суммарная
стоимость перевозки превысит пороговую. Эта вероятность равна:
82

PR R
f
∞
X
(
x
=
=
j
=
(
x
X
≥= =
()
n
R
∫
n
RdR
()
2
u
−
2
(2.63)
Σ
σ
2
Σ
2
Σ
dR e du
Rm
−
n
σ
Σ
Σ
Rm
∞∞
=− =
R
11
∫∫
22
πσ π
n
Σ
⎧⎫
⎪⎪
exp .
⎨⎬
⎪⎪
⎩⎭
−
Rm
()
2
−
n
σ
Σ
Далее исходная задача преобразуется к следующей: найти удов-
летворяющий ограничениям план перевозок
, минимизирующий
(2.63), что эквивалентно максимизации нижнего предела в интеграле
(2.63), который равен
mn
−
Rmx
∑∑
nijij
11
ij
==
mn
⎛⎞
⎜⎟
∑∑
⎜⎟
11
ij
==
⎝⎠
22
σ
ij ij
x
. (2.64)
1
2
⎛⎞
⎜⎟
⎝⎠
,
Критерий
RX
()
−
Rm
n
RX
XR
mn mn mn mn
RR
nn
∑∑ ∑∑ ∑∑ ∑∑
AA
ij ij ij ij
== == == ==
===
11 11 11 11
⎛⎞⎛⎞⎛⎞
⎜⎟⎜⎟⎜⎟
⎜⎟⎜⎟⎜⎟
⎝⎠⎝⎠⎝⎠
ij ij ij ij ij ij ij
mn mn mn
∑∑ ∑∑ ∑∑
ij ij ij
11 11 11
== == ==
==
()
)
преобразуется к виду
−−
22 22 22
σσσ
ij ij ij ij ij ij
mx m x dx
111
222
xxx
σ
Σ
Σ
где
dm
ij ij
Далее понятно, что максимум
ре, что и максимум
n
m
a
∑
i
1
=
i
2
()
XR
FX R X
() ()
, 1, 2, ...,im
−
, 1, 2, ...,
)
XR
достигается на том же набо-
n
,
. В соответствие с этим вводится:
mn
⎛⎞
⎜⎟
∑∑
⎜⎟
11
==
2
==
ij
⎝⎠
mn
∑∑
==
11
ij
dx
σ
2
ij ij
. (2.65)
22
ij ij
R
Таким образом, задача сведена к следующей: найти план
максимизирующий дробно-квадратичную функцию (2.65) и удовле-
83
,

творяющую ограничениям (2.32)–(2.34). При этом показано, что как
()(
)
}
угодно точное решение этой задачи может быть получено в результате решения рекуррентной последовательности задач квадратичного
программирования [67].
Принципиальный недостаток изложенного подхода состоит
в том, что численное значение критерия жестко зависит не только
от плана транспортировки, но и от величины порога.
Исследования показывают, что полученный в результате решения этой задачи (2.65), (2.32)–(2.34) сам план транспортировок существенно зависит от величины порога, что приводит либо к необходимости надлежащего его обоснования, либо к отказу от использования
предложенной в [67] методики.
Рассмотрим альтернативный вариант решения транспортной задачи со стохастически заданной матрицей стоимостей перевозок [66].
Используем вновь описанный выше подход,
приводящий к формированию составного критерия. Аналогично предыдущему, предварительно решается задача отыскания плана транспортировок, оптимального в среднем. Эта задача, с учетом (2.30), формулируется так: найти
план
Xx=
, минимизирующий (2.1) и удовлетворяющий ограни-
{}
ij
чениям (2.2)–(2.4).
00
Пусть
Xx=
{
– решение этой задачи. Далее, с учетом
ij
(2.58), (2.61), (2.62) «компактность» тела неопределенности относительно численного значения случайной суммарной стоимости транс-
1
портировок определяется выражением
⎛⎞
HX e x
()
log 2
=π σ
⎜⎟
⎜⎟
⎝⎠
nm
∑∑
11
==
ji
22
ij ij
2
.
Тогда критерий качества плана транспортировок будет иметь вид
1
RX e x x x
cijijijij
=πσ+ −
()
2
⎛⎞
log 2
⎜⎟
⎜⎟
⎝⎠
mn mn
22
∑∑ ∑∑
ij ij
== ==
11 11
()
2
0
()
. (2.66)
Процедура дальнейшего решения задачи аналогична описанной
в пункте 2.5.1.
Приведенная выше методика моделирования перевозочного
процесса в условиях неопределенности является слишком громоздкой
для перевода ее элементов в линейную зависимость.
Рассмотрим другую нелинейную модель [28], в которой предложена двухуровневая система принятия управленческих решении
в интегрированных предприятиях лесной отрасли в условиях неопре-
84

деленности, вызванной влиянием институциональных факторов риска
и неопределенностей. Здесь на первом уровне решаются традиционные задачи финансового анализа и планирования основных экономических показателей предприятия лесной отрасли (табл. 2.1). На втором уровне дается оценка факторам неопределенности и риска
экспертными методами для определения суммарного эффекта от выбора направления деятельности предприятия. В
качестве критерия
в данной модели предлагается показатель валового дохода (GM),
скорректированный на влияние выделенных факторов неопределенности и риска. Сумма (W
+ В + М + G) представляется, как эффект
L
финансовых потерь (приобретений), вследствие влияния выделенных
факторов неопределенности, где:
W
1. Правило использования (недостаток) сырьевых ресурсов (
).
L
2. Высокий уровень теневых экономических отношений в от-
расли (
В).
3. Деятельность монополий, связанных с холдингом (М).
4. Деятельность государственных законотворческих структур (G).
Влияние каждого из факторов на показатель валового дохода
устанавливается экспертными методами.
В табл. 2.1 приведена схема расчетов основных индикаторов
принятия управленческих решений.
Таблица 2.1
Схема расчетов основных индикаторов принятия
управленческих решений
Фактор Институциональные
факторы
Сцена-
рий
развития
1
2
3
Управленческий
WL B M G EC P TC V S T GK EM GM
1 2 3 4 5 6 7 8 9 10 11 12 13
уровень
Прогнозные планово-
экономические факторы
Планово-экономический
уровень
Индика-
торы
Данная таблица заполняется по трем сценариям развития обо-
значенного события:
1. Благоприятный вариант развития события (оптимальный).
2. Текущий вариант.
3. Неблагоприятный вариант.
85

Ценность полученных экспертных оценок зависит от квалификации эксперта и может быть неодинакова, поэтому иногда применяют
взвешивание мнений экспертов. Таким образом, при подведении итога
исследования и разработки прогноза есть выбор, каким образом сформировать прогноз, особенно если его представляет группа экспертов.
При проведении экспертных опросов необходимо определить достоверность прогноза, можно
только полагаясь на компетентность экспертов. Достоверность у экспертных методов не велика, но, в принципе,
можно повысить вероятность правильного прогноза, если использовать
групповые экспертные опросы в несколько туров. Считается, что вероятность правильных оценок тем выше, чем меньше разброс мнений
(что может быть определено путем вычисления среднеквадратичного
отклонения). Однако нередки
ситуации, когда эксперты, дающие более
достоверные прогнозные оценки, оказываются в меньшинстве. Кроме
того, вероятность правильных оценок косвенно определяется исходя из
устойчивости мнений экспертов. Если оценки экспертов резко меняются от одного тура к другому, то вероятность ошибки прогноза велика. Окончательная оценка прогноза выводится как:
− средняя арифметическая простая или взвешенная (если экс-
перты имеют разную квалификацию);
− правило «диктатора», когда один человек будет выражать
мнение всей группы;
− демократичный подход, когда эксперты должны прийти к об-
щему мнению самостоятельно, за один или несколько туров (метод
Дельфи).
Например, группе экспертов (специалистов лесной отрасли)
предлагается ответить на 4 вопроса, касающихся оценки влияния выделенных факторов неопределенности на результаты деятельности
предприятия, осуществляющего поставки лесного сырья за границу
РФ. В результате проводимого опроса составляется
«шкала ответов»
и каждому из факторов присваивается оценочное значение от 0 до 1.
Это значение используется для расчета коэффициента влияния факторов неопределенности на доход интегрированного предприятия. Промежуточным показателем при принятии управленческих решений
предлагается использовать экономическую маржу (ЕМ), представляющую собой значение дохода от реализуемой партии товаров:
ЕМ
= V
i
· (Рi · ЕС
j
– S
– ТСi + ∆Тi – GCi), (2.67)
i
i
где Vi – объем планируемой партии товара в м3; Pj · ЕСi – планируемая выручка от реализации 1 м
3
лесопродукции i-й партии в нацио-
86

нальной валюте (цена 1 м3 лесопродукции i-й партии в иностранной
валюте умноженная на средневзвешенный месячный валютный курс);
S
– себестоимость (трансфертная цена) 1 м3 лесопродукции i-й партии
i
в национальной валюте; ТСi – транспортные и таможенные затраты на
3
1 м
лесопродукции i-й партии в национальной валюте; ∆Тi – разница
между налогами, уплачиваемыми в бюджет и возмещаемыми из бюджета, с 1 м
общехозяйственные и коммерческие расходы, приходящиеся на 1 м
3
лесопродукции i-й партии в национальной валюте; GC
–
i
3
лесопродукции i-й партии в национальной валюте.
Данное соотношение не является определяющим для принятия
управленческого решения в пользу реализации определенного сценария
ЕМi с наибольшим значением, поскольку не учитывает риски,
связанные с выделенными выше институциональными факторами.
Итоговым показателем двухуровневой системы выбора и принятия
управленческих решений является валовой доход, скорректированный на влияние факторов неопределенности и риска, определяемых
экспертным путем. Формула, выражающая данную зависимость следующая:
GM
– (WL + B
i
× (Vi · (Pi · EC
– S
i
+ M
i
– ТСi + ∆Ti – GCi)). (2.68)
i
+ Gi) ×
i
Степень влияния каждого из факторов на функционал отражает
весовой коэффициент g, определяемый группой экспертов отрасли.
В соответствии с влиянием факторов на показатель GM и значением
показателя экономической маржи, рассчитывается суммарный показатель
GM по каждой из партий продукции. После этого делается
выбор в пользу реализации того экономического решения, где суммарный валовой доход, скорректированный на влияние факторов
неопределенности и риска, принимает максимальное значение [28].
Проведенный анализ методов решения реальных задач транспортной логистики выявил ряд проблем. Первая причина этих проблем связана с необходимостью
учета неопределенности, возникающей при оценивании параметров математических моделей задач
транспортной логистики. Это обстоятельство порождает постановку
ряда задач по следующим направлениям:
− транспортные задачи со стохастическим спросом и предложе-
нием (необходим корректный учет случайного характера спроса и
производства, который позволил бы получить решение, минимизирующее суммарную стоимость транспортировки с учетом затрат на
хранение непроданной части товара и потерь от дефицита);
87

− транспортная задача с неопределенной стоимостью перевозки
A
и неопределенным объемом поставки (решение данной задачи приводит к специфической нетривиальной задаче математического программирования).
Другая группа проблем возникает в связи с возможной неадекватностью линейных моделей транспортной логистики. Следовательно, в экономико-математическую модель должны включаться только
основные факторы. Проверку качества моделей необходимо
производить реальной экономической деятельностью предприятия. Практика
применения моделей показывает, что сложные модели с множеством
параметров оказываются зачастую непригодными для практического
использования. Учитывая вышеизложенное, в целях проверки качества
моделей необходимо производить оценку устойчивости их решений.
2.5.3. Оценка устойчивости решений
транспортных задач
Устойчивость решений оптимизационных задач – это свойство
полученного решения не изменяться в условиях, когда исходные данные задачи изменяются в некотором диапазоне возможных своих значений. Принято, что чем больше этот диапазон, тем выше устойчивость решения. Применительно к транспортным задачам линейного
программирования проблема оценки устойчивости решения задачи
формулируется следующим образом. Пусть
продукт перевозится от
m
поставщиков с объемами поставляемого
некоторый однородный
продукта
12
торых равны
, ,...,
AA
, ,...,
BB B
12
к совокупности потребителей, потребности ко-
m
. Стоимости транспортировок от конкретных
n
поставщиков к конкретным потребителям – случайные величины,
с известными плотностями распределения. Оценка устойчивости
решения транспортной задачи допускает простое аналитическое
решение. Будем считать, что оптимальное решение транспортной задачи искать с использованием метода минимального элемента матрицы. Этот метод конструктивно устроен таким образом, что при его
применении может
быть однозначно получена последовательность
назначений, выстроенная в соответствии с возрастающей последовательностью стоимостей перевозок. Эта последовательность на первом
шаге начинается с процедуры отыскания минимального элемента
матрицы стоимостей и на каждом последующем шаге продолжается
этой же процедурой для матрицы, скорректированной после выполнения очередного шага. В результате получим фиксированную
88

возрастающую последовательность стоимостей, определяющую по-
)
+
+
)
ϕ
+
P
x
+
=
рядок отыскания значений. Перенумеруем эту последовательность:
CC C
<<<
12 1
...
вильного порядка назначений для произвольной пары
стоящих членов последовательности
принимает случайное значение Х. Вычислим вероятность того,
C
1k
+
что случайное число
. Рассчитаем теперь вероятность реализации пра-
mn
++
,1kk+
(
, ,...,
CC C
12 1
примет значение, меньшее Х и затем проин-
C
k
. Пусть число
mn
+
рядом
тегрируем эту вероятность совместно с плотностью распределения
вероятностей случайной величины
Сϕ
и
(
k
()
случайных величин
1kС+
. Плотности распределения
C
1k
C
k
и
естественно считать
C
1k
независимыми. Поэтому искомая вероятность рассчитывается по формуле
∞
=ϕ ϕ
+++
kk k k k k
,1 1 1
CdC CdC
() ()
∫∫
00
X
.
Например:
1
⎧
⎪
ba
С
ϕ=
()
k
⎧
⎪
С
ϕ=
()
1
k
+
ba
⎨
⎪
0, , ;
⎩
−
kk
⎨
⎪
0, , ;
⎩
1
−
11
kk
++
Пусть случайная величина
,,,
,,,
C
k
∈
Cab
Cab
[]
kkk
Cab
∉
[]
kkk
∈
[]
111
kkk
+++
Cab
∉
[]
111
kkk
+++
принимает значения
ab∈
[,]
kk
Тогда
b
1
k
+
PC x C dC
() ()
есть вероятность того, что
>= ϕ
1111
kkkk
++++
Cx
∫
max ,
xa
{}
1
k
+
. Теперь вероятность реализации
>
1k
+
правильной последовательности выбора
b
k
PPCxCdx
CC k k k
kk
=>ϕ
,1
+
1
()()
∫
a
k
+
,
CC
kk
будет равна
1
. (2.69)
Найдем эту вероятность.
.
Если
ab
1kk
+
, то
>
P
CC
kk
,
+
.
1
1
89

Путь теперь
P
+
P
(
)
+
a
kk
+
=>ϕ+>ϕ=+
,1 1 12
CC kk kk
kk
+
JPCxxdx dx
∫∫
1
aa
kk
aa
JPCxxdx dx
=>ϕ= =
11
bbb
kkk
=>ϕ= =
21
∫∫∫
aax
++
kk
ab
1kk
+
1
PC x xdx PC x xdx J J
()
++
11
++
kk
()
+
∫∫
aa
kk
11
kk
()
kk
+
. При этом:
<
()
()
()
a
b
()
1
+
1
ba ba
−−
kk kk
1
−−
ba b a
kk k k
()
aa
−
1
+
kk
+
1
dC
+
k
11
++
,
1
,
1
==
bab a
−−
()( )
kkk k
bbb
kkk
=−=− =
ba bb
()()
++ +
kk kk
==
11 1
bab a bab a
22
()( )
b x dx udu
()
+
1
k
∫∫
aba
+++
111
kkk
−−−
kkk k kkk k
22
−− −−
++ ++
11 11
−=
bxu
+
1
k
−=
dx du
bbaabbbb
bab a
()( )
2222
++++++
kkkkkkkk
11
++
1
−−
kkk k
−+−+−
22
111111
11
++
−
+
1
−
.
()( )
Далее:
−−+−
aabbbaab
11111
+++++
kkkkkkkk
P
CC
kk
ba ba a aa bb ba a b
222222
++ + + + + ++ +
11 1 1 1 1 11 1
kk kk k kk kk kk k k
==
−−++− +−
=+ =
,
+
1
−−−
ba bab a
kk kkk k
222
2
aa bb ba a b
222
++++
1111
kk kk kk k k
=
2
22
2
()( )
bab a
−−
()( )
kkk k
++
11
++
+−−−
bab a
−−
()( )
kkk k
++
11
22
11
22
(2.70)
.
При этом вероятность реализации правильного порядка назна-
чений в алгоритме минимального элемента матрицы равна:
22
mn mn
+− +−
==
PP
∏∏
11 11
kkkkkk
== ++
,
CC
kk
1
+
2
1
mn
+−
2
Вычисленное значение
aa bb ba a b
222
1111
kk kk kk k k
++++
характеризует уровень устойчивости
−−−
−−
bab a
()( )
22
. (2.71)
решения транспортной задачи. При этом выбранное решение будет
абсолютно устойчивым, если для всех пар
90
,1kk
неравенство (3.69)
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]
