Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Применение логистико-математических моделей в лесной отрасли. Монография

.pdf
Скачиваний:
0
Добавлен:
07.09.2026
Размер:
2 Мб
Скачать
☆
∞
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)
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]