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

Шурко Введение в системный анализ

.pdf
Скачиваний:
18
Добавлен:
21.02.2016
Размер:
956.66 Кб
Скачать

1

ДОНБАССКАЯ ГОСУДАРСТВЕННАЯ АКАДЕМИЯ СТРОИТЕЛЬСТВА

И АРХИТЕКТУРЫ

КАФЕДРА “ВЫСШАЯ МАТЕМАТИКА И ЭКОНОМЕТРИЯ”

КОНСПЕКТ ЛЕКЦИЙ по

ВВЕДЕНИЮ В СИСТЕМНЫЙ АНАЛИЗ

для инженерных специальностей

дневного отделения

доц. Шурко Г.К.

УТВЕРЖДЕНО НА ЗАСЕДАНИИ КАФЕДРЫ ВЫСШЕЙ МАТЕМАТИКИ И ЭКОНОМЕТРИИ ПРОТОКОЛ № ___________________

УТВЕРЖДЕНО НА ЗАСЕДАНИИ СОВЕТА ФАКУЛЬТЕТА ФОП ПРОТОКОЛ № ___________________

г. Макеевка 2001 г.

2

СОСТАВЛЕНИЕ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ ОПТИМИЗАЦИОННЫХ ИНЖЕНЕРНОЭКОНОМИЧЕСКИХ ЗАДАЧ

Под математической моделью будем понимать систему (совокупность уравнений, неравенств, функциональных зависимостей), описывающие реальные инженерно-экономические соотношения.

Алгоритм составления математических моделей оптимизационных инженерно-экономических задач опишем на примере задачи.

Задача 1.

Предприятие изготовляет продукцию 2 наименований П1 и П2, используя при этом 4 ресурсаР1, Р2, Р3, Р4.

Норма затраты, каждого из ресурса и прибыль от реализации выпускаемой продукции приведены в таблице №1.

 

Норма затрат на производство ед.

 

Ресурс

продук. вида

 

Запас ресурса

 

 

 

 

 

 

 

 

 

 

П1

 

П2

 

Р1

10

 

12

60

Р2

8

 

10

40

Р3

9

 

0

81

Р4

0

 

8

56

Прибыль от

4

 

5

 

реализации

 

 

 

 

Учитывая наличие ресурсов, нормы затраты и прибыль составить такой план реализации выпуска продукции П1 и П2, которые обеспечил бы наибольшую суммарную прибыль.

Алгоритм.

1 шаг. Введение переменных величин.

Вводим переменные величины количество и смысл (содержание) которых зависит от содержания задачи и того, что требуется определить в ней.

В нашей задаче введем 2 переменные величины Х1 и Х2. Х1- количество единиц продукции вида П1.

Х2- количество единиц продукции вида П2.

2 шаг. Получение системы ограничений.

3

Исходя из введенных на первом шаге переменных величин, анализируя условие задачи, получаем систему ограничений (неравенств, уравнений).

Х1 0,

 

 

Х2 0 т. к. Х1 и Х2

количество единицы продукции.

10

Х1

- затраты Р1 на производство Х1 единицы продукции П1 (аналогично

12 Х2 )

 

 

 

 

10

Х1

+ 12 Х2 60

 

10

Х1

+ 12 Х2 -суммарные затраты Р1 на производство Х1 и Х2.

 

Х1 0,

 

 

Х2 0

 

10

Х1

+ 12 Х2 60

 

8

Х1 + 10 Х2 40

(1)

9

Х1

81

 

8

Х2

56

 

3 шаг. Получение целевой функции.

Под целевой функцией понимается математическое выражение того, что нужно оптимизировать (максимизировать) в инженерно-экономической задаче.

4Х1 -это прибыль от реализации Х1 единиц продукции П1.

5Х2 - это прибыль от реализации Х2 единиц продукции П2.

L = 4 Х1 + 5 Х2 max (2)

- суммарная прибыль от реализации продукта Х1 и Х2.

Полученная математическая модель, т.е. система ограничений (1) и целевая функция (2) представляет собой задачу линейного программирования (планирования) З.Л.П.

Графический способ решения задач

линейного программирования

Графический способ решения задач линейного программирования применяется тогда, когда число переменных входящих ЗЛП не превышает 3.

Наиболее удобен графический способ, когда число переменных равно 2.

Прежде перейдем к описанию алгоритма графического способа решения определим некоторые понятия.

4

Допустимым планом (решением) ЗЛП называется множество значений переменных входящих в ЗЛП и удовлетворяющих системе ограничений.

