юд. вар 1 лаб работа 3. Аппроксимация первой и второй производных через конечные разности
Скачать 52.07 Kb.
|
Аппроксимация первой и второй производныхчерез конечные разностиВспомним определение первой производной. Если f (x) – функция одной переменной и x0 ∈ [a, b], то функцию f ′(x) можно записать , (1) где x0 – фиксированная точка. Геометрическая интерпретация показана на рис. 1. Пусть на графике функции f заданы фиксированная точка P0 [x0, f (x0)] и подвижная точка P [x, f(x)], и секущая, проведенная через эти точки, образует угол β с положительным направлением оси x. , (2) Рис. 1. Геометрическая интерпретация первой производной Разностное отношение функции f в точке x0 равно угловому коэффициенту секущей, проведенной через точки P и P0. Функция f называется дифференцируемой в точке x0∈[a,b], если существует предел разностного отношения функции в точке x0: , (3) Предел (3) называется производной функции f в точке x0 и обозначается (4) Производная функции f в точке x0 – это тангенс угла наклона касательной к графику функции f в точке P0[x0, f(x0)](рис. 1) tgα= f ′(x0 )`. (5) Простейшая формула численного (приближенного) дифференцирования для непрерывной функции в точке x0 через конечные разности имеет вид f ′(x0) = f x( 0 + ∆x) − f x( 0)(6) ∆x или f ′(x0) = f x( 0 + ∆x) − f x( 0 − ∆x) , 2⋅∆x где ∆x=x1– x0или в общем виде ∆x=xi – xi-1– шаг дифференцирования, величина которого должна быть достаточно малой. Если производная функции f / (4) дифференцируема в точке то [f ′( )x ]′ называется второй производной функции f в x0, x x= 0 точке x0 и обозначается одним из приведенных способов ′′(x0),d2 f x(2 0) ,d22f . (7) f dx dx x x= 0 Формула численного нахождения второй производной f ′′(x0) = f ′(x0 + ∆x) − f ′(x0)(8) ∆x При подстановке в (8) выражения для нахождения первой производной получим f ′′(x0) = f x( 0 + 2⋅∆x) − 2⋅ f x2( 0 + ∆x) + f x( 0)(9)∆x 45 или f ′′(x0) = f x( 0 + ∆x) − 2⋅ f x(2 0) + f x( 0 − ∆x) . ∆x При численном дифференцировании исходят из того, что функция f (x) задана конечной последовательностью пар значений (xi , fi) без помехи, и приближенные значения величин f / (xi ) и f // (xi ) находят по формулам (6) и (9). |