
- •Основи чисельних методів математики (з використанням Excel) Передмова
- •Розділ 1. Методи обчислень: предмет, основні поняття та застосування
- •§ 1. Предмет і застосування
- •§ 2. Основні поняття
- •1. Похибки наближень.
- •2. Граничні похибки. Похибки функції.
- •3. Похибки розв'язку.
- •4. Стійкість і коректність.
- •Питання, тести
- •Розділ 2. Інтерполяція функцій
- •§1. Задача інтерполювання
- •§2. Інтерполяційна формула Лагранжа
- •§3. Поділені різниці. Формула Ньютона з поділеними різницями
- •§4. Інтерполяційна формула за допомогою Excel
- •§5. Інтерполювання за схемою Ейткіна
- •§6. Скінчені різниці. Інтерполяційні формули Ньютона для рівновіддалених вузлів
- •§7. Інтерполювання із скінченими різницями за допомогою Excel
- •§8. Інші методи інтерполювання
- •Питання, тести
- •Завдання
- •Розділ 3. Чисельне диференціювання та інтегрування.
- •§ 1. Однобічні формули чисельного диференціювання
- •§ 2. Оцінки похибки чисельного диференціювання
- •§ 3. Чисельне інтегрування. Квадратурні формули
- •§ 4. Квадратурні формули Ньютона – Котеса
- •§ 5. Узагальнені квадратурні формули.
- •§ 6. Метод подвійного перерахунку.
- •1. R2n ( f ) ≈ (правило Рунге) (14)
- •§ 7. Метод кратного перерахунку за допомогою Excel
- •Питання, тести
- •Завдання
- •Розділ 4. Чисельні методи розв‘язування рівнянь з однією змінною
- •§ 1. Відокремлення коренів
- •§ 2. Метод дихотомії (поділу відрізка пополам)
- •§ 3. Ітераційні методи та оператор стиску.
- •§ 4. Похибки ітераційного процесу
- •§ 5. Реалізація методу простої ітерації за допомогою електронних таблиць
- •§ 6. Метод Ньютона. Порядок збіжності ітераційного процесу.
- •§ 7. Метод лінійного інтерполювання.
- •§ 8. Інші приклади ітераційних методів.
- •Питання, тести
- •Завдання
- •Розділ 5. Методи розв’язування систем лінійних рівнянь
- •§ 1. Метод Гаусса
- •Метод Гаусса в матричній формі
- •Елементарні операції над матрицею:
- •§ 2. Метод Гаусса за допомогою Excel
- •§ 3. Матричні операції в Excel
- •3. Множення матриць.
- •§ 4. Метод простої ітерації для слр
- •§ 5. Метод Зейделя
- •Питання, тести
- •Завдання
- •Розділ 6. Методи лінійного програмування
- •§ 1. Оптимізаційні задачі. Математичне програмування
- •§ 2. Геометричний зміст задач лінійного програмування. Графічний метод
- •§3. Канонічна форма задачі лінійного програмування. Опорні розв’язки
- •§4. Симплекс – таблиця
- •§5. Симплекс – метод.
- •§6. Розв’язування задач лінійного програмування за допомогою excel
- •§7 Приклади
- •§8. Пошук початкового опорного розв’язку. Метод штучного базису
- •Властивості допоміжної задачі.
- •Питання, тести
- •Завдання
- •Розділ 7. Чисельні методи розв’язування звичайних диференціальних рівнянь
- •§ 1. Метод Ейлера
- •§ 2. Метод Ейлера за допомогою Excel
- •§ 3. Методи Рунге – Кутта
- •§ 4. Подвійний перерахунок для методів Рунге – Кутта
- •§ 5. Кратний перерахунок для методів Рунге – Кутта за допомогою Excel
- •§ 6. Методи Рунге – Кутта з вищими порядками похибки
- •Питання, тести
- •Завдання
- •Іменний покажчик
- •Предметний покажчик
- •Література
§4. Симплекс – таблиця
Основним об’єктом у симплекс – методі є симплекс – таблиця. Фактично це розширена матриця системи обмежень задачі лінійного програмування у канонічній формі, до якої додано рівняння нульового рівня цільової функції f = 0. Формально кажучи, для задачі
f = (max)
= bi (i = 1, 2, …, m) (3)
xj ≥ 0 (j = 1, 2, …, n)
симплекс – таблицею називають таку таблицю:
x1 |
x2 |
… |
xn |
b |
а11 а21 . . . am1 |
а12 а22 . . . am2 |
… …
… |
a1n a2n . . . amn |
b1 b2 . . . bm |
c1 |
c2 |
… |
cn |
0 |
Нехай
α = (0 ;…;
;…;
;…;
;…;
0) – деякий опорний розв’язок задачі
(3), а відповідні вектори умов
,
,
… ,
утворюють його базис. Як відомо, методом
Гаусса можна привести стовпці i1
, i2 ,
… , ir
до діагонального виду,
а якщо отримані в результаті перетворень
деякі рядки містять лише нулі, то вони
відкидаються. Іншими словами можна
знайти загальний розв’язок системи
лінійних рівнянь, що відповідає симплекс
– таблиці (тобто матриці системи
обмежень), у якому основні змінні – це
,
,
… ,
,
всі решта вільні. Матриця такого розв’язку
на стовпцях основних змінних є
діагональною. Отже маємо:
|
… |
|
… |
|
… |
|
… |
|
b |
. . .
|
… …
… |
0 . . .
|
… …
… |
1 . . .
|
… …
… |
0 . . .
|
… …
… |
. . .
|
. . .
|
|
… |
|
… |
|
… |
|
… |
|
|
У
стовпці вільних членів b системи лінійних
рівнянь у відповідних рядках отримані
значення опорного розв’язку: якщо всі
вільні змінні дорівнюють нулю, то звідси
=
(k = 1,2,…,r). Віднімемо від останнього
рядка отриманої таблиці перший рядок,
помножений на
,
другий рядок, помножений на
,
… , r – й рядок, помножений на
.
В результаті дістанемо нову таблицю,
|
… |
|
… |
|
… |
|
… |
|
b |
. . .
|
… …
… |
0 . . .
|
… …
… |
1 . . .
|
… …
… |
0 . . .
|
… …
… |
. . .
|
. . .
|
|
… |
|
… |
|
… |
|
… |
|
|
де
=
=
… =
=
0, δ0
=
=
–
f (α)
(значення цільової функції на опорному
розв’язку α з протилежним знаком).
Означення 8. Отриману таблицю називають симплекс – таблицею, приведеною до базису , , … , опорного розв’язку α , а числа δj – оцінками цього базису або оцінками її основних (базових) змінних , , … , .
Зміст
оцінок δj
у тому, що це коефіцієнти
цільової функції, вираженої через вільні
змінні опорного розв’язку α. Справді,
δj
= сj
(j=1,2,…,n) ,
δ0
=
.
З другого боку виразимо
через інші змінні з рядка 1
,
з рядка 2, … ,
з рядка k
,
(k = 1,2,…,r).
Підставимо ці вирази замість
(k = 1,2,…,r)
у цільову функцію
f
=
=
)
=
=
=
+
=
=
f (α)
+
=
= – δ0 . Отже, рівняння нульового рівня цільової функції f = 0, виражене через вільні змінні опорного розв’язку α виглядає так: = δ0 . Саме воно і записано у нижньому рядку симплекс – таблиці.
Теорема 1. (Ознака оптимальності). 1) Якщо всі оцінки деякого базису опорного розв’язку α недодатні, то α є оптимальним розв’язком задачі лінійного програмування у канонічній формі (3).
2) Для оптимального опорного розв’язку задачі лінійного програмування у канонічній формі (3) завжди існує базис, всі оцінки якого недодатні.
Справді,
нехай δj
– оцінки деякого базису опорного
розв’язку α, звідки f
=
+
f (α),
δj
≤ 0 . Якщо α1
= (u1
, u2 ,
… , un)
– деякий інший опорний
розв’язок,
uj
≥ 0
(j=1,2,…,n),
то f
(α1)
=
+
f
(α) ≤ f
(α) ,
тобто розв’язок α є максимальним.
Навпаки,
нехай α – невироджений опорний розв’язок
(а у невиродженого опорного розв’язку
базис тільки один). Нехай існує така
оцінка δp
базису опорного розв’язку α,
що δp
> 0. Для
деякого θ > 0
покладемо
(k = 1,2,…,r)
,
= θ, всі інші вільні змінні
=
0. Такий набір значень
змінних задовольняє системі обмежень:
у кожному рядку k
симплекс – таблиці,
приведеної до базису опорного розв’язку
α, маємо
= =
+
= (
)
+
=
.
Оскільки всі
>
0 (
розв’язок є невиродженим), то й всі
>
0 при достатньо малому
θ > 0 . Отже,
такий набір значень змінних є набором
координат деякого допустимого розв’язку
β. Але f
(β) =
+
f (α)
= δpθ
+ f (α)
> f (α)
. Це означає, що не вироджений
опорний розв’язок не є максимальним,
якщо існує оцінка δp
> 0 його
базису. Отже, невироджений опорний
розв’язок є оптимальним тоді і тільки
тоді, коли всі оцінки його базису
недодатні. Але будь – який опорний
розв’язок можна як завгодно точно
наблизити невиродженим. Оцінки базисів
деякої послідовності таких наближень
збігаються до оцінок деякого базису
даного розв’язка і отже всі оцінки
довільного оптимального опорного
розв’язку недодатні.
Теорема 2. (Умова необмеженості цільової функції). Нехай симплекс – таблиця задачі лінійного програмування приведена до базису деякого опорного розв’язку α. Якщо серед оцінок цього базису існує додатна δs > 0, а решта елементів s – го стовпця недодатні, то цільова функція не обмежена зверху. Якщо ж існує оцінка δs < 0, а решта елементів s – го стовпця невід’ємні, то цільова функція не обмежена знизу.
Справді,
для довільного θ > 0
покладемо
(k = 1,2,…,r),
= θ, всі інші вільні змінні
=
0. Такий набір значень
змінних задовольняє системі обмежень:
у кожному рядку k
симплекс – таблиці,
приведеної до базису опорного розв’язку
α, маємо
=
+
= (
)
+
=
.
Оскільки всі
≤ 0, то й
>
0 для довільного θ >
0. Отже, такий набір значень
змінних є набором координат деякого
допустимого розв’язку β.
Але f
(β) =
+
f (α)
= δsθ
+ f (α)
необмежено зростає із
зростанням θ > 0,
оскільки за умовою δs
> 0. Другий
випадок доводиться аналогічно.
Приклад. Розглянемо опорний розв’язок α = (1;0;0;0) задачі
f = – 10x1 + x2 – 9x3 + 6x4 (min) ,
.
Змінивши знак у цільової функції, отримаємо f = 10x1 – x2 + 9x3 – 6x4 (mах). Відповідна симплекс – таблиця:
x1 |
x2 |
x3 |
x4 |
b |
1 1 |
1 – 1 |
3 – 1 |
0 2 |
1 1 |
10 |
– 1 |
9 |
– 6 |
0 |
Приведемо симплекс – таблицю цієї задачі до базису А1 , А2 опорного розв’язку α:
x1 |
x2 |
x3 |
x4 |
b |
1 0 |
0 1 |
1 2 |
1 – 1 |
1 0 |
0 |
0 |
1 |
– 17 |
– 10 |
x1 |
x2 |
x3 |
x4 |
b |
1 0 |
0 1 |
1 2 |
1 – 1 |
1 0 |
10 |
– 1 |
9 |
– 6 |
0 |
x1 |
x2 |
x3 |
x4 |
b |
1 0 |
1 – 2 |
3 – 4 |
0 2 |
1 0 |
10 |
– 1 |
9 |
– 6 |
0 |

Серед оцінок є додатна δ3 > 0. Отже, на цьому етапі не можна стверджувати, що α є оптимальним розв’язком. З іншого боку розглянемо базис А1 , А3 того ж опорного розв’язку α. Так само приведемо симплекс – таблицю до цього базису:
x1 |
x2 |
x3 |
x4 |
b |
1 0 |
– 1/2 1/2 |
0 1 |
3/2 – 1/2 |
1 0 |
0 |
– 1/2 |
0 |
– 33/2 |
– 10 |
Всі оцінки недодатні. Отже, α є оптимальним розв’язком цієї задачі, вектори А1 і А3 складають його базис.