Контрольная работа: Экономико–математические методы в управлении

Название: Экономико–математические методы в управлении
Раздел: Рефераты по экономико-математическому моделированию
Тип: контрольная работа

ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ОБРАЗОВАНИЮ

РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

ИМЕНИ ИММАНУИЛА КАНТА

кафедра экономики

КОНТРОЛЬНАЯ РАБОТА

по дисциплине «Экономико – математические методы в управлении»

вариант №30

КАЛИНИНГРАД

2008


Задание

Задание 1.2.

Смесь можно составить из n продуктов Сj (j=1,n). В каждом из продуктов содержится m компонентов. Минимально допустимый объем содержания i-го компонента в смеси выражается величиной bi (i=1,3). Содержание i-го компонента в единице j-го продукта выражается величиной аij . Цена единицы j-го продукта равна сj . Составить смесь, минимальную по стоимости, выбрав для решения данной задачи наиболее рациональный способ.

C1

C2

C3

bi

cj

9

6

7

a1j

7

5

8

70

a2j

8

2

3

40

a3j

9

6

7

50

Задание 2.2.

Найти графоаналитическим методом оптимальное решение задачи нелинейного программирования.

maxZ = 3.6x1 – 0.2x1 2 + 0.8x2 – 0.2x2 2

2x1 + x2 ≥ 10

x1 2 -10x1 + x2 ≤ 75

x2 ≥ 0

Задание 3.1.

После нескольких лет эксплуатации оборудование может оказаться в одном из трех состояний:

1) требуется профилактический ремонт;

2) требуется замена отдельных деталей и узлов;

3) требуется капитальный ремонт.

В зависимости от ситуации руководство предприятия может принять следующие решения:

1) отремонтировать оборудование своими силами, что потребует затрат а ;

2) вызвать специальную бригаду ремонтников, расходы в этом случае составят b ;

3) заменить оборудование новым, реализовав устаревшее по остаточной стоимости.. Совокупные затраты на это мероприятие составят с .

Требуется найти оптимально решение данной проблемы по критерию минимизации затрат с учетом следующих предположений:

а) на основе обобщения опыта эксплуатации аналогичного оборудования определены вероятности наступления соответствующих состояний – q ;

б) имеющийся опыт свидетельствует о равной вероятности наступления соответствующих состояний;

в) о вероятностях наступления соответствующих состояний ничего определенного сказать нельзя.

П1

П2

П3

a

13

9

15

b

20

12

11

c

18

10

14

q

0.3

0.45

0.25

λ = 0.7

Задание 1.2.

Смесь можно составить из n продуктов Сj (j=1,n). В каждом из продуктов содержится m компонентов. Минимально допустимый объем содержания i-го компонента в смеси выражается величиной bi (i=1,3). Содержание i-го компонента в единице j-го продукта выражается величиной аij . Цена единицы j-го продукта равна сj . Составить смесь, минимальную по стоимости, выбрав для решения данной задачи наиболее рациональный способ.

C1

C2

C3

bi

cj

9

6

7

a1j

7

5

8

70

a2j

8

2

3

40

a3j

9

6

7

50

Смесь, минимальная по стоимости:

7x1 + 5x2 + 8x3 ≥ 70

8x1 + 2x2 + 3x3 ≥ 40

9x1 + 6x2 + 7x3 ≥ 50

x1 ≥ 0; x2 ≥ 0; x3 ≥ 0

F = 9x1 + 6x2 + 7x3 → min

После транспонирования матрицы элементов aij , cсимметричная двойственная задача будет иметь вид:

S(y1 ,y2 ,y3 ) = 70y1 + 40y2 + 50y3 → max , при ограничениях:

7y1 + 8y2 + 9y3 ≥ 9

5y1 + 2y2 + 6y3 ≥ 6

8y1 + 3y2 + 7y3 ≥ 7

y1 ≥ 0; y2 ≥ 0; y3 ≥ 0

Для решения двойственной задачи линейного программирования симплекс – методом, приведём систему неравенств к виду системы уравнений:

7y1 + 8y2 + 9y3 + y4 ≥ 9

5y1 + 2y2 + 6y3 + y5 ≥ 6

8y1 + 3y2 + 7y3 + y6 ≥ 7

y1 ≥0;y2 ≥0;y3 ≥0;y1 ≥0;y2 ≥0;y3 ≥0

S(y1 ,y2 ,y3 ) = 70y1 + 40y2 + 50y3 → max

