Контрольная работа: Симплексный метод

Задача 1.

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

Вариант 3.

Найти наибольшее значение функции f(X) = - x₁ - x₂ + 2x₃ при ограничениях

2x₁ + x₂ + x₃ £ 2

x₁ - x₂ + x₃ £ 1,

x_j ³ 0, j = 1, 2, 3.

Решение.

Приведем задачу к каноническому виду, вводя дополнительные неотрицательные переменные x_4,5 ³ 0.

f(X) = - x₁ - x₂ + 2x₃ ® max

2x₁ + x₂ + x₃ + x₄ = 2

x₁ - x₂ + x₃ + x₅ = 1,

x_j ³ 0, j = 1, 2, 3, 4, 5.

Каноническая задача имеет необходимое число единичных столбцов, т. е. обладает очевидным начальным опорным решением.

Очевидное начальное опорное решение (0; 0; 0; 2; 1).

Решение осуществляется симплекс-методом с естественным базисом. Расчеты оформим в симплекс-таблицах

Начальное опорное решение (0; 0; 0; 1; 1), соответствующее симплекс-таблице 0, неоптимальное, так как в D - строке есть отрицательные значения, наименьшее в столбце А₃. Этот столбец будет направляющим. Минимальное положительное оценочное отношение Q в строке А₅, эта строка направляющая. Направляющий элемент на пересечении направляющих строки и столбца. Столбец А₅ выводим из базиса, а А₃ - вводим в базис. После пересчета получаем симплекс-таблицу 1. Соответствующее опорное решение (0; 0; 1; 1; 0) не оптимально, так как в D - строке есть отрицательные значения, в столбце А₂.Этот столбец будет направляющим. Минимальное положительное оценочное отношение Q в строке А₄. В качестве направляющей строки возьмем А₄. Направляющий элемент на пересечении направляющих строки и столбца. Столбец А₄ выводим из базиса, а А₂ - вводим в базис. Опорное решение, соответствующее симплекс-таблице 2 (0; 1/2; 3/2; 0; 0) - оптимально, так как в D - строке нет отрицательных значений.

Отбрасывая значения дополнительных переменных х₄ и х₅, получаем оптимальное решение исходной задачи:

х₁ = 0, х₂ = 1/2 = 0,5; х₃ = 3/2 = 1,5; f_max = -1×0 - 1×0,5 + 2×1,5 = 2,5.

Задача 2.

Задание 1. Сформулировать экономико-математическую модель исходной экономической задачи.

Задание 2. Решить полученную задачу линейного программирования графическим методом.

Задание 3. Сформулировать двойственную задачу и найти ее оптимальное решение, используя теоремы двойственности.

Вариант 3.

Из 505 м² ткани нужно сшить не более 150 женских и не более 100 детских платьев. На пошив одного женского и детского платья требуется соответственно 3 м² и 1 м² ткани. При реализации каждого женского платья получают 10 ден. единиц прибыли, а детского – 5 ден. единиц. Сколько нужно сшить женских и детских платьев, чтобы получить наибольшую прибыль?

Решение.

Задание 1.

Обозначим x₁ и x₂ количество женских и детских платьев, соответственно (план пошива). Очевидно, x_1,2 ³ 0 и целые. Так как женских платьев должно быть не более 150, то x₁ £ 150, аналогично, для детских платьев получаем x₂ £ 100. Расход ткани на план пошива (x₁, x₂) составит 3x₁ + x₂ м², эта величина не должна превышать запаса ткани 505 м². Следовательно, должно выполняться неравенство 3x₁ + x₂ £ 505.

Прибыль от продажи платьев составит f(X) = 10x₁ + 5x₂ ден. единиц, и она должна быть наибольшей

Получаем экономико-математическую модель задачи:

Найти максимум функции f(X) при заданных ограничениях

f(X) = 10x₁ + 5x₂ ® max

3x₁ + x₂ £ 505

x₁ £ 150

x₂ £ 100

x_1,2 ³ 0, целые.

Задание 2.

Решаем задачу без условия целочисленности решения. Построим множество допустимых решений задачи.

