Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Численные методы. Учебно-методический комплекс.pdf
Скачиваний:
0
Добавлен:
15.08.2026
Размер:
1 Мб
Скачать
f ( x )
F( x )

приближенного числа. Если положительное приближенное число а имеет n верных десятичных знаков, то относительная погрешность d не превосходит (1/10)n-1 , деленную на первую значащую цифру данного числа.

Вопросы

1.Источникипогрешностей результат численного решениязадач.

2.Классификация погрешностей результат численного решения

задач.

3.Абсолютная и относительная погрешности

Тема. Интерполирование функций.

1. Постановка задачи. Простейшая задача интерполирования заключается в следующем. На отрезке [ a,b ] заданы n +1 точки x0 , x1 ,K, xn , которые называются узлами интерполяции, и значения некоторой функции f ( x ) в этих точках

f ( x0 ) = y0 , f ( x1 ) = y1 ,…, f ( xn ) = yn .

(1)

Требуется построить функцию F( x ) (интерполирующая функция), принадлежащую известному классу и принимающую в узлах интерполяции те же значения, что и f ( x ) , т.е. такую, что

F( x0 ) = y0 , F( x1 ) = y1 , …, F( xn ) = yn .

(2)

В такой общей постановке задача может иметь бесчисленное множество решений или совсем не иметь решений. Однако эта задача становится однозначной, если вместо произвольной функции искать полином Pn ( x ) степени не выше n , удовлетворяющий условиям (2), т.е. такой, что

Pn ( x0 ) = y0 , Pn ( x1 ) = y1 ,…, Pn ( xn ) = yn .

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

когда x [ x0 , xn ] , т.е. значение x является промежуточным между x0 и xn , и экстраполирование, когда x [ x0 , xn ] . В дальнейшем под термином интерполирование будем понимать как первую, так и вторую операцию.

КазНПУ 703-10-11 Учебно-методический комплекс дисциплины. Издание третье.

11

2. Интерполяционная формула Лагранжа. Для произвольно заданных узлов интерполирования пользуются так называемой

интерполяционной формулой Лагранжа.

Пусть на отрезке [ a,b ] заданы n +1 различных значений x0 , x1 ,K, xn и известны для функции y = f ( x ) соответствующие значения:

f ( x0 ) = y0 , f ( x1 ) = y1 ,…, f ( xn ) = yn .

Требуется построить полином Ln ( x ) степени не выше n ,

имеющий в заданных узлах x0 , x1 ,K, xn

те же значения,

что и функция

f ( x ) , т.е. такой, что Ln ( xi ) = yi

( i = 0,1,K,n ) .

 

 

 

 

 

 

Полином, удовлетворяющий перечисленным условиям, был

построен. Он имеет следующий вид:

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

( x x0 )( x x1 )K( x xi1

)( x xi+1 )K( x xn )

.

 

 

 

 

 

 

Ln ( x ) = ∑ yi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( xi x0 )( xi x1 )K( xi xi1

)( xi xi+1 )K( xi xn )

 

 

 

 

 

 

 

i=0

 

 

 

Это и есть интерполяционная формула Лагранжа.

 

 

3. Первая и вторая интерполяционные формулы Ньютона.

Пусть

для

равноотстоящих

значений

независимой

переменной

xi = x0 +ih

( i = 0,1,K,n ) , где h - шаг интерполяции, заданы значения f ( xi ) = yi .

Первая

интерполяционная формула

Ньютона

для

функции

y = f ( x ) имеет вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Pn ( x ) = y0 + q y0 +

q( q 1)

2 y0 +K+

q( q 1)K( q n +1)

n y0 ,

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

x x1

 

 

 

2!

 

 

 

n!

 

 

 

где q =

 

x x0

,

q 1 =

,

и т.д.,

h = xi+1 xi

( i = 0,1,2,K,n ) ,

x0

-

начальное

 

 

 

 

h

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

значение аргумента, y0

= f ( x0 ) -

значение функции, соответствующее

начальному

значению

аргумента,

i y0 -

конечная разность i -го

порядка, причем

i y0 = i1 y1

i1 y0 .

Величина q

представляет собой число

шагов, необходимых для достижения точки x , исходя из точки x0 . Формулу (1) выгодно использовать для интерполирования

функции y = f ( x ) в окрестности начального значения x0 , где q мало по абсолютной величине.

Далее, первая интерполяционная формула Ньютона неудобна для интерполирования вблизи конца таблицы. В этом случае применяется вторая интерполяционная формула Ньютона:

 

 

P ( x ) = y

+ q y

 

+

q( q +1)

2 y

 

+K+

q( q +1)K( q + n 1)

n y

,

n1

 

n2

 

 

 

n

n

 

 

2!

 

 

 

n!

0

 

 

x xn

 

x xn1

 

 

 

 

 

 

 

 

где q =

, q +1 =

,

и т.д.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

h

 

 

 

 

 

 

 

 

 

 

4. Интерполяционная формула Ньютона в случае неравно-

отстоящих узлов.

Пусть

на отрезке [ a,b ]

заданы n +1 различных

КазНПУ 703-10-11 Учебно-методический комплекс дисциплины. Издание третье.

12

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]