По правилу соответствия переменных, базисным переменным прямой задачи соответствуют свободные переменные двойственной задачи:

x1 x2 x3 x4 x5 x6

y1 y2 y3 y4 y5 y6

Первая симплексная таблица:

Базис

Сб

А0

y1

70

y2

40

y3

50

y4

0

y5

0

y6

0

y4

0

9

7

8

9

1

0

0

y5

0

6

5

2

6

0

1

0

y6

0

7

8

3

7

0

0

1

0

-70

-40

-50

0

0

0

Вторая симплексная таблица:

Базис

Сб

А0

y1

70

y2

40

y3

50

y4

0

y5

0

y6

0

y4

0

23/8

0

43/8

23/8

1

0

-7/8

y5

0

13/8

0

1/8

13/8

0

1

-5/8

y1

70

7/8

1

3/8

7/8

0

0

1/8

245/4

0

-55/4

45/4

0

0

35/4

Третья симплексная таблица:

Базис

Сб

А0

y1

70

y2

40

y3

50

y4

0

y5

0

y6

0

Y2

40

23/43

0

1

23/43

8/43

0

-7/43

y5

0

67/43

0

0

67/43

-1/43

1

-26/43

y1

70

29/43

1

0

29/43

-3/43

0

8/43

2950/43

0

0

800/43

110/43

0

280/43

В последней таблице в строке Δ нет отрицательных элементов. В соответствии с критерием оптимальности точка максимума Smax = 2950/43 достигнута при значениях: y1 = 29/43; y2 = 23/43; y3 = 0.

По теореме двойственности: Fmin = Smax = 2950/43.

На основании правила соответствия между переменными, оптимальное решение прямой задачи:

y4 x1 = 110/43 y5 x2 = 0 y6 x3 = 280/43

Ответ: В смесь минимальной стоимости 2950/43 целесообразно включить 110/43 единиц продукта C1 , 280/43 единиц продукта C3 , а продукт C2 не включать.

Задание 2.2.

Найти графоаналитическим методом оптимальное решение задачи нелинейного программирования.

maxZ = 3.6x1 – 0.2x1 2 + 0.8x2 – 0.2x2 2

2x1 + x2 ≥ 10

x1 2 -10x1 + x2 ≤ 75

x2 ≥ 0

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

Сначала необходимо преобразовать формулу целевой функции так, чтобы получить её графическое отображение. Воспользуемся методом выделения полного квадрата двучлена относительно x1 и x2 , разделив левую и правую части формулы на -0.2:

-5Z = x1 2 -18x1 + x2 2 – 4x2

Добавим к левой и правой частям уравнения числа, необходимые для выделения полных квадратов двучлена в правой части выражения:

92 и 22 в сумме составляют 85:

85 – 5Z = (x1 – 9)2 + (x2 – 2)2

В результате получилась формула, позволяющая графически изобразить целевую функцию в виде линии уровня на плоскости X1 OX2 . Данные линии уровня представляют собой окружности с общим центром в точке O (9;2). Данная точка является точкой абсолютного экстремума целевой функции.

Для определения характера экстремума нужно провести анализ целевой функции на выпуклость/вогнутость. Для этого необходимо определить вторые частные производные и составить из них матрицу:


Z x1x1 Z x1x2 = -0.4 0

Z x2x1 Z x2x2 0 -0.4

Определим знаки главных миноров данной матрицы.

Главный минор первого порядка -0.4 < 0.

Главный минор второго порядка 0.16 > 0.

Т.к. знаки миноров чередуются, функция Z - строго вогнута. Экстремум вогнутых функций – max, следовательно в точке О у целевой функции находится абсолютный максимум.

Для построения области допустимых значений преобразуем второе неравенство системы ограничений:

x1 2 – 10x1 + x2 ≤ 75

x1 2 – 10x1 + 25 + x2 ≤ 100

(x1 – 5)2 + x2 ≤ 100

(x1 – 5)2 ≤ 100 – x2

Уравнение (x1 – 5)2 = 100 – x2 выразим через переменные x1 * и x2 * :

x1 * = x1 – 5

x2 * = 100 – x2

Уравнение примет вид: x1 *2 = x2 * .

В системе координат X1 * O* X2 * данное уравнение является каноническим уравнением параболы.



На рисунке область допустимых значений – ограниченная часть плоскости ABCD. Из полученного графика видно, что точка абсолютного максимума Z лежит внутри ОДР. Следовательно, целевая функция принимает максимальное значение в этой точке:

max Z = Z(O) = Z(9;2) = 17

