Базовое возможное решение - Википедия - Basic feasible solution

В теории линейное программирование, а базовое возможное решение (BFS) - решение с минимальным набором ненулевых переменных. Геометрически каждая BFS соответствует углу многогранник возможных решений. Если существует оптимальное решение, то существует оптимальная BFS. Следовательно, чтобы найти оптимальное решение, достаточно рассмотреть BFS-ы. Этот факт используется симплексный алгоритм, который, по сути, переходит от одного BFS к другому, пока не будет найден оптимальный.^[1]

Определение

Чтобы определить BFS, мы сначала представляем линейную программу в так называемом эквациональная форма:

максимизировать ${ textstyle mathbf {c ^ {T}} cdot mathbf {x}}$

при условии ${ Displaystyle А mathbf {x} = mathbf {b}}$ и ${ Displaystyle mathbf {х} geq 0}$

куда:

• ${ displaystyle mathbf {c ^ {T}}}$ и ${ displaystyle mathbf {x}}$ являются векторами размера п (количество переменных);
• ${ displaystyle mathbf {b}}$ вектор размера м (количество ограничений);
• ${ displaystyle A}$ является м-к-п матрица;
• ${ Displaystyle mathbf {х} geq 0}$ означает, что все переменные неотрицательны.

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

В качестве предварительного этапа очистки мы проверяем, что:

• Система ${ Displaystyle А mathbf {x} = mathbf {b}}$ имеет хотя бы одно решение (иначе у всей ЛП решения нет и делать больше нечего);
• Все м строки матрицы ${ displaystyle A}$ линейно независимы, т. е. его ранг равен м (в противном случае мы можем просто удалить лишние строки, не меняя LP).

Обычно м < п, поэтому система ${ Displaystyle А mathbf {x} = mathbf {b}}$ есть много решений; каждое такое решение называется возможное решение ЛП.

Позволять B быть подмножеством м индексы из {1, ...,п}. Обозначим через ${ displaystyle A_ {B}}$ квадрат м-к-м матрица из м столбцы ${ displaystyle A}$ проиндексировано B. Если ${ displaystyle A_ {B}}$ является неособый, столбцы, проиндексированные B площадь основа из пространство столбца из ${ displaystyle A}$. В этом случае мы называем B а основа ЛП. С ранга ${ displaystyle A}$ является м, имеет хотя бы одно основание; поскольку ${ displaystyle A}$ имеет п столбцы, он имеет не более ${ Displaystyle { binom {п} {м}}}$ базы.

Учитывая основу B, мы говорим, что допустимое решение ${ displaystyle mathbf {x}}$ это базовое возможное решение с базой B если все его ненулевые переменные проиндексированы B, т.е. для всех ${ displaystyle j not in B: ~~ x_ {j} = 0}$.

Характеристики

1. BFS определяется только ограничениями LP (матрица ${ displaystyle A}$ и вектор ${ displaystyle mathbf {b}}$); это не зависит от цели оптимизации.

2. По определению, BFS имеет не более м ненулевые переменные и не менее п-м нулевые переменные. BFS может иметь менее м ненулевые переменные; в этом случае он может иметь множество различных баз, каждая из которых содержит индексы его ненулевых переменных.

3. Возможное решение ${ displaystyle mathbf {x}}$ является основным, если и только если столбцы матрицы ${ displaystyle A_ {K}}$ линейно независимы, где K - множество индексов ненулевых элементов ${ displaystyle mathbf {x}}$. ^[1]:45

4. BFS однозначно определяется базисом B: для каждой основы B из м индексов, существует не более одного BFS ${ Displaystyle mathbf {x_ {B}}}$ с основанием B. Это потому что ${ Displaystyle mathbf {x_ {B}}}$ должен удовлетворять ограничению ${ Displaystyle A_ {B} mathbf {x_ {B}} = b}$, а по определению базиса матрица ${ displaystyle A_ {B}}$ неособо, поэтому ограничение имеет единственное решение. Обратное неверно: каждая BFS может происходить из разных баз.

