Главная страница
Навигация по странице:

  • Решение Отношение r 1

  • Отношение r 1

  • Эконометрика. Решение задачи с построением области допустимых решений (одр) и целевой функции для каждой итерации


    Скачать 134.5 Kb.
    НазваниеРешение задачи с построением области допустимых решений (одр) и целевой функции для каждой итерации
    Дата09.12.2018
    Размер134.5 Kb.
    Формат файлаdoc
    Имя файлаЭконометрика.doc
    ТипРешение
    #59363


    Задача 1
    Дать графическое решение задачи с построением области допустимых решений (ОДР) и целевой функции для каждой итерации:
    Вариант 3.

    x1 + x2 ≥ 4;

    2x1 - x2 ≥ 1;

    x1 - 2x2 ≥ 1;

    При естественных ограничениях x1 ≥ 0 и x2 ≥ 0;

    L (x1, x2) = 2 x1 - 3x2 + 1→ min.
    Решение
    Построим область допустимых решений, т.е. решим графически систему неравенств. Для этого построим каждую прямую и определим полуплоскости, заданные неравенствами (полуплоскости обозначены штрихом).

    Построим уравнение x1+x2 = 4 по двум точкам. Для нахождения первой точки приравниваем x1 = 0. Находим x2 = 4. Для нахождения второй точки приравниваем x2 = 0. Находим x1 = 4. Соединяем точку (0;4) с (4;0) прямой линией.

    Построим уравнение 2x1-x2 = 1 по двум точкам. Для нахождения первой точки приравниваем x1 = 0. Находим x2 = -1. Для нахождения второй точки приравниваем x2 = 0. Находим x1 = 0.5. Соединяем точку (0;-1) с (0.5;0) прямой линией.

    Построим уравнение x1-2x2 = 1 по двум точкам. Для нахождения первой точки приравниваем x1 = 0. Находим x2 = -0.5. Для нахождения второй точки приравниваем x2 = 0. Находим x1 = 1. Соединяем точку (0;-0.5) с (1;0) прямой линией.


    Рис.1. Изображение ограничений в плоскости переменных.

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

    Обозначим границы области многоугольника решений ( в данном случае он не замкнут).



    Рис.2. Границы области допустимых решений (ОДР).
    Рассмотрим целевую функцию задачи L = 2x1-3x2+1 → min.
    Построим линию уровня - прямую, отвечающую значению функции L = 0: L = 2x1-3x2+1 = 0 (пунктирная линия на рис.3). Вектор-градиент, составленный из коэффициентов целевой функции, указывает направление максимизации L (x1, x2). Начало вектора – точка (0; 0), конец – точка (2; -3). Будем двигать прямую, перпендикулярную этому вектору, параллельным образом в направлении, противоположном направлению градиента – в этом направлении функция наибыстрейшим образом убывает. Поскольку нас интересует минимальное решение, поэтому двигаем прямую до последнего касания ОДР.


    Рис.3. Графическое решение задачи.
    Прямая L (x1, x2) = const пересекает область в точке A. Так как точка A получена в результате пересечения прямых (1) и (3), то ее координаты удовлетворяют уравнениям этих прямых:


    x1+x2=4
    x1-2x2=1


    Решив систему уравнений, получим: x1 = 3, x2 = 1.

    Отсюда найдем минимальное значение целевой функции:
    L (x1, x2) = 2*3 – 3*1 + 1 = 4

    Задача 2.
    Дать аналитическое решение задачи, переходя к канонической форме (добавлением балансных переменных), расчетом допустимых базисных решений для каждой итерации (провести сопоставление с предыдущим графическим решением на каждом шаге решения задачи):
    x1 + x2 ≥ 4;

    2 x1 - x2 ≥ 1;

    x1- 2x2 ≥ 1;

    При естественных ограничениях x1 ≥ 0 и x2 ≥ 0;

    L (x1, x2) = 2 x1 - 3x2 + 1→ min.
    Решение
    Для приведения ЗЛП к канонической форме необходимо ввести неотрицательные балансовые переменные s1, s2, s3.

    x1 + x2 - s1 = 4;    

    2 x1 - x2 - s2 = 1;    

    x1 - 2 x2 - s3 = 1;    

    x1, x2, s1, s2, s3 ≥ 0;

    L (x1, x2) = 2 x1 - 3x2 + 1→ min.

    Ищем в системе ограничений базисные переменные (входящие только в одно уравнение с единичным коэффициентом).

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

    Введем по одной искусственной неотрицательной переменной ri в каждое уравнение системы ограничений.

    Получим следующую систему ограничений,

    x1 + x2 - s1 + r1 = 4    (1)

    2 x1 - x2 - s2 + r2 = 1    (2)

    x1 - 2 x2 - s3 + r3 = 1    (3)

    x1, x2, s1, s2, s3, r1, r2, r3 ≥ 0
    с базисными переменными r1,r2,r3.




    Целью решения вспомогательной задачи является получение допустимого базисного решения не содержащего искусственных переменных (r1,r2,r3). Для этого сформируем вспомогательную целевую функцию :

    G = r1 + r2 + r3

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

    Для решения вспомогательной задачи симплекс-методом выразим функцию G через свободные переменные, для этого последовательно вычтем из функции G все уравнения.

    Функция G примет вид :

    G = - 4 x1 + 2 x2 + s1 + s2 + s3 + 6.

    Теперь мы можем сформировать начальную симплекс-таблицу.


    Начальная симплекс-таблица

    БП

    x1

    x2

    s1

    s2

    s3

    r1

    r2

    r3

    Решение

    Отношение

    r1

    1

    1

    -1

    0

    0

    1

    0

    0

    4

    4/1 = 4

    r2

    2

    -1

    0

    -1

    0

    0

    1

    0

    1

    1/2 = 0.5

    r3

    1

    -2

    0

    0

    -1

    0

    0

    1

    1

    1/1=1

    Q

    2

    -3

    0

    0

    0

    0

    0

    0

    -1

    --

    G

    -4

    2

    1

    1

    1

    0

    0

    0

    -6

    --


    В качестве базовой переменной выбираем x1 (у нее в определяющей строке – пока это строка G – максимальный по модулю отрицательный элемент). Разрешающий элемент, стоящий на главной диагонали, РЭ=1. Строка, соответствующая переменной x1, получена в результате деления всех элементов строки x1 на разрешающий элемент РЭ=1. В остальных клетках столбца x1 записываем нули.

    Все остальные элементы определяются по правилу прямоугольника: выбираются из старого плана четыре числа, расположенные в вершинах прямоугольника и всегда включающие разрешающий элемент РЭ, и новый элемент НЭ находится по формуле НЭ=СЭ-(А*В)/РЭ, где СЭ - элемент старого плана, РЭ - разрешающий элемент, А и В - элементы старого плана, образующие прямоугольник с элементами СЭ и РЭ.
    Итерация 1

    БП

    x1

    x2

    s1

    s2

    s3

    r1

    r3

    Решение

    Отношение

    r1

    0

    1.5

    -1

    0.5

    0

    1

    0

    3.5

    3.5 / 0.5 = 7

    x1

    1

    -0.5

    0

    -0.5

    0

    0

    0

    0.5

    --

    r3

    0

    -1.5

    0

    0.5

    -1

    0

    1

    0.5

    0.5/0.5=1

    Q

    0

    -2

    0

    1

    0

    0

    0

    -2

    --

    G

    0

    0

    1

    -1

    1

    0

    0

    -4

    --


    Итерация 2

    БП

    x1

    x2

    s1

    s2

    s3

    r1

    Решение

    Отношение

    r1

    0

    3

    -1

    0

    1

    1

    3

    3 / 3 = 1

    x1

    1

    -2

    0

    0

    -1

    0

    1

    --

    s2

    0

    -3

    0

    1

    -2

    0

    1

    --

    Q

    0

    1

    0

    0

    2

    0

    -3

    --

    G

    0

    -3

    1

    0

    -1

    0

    -3

    --

    Итерация 3

    БП

    x1

    x2

    s1

    s2

    s3

    Решение

    Отношение

    x2

    0

    1

    -1/3

    0

    1/3

    1

    --

    x1

    1

    0

    -2/3

    0

    -1/3

    3

    --

    s2

    0

    0

    -1

    1

    -1

    4

    --

    Q

    0

    0

    1/3

    0

    5/3

    -4

    --

    G

    0

    0

    0

    0

    0

    0

    --

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


    Итерация 4

    БП

    x1

    x2

    s1

    s2

    s3

    Решение

    Отношение

    x2

    0

    1

    - 1/3

    0

    1/3

    1

    --

    x1

    1

    0

    - 2/3

    0

    -1/3

    3

    --

    s2

    0

    0

    -1

    1

    -1

    4

    --

    Q

    0

    0

    1/3

    0

    5/3

    -4

    --


    Достигнуто оптимальное решение, т.к. в строке целевой функции нет отрицательных коэффициентов.
    Ответ: Оптимальное значение функции Q(x)= 4 достигается в точке с координатами: x1= 3 и x2= 1 (s1= 0, s2= 4, s3= 0).



    написать администратору сайта