1.8.5. Частные производные
Рассмотрим функцию двух переменных U=f(x, y), заданную в виде таблицы Uij=f(x, y), где xi=x0+ihx (i=0,1,...,i), yi=y0+ihy (j=0, 1,...,j). Используя понятие частной производной, можно приближенно записать для малых значений шагов hy, hy
Введем обозначения по правилам понятным из таблицы
y x |
xi-2 |
xi-1 |
xi |
xi+1 |
xi+2 |
yj-2 |
Ui-2,j-2 |
Ui-1,j-2 |
Ui,j-2 |
Ui+1,j-2 |
Ui+2,j-2 |
yj-1 |
|
Ui-1,j-1 |
Ui,j-1 |
Ui+1,j-1 |
Ui+2,j-1 |
yj |
|
Ui-1,j |
Ui,j |
Ui+1,j |
Ui+2,j |
yj+1 |
|
Ui-1,j+1 |
Ui,j+1 |
Ui+1,j+1 |
Ui+2,j+1 |
yj+2 |
|
|
|
Ui+1,j+2 |
Ui+2,j+2 |
Тогда получим следующие приближенные выражения (аппроксимации) для частных производных в узлах хi, yj с помощью отношений конечных разностей:
.
Для численного дифференцирования функций многих переменных воспользуемся разложением в ряд Тейлора функции многих переменных
используем эту формулу дважды:
найдем Ui+1,j=f(xi+hx, yj) при x=hx, y=0;
найдем Ui-1,j=f(xi-hx, yj) при x=-hx, y=0;
Вычтем из первого равенства второе:
Отсюда найдем аппроксимацию производной с помощью центральных разностей
Она имеет второй порядок.
Аналогично могут быть получены аппроксимации производной , а также старших производных. В частности для второй производной можно получить
Записывая разложения в ряд при разных значениях x и y, можно вывести различные формулы численного дифференцирования с необходимым порядком аппроксимации.
Данные аппроксимации производных могут быть использованы при построении разностных схем для решения уравнений с частными производными.