Если единственное решение ${ Displaystyle A_ {B} mathbf {x_ {B}} = b}$ удовлетворяет ограничениям неотрицательности, то базис называется возможная основа.

5. Если линейная программа имеет оптимальное решение (т. Е. Допустимое решение, а набор допустимых решений ограничен), то она имеет оптимальную BFS. Это следствие Принцип максимума Бауэра: цель линейной программы - выпуклая; множество допустимых решений выпукло (это пересечение гиперпространств); поэтому цель достигает своего максимума в крайней точке множества возможных решений.

Поскольку количество BFS-s конечно и ограничено ${ Displaystyle { binom {п} {м}}}$, оптимальное решение любой ЛП может быть найдено за конечное время, просто оценив целевую функцию во всех ${ Displaystyle { binom {п} {м}}}$БФС-ы. Это не самый эффективный способ решить ЛП; то симплексный алгоритм исследует BFS-ы более эффективным способом.

Примеры

Рассмотрим линейную программу со следующими ограничениями:

${ displaystyle { begin {align} x_ {1} + 5x_ {2} + 3x_ {3} + 4x_ {4} + 6x_ {5} & = 14 x_ {2} + 3x_ {3} + 5x_ { 4} + 6x_ {5} & = 7 forall i in {1, ldots, 5 }: x_ {i} & geq 0 end {выровнено}}}$

Матрица А является:

${ displaystyle A = { begin {pmatrix} 1 & 5 & 3 & 4 & 6 0 & 1 & 3 & 5 & 6 end {pmatrix}} ~~~~~ mathbf {b} = (14 ~~ 7)}$

Здесь, м= 2 и имеется 10 подмножеств по 2 индексов, однако не все из них являются базисами: набор {3,5} не является базисом, поскольку столбцы 3 и 5 линейно зависимы.

Набор B= {2,4} является базисом, поскольку матрица ${ displaystyle A_ {B} = { begin {pmatrix} 5 & 4 1 & 5 end {pmatrix}}}$ неособен.

Единственная BFS, соответствующая этому базису, есть ${ displaystyle x_ {B} = (0 ~~ 2 ~~ 0 ~~ 1 ~~ 0)}$.

Геометрическая интерпретация

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

Каждой BFS соответствует вершина этого многогранника. ^[1]:53–56

Применение в симплексном алгоритме

В симплексный алгоритм сохраняет в каждой точке своего исполнения «текущую основу» B (подмножество м снаружи п переменные), "текущая BFS" и "текущая таблица". Таблица представляет собой представление линейной программы, где основные переменные выражены через неосновные: ^[1]:65

куда ${ displaystyle x_ {B}}$вектор м основные переменные, ${ displaystyle x_ {N}}$вектор п неосновные переменные и ${ displaystyle z}$ является целью максимизации. Поскольку неосновные переменные равны 0, текущая BFS равна ${ displaystyle p}$, а текущая цель максимизации - ${ displaystyle z_ {0}}$.

Если все коэффициенты в ${ displaystyle r}$отрицательны, то ${ displaystyle z_ {0}}$является оптимальным решением, поскольку все переменные (включая все неосновные переменные) должны быть не меньше 0, поэтому вторая строка подразумевает ${ Displaystyle Z Leq Z_ {0}}$.

Если некоторые коэффициенты в ${ displaystyle r}$положительны, тогда можно будет увеличить цель максимизации. Например, если ${ displaystyle x_ {5}}$неосновной и ее коэффициент в ${ displaystyle r}$положительно, то увеличение его выше 0 может привести к ${ displaystyle z}$больше. Если это возможно сделать без нарушения других ограничений, тогда увеличенная переменная становится базовой (она «входит в базовую»), в то время как другая небазовая переменная уменьшается до 0, чтобы сохранить ограничения равенства и, таким образом, становится не базовой (она «выходит с базы»).

Если этот процесс проделан аккуратно, то можно гарантировать, что ${ displaystyle z}$ увеличивается, пока не достигнет оптимального BFS.

Рекомендации