12467

Прямі методи розв’язання систем лінійних алгебраїчних рівнянь. Метод Гаусса та LU-розкладу

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

Информатика, кибернетика и программирование

Лабораторна робота №1 Прямі методи розвязання систем лінійних алгебраїчних рівнянь. Метод Гаусса та LUрозкладу. Мета роботи: ознайомитися з методами розв‘язання систем лінійних алгебраїчних рівнянь. Розглянути особливості реалізації прямих методів розв‘язання ...

Украинкский

2013-04-27

56.5 KB

19 чел.

Лабораторна робота №1

Прямі методи розв’язання систем

лінійних алгебраїчних рівнянь. Метод Гаусса та LU-розкладу.

Мета роботи: ознайомитися з методами розв‘язання систем лінійних алгебраїчних рівнянь. Розглянути особливості реалізації прямих методів розв‘язання систем лінійних алгебраїчних рівнянь у вигляді m-файлу функції у середовищі MatLab.

Задачі лабораторної роботи: реалізувати один з прямих методів розв‘язання систем лінійних алгебраїчних рівнянь у вигляді m-файлу функції у середовищі MatLab у відповідності з варіантом. Продемонструвати його роботу, оцінити точність отриманих розв‘язків та визначити джерела виникнення похибок на прикладі розв‘язання конкретної систем лінійних алгебраїчних рівнянь.

Теоретичні відомості.

Математичні моделі багатьох технічних задач представлені системами лінійних рівнянь. Багато методів розв’язання нелінійних задач також зводяться до розв’язання деякої послідовності систем лінійних алгебраїчних рівнянь (систем ЛАР) на ЕОМ. Для багатьох методів розроблений математичний апарат, що дозволяє оцінити точність отриманого розв’язку. Чисельні методи розв’язку систем ЛАР поділяються на прямі та ітераційні (наближені).

Прямі (точні) методи дозволяють розв’язати систему рівнянь за скінчене число арифметичних операцій. Якщо всі операції виконуються точно (без помилок округлення), то розв’язок заданої системи також отримуємо точним. До прямих методів належать: метод послідовного виключення невідомих (метод Гаусса та його модифікації: метод головного елемента, метод квадратного кореня, метод відображень та ін.), метод ортогоналізації, метод LU-розкладу. Прямі методи застосовують на практиці для розв’язання систем ЛАР за допомогою обчислювальної техніки, як правило, з числами порядку не вище 103 .

Ітераційні методи є наближеними. Вони дозволяють знайти розв’язок системи, як межу послідовних наближень, що обчислюються по однаковому алгоритму. Для застосування ітераційних методів у початкових умовах необхідно задати точність обчислень ε і початкове наближення х0 чи (х01, х02, х03, …). До ітераційних методів належать: метод Зейделя, метод простої ітерації, метод релаксації, градієнтні методи та їх модифікації. На практиці ітераційні методи застосовують для розв’язання систем ЛАР з числами порядку 106 і вище.

Розглянемо систему m лінійних алгебраїчних рівнянь з n невідомими:

              (1.1)

Коротко її можна записати в матричному вигляді:

     (1.2)

де матриця m×n;  - вектор n-го порядку; - вектор m-го порядку.

Розв’язком системи ЛАР називається така впорядкована сукупність чисел х1= с1, х2 = с2, … хn = сn, яка обертає всі рівняння системи у істинні рівняння.

Система ЛАР називається сумісною, якщо вона має хоча б один розв’язок, і несумісною – якщо вона не має розв’язків. Сумісна система називається визначеною, якщо вона має один розв’язок, і невизначеною, якщо має більше одного розв’язку.

Метод Гаусса

Розглянемо систему ЛАР:

    а11х1 + а12х2 + …+а1nхn = b1

    а21х1 + а22х2 + …+а2nхn = b2

    ……………………………………………                                          (1.3)

    аn1х1 + аn2х2 + …+аnnхn = bn

у якій матриця А = (аij) не вироджена. Метод Гаусса полягає у послідовному виключенні невідомих. Суть його у перетворенні системи (1.3) в систему з трикутною матрицею, з якої послідовно (при зворотному розв‘язку, тобто від останнього рівняння до першого) отримують значення всіх невідомих.

Алгоритм послідовного виключення невідомих може бути побудований за різними обчислювальними схемами. Побудуємо його по схемі єдиного ділення.

