– найти точку в «середине» области ограничений (суммарное расстояние от выбранной точки до границы области минимально по квадратичному критерию) и координаты этой точки берутся за оптимальные значения.
Пусть дана система неравенств:
(1.7.1)
Графически система ограничений представлена на рисунке 1.4.4.
Система ограничений (1.7.1) записана в виде:
(1.7.2)
Или в матричном виде: .(1.7.3)
Сделаем замену в левой части системы:
; (1.7.4)
или в матричном виде: .(1.7.5)
С учетом (1.7.4) имеем:
(1.7.6)
Решим задачу, используя предварительное решение. Для исходной системы ограничений (1.7.1) получаем одно ограничение:
10*P1 + 12,5*P2 + 17,5P3 + 21,875P4 = 1(1.7.7)
и , выраженные через :
X1 = 70*Р1 + 87,5*Р2 + 70*Р3 + 87,5*Р4;
Х2 = 50*Р1 + 50*Р2 + 87,5*Р3 + 87,5*Р4.(1.7.8)
Учитывая замену (1.7.4), имеем:
Для данной задачи
Заменим переменные на :
и ограничение (1.7.7):
10*X1 + 12,5*X2 + 17,5X3 + 21,875X4 = 1.(1.7.9)
Находим максимальное и минимальное значение на области (1.7.9). Находим (табл. 1.7.1, 1.7.3).
Таблица 1.7.1
-x1 -х2 -х3 -х4 В
0 10 12,5 17,5 21,9 1
Fmax 70 87,5 70 87,5 595630-2540Разрешающий элемент А(1,1).
00Разрешающий элемент А(1,1).
0
Таблица 1.7.3
0 -х2 -х3 -х1 В
-x4 0,05 0,57 0,8 0,46 0,05
Fmax -4 37,5 0 30 -4
Максимальное значение . Количество пересчетов элементов матрицы:
Находим (табл. 1.7.4 и 1.7.5).
Таблица 1.7.4
-x1 -х2 -х3 -х4 В
0 10 12,5 17,5 21,9 618490189865Разрешающий элемент А(1,1).
00Разрешающий элемент А(1,1).
1
Fmin -70 -87,5 -70 -87,5 0
Таблица 1.7.5
0 -х2 -х3 -х4 В
-x1 0,1 1,3 1,8 2,2 0,1
Fmin 7 0 52,5 65,6 7
Минимальное значение . Количество пересчетов элементов матрицы: .
Находим (табл. 1.7.6 и 1.7.7).
Таблица 1.7.6
-x1 -х2 -х3 -х4 В
0 10 12,5 17,5 21,9 584200187325Разрешающий элемент А(1,1).
00Разрешающий элемент А(1,1).
1
Fmax -70 -87,5 -70 -87,5 0
Таблица 1.7.7
0 -х2 -х3 -х4 В
-x1 0,1 1,3 1,8 2,2 0,1
Fmax 7 0 52,5 65,6 7
Максимальное значение .
Количество пересчетов элементов матрицы:
.
Находим (табл. 1.7.8, 1.7.10).
Таблица 1.7.8
-x1 -х2 -х3 -х4 В
0 10 12,5 17,5 21,9 1
Fmin 70 87,5 70 87,5 549910-3810Разрешающий элемент А(1,1).
00Разрешающий элемент А(1,1).
0
Таблица 1.7.10
0 -х2 -х3 -х1 В
-x4 0,1 1,3 1,8 2,2 0,1
Fmin -4 37,5 0 30 -4
Минимальное значение .
Количество пересчетов элементов матрицы:
Находим (табл. 1.7.11, 1.7.13).
Таблица 1.7.11
-x1 -х2 -х3 -х4 В
0 10 12,5 17,5 21,9 1
Fmax 50 50 87,5 87,5 576580-10795Разрешающий элемент А(1,1).
00Разрешающий элемент А(1,1).
0
Таблица 1.7.13
0 -х2 -х3 -х1 В
-x4 0,05 0,57 0,8 0,46 0,05
Fmax -4 0 17,5 10 -4
Максимальное значение .
Количество пересчетов элементов матрицы:
Находим (табл. 1.7.14 и 1.7.15).
Таблица 1.7.14
-x1 -х2 -х3 -х4 В
0 10 12,5 17,5 21,9 1
Fmin -50 -50 -87,5 -87,5 6070606985Разрешающий элемент А(1,1).
00Разрешающий элемент А(1,1).
0
Таблица 1.7.15
0 -х2 -х3 -х4 В
-x1 0,1 1,3 1,8 2,2 0,1
Fmin 5 12,5 0 21,9 5
Минимальное значение .
Количество пересчетов элементов матрицы:
Находим (табл. 1.7.16 и 1.7.17).
Таблица 1.7.16
-x1 -х2 -х3 -х4 В
0 10 12,5 17,5 21,9 1
Fmax -50 -50 -87,5 -87,5 626110-8255Разрешающий элемент А(1,1).
00Разрешающий элемент А(1,1).
0
Таблица 1.7.17
0 -х2 -хЗ -х4 В
-x1 0,1 1,3 1,8 2,2 0,1
Fmax 5 12,5 0 21,9 5
Максимальное значение .
Количество пересчетов элементов матрицы:
Находим (табл. 1.7.18, 1.7.20).
Таблица 1.7.18
-x1 -х2 -х3 -х4 В
0 10 12,5 17,5 21,9 1
Fmin 50 50 87,5 87,5 610870-16510Разрешающий элемент А(1,1).
00Разрешающий элемент А(1,1).
0
Таблица 1.7.20
0 -х2 -х3 -х1 В
-x4 0,05 0,57 0,8 0,46 0,05
Fmin -4 0 17,5 10 -4
Минимальное значение .
Количество пересчетов элементов матрицы:
Так как в данной задаче , то значения их максимумов и минимумов совпадают:
Аналогично:
Общее количество пересчетов:
Среднее значение находим по формуле:
.
Следовательно:
Учитывая равенство максимумов и минимумов у соответствующих функций, получаем:
С учетом (1.7.5) можно записать:
. Отсюда :
.
Таким образом, решение данной задачи:
При решении задачи без предварительного решения было выполнено пересчетов элементов таблицы, а при использовании предварительного решения – пересчетов.
Коэффициенты целевой функции и правых частей ограничений могут изменяться с течением времени, поэтому желательно их оптимальное прогнозирование.
Метод повышения оперативности анализа решения задачи линейного программирования
Оперативный анализ при изменении коэффициентов С целевой функции
Пример 1.10.1. На этапе планирования решить задачу линейного программирования
I = x1 + 5×2 max(1.10.1)
при наличии ограничений:
(1.10.2)
Пример 1.10.2. Найти решение задачи из примера 1.10.1 при различных значениях коэффициентов целевой функции (с1 = 1; с2 = 5) и(с1 = 1; с2 = 1).
Применим метод предварительного решения системы ограничений.
Решение системы линейных неравенств (1.10.2) имеет вид:
(1.10.3)
Подставляя х1, х2 в выражение для целевой функции при заданных значениях c1 и c2 в первом варианте, приходим к задаче линейного программирования:
F = 6P1 + 10P2 + 26P3 + 30P4 max(1.10.4)
при наличии ограничений
(1.10.5)
Решение задачи (1.10.4), (1.10.5) имеет вид:
P1 = 0; P2 = 0; P3 = 0; P4 = 1.
При этом х1 = 5, х2 = 5. Число вычисляемых в симплекс-методе элементов матриц равно 20 (вместо 60 при известном подходе), F = 30.
Для второго варианта число вычисляемых в симплекс-методе элементов матрицтакже уменьшилось иравно 20, F = 10.
Оперативный анализ при изменении величины ресурсов B
Для решения этой задачи составим двойственную задачу линейного программирования и к ней применим метод предварительного решения.
Пример 1.10.3. Оперативный анализ при изменении величины ресурсов Bиз примера 1.10.1.
Двойственная задача линейного программирования имеет вид:
найти минимум целевой функции
f = –U1+5U2 – U3+5U4
при наличии ограничений
(1.10.6)
Число вычисляемых в симплекс-методе элементов матрицравно 30, f = 30.
Анализ влияния изменения величины ресурсов B производится путем изменения коэффициентов целевой функции.
Воспользуемся теперь методом предварительного решения системы ограничений.
Решение системы линейных неравенств (1.10.6) имеет вид
U1 = P3 ; U2 = P1 + P3 + P5; U3 = P4; U4 = P2 + P4 + 5P5; U5 = P5.(1.10.7)
Ui≥ 0, i=1,5.
Подставляя Ui в выражение для целевой функции (1.10.5), приходим к задаче линейного программирования:
F = 5P1 + 5P2 + 4P3 +4P4 + 30P5 max .(1.10.8)
При наличии ограничений
(1.10.9)
Число вычисляемых элементов матриц равно 12 (вместо 30 при использовании симплекс-метода без предварительного решения),f = 30.
Рассматриваемый подход к повышению оперативности анализа может быть применен к различным задачам, сводящимся к задачам линейного программирования. Рассмотрим такую возможность на примере оперативного анализа решения задачи дробно-линейного программирования, имеющей большое практическоезначение.
Оперативный анализ решения задачи дробно-линейного программирования
Пример 1.10.4. Найти максимум целевой функции
Z = (x1 + 2×2) / (3×1 + x2)(1.10.10)
при выполнении ограничений
(1.10.11)
Преобразованная известным методом задача имеет вид:
найти максимум целевой функции
Z = y1 + 2y2(1.10.12)
при наличии ограничений
(1.10.13)
Решение задачи осуществляется стандартными методами линейного программирования на основе симплекс-метода. Число вычисляемых элементов матриц равно 192. При этом величина критерия Z = 1,3.
Пример 1.10.5. Найти решение задачи из примера 1.10.4 на основе применения метода предварительного решения системы ограничений.
Решение системы линейных ограничений(1.10.11) имеет вид
(1.10.14)
Подставляя y1, y2 в выражение для целевой функции (1.10.12) приходим к задаче линейного программирования:
F = 6p1 + 9p2 + 24p3 + 30p4 +12p5+15p6 max . (1.10.15)
При наличии ограничений
(1.10.16)
Решение задачи (1.10.15), (1.10.16) имеет вид:
p1 = 0; p2 = 0,14; p3 = 0; p4 = 0; p5=0; p6 =0.
Число вычисляемых в симплекс-методе элементов матриц равно 42 (вместо 192 приподходе без предварительного решения), F = 1,3.
Получение математической модели динамического промышленного объекта
Моделирование динамического объекта начинается с установления его типа: стационарной или нестационарной, линейный или нелинейный.
Линейные стационарные объекты описываются линейными дифференциальными уравнениями с постоянными коэффициентами. Если коэффициенты линейных дифференциальных уравнений являются функциями независимых переменных, то объект относится к классу линейных нестационарных.
Нелинейные стационарные объекты описываются линейными уравнениями с постоянными коэффициентами, а нелинейные нестационарные — нелинейными уравнениями с переменными коэффициентами.
В изучаемом курсе рассматриваются модели линейных объектов.
Аналитическое представление модели динамического объекта в виде дифференциального уравнения не является единственно возможным. Для систем автоматического регулирования принято представление модели в виде типовых линейных и нелинейных звеньев и их передаточных функций.
Примером линейного стационарного динамического объекта является электрическая цепь, содержащая активные и реактивные элементы, (рисунок 1).
Рисунок 1 – Схема электрической цепи
Переходный процесс при замыкании ключа в такой цепи описывается дифференциальным уравнением
, (8)
в котором i и Е являются функциями времени, а параметры цепи L и R — постоянными коэффициентами.
В качестве другого примера рассмотрим движение механизма, имеющего приведенный момент инерции I и момент нагрузки Мнагр, в общем случае переменный. Механизм приводится в движение моментом двигателя М, (рисунок 2)
Рисунок 2- Расчётная схема механизма
Изменение угловой скорости механизма w описывается дифференциальными уравнениями, называемыми уравнениями движения
Математическими моделями объектов в приведенных примерах являются дифференциальные уравнения первого порядка. Такие уравнения имеют семейства решений. Чтобы выбрать одно решение из многих, необходимо знать начальное значение функции, то есть ее значение в начальный момент времени.
В общем виде можно записать
y¢ = ¦ (y, t)
y (t0) = y0. (10)…
0104OLV 4.7
Преподаватель высшей школы менеджмента Преподаватель английского языка высшей категории
На странице представлен фрагмент
Уникализируй или напиши новое задание с помощью нейросети
Похожие работы
Руководитель Комитета по земельным ресурсам представил на рассмотрение главе администрации N-ского муниципального района предложения относительно мест
Руководитель Комитета по земельным ресурсам представил на рассмотрение главе администрации N-ского муниципального района предложения относительно местоположения земельных участков, которые следует...
Токсичные элементы как загрязняющие вещества пищевых продуктов предельно допустимые концентрации в пищевых продуктах
Токсичные элементы как загрязняющие вещества пищевых продуктов, предельно допустимые концентрации в пищевых продуктах Часть выполненной работыВ результате воздействия загрязненной окружающей среды, а...