![](/user_photo/2706_HbeT2.jpg)
- •Глава 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.
В дальнейших параграфах будут рассмотрены примеры каждого из трех типов уравнений и методы их приближенного решения.