Задание 3.1

После нескольких лет эксплуатации оборудование может оказаться в одном из трех состояний:

1) требуется профилактический ремонт;

2) требуется замена отдельных деталей и узлов;

3) требуется капитальный ремонт.

В зависимости от ситуации руководство предприятия может принять следующие решения:

1) отремонтировать оборудование своими силами, что потребует затрат а ;

2) вызвать специальную бригаду ремонтников, расходы в этом случае составят b ;

3) заменить оборудование новым, реализовав устаревшее по остаточной стоимости.. Совокупные затраты на это мероприятие составят с .

Требуется найти оптимально решение данной проблемы по критерию минимизации затрат с учетом следующих предположений:

а) на основе обобщения опыта эксплуатации аналогичного оборудования определены вероятности наступления соответствующих состояний – q ;

б) имеющийся опыт свидетельствует о равной вероятности наступления соответствующих состояний;

в) о вероятностях наступления соответствующих состояний ничего определенного сказать нельзя.

П1

П2

П3

a

13

9

15

b

20

12

11

c

18

10

14

q

0.3

0.45

0.25

λ = 0.7

Составим платёжную матрицу, в которой Пj – состояния оборудования, Аi – альтернативы принятия решений:

П1

П2

П3

А1

-13

-9

-15

А2

-20

-12

-11

А3

-18

-10

-14

Для принятия оптимального решения в случае а). воспользуемся критерием Байеса; в случае б). критерием Лапласа; в случае в). критериями Вальда, Сэвиджа, Гурвица.

а). на основе обобщения опыта эксплуатации аналогичного оборудования определены вероятности наступления соответствующих состояний: q1 = 0.3; q 2 = 0.45; q 3 = 0.25

Критерий Байеса.

Для каждой альтернативы найдём средний выигрыш: ` ai = ∑ aij × qj

`a1 = -11.7 `a2 = -14.15 `a3 = -13.4

П1

П2

П3

`ai

А1

-13

-9

-15

-11.7

А2

-20

-12

-11

-14.15

А3

-18

-10

-14

-13.4

qj

0.3

0.45

0.25

Из средних выигрышей выбираем максимальный: max ai = ` a 1 = -11.7 – первая альтернатива оптимальна в случае известных вероятностей наступления событий при выборе решения по критерию Байеса.

б). имеющийся опыт свидетельствует о равной вероятности наступления соответствующих состояний;

Критерий Лапласа.

Для каждой альтернативы найдём средний выигрыш: ` ai = 1/3∑ aij

`a1 = -12.3 `a2 = -14.3 `a3 = -14

П1

П2

П3

`ai

А1

-13

-9

-15

-12.3

А2

-20

-12

-11

-14.3

А3

-18

-10

-14

-14

Из средних выигрышей выбираем максимальный: max ai = ` a 1 = -12.3 – первая альтернатива оптимальна в случае равной вероятности наступления событий при выборе решения по критерию Лапласа.

в). о вероятностях наступления соответствующих состояний ничего определенного сказать нельзя.

Критерий Вальда.

Для каждой альтернативы определим наихудший исход. di – минимальный элемент строки. Из наихудших исходов выбираем наилучший, т.е. максимальный di .

П1

П2

П3

di

А1

-13

-9

-15

-15

А2

-20

-12

-11

-20

А3

-18

-10

-14

-18

max di = d 1 = -15 – первая альтернатива оптимальна по критерию Вальда.

Критерий Сэвиджа.

Для каждого столбца находим максимальный элемент βj .

П1

П2

П3

А1

-13

-9

-15

А2

-20

-12

-11

А3

-18

-10

-14

βj

-13

-9

-11

Построим матрицу рисков, элементы которой: rij = βj - aij

max ri

0

0

4

4

7

3

0

7

5

1

3

5

В матрице рисков в каждой строке найдём максимальный риск, и из них выберем минимальный: min r = r 1 = 4 – первая альтернатива оптимальна по критерию Сэвиджа.

Критерий Гурвица.

Для каждой строки находим минимальный di и максимальный βj .

П1

П2

П3

di

βj

χi

А1

-13

-9

-15

-15

-9

-13.2

А2

-20

-12

-11

-20

-11

-17.3

А3

-18

-10

-14

-18

-10

-15.6

χ i = λ × di + (1 – λ) × βj λ = 0.7

Максимальный из элементов последнего столбца: max χ i = χ1 = -13.2 – первая альтернатива оптимальна по критерию Гурвица.