Перетворимо систему (1.3): припустимо а11  0, розділимо всі коефіцієнти першого рівняння на а11. Віднімемо з кожного рівняння системи перше рівняння, що помножене на коефіцієнт при х1. На перетворюючи першого рівняння, зробимо аналогічні перетворення над всіма іншими рівняннями системи і т.д. Отримуємо систему з трикутною матрицею:

   х1 + 12х2 + …+1nхn = 1

          22х2 + …+2nхn = 2

   ……………………………………………                            (1.4)

                                   хn = n

Із системи рівнянь (1.4) послідовно знаходимо значення всіх невідомих  хn, хn-1,…, х1.

Таким чином процес розв‘язання системи (1.3) по методу Гаусса можна розділити на два етапи.

1) Перший етап полягає в послідовному виключенні невідомих (приведенні матриці коефіцієнтів А = (аij) до верхньої трикутної матриці). Його називають прямим ходом.

2) Другий етап – безпосереднє отримання значень невідомих, як результату розв’язання низки рівнянь однією змінною – називають зворотнім ходом. Таку назву етап отримав, оскільки розв‘язання починають із n-го останнього рівняння, послідовно проводячи розрахунки з n-1, n-2, n-3…3, 2, 1 рівняннями.

В ході виконання операцій ділення в методі Гауса виникає обчислювальна похибка, яка при збільшенні кількості операцій (≥40) може значно спотворювати результат. Для зменшення обчислювальної похибка застосовують перестановки рівнянь у системі в залежності від головного елементу. Головний елемент – це максимальний коефіцієнт рівняння. Сукупність рівнянь в системі потрібно переставити таким чином, щоб на головній діагоналі матриці коефіцієнтів стояли максимальні елементи.

Метод LU-розкладу

 При розв‘язанні системи лінійних алгебраїчних рівнянь цим методом матрицю коефіцієнтів А розкладають на добуток двох матриць нижньої трикутної матриці L, на головній діагоналі якої стоять одиниці, та верхньої трикутної матриці U, елементи головної діагоналі якої не дорівнюють нулю. У матричному вигляді при n = 4 це можна записати таким чином:

 (1.5)

Для розкладання матриці коефіцієнтів А на трикутні матриці використаємо метод виключення Гауса. Отримаємо матрицю L з допомогою одиничної матриці. Для цього помножимо зліва матрицю А на одиничну матрицю. Наприклад:

   (1.6)

Перший рядок матриці А використовуємо для обнулення елементів першого стовпчика цієї ж матриці. Помножимо перший рядок на –0,5 і віднімемо від другого, помножимо перший рядок на 0,25 і віднімемо від третього рядка. Вказану операцію проводимо не тільки з матрицею А, а й з одиничною матрицею, але замість віднімання першого рядка від чергового здійснюємо додавання першого рядка до чергового. В результаті отримаємо:

   (1.7)

Другий рядок матриці А використовуємо для обнулення елементів другого стовпчика цієї ж матриці, що знаходяться нижче другого рядка (у матриці, що утворена з одиничної, проводимо додавання відповідного елементу):

  (1.8)

В результаті отримуємо LU - розклад матриці А і розв’язуємо систему рівнянь у два етапи:

   (1.9)

На першому етапі знаходимо проміжний вектор Y, використовуючи пряму підстановку, на другому етапі знаходимо безпосередньо вектор розв’язків Х, застосовуючи зворотну підстановку. Зокрема, приймаючи в даному прикладі з проміжної системи LY = B матимемо . Здійснюючи другий етап розв‘язання системи, тобто розв‘язуючи рівняння UX = Y , отримаємо .

 Завдання для виконання лабораторної роботи:

 Створити програму на внутрішній мові середовища МatLAB, що реалізує метод Гауса з вибором головного елементу (для непарних варіантів) чи метод LU-розкладу (для парних варіантів). Провести тестування створеної програми на прикладі, вибраному за варіантом. Здійснити перевірку отриманих результатів, в разі виникнення похибки, пояснити джерело її виникнення та накопичення (компенсації).

Варіанти завдань:

1 варіант.  2 варіант.

3 варіант.  4 варіант.

5 варіант.  6 варіант.

7 варіант.  8 варіант.

9 варіант.  10 варіант.


 

А также другие работы, которые могут Вас заинтересовать

