Задачи линейного программирования и методы их решения

МИНИСТЕРСТВА ОБРАЗОВАНИЯ И НАУКИ РЕСПУБЛИКИ КАЗАХСТАН
ДЕПАРТАМЕНТ ОБРАЗОВАНИЯ ГОРОДА АЛМАТЫ
АЛМАТИНСКИЙ КОЛЛЕДЖ УПРАВЛЕНИЯ И РЫНКА
НА ТЕМУ: “ ЗАДАЧИ ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ И МЕТОДЫ ИХ РЕШЕНИЯ”
ПРОВЕРИЛА: АЙМУХАНБЕТОВА Ш. Е
ВЫПОЛНИЛА: УЧАЩИЕСЯ III КУРСА СПЕЦИАЛЬНОСТИ ПРОГРАММНОЕ
ОБЕСПЕЧЕНИЕ ВТ И АС
НАБИЕВА А. Б
АЛМАТЫ 2007 Г
СОДЕРЖАНИЕ
1. Общая задача линейного программирования
2. Элементы линейной алгебры и геометрии выпуклых множеств
- Свойства задачи линейного программирования
- Теоретические основы методов линейного программирования
- Геометрические методы решения задач линейного программирования
- Применение метода Жордана - Гаусса к решению системных линейных уравнений
- Библиографический список
Примеры задач линейного программирования позволяют сформировать общую задачу линейного программирования. Дана система m линейных уравнений и неравенств с n переменными
α 11 x 1 + α 12 x 2 + …+ α 1n x n ≤ b 1 ,
α 21 x 1 + α 22 x 2 + …+ α 2n x n ≤ b 2 ,
. .
α k1 x 1 + α k2 x 2 + …+ α kn x n ≤ b k ,
α k +l, l + α k+1, 2 x 2 + …+ α k+1, n x n = b k+1
α k+2, 1 x 1 + α k+2, 2 x 2 + …+ α k+2, n x n = b k+2
α m1 x 1 + α m2 x 2 + …+ α mn x n = b m
и линейная функция
F = c 1 x 1 + c 2 x 2 + …+ c n x n
Необходимо найти такое решение системы Х=(x 1 , x 2 , …, х j , …, x n ), где х j ≥0 (j = 1, 2, …, l:l ≤ n), при котором линейная функция F принимает оптимальное значение. Система называется системой ограничений, а функция F - линейной функцией цели. Более кратко общую задачу линейного программирования можно представить в виде:
F = Σc j x j → max ( или → min)
при ограничениях:
Σ α ij x j ≤ b i (i = 1, 2, …, k),
Σ α ij x j ≤ b i (i =k + 1, k + 2, …, m),
X j ≥0 (j = 1, 2, …, l; l ≤ n) .
Оптимальным решением задачи линейного программирования называется решения Х=(x 1 , x 2 , …, х j , …, x n ) система ограничений удовлетворяющее условию, при котором линейная функция принимает оптимальное значение. Термины «решение» и «план» - синонимы, однако первый используется чаще, когда речь идет о формальной стороне задачи (ее математическом решении), а второй - о содержательной стороне (экономические интерпретации) . При условии, что все переменные неотрицательны ( X j ≥0, j = 1, 2, …, n), система ограничений состоит лишь из одних неравенств, - такая задача линейного программирования называется стандартной, если система ограничений из одних уравнений, та задача называется канонической. Так, в приведенном выше примерах задач линейного программирования задачи 1 и 2 - стандартные, задача 4 - каноническая, а задача, а задача 3 - общая.
Любая задача линейного программирования может быть сведена к канонической, стандартной или общей задаче.
Ранг матрицы. Эквивалентные матрицы. Дана прямоугольная матрица.
α 11 α 12 … α 1n
A = α 21 α 22… α 2n
α n1 α n2…. … α nn
выделим в этой матрице к произведенных столбцов и к произведенных строк (k ≤ m, k ≤ n) . Определитель которого порядка, составлений из элементов матрицы А, расположенных на пересечении выделенных строк и столбцов, называются минором которого порядка матрицы А . матрица А имеет
С · С миноров которого порядка.
Рассмотрим всевозможные миноры матрицы А, отличные от нуля. Рангом матрицы А называется наибольший порядок минора этой матрицы принимают равным нулю. Всякий отличный от нуля минор матрицы, порядок которого равен рангу этой матрицы называется базисным минором матрицы. Ранг матрицы А будем обозначать через r(A) . r(A) = r(B) , матрицы А и В называется эквивалентными. В этом случае записывают А ~ В. ранг матрицы не изменяется от элементарных приобретений. Под элементами преобразованиями понимается:
- замена строк столбцами, αстолбцов соответственно строками;
- перестановка строк матрицы;
- вычеркните строки, все элементы которой равны нулю;
- умножение какой - либо строки на число не равно нулю;
- прибавление к матрицам одной строки соответствующих элементов другой строки.
ПРИМЕР 1
Определить ранг матрицы
2 3 4
4 6 8
6 9 12
определитель матрицы второго и третьего порядка.
4 6 1 2
М = = 0; М =
6 9 2 4
все миноры второго и третьего порядка равны нулю, следовательно, ранг матрицы r (A) = 1
ПРИМЕР 2
1 0 0 0 5
0 0 0 0 0
2 0 0 0 11
Вычеркним строки и столбцы с нулевыми элементами, получим матрицу эквивалентностью данной.
1 5
2 11
так как 1 5
= 0 = 1, то ранг матрицы r (A) = 2
2 11
ПРИМЕР 3. Сколько миноров второго порядка имеет матрица
α 11 α 12 α 13
А = α 21 α 22 α 23
α 31 α 32 α 33
С • С = 3·3 = 9 миноров второго порядка
Исследования систем m линейных уравнений с n неизвестными. Дана система m линейной уравнений с n неизвестными.
α 11 x 1 + α 12 x 2 +…+ α 1n x n = b 1
α 21 x 1 + α 22 x 2 +…+ α 2n x n = b 2
. .
α m1 x 1 + α m2 x 2 +…+ α mn x n = b m
Решением этой системы называется совокупность n чисел
( x 1 , x 2 , …, x n ), которое при подстановки в уравнения, обращают их в тождества. Система уравнений называется совместный, если она имеет хотя бы одно решение. Если же система не имеет ни одного решения, то она называется несовместной. Совместная система называется определенной, если она имеет больше одного решения.
Матрица
α 11 α 12 … α 1n
А = α 21 α 22 … α 2n
.
α 31 α 32 …α 3n
называется матрицей системы. Матрица
α 11 α 12 … α 1n b 1
А = α 21 α 22 … α 2n b 2
. .
α m1 α m2 …α mn b m
называется расширенной матрицей этой системы.
Для совместной системы необходимо и достаточно, чтобы ранг матрицы система равнялся рангу расширенной матрицы этой системы.
r(A) = r(A 1 ) = r
Если b = b 2 = … b n = 0, то система линейного уравнения называется однородной. Однородная система уровней всегда совместна. Если ранг совместной системы равен числу неизвестных (т. е. r = n), то система будет определенной. Если же ранг совместимой системы меньших числа неизвестных, то система неопределенная.
Задачи линейного программирования показывает, что любая задача линейного программирования может быть представлена в виде общей, канонической или стандартной задачи. В данном разделе будем рассматривать каноническую задачу, в которой система ограничений есть система уравнений. Для обоснования свойства задачи линейного программирования и методов ее решения целесообразно рассматривать еще два вида записи канонической задачи.
Матричная форма записи:
F = CX → max (min)
при ограничениях
АX = В,
X ≥ 0
где
α 11 α 12… α 1n x 1 b 1
α 11 α 11 … α 11 x 2 b 2
С = (с 1 , с 2 , …, с n ) ; A = ; X = … B = . . .
α 11 α 11 … α 11 x n b m
здесь С - матрица - строка, А - матрица систем, Х - матрица - столбец переменных, В - матрица - столбец свободных членов.
Векторная форма записи:
F = CX → max (min)
при ограничениях
p 1 x 1 + p 2 x 2 + …+ p n x n = p
x ≥ 0
где CX - скалярное произведение векторов С = (с 1 , с 2 , …, с n ) и
X = ( x 1 , x 2 , …, x n ), векторы
α 11 α 11 α 1n b 1
P 1 = α 11 , P 2 = α 11 , P 1 = α 2n , P 1 = b 1
… … . . . …
α 11 α 11 α mn b 1
состоят соответственно из коэффициентов при переменных и свободных членов. Векторное неравенство x ≥ 0 означает, что все компоненты вектора неотрицательны, т. е. x j ≥ 0, j = 1, 2, …, n.
Множества всех допустимых решений систем ограничений задачи линейного программирования является выпуклым. Пусть
Х 1 =(x 1 , x 2 , …, x n ) и Х 2 =(x 1 , x 2 , …, x n ) - два допустимых решения задачи, заданной в матричной форме. Тогда АХ 1 = В и АХ 2 =В
Рассмотрим выпуклую линейную комбинацию решений Х 1 и Х 2 , т. е. х = α 1 х 1 + α 2 х 2 при α 1 ≥0, α 2 ≥0 и α 1 + α 2 = 1 и покажем, что она также является допустимым решением системы самом деле
AX = A(α 1 х 1 + α 2 х 2 ) = α 1 AX 1 + (1 - α 1 ) AX 2 = α 1 B + (1 - α 1 ) B = B, т. е. решение Х удовлетворяет систему. Но так как х 1 ≥0, х 2 ≥0, α 1 ≥0,
Α 2 ≥0, то и х ≥0, т. е решение Х удовлетворяет и условно. Итак, доказано, что множество всех допустимых решений задачи линейного программирования выпуклым, а точнее, представляет выпуклый многогранник или выпуклую многогранную область, которое в дальнейшем будет называть одним термином - многогранником решений. Ответ на вопрос, в какой точке многогранника решений возможно оптимальное решение задачи линейного программирования, дается в следующей фундаментальной теореме. В этом разложении среди значений
F(x j ) (j = 1, 2, …, p) выберем максимальную. Пусть оно соответствует угловой точке x к (1≤ k ≤ p) ; обозначим его через M, т. е. F(x к ) = М.
Заменим в выражении каждое значение этим максимальным значением М. тогда, учитывая, что α j ≥0, Σ α j = 1,
F(x*) ≤ α 1 M + α 2 M + …+ α p M = M Σ α j =M. Под предположению оптимальное решение, поэтому, с одной стороны, F(x*) ≥ F(x к ) = М ,
то с другой стороны, доказано, что F(x*) ≤ М , следовательно,
F(x*) = М = F(x к ), где x к - угловая точка. Итак, существует угловая точка x к , линейная функция принимает максимальное значение более чем в одной угловой точке, например, в точках
X 1 , X 2 , …, X q 1≤ q ≤p ; F(X 1 , ) = F(X 2 ) = …= F(X q ) = M.
Пусть Х - выпуклая линейная комбинация этих угловых точек, т. е
F(X) = F(α, X 1 + α 2 X 2 + …+ α q X q ) = α 1 F(X 1 ) + α 2 F(X 2 ) + …+ α q F(X q ) = = α 1 M + α 2 M + …+ α q M = MΣα j = M,
т. е. линейная функция. F принимает максимальное значение в произвольной точке Х, является выпуклой линейной комбинацией угловых точек Х 1 , Х 2 , …, Х q .
Каждому допустимому базисному решению задачи линейного программирования соответствует угловая точке многогранника решений соответствует допустимое базисное решение.
Пусть Х 1 =(x 1 , x 2 , …, x m ; 0, 0, …, 0) - допустимое базисное решение системных ограничений задачи, в котором первые m компонент - основные переменные, а остальные n - m компонент - не основанные переменные, равные нулю в базисном решении (если это не так, то соответствующие переменные можно перенумеровать) . Покажем, что Х - угловая точка многогранника решений. Предположим противное, т. е. что Х не является угловой точкой. Точку Х можно представить внутренней точкой отрезка, соединяющего две различные не совпадающие с Х , точки
Х 1 =(x 1 , x 2 , …, x m ; 0, 0, …, 0) и
Х 2 =(x 1 , x 2 , …, x m ; 0, 0, …, 0)
другими словами, - выпуклой линейной комбинацией точек x 1 и x 2 многогранника решений, т. е
X = α 1 X 1 + α 2 X 2
где α 1 ≥0, α 2 ≥0, α 1 + α 2 = 1 ( полагаем, что α 1 = 0, α 2 = 0, ибо в противном случае точка Х совпадает с точкой X 1 или X 2 )
Запишим векторное равенство в координатной форме:
Х 1 = α 1 X 1 + α 2 X 2 ,
. .
Х 1 = α 1 X m + α 2 X m ,
... продолжение- Информатика
- Банковское дело
- Оценка бизнеса
- Бухгалтерское дело
- Валеология
- География
- Геология, Геофизика, Геодезия
- Религия
- Общая история
- Журналистика
- Таможенное дело
- История Казахстана
- Финансы
- Законодательство и Право, Криминалистика
- Маркетинг
- Культурология
- Медицина
- Менеджмент
- Нефть, Газ
- Искуство, музыка
- Педагогика
- Психология
- Страхование
- Налоги
- Политология
- Сертификация, стандартизация
- Социология, Демография
- Статистика
- Туризм
- Физика
- Философия
- Химия
- Делопроизводсто
- Экология, Охрана природы, Природопользование
- Экономика
- Литература
- Биология
- Мясо, молочно, вино-водочные продукты
- Земельный кадастр, Недвижимость
- Математика, Геометрия
- Государственное управление
- Архивное дело
- Полиграфия
- Горное дело
- Языковедение, Филология
- Исторические личности
- Автоматизация, Техника
- Экономическая география
- Международные отношения
- ОБЖ (Основы безопасности жизнедеятельности), Защита труда
