КАТЕГОРИИ: Архитектура-(3434)Астрономия-(809)Биология-(7483)Биотехнологии-(1457)Военное дело-(14632)Высокие технологии-(1363)География-(913)Геология-(1438)Государство-(451)Демография-(1065)Дом-(47672)Журналистика и СМИ-(912)Изобретательство-(14524)Иностранные языки-(4268)Информатика-(17799)Искусство-(1338)История-(13644)Компьютеры-(11121)Косметика-(55)Кулинария-(373)Культура-(8427)Лингвистика-(374)Литература-(1642)Маркетинг-(23702)Математика-(16968)Машиностроение-(1700)Медицина-(12668)Менеджмент-(24684)Механика-(15423)Науковедение-(506)Образование-(11852)Охрана труда-(3308)Педагогика-(5571)Полиграфия-(1312)Политика-(7869)Право-(5454)Приборостроение-(1369)Программирование-(2801)Производство-(97182)Промышленность-(8706)Психология-(18388)Религия-(3217)Связь-(10668)Сельское хозяйство-(299)Социология-(6455)Спорт-(42831)Строительство-(4793)Торговля-(5050)Транспорт-(2929)Туризм-(1568)Физика-(3942)Философия-(17015)Финансы-(26596)Химия-(22929)Экология-(12095)Экономика-(9961)Электроника-(8441)Электротехника-(4623)Энергетика-(12629)Юриспруденция-(1492)Ядерная техника-(1748) |
Метод прогонки
Якщо звичайне диференціальне рівняння (ЗДР) 2-го порядку розв’язати методом скінченних різниць, то початкове ЗДР зводиться до системи лінійних алгебраїчних рівнянь з тридіагональною матрицею коефіцієнтів, тобто такою, що кожне рівняння системи має три сусідні невідомі: (10.11) Для розв’язання такої системи розроблено спеціальний метод – прогонки. Розглянемо лінійне диференціальне рівняння другого порядку (10.12) з двоточковими граничними умовами: (10.13) ; , (10.14) де функції , , неперервні на . Прийдемо до скінченнорізницевого рівняння ; ; ; ; . (10.15) Підставляючи вираз (10.15) у початкове ЗДР (10.12), отримаємо для внутрішніх точок систему скінченнорізницевих рівнянь: . (10.16) Після деяких перетворень (10.16) матимемо ; (10.17) де ; ; . (10.18) Для похідних на кінцях відрізка інтегрування і знаходимо скінченнорізницеві рівняння виду ; (10.19) і, підставляючи їх у граничні умови (10.13) і (10.14), отримаємо два рівняння: (10.20) Система рівнянь (10.17), (10.20) відносно невідомих являє собою систему лінійних алгебраїчних рівнянь з три діагональною матрицею коефіцієнтів. Цю систему зручно розв’язати методом прогонки. Розглянемо суть цього методу. 1. Представимо рівняння (10.17) відносно : . (10.21) 2. Допустимо, що за допомогою рівнянь (10.20) і (10.21) виключена складова , тоді рівняння будуть мати вигляд (10.22) де – деякі коефіцієнти. 3. За аналогією з (10.22) представимо у вигляді (10.23) і, підставляючи його в (10.17), матимемо (10.24) і, отже, . (10.25) 4. Порівнюючи вирази (10.22) та (10.25), отримуємо рекурентні формули для визначення коефіцієнтів , вигляду ; , . (10.26) 5. Для визначення виразу для та використовуємо рівняння (10.20) та (10.22): ; (10.27) та . (10.28) Порівнюючи останні дві нерівності, знаходимо:
; . (10.29) На основі формул (10.17) та (10.29) послідовно визначаємо коефіцієнти до і включно. Процес знаходження коефіцієнтів називають прямим ходом методу прогонки. Зворотній хід методу прогонки починається з визначення . Використовуючи другу граничну умову (10.20) та формулу (10.22), отримуємо систему з двох рівнянь: (10.30) Розв’язуючи цю систему відносно , отримаємо: . (10.31) Тепер за формулою (10.22) послідовно знаходимо розв’язок початкового ЗДР: Таким чином, метод прогонки достатньо простий, легко алгоритмізується та включає прямий та зворотній хід. Причому: 1) прямий хід полягає в знаходженні коефіцієнтів за рекурентними формулами (10.26); 2) зворотній хід полягає у знаходженні розв’язку початкового ЗДР за формулою (10.22). Схему алгоритму методу прогонки показано на рис. 10.1. Рисунок 10.1- Схему алгоритму методу прогонки
Дата добавления: 2014-11-08; Просмотров: 470; Нарушение авторских прав?; Мы поможем в написании вашей работы! Нам важно ваше мнение! Был ли полезен опубликованный материал? Да | Нет |