Прямые ограничения x_1,2 ³ 0 выделяют первую четверть плоскости.

Проведем прямую 3x₁ + x₂ = 505 через точки (110; 175) и (175; -5). Подставим в первое неравенство координаты точки (0; 0): 3×0 +1×0 = 0 < 505, так как неравенство выполняется, то выбираем полуплоскость, содержащую эту точку.

Проведем прямую x1 = 150 и выберем левую полуплоскость.

Проведем прямую x2 = 100 и выберем нижнюю полуплоскость.

Множество допустимых решений – это многоугольник ABCDO.

Построим линию уровня целевой функции f(X) = 10x₁ + 5x₂ 10x₁ + 5x₂ = 0 через точки (0; 0 ) и (-10; 20). Вектор-градиент {10; 5} задает направление, перемещаясь вдоль которого, можно увеличить значение целевой функции; перемещаясь в противоположном направлении, можно уменьшить ее значение.

Из чертежа видно, что наибольшее значение целевой функции будет на линии уровня, проходящей через точку В.

Координаты этой точки найдем из системы

3x₁ + x₂ = 505,

x₂ = 100.

x₁ = 135,

x₂ = 100.

f_mах = 10 ×135 + 5 ×100 = 1850 ден. единиц.

Полученное оптимальное решение оказалось целым, следовательно, это решение поставленной задачи. Получили: в оптимальном плане пошива следует сшить женских платьев 135 шт., детских – 100 шт. При этом прибыль составит 1850 ден. единиц и будет наибольшей.

Задание 3.

Двойственная задача.

Найти минимум функции g(Y) при ограничениях:

g(Y) = 505y₁ + 150y₂ + 100y₃ ® min

3y₁ + y₂ ³ 10

y₁ + y₃ ³ 5

y_1,2,3 ³ 0.

Оптимальное решение прямой задачи Х = (135; 100). Подставим его в ограничения этой задачи

3×135 + 1×100 = 505

135 < 150

100 = 100.

Условия дополняющей нежесткости (вторая теорема двойственности): для оптимальных планов двойственных задач имеют место соотношения:

Так как для оптимального решения прямой задачи второе ограничение выполняется как неравенство, то в оптимальном решении двойственной задачи y₂ = 0.

Так как для оптимального решения прямой задачи х₁ > 0и х₂ > 0, то оба ограничения двойственной задачи выполняются как равенство. Для нахождения решения двойственной задачи получаем систему

y₂ = 0

3y₁ + y₂ = 10

y₁ + y₃ = 5.

Получаем решение: y₁ = 10/3, y₂ = 0, y₃ = 5/3.

Найдем значение целевой функции двойственной задачи:

g(Y) = 505×10/3 + 150×0 + 100×5/3 = 5550/3 = 1850.

Получили g_min = f_max = 1850 ден. единиц.

Так как значения прямой и двойственной функций равны, то Y = (10/3; 0; 5/3) является оптимальным решением двойственной задачи (по первой теореме двойственности).

Задача 3.

Задание 1. Записать исходные данные задачи в виде транспортной таблицы, определить, открытой или закрытой является транспортная задача.

Задание 2. Сформулировать экономико-математическую модель исходной транспортной задачи.

Задание 3. Найти оптимальный план перевозок, отметив при этом единственность или неединственность оптимального плана.

Вариант 3.

Картофель из четырех районов должен быть перевезен в три хранилища. Запасы картофеля в районах соответственно равны 400 т, 500 т, 800 т и 500 т. Возможности хранилищ соответственно равны 700 т, 800 т и 700 т. Затраты на перевозку одной тонны картофеля из первого района в каждое из хранилищ равны соответственно 1, 4 и 3 ден. единиц; аналогичные затраты на перевозку из второго района составляют 7, 1 и 5 ден. единиц, из третьего - 4, 8 и 3 ден. единиц, из четвертого - 6, 2 и 8 ден. единиц. Найти план перевозок картофеля из районов в хранилища, при котором транспортные расходы были бы минимальными.

