Автор работы: Пользователь скрыл имя, 25 Августа 2013 в 23:46, контрольная работа
Постановка задачи, критерий оптимальности и матрица ЭММ распределения и использования удобрений
Все модели по расчетам оптимальных рационов кормления скота и птицы, применяемые на практике с использованием экономико-математических методов и ЭВМ, можно свести к трем основным:
- модель оптимальных рационов кормления скота;
- модель планирования оптимальных кормовых смесей с учетом всех ингредиентов питания;
- модель оптимального плана использования (распре¬деления) заготовленных кормов в сельскохозяйственном предприятии.
Задание №1……………………………………...……………………….…..3
Постановка задачи, критерий оптимальности и матрица ЭММ распределения и использования удобрений
Задание №2…………………………………………………………………..8
Задание №3………………………………………………………………….12
Задание №4………………………………………………………………….16
Задание №5………………………………………………………………….18
Список литературы………………………………………………………….
Задание №4
Разработать рацион кормления
коров с минимальной
Вид питательного вещества |
Содержание питательных веществ в 1 кг |
Минимальная потребность | |
сена |
картофеля | ||
Кормовые единицы, кг. |
0,45 |
0,3 |
24 |
Переваримый протеин, гр. |
120 |
10 |
3000 |
Каротин, мг. |
30 |
2 |
1200 |
Себестоимость, руб. |
1,5 |
0,9 |
Содержание картофеля в рационе не должно быть менее 20% его веса.
Содержание сена в рационе не должно быть менее 50% питательного рациона.
Решение:
х1- сено
х2 – картофель
Ограничения по потребности
0,45х1 + 0,3х2 ³ 24
120х1 + 10х2 ³ 3000
30х1 + 2х2 ³ 1200
х1³ 0, х2 ³ 0
Ограничение по составу
х2 ³ 0,2(х1 + х2) или 0,2х1 - 0,8х2 £ 0
х1 ³ 0,5(х1 + х2) или -0,5х1 + 0,5х2 £ 0
Целевая функция
1,5х1 + 0,9х2 → min
Наносим на график уравнения ограничения.
После этого определяем область допустимых значения.
Чертим вектор с координатами (1,5; 0,9) и линии уровня, перпендикулярные ему. Видим, что линия уровня пересекает область в точке (1).
Найдем координаты точки (1). Это точка пересечения прямых
0,45х1 + 0,3х2 = 24 и 30х1 + 2х2 = 1200
х1 = (24-0,3*x2)/0,3
Подставим во второе уравнение.
30*(24-0,3*x2)/0,3 + 2x2 = 1200
Откуда х2 = 22,22 кг
х1 = 38,52 кг
Себестоимость: Z = 1,5*38,52+ 0,9*22,22 = 77,78 руб.
Задание №5
Дана математическая запись модели:
-x1 - 5x2 + 3х3 = 4;
2x1 + 5x2 - 3х3 ≥ 2;
2х1 + 4х2 ≥ -4;
F(x)= -5x1 + х2 - 2х3 → max.
Решить задачу оптимизации модели модифицированным симплексным методом.
Решение:
Решим прямую задачу линейного программирования модифицированным симплексным методом.
Определим максимальное значение целевой функции F(X) = -5x1+x2-2x3 при следующих условиях-ограничений.
-x1-5x2+3x3=4
2x1+5x2-3x3≥2
2x1+4x2≥-4
Для построения
первого опорного плана
-1x1-5x2 + 3x3 + 0x4 + 0x5 = 4
2x1 + 5x2-3x3-1x4 + 0x5 = 2
2x1 + 4x2 + 0x3 + 0x4-1x5 = -4
Введем искусственные переменные x.
-1x1-5x2 + 3x3 + 0x4 + 0x5 + 1x6 + 0x7 + 0x8 = 4
2x1 + 5x2-3x3-1x4 + 0x5 + 0x6 + 1x7 + 0x8 = 2
2x1 + 4x2 + 0x3 + 0x4-1x5 + 0x6 + 0x7 + 1x8 = -4
Поскольку в начальном плане присутствуют отрицательные значения bi<0, то с помощью двойственного симплекс-метода устраняем отрицательные значения.
Полагая,
что свободные переменные
X1 = (0,0,0,0,0,4,2,-4)
x0 = 0
x6 = -1+x1+5x2-3x3
x7 = 2-2x1-5x2+3x3+x4
x8 = 2-2x1-4x2+x5
Среди свободных
членов в системе уравнений
есть отрицательные элементы. Используем
двойственный симплекс-метод.
Чтобы теперь выразить все переменные через небазисные, в выражении для x8 выразим x5 и подставим полученное выражение во все остальные равенства.
x0 = 0-5x1+x2-2x3
x6 = 4+x1+5x2-3x3
x7 = 2-2x1-5x2+3x3+x4
x5 = 4+2x1+4x2+x8
Переходим
к первому этапу
Первый этап. Для нахождения начальной допустимой базы воспользуемся методом искусственного базиса.
Имеем:
Матрица коэффициентов A = aij
-1 |
-5 |
3 |
0 |
0 |
1 |
0 |
0 |
2 |
5 |
-3 |
-1 |
0 |
0 |
1 |
0 |
-2 |
-4 |
0 |
0 |
1 |
0 |
0 |
-1 |
Матрица b.
Итерация №1.
<X> = (6, 7, 5)
Матрица c.
c = (0, 0, 0, 0, 0, 1, 1, 0)
cB = (1, 1, 0)
cN = (0, 0, 0, 0, 0)
Вычисляем:
Матрицу B-1 вычисляем через алгебраические дополнения.
u = cBB-1 = (1, 1, 0)
c* = cN - uN = (-1, 0, 0, 1, 0)
Откуда s = 1
Откуда r = 2
Итерация №2.
<X> = (6, 1, 5)
Матрица c.
c = (-1, 0, 0, 1, 0, 0, 0, 0)
cB = (0, -1, 0)
cN = (0, 0, 1, 0, 0)
Вычисляем:
u = cBB-1 = (0, -0.5, 0)
c* = cN - uN = (2.5, -1.5, 0.5, 0.5, 0)
Откуда s = 2
Откуда r = 1
Итерация №3.
<X> = (3, 1, 5)
Матрица c.
c = (0, 2.5, -1.5, 0.5, 0, 0, 0.5, 0)
cB = (-1.5, 0, 0)
cN = (2.5, 0.5, 0, 0.5, 0)
Вычисляем:
u = cBB-1 = (-1, -0.5, 0)
c* = cN - uN = (-0, 0, 1, 1, 0)
Нулевая строка
симплексной таблицы
Второй этап. Удаляем столбцы с искусственными переменными. Заменим вектор оценок С на целевую функцию.
Выразим базисные переменные:
x3 = 3.33-1.67x2-0.3333x4
x1 = 6-1x4
которые подставим в целевую функцию:
F(X) = 5(6-1x4)-x2 + 2(3.33-1.67x2-0.3333x4)
или
F(X) = 36.67+2.33x2+5.67x4
Имеем:
Матрица коэффициентов A = aij
Матрица b.
Итерация №1.
<X> = (3, 1, 5)
Матрица c.
c = (0, -2.3333, 0, -5.6667, 0)
cB = (0, 0, 0)
cN = (-2.3333, -5.6667, 0, 0, 0)
Вычисляем:
Матрицу B-1 вычисляем через алгебраические дополнения.
u = cBB-1 = (0, 0, 0)
c* = cN - uN = (-2.3333, -5.6667, 0, 0, 0)
Откуда s = 2
Выводимую переменную r найти невозможно. Прерываем процесс поиска первого опорного плана.
Вектор результатов X = (6, 0, 3.33)T
Значение целевой функции F(X) = bc = -36.67
Список литературы