Курс лекций по основам алгоритмизации и программирования задач машиностроения. Кравченко Д.В. - 93 стр.

UptoLike

Составители: 

Рубрика: 

91
*
2
1i1ii
y
6
h
)yy(
h
2
1
y
=
+
;
n
= 2, (42)
V
*
4
2i1i1i2ii
y
30
h
)yy8y8y(
h
12
1
y ++=
++
;
n
= 4
Выражения (42) называются аппроксимациями производных с помо-
щью центральных разностей и широко применяются на практике.
Остаточный член многочлена Лагранжа:
)x(f
)!1n(
)xx)...(xx)(xx(
)x(R
*
)1n(
n10
L
+
+
=
,
где )x(f
*
)1n( +
производная
)
1
n
(
+
порядка функции
)
x
(
f
в некоторой точке
*
xx
,
]x,x[x
n0*
.
Если использовать многочлен Лагранжа для случая произвольно располо-
женных узлов, то придется вычислять громоздкие выражения. Поэтому здесь
удобнее применять метод неопределенных коэффициентов.
Сущность его заключается в следующем.
Искомое выражение для производной
k
-го порядка в некоторой точке
i
xx
=
представляют в виде линейной комбинации заданных значений функции
в узлах
0
x ,
1
x , … ,
n
x :
nn1100
)k(
i
yc...ycycy +++=
. (43)
При этом предполагают, что формула (43) имеет место для многочленов
1
y
=
,
i
xxy
=
,…,
n
i
)xx(y =
.
Подставляя последовательно эти выражения в (40), получают систему
(
1
n
+
) линейных алгебраических уравнений для определения неизвестных ко-
эффициентов
0
c ,
1
c , …,
n
c .
Пример 16. Найти выражение для производной
1
y
в случае четырех рав-
ноотстоящих узлов (
3
n
=
).
Равенство (43) запишем в виде:
332211001
ycycycycy
+
+
+
=
. (44)
Используем следующие многочлены:
1
y
=
;
0
xxy
=
;
2
0
)xx(y = ;
3
0
)xx(y = . (45)
Вычислим их производные:
0y
=
; 1y
=
; )xx(2y
0
=
;
2
0
)xx(3y =
. (46)
Подставим последовательно соотношения (45) и (46) соответственно в пра-
вую и левую части равенства (44) при
1
xx
=
:
1c1c1c1c0
3210
+
+
+
=
;