- •Глава 1. Интерполирование
- •Постановка задачи интерполирования
- •1.2. Единственность интерполяционного многочлена
- •. Интерполяционный многочлен Лагранжа
- •1.4. Погрешность интерполяционной формулы Лагранжа
- •. Конечные разности и их свойства
- •1.6. Первая интерполяционная формула Ньютона
- •1.7. Вторая интерполяционная формула Ньютона
- •1.8. Интерполяционная формула Гаусса
- •1.10. Численное дифференцирование
- •Глава 2. Численное интегрирование
- •2.1. Общие замечания
- •Глава 3. Численные методы решения алгебраических и трансцендентных уравнений
- •3.1. Вводные замечания
- •3.2. Отделение корней
- •3.3. Метод половинного деления
- •3.5. Метод Ньютона (метод касательных)
- •3.6. Комбинированный метод
- •3.7. Метод итерации (метод последовательных приближений)
- •Глава 4. Решение систем линейных уравнений
- •4.1. Метод Гаусса (схема единственного деления)
- •4.2. Вычисление определителей
- •4.3. Вычисление обратной матрицы
- •4.4. Некоторые сведения из линейной алгебры
- •Рассмотрим систему линейных алгебраических уравнений
- •Вычитая второе равенство из первого, получаем
- •4.6. Приведение линейной системы к виду,
- •4.7. Метод Зейделя Пусть система линейных алгебраических уравнений приведена к виду
- •Глава 5. Обработка результатов наблюдений
- •5.1. Постановка задачи
- •5.2. Приближение функции, заданной таблично,
- •Глава 6. Численные методы решения
- •6.1. Методы решения задачи Коши. Вводные замечания
- •6.2. Решение дифференциальных уравнений
- •6.3. Метод Эйлера Пусть дано дифференциальное уравнение
- •6.4. Уточненный метод Эйлера Рассмотрим дифференциальное уравнение
- •6.6. Методы Рунге – Кутта Пусть дано дифференциальное уравнение
- •6.7. Метод Милна
- •6.8. Приближенное решение систем дифференциальных
- •6.9.Краевые задачи для обыкновенных дифференциальных
- •6.10. Решение краевых задач
- •6.11. Метод прогонки Рассмотрим линейную систему уравнений (6.54). Преобразуя первые (n-1) уравнения:
- •Глава 7. Численное решение дифференциальных уравнений в частных производных
- •7.1. Классификация дифференциальных уравнений
- •7.2. Метод сеток решения краевых задач
- •7.3. Погрешность аппроксимации дифференциальных уравнений
- •7.4. Аппроксимация граничных условий
- •7.6. Метод сеток для уравнений
6.11. Метод прогонки Рассмотрим линейную систему уравнений (6.54). Преобразуя первые (n-1) уравнения:
или
.
Отсюда
, (6.55)
где
(6.56)
Предположим, что с помощью полной системы (6.55) из i-го уравнения системы исключен член, содержащий . Тогда это уравнение можно записать в виде
(6.57)
Аналогично
.
Подставляя выражение для в i-е уравнение системы (6.55), получим
.
Отсюда
. (6.58)
Сравнивая (6.57) с (6.58), получим
. (6.59)
Получены рекуррентные соотношения для определения и . Для начала подсчетов необходимо знать и . Рассмотрим систему уравнений:
(получено из (6.55) при i=1);
(получено из первого краевого условия).
Находим у0 из второго уравнения системы:
(6.60)
и, подставляя его в первое уравнение системы, получаем
.
Сравнивая полученное выражение с выражением
,
получаем
; (6.61)
Итак, зная и , по формулам (6.59) можно последовательно вычислить коэффициенты и для . Вычисление коэффициентов и называется прямым ходом.
О пределим yn из системы уравнений:
( получено из (6.57) при i=n-1);
(получено из второго краевого условия);
. (6.62)
Теперь, используя формулы (6.57) и (6.60), последовательно находим:
, , …, ;
……………………………
.
Вычисление yi называется обратным ходом.
Вычисления удобно оформить в виде табл. 6.4.
Таблица 6.4
|
|
|
|
|
Прямой ход |
Обратный ход |
|
|
|
|
|||||
0 |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
… |
… |
… |
… |
… |
… |
… |
… |
n-1 |
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
Глава 7. Численное решение дифференциальных уравнений в частных производных
7.1. Классификация дифференциальных уравнений
в частных производных
Дифференциальные уравнения в частных производных имеют широкое применение в различных областях науки и техники. Получить их решение в явном виде удается только в самых простейших случаях. В связи с этим важное значение приобретают приближенные методы решения различных задач для дифференциальных уравнений в частных производных. На практике чаще всего используются линейные дифференциальные уравнения второго порядка.
Пусть D – некоторая область изменения независимых переменных х и у, ограниченная контуром Г. в области задано линейное дифференциальное уравнение второго порядка для функции U(x, y), если для любой точки из области D имеет место соотношение
. (7.1)
где - коэффициенты;
f(x,y) – свободный член уравнения.
Если коэффициенты a, b, c, d, e, g не зависят от х и у, то уравнение (7.1) является линейным дифференциальным уравнением с постоянными коэффициентами.
Обозначим
.
В зависимости от знака функции линейное дифференциальное уравнение (7.1) относится в данной области D к одному из следующих типов:
К эллиптическому - <0;
К параболическому - =0;
К гиперболическому - >0.
В дальнейших параграфах будут рассмотрены примеры каждого из трех типов уравнений и методы их приближенного решения.