ВУЗ:
Составители:
50
Теорема доказана.
Перейдем к изучению остатка квадратуры.
Теорема 4. Если f имеет производную порядка n2 на ],[ ba , то
существует такая точка ],[ ba
∈
η
,что для остатка квадратуры наи-
высшей степени точности верно равенство
∫
≤η≤ω
η
=
b
a
n
badxxxp
n
f
fR .,)()(
)!2(
)(
)(
2
)2(
(6)
Доказательство
. Построим интерполяционный многочлен )(xH
степени
12 −≤ n , удовлетворяющий условиям
).()(),()(
kkkk
xfxHxfxH
′
=
′
=
Остаток интерполирования )()()( xHxfxr
−
=
имеет следующее
представление:
()
(
)
(
)
()
()
x
n
f
xr
n
2
2
!2
ω
ξ
= ,
где ξ – некоторая точка, лежащая между
x
и узлами
k
x .
Таким образом,
∫∫ ∫
ωξ+=
b
a
b
a
b
a
n
dxxfxp
n
dxxHxpdxxfxp .)()()(
)!2(
1
)()()()(
2)2(
Поскольку квадратурная формула верна для всех многочленов
степени
12 −≤ n
, а степень многочлена )(xH также
12
−
≤
n
, то
∫
∑∑
==
==
b
a
n
k
n
k
kkkk
xfAxHAdxxHxp
11
)()()()(
и для остатка квадратуры имеет место равенство
∫
ωξ=
b
a
n
dxxxpf
n
fR .)()()(
)!2(
1
)(
2)2(
Теорема доказана.
Теорема доказана.
Перейдем к изучению остатка квадратуры.
Теорема 4. Если f имеет производную порядка 2n на [a, b] , то
существует такая точка η ∈ [ a, b] ,что для остатка квадратуры наи-
высшей степени точности верно равенство
f ( 2n) (η) b 2
R( f ) = ∫ p( x)ω ( x)dx, a ≤ η ≤ b. (6)
( 2n)! a
Доказательство. Построим интерполяционный многочлен H (x)
степени ≤ 2n − 1 , удовлетворяющий условиям
H ( xk ) = f ( xk ), H ′( xk ) = f ′( xk ).
Остаток интерполирования r ( x) = f ( x) − H ( x) имеет следующее
представление:
f (2n ) (ξ ) 2
r (x ) = ω (x ) ,
(2n ) !
где ξ – некоторая точка, лежащая между x и узлами xk .
Таким образом,
b b
1 b ( 2n)
∫ p( x) f ( x)dx = ∫ p( x) H ( x)dx + ∫ p( x) f (ξ)ω2 ( x)dx.
a a (2n)! a
Поскольку квадратурная формула верна для всех многочленов
степени ≤ 2n − 1 , а степень многочлена H (x) также ≤ 2n − 1 , то
b n n
∫ p( x) H ( x)dx = ∑ Ak H ( xk ) = ∑ Ak f ( xk )
a k =1 k =1
и для остатка квадратуры имеет место равенство
1 b ( 2n)
R( f ) = ∫ f (ξ) p ( x)ω2 ( x)dx.
(2n)! a
Теорема доказана.
50
Страницы
- « первая
- ‹ предыдущая
- …
- 48
- 49
- 50
- 51
- 52
- …
- следующая ›
- последняя »
