і отримати безкоштовне
свідоцтво про публікацію
! В а ж л и в о
Предмети »

Метод послідовного виключення змінних

Перегляд
матеріалу
Отримати код Поділитися

Метод послідовного виключення змінних

Найпростішим методом розв’язання систем лінійних алгебраїчних рівнянь є метод послідовного виключення змінних, або метод Гауса. Є кілька модифікацій цього методу. Розглянемо схему єдиного ділення, за якою систему розв’язують в два етапи. На першому етапі вихідну систему рівнянь зводять до рівносильної їй системи трикутної форми. Цей процес перетворення називають прямим ходом. На другому етапі, який називають зворотним ходом, знаходять розв’язок лінійної системи рівнянь трикутної форми.

Обмежимося розглядом системи трьох рівнянь з трьома змінними

(1)

визначник якої не дорівнює нулю.

Нехай a (цього завжди можна досягти перестановкою рівнянь системи, бо визначник не дорівнює нулю, і тому завжди є таке рівняння, в якому a). Поділимо коефіцієнти першого рівняння системи (1), включаючи й вільний член, на коефіцієнт a. Дістанемо нове рівняння

++=, (2)

де

=/, (=) (3).

Виключимо тепер змінну з другого і третього рівнянь системи(1). Для цього рівняння (2) помножимо послідовно спочатку на коефіцієнт і віднімемо його від другого рівняння системи (1), а потім на a і віднімемо від третього рівняння системи (1). Дістанемо систему двох рівнянь з двома змінними і

(4)

де коефіцієнти обчислюють за формулами

=- (). (5)

Верхній індекс (1) вказує на те, що над коефіцієнтами системи (4) виконано перше перетворення.

Далі поділимо коефіцієнти першого рівняння системи (4) на (якщо =0, то переставляємо рівняння місцями). Дістанемо рівняння

+=, (6)

де

=/, (). (7)

Із системи (4) виключимо змінну так само, як і з системи (1). Дістанемо рівняння

x=, (8)

де

=- (). (9)

З рівняння (8) маємо

= , (10)

де

=/. (11)

Після трьох кроків перетворення дістанемо систему рівнянь трикутної форми

(12)

яка еквівалентна системі (1).

На цьому прямий хід методу Гаусса завершено. Описаний процес перетворень системи (1) до рівносильної ій системи (12) можна здійснити, якщо виконуються умови 0, 0, 0. Можна довести, що коли 0, то для цього достатньо, щоб поряд з умовою визначник 0 виконувалась нерівність

Близькість діагональних коефіцієнтів системи (1) до нуля може призвести до значної втрати точності.

Оскільки системи (1) і (12) еквівалентні, то розв’язком системи (1) буде розв’язок системи(12), який можна записати формулами

(13)

Цим завершено зворотний хід методу Гаусса.

Оскільки перетворення систем рівнянь є фактично перетвореннями їх коефіцієнтів при змінних і вільних членів, то для виконання перетворень немає потреби виписувати системи. Досить виписати лише матрицю коефіцієнтів і вільні члени і над ними виконати зазначені перетворення. Всі записи доцільно розміщувати в окремій таблиці, яка складається з кількох частин, що відповідають певним крокам перетворення вихідної системи рівнянь.

Кроки перетв-

орень

рядок

Коефіцієнти при змінних

Вільний член

Контроль

Х1

Х2

Х3

Контрольна сума

Рядкова сума

1

2

3

4

5

6

7

8

1

1

2

3

2

4

5

6

1

3

7

8

1

4

9

10

11

1

1

1

У процесі розв’язування системи рівнянь треба організувати контроль правильності обчислень. Щоб вчасно виявити (і виправити) випадкові обчислювальні помилки, доцільно забезпечити контроль правильності обчислень у кожному рядку таблиці (так званий поточний контроль). Для цього до схеми обчислень введено два додаткових стовпців: 7-й- контрольна сума і 8-й- рядкова сума.

Контрольна сума () – це сума коефіцієнтів при змінних і вільного члена для кожного рівняння системи (1)

= (). (14)

У процесі розв’язування системи (1) за схемою єдиного ділення над контрольними сумами виконують ті самі перетворення, що й над відповідними елементами рядка, а знайдені в результаті цих перетворень значення контрольної суми записують у тому рядку, що й нові елементи відповідного рядка.

