Файл: Пивоваров, С. Э. Моделирование процессов прогнозирования в приборостроении.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 23.10.2024
Просмотров: 82
Скачиваний: 0
Г л а в а 5
РАЗРАБОТКА МЕТОДОВ И АЛГОРИТМОВ
РЕШЕНИЯ ЗАДАЧИ ОПТИМИЗАЦИИ ПРОГНОЗА
РАЗВИТИЯ ОТРАСЛИ
5.1. Вывод математической модели выбора оптимального варианта прогноза
развития отрасли
Задачу выбора (определения) оптимального варианта прогноза развития отрасли будем формулировать с помощью модели, описан ной в п. 4.2, несколько упростив ее и приняв все коэффициенты мо дели за функции времени. Причем время характеризует весь период прогнозирования. Таким образом модель примет следующий вид: требуется найти минимум функционала
т р т п ft
2 ?и„т. г (/) X? + |
2 I ] 2 Z j r s X j r s (0 (5.1.1) |
г= 1 q—\ |
г = 1 s = 1 / = 1 |
при условиях:
i j |
( 0 - i ] |
$ ( * ) * ? < 0; |
(5Л .2) |
S = 1 |
(7 = |
1 |
|
т |
|
|
(5.1.3) |
Z |
X jrs(t) = PJS(ty, |
||
г = |
1 |
|
|
I ; 2 |
bJrl(f)XJr(t)^Q,(f); |
(5.1.4) |
1/—1
тр
2 |
2 |
(/) X? < |
Ц (0 + Q* {*); |
(5.1.5) |
г = 1 ? = 1 |
|
|
|
|
2 |
2 ^ |
(0 X? ^ |
(0 + Qx (0; |
(5.1.6) |
Г = 1 |
( 7 = 1 |
|
|
|
m |
р |
|
|
|
2 |
21 tf(0 * !< L (0 + Q (О; |
(5.1.7) |
%г = 1 <7 = 1
|
ХуЛО, |
Ху„ (0 . |
Х чг ^ |
0. |
(5.1.8) |
Обозначим |
переменные |
модели |
(5.1.1) —(5.1.8) соответственно |
||
X qr, (л = 1, т; |
9 = 1 ,р)\ Xj,s, (/'= 1 , h\ |
г = 1 , т\ |
s = 1, п); Ху,, |
4 С, Э, Пивоваров |
97 |
Обозначим коэффициенты модели (5.1.1) — (5.1.8) соответственно
Обозначим |
индексы |
hm + hn + I g + Л + 1, |
rnp -f- hmn + |
+ hm, mp + |
hmn через |
M, N и N1 соответственно. |
Тогда задачу |
(5.1.1) — (5.1.8) сформулируем следующим образом: требуется обра тить в максимум линейную форму
N |
|
|
|
|
Z C j ( t ) X j ( t ) |
|
(5.1.9) |
||
при „условиях: |
|
|
|
|
N |
|
|
|
|
£ а/у (/)*,(/) |
МО. |
(« = |
1, 2, .... М); |
(5.1.10) |
Xy(/)2 s 0, |
(/= 1 , |
2, |
.... N). |
(5.1.11) |
Предположив, что фактор времени выступает как параметр коэф фициентов линейной формы, составляющих векторов условий и век торов ограничений и линейно зависит от них, задачу (5.1.9) — (5.1.11) сформулируем как общую задачу линейного параметриче ского программирования.
Требуется обратить в максимум линейную форму
|
N |
|
|
|
Ц (C'i + tel) Xj |
(5.1.12) |
|
при условиях: |
|
|
|
|
N |
(i = 1, М); |
|
|
^ ( a 'u + t a V X j ^ b i + tb!. |
(5.1.13) |
|
<• |
X j ^ O , (/ = 1, |
N). |
(5.1.14j |
98
5.2. Анализ методов решения задачи оптимизации
Прежде, чем начать анализировать метод решения общей задачи линейного параметрического программирования, сформулирован ной в п. 5.1, разберем ряд более простых задач параметрического программирования [12], [67, 68], которые будут способствовать вы воду алгоритма общей задачи.
1-я частная задача линейного параметрического программиро
вания. Требуется для каждого |
t из некоторого множества опре |
||||||
делить вектор Xt, обращающий в максимум линейную форму: |
|||||||
L = |
2 |
q x j = |
2 |
(С/ + tcj) Xj |
(5.2.1) |
||
|
/=1 |
/=1 |
|
|
|
||
при условиях: |
|
|
|
|
|
___ |
|
N |
|
|
|
|
|
|
|
^ |
a |
ijXj = bi, |
0 |
= |
1, М); |
(5.2.2) |
|
/=1 |
|
|
|
|
|
|
|
|
X j ^ O , |
(j = |
T7JV). |
(5.2.3). |
|||
Набор переменных |
X j (j = |
1, |
N), |
характеризующих |
решение |
||
задачи (5.2.1) — (5.2.3), |
будем называть планом. Базис |
опорного |
плана данной задачи является оптимальным для некоторого t, если оценки относительно этого базиса всех векторов условий, вычислен ные при данном t, неотрицательны. Полная совокупность значений параметра t, при которых базис оптимален, называется множеством оптимальности этого базиса.
Проанализируем изменение решения задачи (5.2.1) — (5.2.3) с изменением параметра t. Пользуясь любым алгоритмом метода последовательного улучшения плана (метод подробно будет описан ниже), решим задачу для некоторого t = t0. После анализа конеч ного числа шагов получим оптимальный план задачи (случай 1°) или убедимся в неограниченности линейной формы (5.2.1) на множестве планов задачи (случай 2°). Случай несовместности условий (5.2.2) — (5 .2 .3) не представляет интереса, так как при этом задача для любого t не имеет решения. Рассмотрим первые два случая отдельно.
С л у ч а й 1°. Вычислим оценки Дj (t) векторов Aj относительно найденного оптимального (для t = t0) базиса A h, A iit ..., A iM задачи
м
Д ; « ) = Е С‘ ( 0 * у - С у (* ) =
S = 1 |
М |
М |
|
= X c i X i . - q + t |
Z c l x i j - Q |
s = l |
s= 1 |
4* |
99 |
Обозначим:
м
м
Тогда
Aj (t) = A(l -f- 1А1г.
Так как исследуется оптимальный для i = t0 базис, то А Э * О, (/ = 1, N), т. е. система неравенств
Au + t A , , ^ 0 (/= 1 , N) |
(5.2.4) |
совместна. Для всех А/, < 0 неравенства системы могут быть запи саны в виде
адля всех Ду, > 0 — в виде
Обозначим:
(5.2.5)
Рассматриваемый базис Ailt Л/2, |
А{м по определению явля |
ется оптимальным только для тех значений t, которые удовлет воряют системе неравенств (5.2.4). Множество оптимальности дан ного базиса состоит из всех значений t, для которых справедливо
условие / «s t ^ /.
Таким образом, числа / и 1 полностью характеризуют множество
оптимальности Е (t, t) данного базиса:
а) если t_> то Е (t, t) — пустое множество;
б) если |
t = |
t, то |
Е (J, 1) — точка |
t = t = ~t\ |
|
|
|
в) если |
—o o < t < t < c o , |
то Е (i, |
^ — отрезок |
[/, |
Г]; |
||
г) если |
t = |
— оо, |
I < оо, |
то Е (t, |
t ) — луч (— оо, |
<]; |
100
д) |
если |
— о о < t, |
7 = оо, |
то |
E (t, |
7) — луч |
[*, оо); |
е) |
если |
t = — оо, |
/ = оо, |
то |
Е (t, |
7) —вся |
прямая. |
Исследуем далее задачу (5.2.1) — (5.2.3) для значений t > t .
При этом будем считать 7 •< оо, т. е. среди чисел А/, имеются отри цательные. Пусть k один из индексов, определяемых соотношением
|
|
7 = min |
|
|
ДА, < 0 . |
(5.2.6) |
|
|
|
Л/,<0\ Д/,/ |
|
|
|
||
Очевидно, Аа(7) — 0. |
Тогда возникают следующие |
две возмож |
|||||
ности: |
все компоненты xik разложения вектора условий Ак по векто |
||||||
1) |
|||||||
рам базиса неположительны. При этом условии линейная форма |
|||||||
задачи |
для / > |
7 не |
ограничена в. области своего |
определения. |
|||
Действительно, |
в этом случае |
|
|
|
|||
|
Да (t) = Да (7) + (t - |
7) Да, = (t - 7) Да, < 0. |
|||||
Это неравенство вместе с условием xik |
0 для всех i |
представляет |
|||||
собой признак неразрешимости задачи для метода последовательного |
|||||||
улучшения плана; |
|
|
|
|
|
||
2) |
по крайней мере одна из составляющих xik положительна. |
||||||
Продолжая решение задачи методом последовательного улучшения |
|||||||
плана, |
вводим в базис вектор А к, для которого ДА(7) |
= 0 и Да, < 0 , |
|||||
и исключаем из базиса вектор Аг, для |
которого |
|
|||||
|
|
|
х го |
min |
^ |
X r k > 0. |
|
|
|
|
x r k |
|
|||
|
|
|
* i k > |
° X i k |
|
|
|
|
|
|
|
|
|
Новый базис является оптимальным по крайней мере для t = t. Отсюда по формулам, связывающим параметры двух последователь ных итераций метода последовательного улучшения плана,
Д/(0 = Ду( 0 - ^ А а(7), |
0 = Т7а0 |
|
||
имеем |
|
(/ = 1, ДО. |
|
|
Д / ( 0 = Д у (7 )2 г 0 , |
|
|||
Следовательно, при t = t |
неравенства |
|
||
Д/. + Щ, |
0, |
0 = |
П~Д7) |
(5.2.7) |
совместны. |
|
|
|
|
Неравенства (5.2.7) могут быть совместны и при t > |
t, однако |
любое i < 7 не удовлетворяет системе (5.2.7).
101