Решение.

Задание 1.

Сумма мощностей поставщиков (запасы картофеля в всех районах) 400+500+800+500 = 2200, сумма мощностей потребителей (возможности всех хранилищ) 700+800+700 = 2200. Суммы равны, данная задача является транспортной задачей закрытого типа.

Задание 2.

Обозначим x_ij объем поставок картофеля от i – го поставщика (района) j – му потребителю (хранилищу), i = 1, 2, 3, 4; j = 1, 2, 3. Очевидно, x_ij ³ 0. В закрытой транспортной задаче все ограничения являются равенствами.

Так как потребности должны быть удовлетворены, то выполняются условия:

х₁₁ + х₂₁ + х₃₁ + х₄₁ = 700

х₁₂ + х₂₂ + х₃₂ + х₄₂ = 800                (1)

х₁₃ + х₂₃ + х₃₃ + х₄₃ = 700

Так как поставки от поставщика всем потребителям не могут быть больше его возможностей, то выполняются условия:

х₁₁ + х₁₂ + х₁₃ = 400

х₂₁ + х₂₂ + х₂₃ = 500       (2)

х₃₁ + х₃₂ + х₃₃ = 800

х₄₁ + х₄₂ + х₄₃ = 500

Затраты на транспортировку составят

F(X) = х₁₁ + 4х₁₂ + 3х₁₃ +

+ 7х₂₁ + х₂₂ + 5x₂₃ +

+ 4х₃₁ + 8х₃₂ + 3х₃₃ +

+ 6х₄₁ + 2х₄₂ + 8х₄₃ .

Требуется найти неотрицательное решение системы уравнений (1) – (2), на котором целевая функция затрат F(X) принимает минимальное значение.

Задание 3.

Начальный план перевозок находим методом минимальной стоимости:

Заполняем клетку (1; 1) х1₁ = min {700, 400} = 400, от поставщика 1 вывезено все, в строке 1 больше поставок нет. Заполняем клетку (2; 2) х₂₂ = min {800, 500} = 500, от поставщика 2 вывезено все, в строке 2 больше поставок нет. Клетка (4; 2) х₄₂ = min {800 - 500, 500} = 300, потребителю 2 все завезено, в столбец 2 больше поставок нет. Клетка (3; 3) х₃₃ = min {700, 800} = 700, потребителю 3 все завезено, в столбец 3 больше поставок нет. Далее клетка (3; 1) х₃₁ = 100. Клетка (4; 1) х₄₁ = 200. Все клетки, в которые даны поставки, считаем занятыми, остальные – свободными. Первоначальный план перевозок задается таблицей 1.

Таблица 1.

Исследуем этот план перевозок на оптимальность методом потенциалов. Потенциалы для занятых клеток удовлетворяют уравнениям: v_j = c_{ij +} u_i.

Пусть u₁ = 0; по клетке (1; 1) находим v₁ = 1; по клетке (3; 1) находим u₃ = -3; по клетке (4; 1) находим u₄ = -5; по клетке (4; 2) находим v₂ = -3; по клетке (3; 3) находим v₃ = -0; по клетке (2; 2) находим u₂ = -4.

Для всех клеток матрицы перевозок найдем оценки клеток d_ij = (u_i + c_ij) - v_j :

Среди оценок нет отрицательных, следовательно план перевозок Х₀ (таблица 1) оптимальный.

Так как среди оценок свободных клеток есть нулевые (клетка (1; 3)), то оптимальный план перевозок не единственный.

Общие затраты на перевозки

F(X₁) = 1*400 + 1*500 + 4*100 + 3*700 + 6*200 + 2*300 = 5200 ден. единиц будут минимальными при:

x₁₁ = 400, x₂₂ = 500, x₃₁ = 100, х₃₃ = 700, x₄₁ = 200, x₄₂ = 300, остальные x_ij = 0.

По оптимальному плану перевозок следует перевезти картофеля:

из первого района в первое хранилище                   - 400 т;

из второго района во второе хранилище                 - 500 т;

из третьего района в первое хранилище                   - 100 т,