24763. Особенности правовой защиты персональных данных 134.5 KB
  Особенности правовой защиты персональных данных. Эти процессы стимулируют создание системы правовой защиты персональных данных. Персональные данные любая информация относящаяся к прямо или косвенно определенному или определяемому физическому лицу субъекту персональных данных;ФЗ 152 Государственный надзор за выполнением требований законодательства в области защиты ПДн распределен между тремя ведомствами: 1 Роскомнадзор основной исполнительный и надзорный орган по защите прав физических лиц чьи персональные данные обрабатываются; 2...
24764. Конституционные основы организации публичной власти в России 62 KB
  Государственную власть в РФ осуществляют Президент РФ Федеральное Собрание парламент состоит из двух палат: Совета Федерации и Государственной Думы Правительство РФ суды РФ. Признание человека его прав и свобод высшей ценностью является фундаментальной нормой конституционного строя Российской Федерации образующей основу не только конституционно организованного общества но и правовой защиты этого общества от попыток подавления личности ущемления ее прав игнорирования индивидуальных интересов и потребностей людей. Экономической...
24765. Социальное управление: специфика, методы, подсистемы, функции 111.5 KB
  Содержание системы управления и внутренние принципы ее организации и самоорганизации определяет целеполагание связанное с прогнозированием и моделированием социальных систем. Организационноадминистративное воздействие осуществляется в следующих основных видах: прямое административное указание; установление правил регулирующих деятельность подчиненных; разработка и внедрение рекомендаций по организации и совершенствованию тех или иных процессов; контроль и надзор за деятельностью организаций и отдельных работников. Функция целеполагания...
24766. Понятие и содержание трудового договора 123.5 KB
  15 Трудового кодекса РФ под трудовыми отношениями понимаются отношения основанные на соглашении между работником и работодателем о личном выполнении работником за плату трудовой функции то есть работы по определенной специальности квалификации или должности подчинении работника правилам внутреннего трудового распорядка при обеспечении надлежащих условий труда а также своевременной и в полном объеме выплате ему заработной платы. Трудовой договор это соглашение между работодателем и работником в соответствии с которым работодатель...
24767. Государственное управление как объект АП регулирования 86.5 KB
  Государственное управление это один из видов деятельности государственных органов РФ по осуществлению государственной власти реализации государственновластных полномочий. Ветви власти взаимосвязаны но в то же время характеризуются определенной самостоятельностью. В широком смысле государственное управление осуществляется всеми ветвями власти. Система органов исполнительной власти в РФ.
24768. Сущность, содержание и структура системы национальной безопасности государства 42.5 KB
  Под национальной безопасностью страны необходимо понимать систему элементов связей и отношений обеспечивающую реализацию жизненно важных политических экономических военных гуманитарных экологических информационных и других интересов личности общества и государства. Национальные интересы объективные потребности гражданина общества и государства вытекающие из особенностей социальноэкономического и политического устройства государства уровня его экономического развития исторически сложившегося места в международном разделении...
24769. Гражданские правоотношения: понятие, элементы, виды 109.5 KB
  Гражданские правоотношения: понятие элементы виды Гражданские правоотношения общественные отношения урегулированные нормами ГП это связь субъектов наделенных взаимными правами и обязанностями.Субъекты: лица обладающие гражданскими правами и несущие гражданские обязанности в связи с участием в конкретном гражданском правоотношении.Объекты определенные цели на достижение которых направлены те или иные права. характер санкций применяемый в гражданском праве: Меры принуждения имеют имущественный характер и санкции являются...
24770. Экономическое понятие собственности и ее основные формы 213.5 KB
  В каждой стране в экономике имеется государственный сектор экономики что представляет общую черту современного рыночного хозяйства. В России предприятия государственного сектора экономики подразделяются на федеральные и муниципальные. Можно выделить следующие типы экономических систем: система современной рыночной экономики; система смешанной экономики; система традиционной нерыночной экономики; система административнокомандной экономики. Система современной рыночной экономики сегодня является доминирующей по причине высокой...
24771. Финансово-кредитная система РФ, ее структура и принципы формирования 197.5 KB
  Распределительная функция государственных и муниципальных финансов заключается в том что через распределение и перераспределение вновь созданной стоимости обеспечиваются общегосударственные потребности формируются источники финансирования общественного сектора экономики достигается сбалансированность бюджетов и внебюджетных фондов в рамках единой бюджетной системы РФ. Особую роль играет процесс перераспределения доходов между различными уровнями бюджетов. Бюджетная система Российской Федерации состоит из бюджетов трех уровней:...