Метод конечных разностей

В соответствии с идеей метода конечных разностей строится локальная аппроксимация решения (для трех соседних узлов разностной сетки с номерами i, j и k, соответственно) в виде:

,

где – квадратичные функции, принимающие значения 1 в своем (одноименном) узле и 0 в соседних узлах. Невязка дифференциального уравнения ( 1 .11) на таком приближении решения для всего отрезка [x_i, x_j] длиной h имеет вид:

.

В качестве взвешивающей возьмем –функцию Дирака, . Тогда выражение ( 1 .16) метода взвешенных невязок приводится к виду

,

,

где – расстояние между двумя соседними узлами, то есть шаг сетки. Последнее выражение является конечно-разностной аппроксимацией уравнения ( 1 .11), применяемой в сеточных методах.

-функция Дирака

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

Например, распределение масс вдоль прямой удобно задавать плотностью этого распределения. Однако, если на прямой имеются отдельные точки, несущие сосредоточенную положительную массу, то плотность такого распределения не может быть описана никакой обычной функцией. Применение классического аппарата математического анализа для решения целого ряда задач приводит к невозможности выполнения некоторых операций, когда функцию, не имеющую производную, невозможно продифференцировать, если эту производную понимать в обычном смысле. Оказывается, что подобные затруднения можно преодолеть введением понятия обобщенной функции. В физике интуитивное понятие обобщенной -функции введено и используется достаточно давно, значительно раньше, чем была построена строгая математическая теория обобщенных функций.