в третье хранилище              - 700 т;

из четвертого района в первое хранилище     - 200 т,

во второе хранилище - 300 т.

Задача 4

В таблице приведены годовые данные о трудоемкости производства I т цемента (нормо-смен) (N —последняя цифра зачетной книжки студента):

Задание 1. Сгладить временной ряд методом простой скользящей средней, выбрав длину интервала сглаживания m = 3; результаты отра­зить на графике.

Задание 2. Определить наличие тренда во временном ряду методом Фостера - Стьюарта. Табличные значения статистики Стьюдента t_a принять равными при уровне значимости a = 0.05 t_a = 2,23 , а при a = 0,30 - t_a = 1,09; другие необходимые табличные данные приведены в таблице 4.5 учебника на с.153 (описание метода Фостера - Стьюар­та см. учебник с. 151- 153).

Задание 3. Для исходного временного ряда построить линейную трендовую модель , определив ее параметры на основе метода наименьших квадратов (соответствующую систему нормальных уравнений см. в учебнике на с. 196 формула (5.5)).

Задание 4. Оценить адекватность построенной модели на основе исследования

а) близости математического ожидания остаточной компоненты (ряда остатков) нулю; критические значения r-критерия принять равным тому числу, как указанно в задании 2;

б) случайности отклонений остаточной компоненты по критерию пиков (поворотных точек); Расчеты выполнить на основе соотношения 5.9. учебника на с. 200;

в) независимости уровней ряда остатков (отсутствие автокорреляции) на основе критерия Дарбина — Уотсона (см. учебник с. 203— 204), используя в качестве критических значений d_l = 1.08 и d₂ = 1,36; если критерий Дарбина — Уотсона ответа не дает, исследование независимости провести по первому коэффициенту автокорреляции:

,

где e_i -- уровни остаточной компоненты;

Модуль первого коэффициента автокорреляции сравнить с критическим уровнем этого коэффициента, значение которого принять равным 0,36;

г) нормальности закона распределения уровней остаточной компоненты на основе RS-критерия;

В качестве критических значений принять интервал от 2,7 до 3,7 (см. учебник, стр. 201—-202).

Задание 5. Оценить точность построенной трендовой линейной модели, используя показатели среднего квадратического отклонения от линии тренда (формула (5,17) учебника на с. 210, k = 1) и средней относительной ошибки аппроксимации (формула (5.14) учебника на с. 204).

Задание 6. Построить точечный и интервальный прогноз трудоемкости производства 1 т цемента на два шага вперед (формула (5.18) учебника на с. 210). Результаты моделирования и прогнозирования отразить на графике.

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

Вариант 3. Условия при N = 3

Решение.

Задание 1. Сглаживание ряда Y(t) произведем по простой скользящей средней

Результаты в таблице 1.

Задание 2.

Этап 1. Строим две числовые последовательности k_t и l_t

Этап 2. Находим величины

7; 1 – 6 = -5.

Этап 3. Для n = 10 выпишем табличные значения m = 3,858; s₁ = 1,288; s₂ =         1,964.

Вычисляем

2,44; 2,55.

Этап 4.

Так как расчетные значения t_s = 2,44 и t_d = 2,55 больше табличного значения t_a = 2,23, то в данном временном ряду присутствуют тренд и тенденция в дисперсии ряда.

Из таблицы 1 видно, что ряд Y(t) имеет тенденцию к снижению.

Задание 3. Линейную трендовую модель ищем в виде . Параметры модели а₀, а₁ найдем, решив систему уравнений

.

n = 10.

Составим расчетную таблицу 2.

Получаем систему

; .

Получили 1,5а₁ = -0,64, а₁ = -0,64:1,5 = -0,43; а₀ = 6,75    - 5,5а₁ = 6,75 - 5,5×(-0,43) =         9,12.

Получили трендовую модель: .

Задание 4.

Оценим качество модели. Для этого найдем расчетные значения Y_p(t), подставляя t =1, …, 10 в трендовую модель, найдем отклонения расчетных значений от исходных E(t) = Y(t) - Y_p(t). Для исследования модели на адекватность составим таблицу 3.