Оптимальным планом (решением) ЗЛП называется такой допустимый план при котором достигается оптимальное значение целевой функции.

Алгоритм графического способа.

Алгоритм опишем на примере задачи линейного программирования (1) и

(2).

Х1 0,

 

Х2 0

 

1+7х2 21

 

10Х1 + 8 Х2 40

(1)

L = - Х1 - Х2 min

(2)

1 шаг. Получение выпуклого многоугольника допустимых решений.

Выпуклый многоугольник допустимых решений получается как пересечение множеств, каждое из которых описуется неравенством из системы ограничений.

Для построения многоугольника введем прямоугольную систему координат с осями Х1 и Х2.

2 шаг. Строим область допустимых решений для системы

ограничений.

Х1 0,

Х2 0

1 + 7х2 21

10Х1 + 8 Х2 40

5

3шаг. Строим вектор нормали n, выходящий из начала координат

иимеющий координаты равные числовым коэффициентам при соответственных переменных в выражении для целевой функции.

4шаг. Перпендикулярно направлению вектора n проводим прямую

l, так чтобы полученный многоугольник допустимых решений оказался по ту сторону от прямой l, куда указывает вектор n. Прямая l называется опорной прямой.

5 шаг. Перемещаем прямую l вдоль направления вектора n

параллельно самой себе. Тогда в первой точке принадлежащей многоугольнику допустимых решений, которую встретит прямая l

достигается min целевой функции.

В последней же точке - max.

Если при решении задачи одна из сторон многоугольника допустимых решений окажется параллельной прямой l, то max (min) достигается в любой точке этой стороны.

Проверка.

Решим систему уравнений:

1 + 7х2 = 21

10Х1 + 8 Х2 = 40

х1=(21-7*х2)/3=2,43 10*(21-7*х2)/3+8*x2=40 70-23,3*х2+8*х2=40 -15,3*х2=-30

х2=1,96

6

СИМПЛЕКС-МЕТОД

("МЕТОД ПОСЛЕДОВАТЕЛЬНОГО УЛУЧШЕНИЯ ПЛАНА")

1. Понятие о симплекс-методе.

Впервые был разработан Данцигом в 1947г. Позволяет переходить от одного

допустимого базисного решения к другому, причем так, что значение целевой

функции возрастают. В результате оптимальное решения находят за конечное

число шагов.

Имеем О.З.Л.П.:

 

n

 

L=Sci*xi + c0 max

 

i=1

 

при выполнении условий:

 

a11x1 + a12x2 + a13x3 + ... +a1nxn = b1

 

a21x1 + a22x2 + a23x3 + ... +a2nxn = b2

(1)

......................................................

 

am1x1 + am2x2 + am3x3 + ... +amnxn = bm

 

Или, в векторной форме:

L=( c * x ) + c0 max

При выполнении условий:

A * X = B,

где:

c = (c1,c2,...,cn) ,

x

 

 

a

 

1

 

11

x2

 

a21

.

 

.

x =

 

,

A =

.

 

.

.

 

.

 

 

 

 

xn

 

am1

7

a

. . .

a

 

b

 

12

. . .

1n

 

1

 

a22

a2n

b2

.

. . .

.

 

.

 

.

. . .

.

,

B =

 

 

.

 

.

. . .

.

 

.

 

am2

 

 

 

 

 

. . .

amn

bm

I шаг.

Приведение З.Л.П. к каноническому виду, т.е. преобразование неравенств в равенства путем введения новых неотрицательных переменных Xn+1, …, Xn+m,

которые называются балансовыми (выравнивающими).

Тогда вместо системы неравенств (1) получим систему уравнений (2):

a11x1 + a12x2 + a13x3 + ... +a1nxn = b1 a21x1 + a22x2 + a23x3 + ... +a2nxn = b2 (2)

....................................................

am1x1 + am2x2 + am3x3 + ... +amnxn = bm

Матрица коэффициентов Асистемы (2) имеем вид:

 

a

a

. . .

 

11

12

 

A =

a21

a22 . . .

.

. . . .

 

 

am2 . . .

 

am1

a

1

0

. . . 0

.

0

1n

0

1

0 . . .0

 

 

a2n

. 0

.

. .

. . . .

.

. .

amn 0

0

. . . 0

 

 

. 1

Поэтому задачу Л.П. можно сформулировать в каноническом виде:

максимизировать функцию

n

L(x)=S ci*xi + c0 (3)

i=1

при ограничениях:

А * х = В, х 0

(4)

8

Условие (4) аналогично условию:

n

