ЧИСЕЛЬНІ МЕТОДИ РОЗВ\'ЯЗУВАННЯ СИСТЕМ ЛІНІЙНИХ АЛГЕБРАЇЧНИХ РІВНЯНЬ.. Лабораторна робота з Чисельні методи в інформатиці. Робота № 201001

Перехід до торгівельного партнера Binance

ЧИСЕЛЬНІ МЕТОДИ РОЗВ'ЯЗУВАННЯ СИСТЕМ ЛІНІЙНИХ АЛГЕБРАЇЧНИХ РІВНЯНЬ.

Інформація про навчальний заклад

ВУЗ:

Національний університет Львівська політехніка

Інститут:

Не вказано

Факультет:

Не вказано

Кафедра:

Системи автоматизованого проектуваня

Інформація про роботу

Рік:

2008

Тип роботи:

Лабораторна робота

Предмет:

Чисельні методи в інформатиці

Група:

КН-3

Завантажити

Частина тексту файла (без зображень, графіків і формул):

Міністерство освіти і науки України Національний університет «Львівська політехніка» кафедра САПР Лабораторна робота №4 з курсу "Чисельні методи в інформатиці" на тему: ЧИСЕЛЬНІ МЕТОДИ РОЗВ'ЯЗУВАННЯ СИСТЕМ ЛІНІЙНИХ АЛГЕБРАЇЧНИХ РІВНЯНЬ Виконав: cт. гр. КН-3 Львів-2008 1. МЕТА РОБОТИ Мета роботи – ознайомитись із чисельними методами розв'язування систем лінійних алгебраїчних рівнянь та їх практичним застосуванням. 2. ТЕОРЕТИЧНА ЧАСТИНА Методи розв'язування систем лінійних рівнянь можна поділити на два типи: прямі, або точні, та ітераційні. Прямі методи дають можливість дістати розв'язок, виконавши скінченну апріорі відому кількість операцій. Якщо всі проміжні обчислення виконувати точно (без заокруглень), то отримаємо точний розв'язок. Ітераційні методи дають нескінченну послідовність наближених розв'язків, границі яких є розв'язком системи. 2.1. Прямі методи розв'язування систем лінійних алгебраїчних рівнянь Із точних методів розв'язування систем лінійних рівнянь відомий з курсу лінійної алгебри метод Гаусса, реалізація якого за схемою єдиного ділення вимагає виконання n(4n2-3n-4)/6 арифметичних операцій. Розглянемо розв'язування систем лінійних рівнянь методом Гаусса за схемою Халецького, для реалізації якої необхідно n(2n-1) арифметичних операцій. Система лінійних алгебраїчних рівнянь в матричному представленні має вигляд: (1) Матрицю коефіцієнтів А подаємо у вигляді добутку двох трикутних матриць, тобто: А=СВ. (2) Матриця С – нижня трикутна, В – верхня трикутна, причому діагональні коефіцієнти матриці В дорівнюють одиниці. Підставимо (2) в (1): CBx=d. (3) Позначимо Bx=y і перепишемо (3) у вигляді: Cy=d (4) Розв'язок (4) називається прямим ходом методу Гаусса, розв'язок (3) – зворотним ходом. Отже, першим кроком при розв'язуванні сформульованої задачі є розбиття матриці А на дві трикутні матриці С і В. Коефіцієнти матриць С і В обчислюються послідовно. Спочатку обчислюється перший стовпець матриці С: , тоді обчислюється перший рядок матриці В: і елемент , після чого – елементи другого стовпця матриці С: тоді обчислюються елементи другого рядка матриці В: і елемент: і т.д. Після того, як будуть обчислені елементи матриць С і В, а також вектор y, знаходять невідомі : . Алгоритм методу Гаусса за схемою Халецького має такий компактний вигляд: 2.2.1. Метод простої ітерації Найпростішим ітераційним методом розв'язування систем лінійних рівнянь ах=d (5) є метод простої ітерації. Система рівнянь (5) перетворюється до вигляду: x=вx+e, (6) де , . (7) Ітераційний процес записується у вигляді: . (8) Для збіжності методу простих ітерацій при довільному початковому наближенні ітераційного процесу (8) необхідно і достатньо, щоб всі власні значення матриці B були за модулем менші від одиниці. Обчислювальна схема методу простої ітерації: 1. Обчислюємо елементи матриці B і вектора e за (7). 2. Обчислюємо . Позначаємо: , . Якщо C<1, то ітераційний процес збігається. 3. Як початкове наближення вибираємо вектор e:=e. Позначимо . Розв'язуємо систему рівнянь:. Тоді:. 4. Обчислюємо необхідну кількість ітерацій для досягнення потрібної точності розв'язку: . Звідси . (9) 2.2.2. Метод Зейделя Ітераційний метод Зейделя відрізняється від методу простої ітерації тим, що на (k+1)-й ітерації при обчисленні і-ї компоненти вектора використовуються значення , обчислені на цій ітерації. Обчислювальна схема методу Зейделя: Для реалізації методу Зейделя необхідно: 1. Перевірити достатню умову збіжності:. (10) Якщо умова (10) виконується для всіх і, то обчислюємо елементи матриці B і вектора e за (7). 2. Обчислити норму матриці В: . 3. Вибрати як початкове наближення . Визначити необхідну для досягнення потрібної точності кількість ітерацій К за (9). 4. Представити ітераційний процес у вигляді: ; і реалізувати його К разів. Лабораторне завдання (Варіант №11) ---- Метод Гаус ---- Вхідні дані: ---- Матриця А ----- ,110 1,13 -,170 ,180 ,130 -1,17 -,180 ,140 ,110 -1,05 -1,70 - ,150 ,150 -,500Е-01 ,180 -,110 --- Матриця Dt ------ 1.00 .130 .110 1.00 Результати: X( 1)= 5.7217 X( 2)= .51589 X( 3)= .49623 X( 4)= -.71106 ---- Метод Гаус (Схема Халецького)---- Вхідні дані: ---- Матриця А ----- ,110 1,13 -,170 ,180 ,130 -1,17 -,180 ,140 ,110 -1,05 -1,70 - ,150 ,150 -,500Е-01 ,180 -,110 --- Матриця Dt ------ 1.00 .130 .110 1.00 Результати: X( 1)= 5.7217 X( 2)= .51589 X( 3)= .49623 X( 4)= -.71106 ---- Метод простої ітерації розв’язку СЛАР ---- Вхідні дані: ---- Матриця А ----- ,370 ,300Е-01 -,160 ,800Е-01 ,110 1,13 -,170 ,180 ,220 ,800Е-01 -,340 - ,300Е-01 ,240 -,400Е-01 ,100Е-01 -,320 --- Матриця Dt ------ 2.11 1.00 1.11 1.13 Результати: Кількість ітерацій - 253 Корені: X( 1)= 5.7217 X( 2)= .51589 X( 3)= .49623 X( 4)= -.71106 Вхідні дані: ---- Матриця А ----- ,370 ,300Е-01 -,160 ,800Е-01 ,110 1,13 -,170 ,180 ,220 ,800Е-01 -,340 - ,300Е-01 ,240 -,400Е-01 ,100Е-01 -,320 --- Матриця Dt ------ 2.11 1.00 1.11 1.13 Результати: Кількість ітерацій - 407 Корені: X( 1)= 5.7217 X( 2)= .51589 X( 3)= .49623 X( 4)= -.71106 Вхідні дані: ---- Матриця А ----- ,370 ,300Е-01 -,160 ,800Е-01 ,110 1,13 -,170 ,180 ,220 ,800Е-01 -,340 - ,300Е-01 ,240 -,400Е-01 ,100Е-01 -,320 --- Матриця Dt ------ 2.11 1.00 1.11 1.13 Результати: Кількість ітерацій - 562 Корені: X( 1)= 5.7217 X( 2)= .51589 X( 3)= .49623 X( 4)= -.71106 ---- Метод Зейделя розв’язку СЛАР ---- Вхідні дані: ---- Матриця А ----- ,370 ,300Е-01 -,160 ,800Е-01 ,110 1,13 -,170 ,180 ,220 ,800Е-01 -,340 - ,300Е-01 ,240 -,400Е-01 ,100Е-01 -,320 --- Матриця Dt ------ 2.11 1.00 1.11 1.13 Результати: Кількість ітерацій - 253 Корені: X( 1)= 5.7217 X( 2)= .51589 X( 3)= .49623 X( 4)= -.71106 Вхідні дані: ---- Матриця А ----- ,370 ,300Е-01 -,160 ,800Е-01 ,110 1,13 -,170 ,180 ,220 ,800Е-01 -,340 - ,300Е-01 ,240 -,400Е-01 ,100Е-01 -,320 --- Матриця Dt ------ 2.11 1.00 1.11 1.13 Результати: Кількість ітерацій - 407 Корені: X( 1)= 5.7217 X( 2)= .51589 X( 3)= .49623 X( 4)= -.71106 Вхідні дані: ---- Матриця А ----- ,370 ,300Е-01 -,160 ,800Е-01 ,110 1,13 -,170 ,180 ,220 ,800Е-01 -,340 - ,300Е-01 ,240 -,400Е-01 ,100Е-01 -,320 --- Матриця Dt ------ 2.11 1.00 1.11 1.13 Результати: Кількість ітерацій - 562 Корені: X( 1)= 5.7217 X( 2)= .51589 X( 3)= .49623 X( 4)= -.71106 Висновок: У даній лабораторній роботі ми ознайомились із чисельними методами розв'язування систем лінійних алгебраїчних рівнянь та їх практичним застосуванням.

Лабораторна робота Чисельні методи в інформатиці

01.01.1970 03:01-

Коментарі

Ви не можете залишити коментар. Для цього, будь ласка, або зареєструйтесь.

Ділись своїми роботами та отримуй миттєві бонуси!

Маєш корисні навчальні матеріали, які припадають пилом на твоєму комп'ютері? Розрахункові, лабораторні, практичні чи контрольні роботи — завантажуй їх прямо зараз і одразу отримуй бали на свій рахунок! Заархівуй всі файли в один .zip (до 100 МБ) або завантажуй кожен файл окремо. Внесок у спільноту – це легкий спосіб допомогти іншим та отримати додаткові можливості на сайті. Твої старі роботи можуть приносити тобі нові нагороди!

поділитись