а) Близость математического ожидания остаточной компоненты нулю.

Сумма остатков равна 0. Расчетное значение критерия Стьюдента

0.

Критическое значение t_a = 2,23 больше расчетного, следовательно, математическое ожидание остаточной компоненты равно нулю.

б) Проверка остатков E(t) на случайность.

Критическое количество поворотных точек для n =10 равно 2.

2.

Для данного ряда количество таких точек k = 5. Это больше 2, поэтому остатки E(t) случайные.

в) Проверка остатков E(t) на независимость.

Независимость (отсутствие автокорреляции) проверим, используя критерий Дарбина-Уотсона:

, 2,20.

d > 2, преобразуем d' = 4 - d = 4 - 2,20 = 1,80, получили 1,36 < d' = 1,80 < 2. Это означает, что остатки не зависимы.

г) Проверка остатков на соответствие нормальному закону распределения.

Используется RS - критерий:

, где 0,36.

2,87,

RS_т = 2,7 - 3,7; так как расчетное значение RS - критерия RS_расч = 2,87 попадает внутрь интервала от 2,7 до 3,7, то остатки E(t) подчиняются по нормальному закону распределения.

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

Задание 5. 

Определим точность модели.       

Среднее квадратическое отклонение от линии тренда

0,36.

Средняя относительная ошибка

.

Так как 3,85% < 5%, то точность модели высокая.

Задание 6.

Точечный прогноз для Y получим, подставляя в трендовую модель t =11 и t = 12.

4,385; 3,955.

Для интервального прогноза найдем ширину интервала

.

Для числа степеней свободы k = n -2 = 10 - 2 = 8 и уровня значимости a = 0,05 t_a = 2,31.

1,00; 1,04.

Границы интервалов прогноза: НГ = Y_n+k - U(k), ВГ = Y_n+k + U(k).

Результаты прогноза представлены таблицей 4.

Построим график.

Задача 5.

В таблице представлены первый (х_ij) и второй (Y_i) квадранты схемы межотраслевого баланса производства и распределения продукции для трехотраслевой экономической системы (N — последняя цифра зачет­ной книжки студента):

Задание 1. Рассчитать объемы валовой продукции отраслей (форму­ла (6.2) учебника на с. 237).

Задание 2. Рассчитать матрицу коэффициентов прямых затрат А = (a_ij) (формула (6.4) учебника на с. 238).

Задание 3. Найти матрицу коэффициентов полных затрат B = (b_ij), используя формулу (6.16) учебника на с. 244.

Задание 4. Рассчитать объемы условно чистой продукции отраслей Zj, используя формулу (6.1) учебника на с. 236.

Задание 5. Представить в таблице полную схему межотраслевого баланса (в соответствии с принципиальной схемой МОБ; табл. 6.1 учебника на с.234).

Вариант 3. Условия при N = 3.

Решение.

Задание 1.

Объем валовой продукции находим по формуле

.

Результаты в таблице 2.

Задание 2.

Коэффициенты матрицы прямых затрат находим по формуле

Получаем матрицу А.

Задание 3. Чтобы найти матрицу коэффициентов полных затрат В, запишем матрицу Е - А, где Е - единичная матрица.

Матрица В находится по формуле

Задание 4. Условно чистую продукцию найдем по формуле

Результаты в таблице 3.

Задание 5. Межотраслевой баланс производства и распределения продукции отраслей представлен таблицей 4.

Литература

1. Экономико-математические методы и прикладные модели: Учеб. пособие для вузов /В.В.Федосеев, А.Н.Гармаш, Д.М. Дайитбегов и др.; Под ред. В.В. Федосеева. М.: ЮНИТИ, 1999.

2. Экономико-математические методы и прикладные модели. Методические указания по изучению дисциплины и задания к контрольной работе для студентов III курса специальностей 061000 «Государственное и муниципальное управление», 061100 «Менеджмент организации», 061500 «Маркетинг». – М.: ВЗФЭИ, 2002