S аi*xк = bi , i=1,2,...,m; n m.

к=1

хк 0 , к=1,2,...,n.

При этом можем считать, что bi 0, i=1,2,...,m.

Предположим, что ранг матрицы А равен r; r m. Выразим переменные х1, х2, ... , хr (которые называются базисными переменными) через остальные

переменные (которые называются свободными):

x1 = a’1r+1xr+1 + a’1r+2xr+2 + a’1r+3xr+3 + ...

+a’1nxn +b’1

 

x2 = a’2r+1xr+1 + a’2r+2xr+2 + a’2r+3xr+3 + ... +a’2nxn +b’2

(5)

.....................................................

 

xr = a’n+1xr+1 + a’n+2xr+2 + a’n+3xr+3 + ...

+a’mxn +b’r

 

Система уравнений (5) называется системой, приведенной к единичному базису, а

набор (х1, х2,…, хк) – базисом. Подставляя в функцию L вместо базисных переменных их выражения через свободные из (5), получим:

n

L = S ixi + 0

= 0 + r+1xr+1 + … + nxn

(6)

 

i=1

 

 

Положим все свободные переменные в (5)

равными нулю (x r+1 =0,…, x n

=0), найдем значения базисных переменных.

 

 

X1=b`1, ... , Xr=b`r.

 

 

Таким образом, решение

системы (b`1, b`2, ...,

b`r,0, ...., 0) будет

допустимым. Оно называется базисным.

Для полученного базисного решения значения целевой функции будет равно:

LБ = 0.

Решение задачи при помощи симплекс–метода распадается на ряд шагов,

заключающихся в том, что от данного базиса Б мы переходим к другому базису Б` с таким расчетом, чтобы

LБ LБ`.

9

2. Симплексные таблицы.

В виде таблицы эти данные можно представить так:

БАЗИСЫЕ СВОБОД-

ПЕРЕМЕН-

НЫЕ

X1 ...

XI XR

XR+1

...

XJ

...

...

XN

НЫЕ

ЧЛЕНЫ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Х1

B1

1 ...

0

0

A1,R+1

...

A1,J

...

...

A1,N

...

...

... ... ... ...

...

... ... ... ... ...

 

 

 

 

 

 

 

 

 

 

 

ХI

BI

0 ...

1

0

AI,R+1

...

AI,J

...

...

AIN

...

...

... ... ... ...

...

... ... ... ... ...

 

 

 

 

 

 

 

 

 

 

 

ХR

BR

0 ...

0

1

AR,R+1

...

AR,J

...

...

AR,N

L

0

0 ...

0

0

R+1

...

J

...

...

N

Равенство (6) будем называть приведенным (к свободным переменным)

выражением

для L,

а коэффициенты

о

оценками

(индексами)

соответствующих свободных переменных Xj.

АЛГОРИТМ СИМПЛЕКС-МЕТОДА.

1 шаг. Выбирается р-ый разрешающий столбец из условия: оценка

p <0 и хотя бы один элемент aip >0.

2 шаг. Выбирается q-я разрешающая строка из условия:

 

bq

 

 

 

 

 

bi

 

 

 

 

 

min

 

 

для

0

 

 

 

 

aqp

 

 

 

aip

 

aip

 

 

3 ШАГ. Производится

перерасчет элементов разрешающей q -й

строки по формуле.

10

a aqk (k 1,....., n)

qk aqk

4 шаг. Вычисляются элементы всех остальных строк (при к p) по формуле:

a`ik = aik - a`qk * aip ( i=1,...,q-1,q+1,...,n ).

Вычитаем из каждой строки ведущую строку, умноженную на такой коэффициент, чтобы в ведущем столбце получились нули.

5 шаг. В базис вводится переменная хp и выводится переменная xq.

Следует иметь в виду основную теорему симплекс-метода:

Теорема. Если после выполнения очередной операции;

1)найдется хотя бы одна отрицательная оценка и в каждом столбце с такой оценкой будет хотя бы один положительный элемент, т.е. k<0 для некоторого k и aik >0, для тех же к и некоторого i, то можно улучшить решение, выполнив следующую операцию;

2)найдется хотя бы одна отрицательная оценка, столбец которой не содержит положительных элементов, т.е.

k < 0,

aik < 0.

для какого-то k и всех i, то функция L неограничена в области допустимых решений;

(Lmax )

3) все оценки окажутся неотрицательными, т.е. k

0 для всех к, то

достигнуто оптимальное решение.

 

 

Пример 1.

2x1+3x2 19

 

2x1+ x2 13

3x2 15

3x1 18