4.2.1. Конечно-разностная аппроксимация

МКР включает следующие основные шаги:

1. Для получения основного уравнения производные в решаемом дифференциальном уравнении выражаются в разностной форме.

2. Область дискретизируется в форе ячеек, содержащих конечное число узлов.

3. Конечно-разностная форма дифференциальных уравнений применяется для каждого узла, тем самым, формируя систему линейных уравнений.

4. Решается система линейных уравнений, и находятся значения неизвестной функции в узлах ячеек.

5. Значения неизвестной функции в узлах ячеек используются для определения характерных параметров устройства (характеристического импеданса, критической длины волны, резонансной частоты и т.п.).

Пусть задана функция – рис. 4.2. Производная функции в точке определяется как

. (4.13)

Рис. 4.2. К разностным формам аппроксимации первой производной.

Выражение (4.13) является точным, однако неудобным для численной реализации, поскольку значения функция могут быть заданы только в дискретных точках. Поэтому аппроксимируем (4.13) в виде:

, (4.14 а)

где - шаг дискретизации, а представление (4.14а) носит название разностная аппроксимация вперёд. Возможно аналогичное представление, называемое разностной аппроксимацией назад

. (4.14 б)

Существует ещё одно возможное представление, называемое центрально разностной аппроксимацией

. (4.15)

Запишем разложение функции в ряд Тейлора в окрестности точки

(4.16)

Используя разностную аппроксимацию вперёд, получаем

(4.17а)

Второе, третье и последующие слагаемые в правой части (3.5) представляют погрешности аппроксимации для . Аналогично для аппроксимации назад

(4.17б)

Складывая (3.5а) с (3.5б) получаем выражение для центрально разностной аппроксимации

(4.17в)

Видно, что погрешность конечно-разностной аппроксимации (4.17а), (4.17б) порядка , а для центрально разностной аппроксимации, которая существенно меньше, поскольку . Эту погрешность называют погрешностью дискретизации.

При решении методом конечных разностей волновых уравнений, уравнений Лапласа, Пуассона необходима конечно-разностная аппроксимация вторых производных . Обратимся к рис. 4.2. и используем разложение в ряд Тейлора (4.16), располагая тремя точками и

,

(4.18а)

то есть погрешность аппроксимации порядка . Аналогично для функции двух переменных

. (4.18б)