4.2.2. Конечно-разностная аппроксимация уравнений Лапласа и Пуассона
Проиллюстрируем процедуру конечно-разностной аппроксимации на примере решения основных двумерных уравнений типа (4.3), (4.5) конечно-разностным методом
. (4.19)
Запишем конечно-разностную аппроксимацию (4.19), используя формулы (4.18) и ориентируясь на рис. 4.3.
(4.20)
Рис. 4.3. Геометрия узлов в квадратной сетке.
Складывая выражения в (4.20), используя уравнение (4.19) и умножая на , получим
(4.21)
Таким образом, мы аппроксимировали значение функции в центральном узле первым слагаемым в (4.21); второе слагаемое – погрешность аппроксимации
. (4.22)
Уравнение (3.10) может быть записано для любого узла анализируемой области как
. (4.23)
Уравнение (4.33) – основное конечно разностное уравнение для решения уравнения Пуассона (Лапласа); его также называют пятиточечной конечно разностной аппроксимацией обсуждаемых уравнений. Записывая (4.33) для каждого узла анализируемой области и преобразуя эти уравнения можно получить совместную систему уравнений относительно значений функции в узлах сетки, решая которую далее найти эти значения.
- Оглавление
- Предисловие
- Раздел 1. Основные сведения о среде программирования «matlab»
- Глава 1. Вычисления в командном режиме
- 1.1. Простейшие математические операции в matlab
- 1.2. Переменные
- 1.3. Создание матриц
- 1.4. Доступ к элементам матриц
- 1.5. Операции с матрицами
- 1.6. Ввод, вывод и работа со строками
- Глава 2. Построение графиков в matlab
- 2.1. Построение графика в виде двумерной линии
- 2.2. Оформление графиков
- 2.3. Построение трехмерных графиков
- 2.4. Построение линий уровня
- 2.5. Построение векторного поля
- 2.6. Отображение нескольких графиков в одном окне
- Глава 3. Скрипты в matlab и управляющие конструкции
- 3.1. Создание и выполнение скриптов в matlab
- 3.2. Оператор for
- 3.3. Логические операции
- 3.4. Оператор if / elseif / else
- 3.5. Оператор while
- 3.6. Операторы break / continue
- 3.7. Оператор switch
- 3.8. Создание функций
- Раздел 2. Краткие теоретические сведения и задания Тема 1. Векторный анализ
- 1.1. Элементы векторного анализа
- Задания
- Тема 2. Уравнения Максвелла (произвольная и гармоническая временная зависимость, статические, стационарные и квазистационарные поля)
- 2.1. Система уравнений электродинамики – уравнения Максвелла
- 2.2. Граничные условия. Принцип эквивалентности
- Задания
- Тема 3. Плоские волны
- 3.1. Явление дисперсии и групповая скорость
- Задания
- Тема 4. Граничные задачи, уравнения и методы
- 4.1. К классификации электромагнитных явлений
- Задания
- 4.2. Метод конечных разностей
- 4.2.1. Конечно-разностная аппроксимация
- 4.2.2. Конечно-разностная аппроксимация уравнений Лапласа и Пуассона
- 4.2.3. Конечно-разностная аппроксимация для граничных узлов
- Задания
- Литература