25. Особенности задач численного дифференцирования. Формулы численного дифференцирования

Задача численного дифференцирования состоит в приближенном вычислении производных функции f(x) по заданным в конечном числе точек значениям этой функции. Один из универсальных способов построения формул численного дифференцирования состоит в том, что по значениям функции f(x) в некоторых узлах x0 , x1 , ... , xN строят интерполяционный полином PN(x) (обычно в форме Лагранжа) и приближенно полагают f(r)(x) P(r)N(x), 0 ≤ r ≤ N В ряде случаев наряду с приближенным равенством удается (например, используя формулу Тейлора) получить точное равенство, содержащее остаточный член R (погрешность численного дифференцирования): f(r)(x) = P(r)N(x) + R, 0 ≤ r ≤ N Такие формулы называются формулами численного дифференцирования с остаточными членами. Степень, с которой входит величина в остаточный член, называется порядком погрешности формулы численного дифференцирования. Формулы с отброшенными остаточными членами называются просто формулами численного дифференцирования. 26. Численные методы решения обыкновенных дифференциальных уравнений

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

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

Приближенные методы основаны на различных упрощениях самих уравнений путем обоснованного отбрасывания (пренебрежения) некоторых содержащихся в них членов. В некоторых случаях сначала находят точное решение упрощенной задачи, а затем приближенно вычисляют поправки, обусловленные малыми членами, отброшенными на первом этапе. На этом основаны методы теории возмущений. Другой подход связан с представлением решения в виде разложения по малому параметру, который содержится в задаче. К данной группе относятся асимптотические методы, с помощью которых получают решения, описывающие некоторую предельную картину рассматриваемого явления.

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

Наиболее распространенным и универсальным подходом к численному решению дифференциальных уравнений является метод конечных разностей. Суть этого подхода состоит в следующем. Область непрерывного изменения независимой переменной (аргумента) заменяется дискретным множеством точек, называемых узлами. Эти узлы образуют расчетную сетку. Искомая функция непрерывного аргумента приближенно заменяется функцией дискретного аргумента на заданной сетке. Эта функция называется сеточной. Исходное дифференциальное уравнение заменяется разностным соотношением относительно сеточной функции. При этом входящие в уравнения производные заменяются разностными отношениями (см. гл. 4, §1). Такая замена дифференциального уравнения разностным уравнением называется его аппроксимацией на сетке. Таким образом, решение дифференциального уравнения сводится к отысканию значений сеточной функции в узлах сетки.

 27. Численные методы решения обыкновенных дифференциальных уравнений. Одношаговые методы-метод Эйлера (явный).

Метод Эйлера — наиболее простой численный метод решения (систем) обыкновенных дифференциальных уравнений.

Исторически первым и наиболее простым способом численного решения задачи Коши для ОДУ первого порядка является метод Эйлера. В его основе лежит аппроксимация производной отношением конечных приращений зависимой (y) и независимой (x) переменных между узлами равномерной сетки:

где y_i+1  это искомое значение функции в точке x_i+1.

Если теперь преобразовать это уравнение, и учесть равномерность сетки интегрирования, то получится итерационная формула, по которой можно вычислить  y_i+1 , если известно y_i  в точке х_i:

                                                                                             (6.4)

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

Графическая интерпретация метода Эйлера также не представляет затруднений (см. рисунок ниже). Действительно, исходя из вида решаемого уравнения (6.2) следует, что значение  есть значение производной функции y(x) в точке x=x_i -    , и, таким образом, равно тангенсу угла наклона каcательной, проведенной к графику функции y(x) в точке x=x_i.

 Из прямоугольного треугольника на рисунке можно найти

  ,

откуда и получается формула Эйлера. Таким образом, суть метода Эйлера заключается в замене функции y(x) на отрезке интегрирования прямой линией, касательной к графику в точке x=x_i. Если искомая функция сильно отличается от линейной на отрезке интегрирования, то погрешность вычисления будет значительной. Ошибка метода Эйлера прямо пропорциональна шагу интегрирования:

Ошибка ~ h

Процесс вычислений строится следующим образом. При заданных начальных условиях x₀ и y₀ можно вычислить

Таким образом, строится таблица значений функции y(x) с определенным шагом (h) по x на отрезке [x₀, x_N]. Ошибка в определении значения y(x_i) при этом будет тем меньше, чем меньше выбрана длина шага h (что определяется точностью формулы интегрирования).

При больших h метод Эйлера весьма неточен. Он дает все более точное приближение при уменьшении шага интегрирования. Если отрезок [x_i, x_i+1] слишком велик, то каждый участок [x_i, x_i+1] разбивается на N отрезков интегрирования и к каждому их них применяется формула Эйлера с шагом , то есть шаг интегрирования h берется меньше шага сетки, на которой определяется решение.