Рядкові суми – це суми всіх елементів відповідного рядка, що містяться в стовпцях з 3-го по 6-й включно. Поточний контроль обчислень за схемою єдиного ділення полягає в тому, що в кожному рядку обчислюють контрольну і рядкові суми та порівнюють їх між собою. Якщо вони збігаються або відрізняються на 1-2 одиниці нижчого розряду, що обумовлено впливом похибок округлення, то обчислення виконано правильно і можна переходити до обробки наступного рядка. Якщо ці суми в будь –якому рядку значно відрізняються між собою, то обчислення треба припинити і зясувати причину розбіжності. Значна розбіжність між контрольною і рядковими сумами може свідчити про випадкові помилки в обчисленнях або про нестійкість алгоритму обчислень відносно похибок округлень. Якщо причиною розбіжності є випадкова помилка, то її виправляють і продовжують обчислення. А якщо причиною розбіжності сум є нестійкість алгоритму обчислень відносно похибок округлення, то від нього слід відмовитися.

Зазначений спосіб організації контролю за обчисленнями по суті означає, що одночасно розв’язується дві системи лінійних алгебраїчних рівнянь з однаковою матрицею коефіцієнтів при змінних, але з різними вільними членами. Вільними членами вихідної системи (1) є числа , , , а допоміжної системи - , , .

Розв’язок , , системи (1) зв’язаний з розв’язком , , допоміжної системи

(). (15)

простою залежністю

(). (16)

Справді, підставивши (16) і врахувавши формули (14), дістанемо тотожність

().

Під час виконання зворотного ходу методу Гауса одночасно обчислюють як розв’язок системи (1), так і розв’язок допоміжної системи (15). Ці розв’язки в межах заданої точності повинні задовольняти рівність (16). Саме в цьому й полягає суть заключного контролю.

Другою формою заключного контролю є безпосередня перевірка знайденого розв’язку підстановкою його і систему (1). Якби всі обчислення виконувалися точно, то в результаті підстановки дістали б правильну числову рівність. Але в процесі обчислень виконувалися округлення, тому значення лівих частин рівнянь системи (1), взагалі кажучи, не збігатимуться із значенням їх правих частин. Значення різниць між вільними членами вихідної системи лінійних рівнянь і результатами підстановки в ці рівняння обчислених значень змінних називають нев’язками. Якщо невязки досить малі, то можна стверджувати, що розв’язок системи (1) знайдено з малими похибками. Якщо невязки досить значні, то це означає, що значення шуканих змінних обчислено з недостатньою точністю і їх треба уточнити. Це буває здебільшого тоді, коли проміжні обчислення виконують з недостатньою точністю. Зменшити невязки можна так: розв’язують систему повторно, залишаючи в проміжних результатах більшу кількість десяткових розрядів, ніж при попередньому розв’язуванні, або обчислюють значення поправок до знайденого раніше розв’язку системи. Перший шлях досить громіздкий, потребує тим більших затрат обчислювальної роботи, чим вищий порядок системи. Тому перевагу доцільно надати другому шляху.

Нехай - наближений розв’язок системи (1). Підставивши його в систему (1), дістанемо невязка

(17)

Шукатимемо тепер новий розв’язок системи (1) в такому вигляді

,

, (18)

,

де , , - шукані похибки розв’язку. Підставивши (18) в систему (1) і взявши до уваги невязки (17), для знаходження поправок (j=1,2,3) дістанемо систему рівнянь

(). (19)

Система (19) відрізняється від системи (1) лише вільними членами; її розвязок можна знайти за схемою єдиного ділення. Для цього в нашій таблиці досить добавити новий стовпець 9, елементами якого будуть невязки (17), тобто стовпець вільних членів системи (19). Значення поправок , , будуть знайдені, якщо над невязками виконані ті самі перетворення, що й над вільними членами системи (1). Обсяг обчислювальної роботи при цьому значно зменшується, бо перетворення виконуються лише над елементами одного (9-го) стовпця таблиці.

Відображення документу є орієнтовним і призначене для ознайомлення із змістом, та може відрізнятися від вигляду завантаженого документу

Опис документу:
У цьому документі йде мова про метод послідовного виключення змінних.

Зверніть увагу, свідоцтва знаходяться в Вашому особистому кабінеті в розділі «Досягнення»

Курс:«Інтегрований курс «Мистецтво»: теоретико-прикладний аспект»
Просіна Ольга Володимирівна
36 годин
590 грн
295 грн

Бажаєте дізнаватись більше цікавого?


Долучайтесь до спільноти