Теми рефератів
> Реферати > Курсові роботи > Звіти з практики > Курсові проекти > Питання та відповіді > Ессе > Доклади > Учбові матеріали > Контрольні роботи > Методички > Лекції > Твори > Підручники > Статті Контакти
Реферати, твори, дипломи, практика » Курсовые проекты » Чисельне рішення звичайних диференціальних рівнянь

Реферат Чисельне рішення звичайних диференціальних рівнянь





3.3)

Ця система має єдине рішення х =1; у =1. Тепер розглянемо пару значень невідомих х =2,415; у =0. При підстановці цих значень в вихідні рівняння отримуємо


(3.4)

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


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


Справа в тому, що дві прямі лінії, описувані двома рівняннями цієї системи, майже паралельні, як це показано на рис. 8 .. Точка х =2,415; у =0 хоча і не лежить на жодній з цих прямих ліній, але дуже близька до них.

Системи типу (3.3) називаються погано зумовленими . У будь-якому випадку, коли дві лінії (або площини і гіперплощини) майже паралельні, система рівнянь стає погано обумовленою. У цьому випадку знайти чисельне рішення системи важко, а точність його вельми сумнівна. Більш того, система з трьох і більше рівнянь може виявитися погано обумовленою, навіть якщо ніякі площини не є паралельними або майже паралельними 1). Незалежно від геометричної інтерпретації обумовленість системи можна досліджувати при її вирішенні методом Гаусса або за допомогою спеціальної перевірки.

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

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

Нехай задана система п лінійних алгебраїчних рівнянь з п невідомими:


. (3.5)


Система рівнянь (3.5) в матричній формі представляється наступним чином:

АХ=В, (3.6)

де А - квадратна матриця коефіцієнтів, розміром п ? п рядків і стовпців;

Х - вектор-стовпець невідомих;

В - вектор-стовпець правих частин.

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

Недоліком методу є накопичення похибок в процесі округлення, тому метод Гаусса без вибору головних елементів використовується зазвичай для вирішення порівняно невеликих ( п ? 100) систем рівнянь з щільно заповненою матрицею і не близьким до нуля визначником. Якщо матриця А сильно розріджена, а її визначник неблизький до нуля, то метод Гаусса придатний для вирішення великих систем рівнянь. У MATLAB є обширний арсенал методів розв'язання систем рівнянь (3.5) методом виключення Гауса. Для цього застосовуються такі оператори дають вирішення ряду систем лінійних рівнянь АХ=В, де А - матриця розміром m ? n , В - матриця розміром п ? до .

Вирішити систему лінійних рівнянь:



Програма обчислення в системі Matlab наведена в Додатку 3. Результат рішення:


x=1.000000000000000

2.000000000000001

3.000000000000001

4.000000000000000=1.000000000...


Назад | сторінка 6 з 9 | Наступна сторінка





Схожі реферати:

  • Реферат на тему: Рішення систем лінійних алгебраїчних рівнянь методом Гаусса
  • Реферат на тему: Рішення систем лінійних рівнянь за методом Гаусса
  • Реферат на тему: Рішення систем лінійних алгебраїчних рівнянь методом Гауса
  • Реферат на тему: Рішення системи двох лінійних рівнянь з поданням про вирішення в числовому ...
  • Реферат на тему: Метод Гаусса розв'язання систем лінійних алгебраїчних рівнянь