mathserfer.narod.ru

Теория по высшей математике

‹-- Назад

Приближённое вычисление производных

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

Для приближённого нахождения в заданной точке часто поступают следующим образом. Исходя из того, что при достаточно малых приращениях разностное отношение мало отличается от своего предельного значения, равного производной ,¹⁵ мы можем приближённо заменить этим разностным отношением с малым , полагая , например, равным или . Таким образом, получаем приближённую формулу

Правая часть этой формулы при называется разностной производной вправо (или вперёд) с шагом .

Если же взять отрицательное приращение , , то аналогично получаем, что

Правая часть этой формулы при называется разностной производной влево (или назад) с шагом .

Согласно геометрическому смыслу производной, при замене производной разностной производной вправо или разностной производной влево, мы заменяем угол наклона касательной к графику углом наклона секущей , равным , или углом наклона секущей , равным , соответственно (см. следующий чертёж).

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

Правая часть полученной формулы называется центральной разностной производной с шагом . Эта формула применяется чаще других для практического нахождения .

Имеются и ещё более точные формулы для нахождения первой производной; приведём, например, гораздо более точную, чем (4.18), формулу

которая, правда, требует для своего применения не двух, а четырёх вычислений значения функции . Однако выигрыш в точности с лихвой перекрывает увеличение количества вычислений. По поводу методов получения приближённых формул вычисления производной см. [Бахвалов Н. С., Жидков Н. П., Кобельков Г. М., Численные методы. -- М.: Наука, 1987. -- Гл. II] или [Амосов А. А., Дубинский Ю. А., Копченова Н. В., Вычислительные методы для инженеров. -- М.: Высш. шк., 1994. -- Гл. 12].

Для нахождения способа приближённого вычисления второй производной введём такие обозначения. Разностную производную вперёд с шагом в точке обозначим как

разностную производную назад -- как

а центральную разностную производную с шагом  -- как

Поскольку вторая производная  -- это производная от первой производной , то естественно для получения приближённой формулы для заменить первую производную на какое-нибудь её приближение, а затем применить тот же способ приближённого вычисления производной. Например, если применять оба раза разностную производную вперёд, получим:

Точно так же, применяя два раза разностную производную назад, получим формулу

а применяя два раза центральную разностную производную с шагом  -- формулу

Последняя из трёх полученных формул предпочтительнее, поскольку основывается на более точной из трёх приближённых формул для первой производной.

Применяя тот же приём ещё раз, мы можем получить приближённые формулы для третьей производной. Например, основываясь на формуле разностной производной вправо, получим

а основываясь на центральной разностной производной --

Имеются и формулы для старших производных, дающие большую точность, чем приведённые выше. Например, для второй производной приведём формулу

О методах получения таких формул можно прочитать в книгах [Бахвалов Н. С., Жидков Н. П., Кобельков Г. М., Численные методы. -- М.: Наука, 1987. -- Гл. II] и [Амосов А. А., Дубинский Ю. А., Копченова Н. В., Вычислительные методы для инженеров. -- М.: Высш. шк., 1994. -- Гл. 12].