Файл: Пивоваров, С. Э. Моделирование процессов прогнозирования в приборостроении.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