By | 25/10/2013

Нехай для функції Перша інтерполяційна формула Нютона задані значення i_interpolacijna_formyla_nytona2 для рівновіддалених вузлів, тобто i_interpolacijna_formyla_nytona31, де h крок інтерполяції. Потрібно знайти поліном i_interpolacijna_formyla_nytona4, степінь якого не перевищує n, і який в точках i_interpolacijna_formyla_nytona5 набуває значень i_interpolacijna_formyla_nytona61.

Даний поліном будемо шукати у наступному вигляді:

i_interpolacijna_formyla_nytona71

Використовуючи узагальнену степінь числа, вираз (2) запишемо у наступному вигляді:

i_interpolacijna_formyla_nytona8

Задача полягає у знаходженні коефіцієнтів i_interpolacijna_formyla_nytona91. У виразі (2) покладемо i_interpolacijna_formyla_nytona10. В результаті отримаємо i_interpolacijna_formyla_nytona11.

Для того, щоб знайти коефіцієнт i_interpolacijna_formyla_nytona12 запишемо скінченну різницю першого порядку (скінченною різницею першого порядку називають різницю між значеннями функції у сусідніх вузлах інтерполяції, тобто, i_interpolacijna_formyla_nytona26):

i_interpolacijna_formyla_nytona13

Покладаючи в останньому виразі i_interpolacijna_formyla_nytona10, отримаємо i_interpolacijna_formyla_nytona14. Звідси i_interpolacijna_formyla_nytona15.

Для знаходження коефіцієнта i_interpolacijna_formyla_nytona16 запишемо скінченну різницю другого порядку. Для знаходження яких використовуються скінченні різниці першого порядку (i_interpolacijna_formyla_nytona27).

i_interpolacijna_formyla_nytona17

Знову, покладаючи в останній вираз i_interpolacijna_formyla_nytona10, отримаємо i_interpolacijna_formyla_nytona18. Звідки отримуємо i_interpolacijna_formyla_nytona19.

Продовжуючи даний процес, отримаємо загальну формулу для обчислення коефіцієнтів: i_interpolacijna_formyla_nytona20, де скінченні різниці i-го порядку визначаються наступним чином: i_interpolacijna_formyla_nytona28

Підставляючи знайдені коефіцієнти в формулу (2), отримаємо інтерполяційний поліном Ньютона:

i_interpolacijna_formyla_nytona21

Для практичного використання інтерполяційну формулу Ньютона (3), зазвичай записують у дещо іншому вигляді. Для того вводять нову зміннуi_interpolacijna_formyla_nytona22. Підставляючи дану змінну в (3), перша інтерполяційна формула Ньютона набуде наступного вигляду:

i_interpolacijna_formyla_nytona23

де q представляє собою число кроків необхідних для досягнення точки x, виходячи з точки i_interpolacijna_formyla_nytona24. Таким чином, ми отримали кінцевий варіант першої інтерполяційної формули Ньютона.

Перша інтерполяційна формула Ньютона – приклад:

Нехай функція Перша інтерполяційна формула Нютона задана таблично:

i_interpolacijna_formyla_nytona30

Таблиця фіксованих значень функції

Необхідно, скориставшись першою інтерполяційною формулою Ньютона, обчислити значення функції в точці i_interpolacijna_formyla_nytona311, яка являється відмінною від заданих. Для цього, на першому кроці обчислимо скінченні різниці до шостого порядку включно:

    1. Скінченні різниці першого порядку:

i_interpolacijna_formyla_nytona38

    1. Скінченні різниці другого порядку:

i_interpolacijna_formyla_nytona39

    1. Скінченні різниці третього порядку:

i_interpolacijna_formyla_nytona34

    1. Скінченні різниці четвертого порядку:

i_interpolacijna_formyla_nytona351

    1. Скінченні різниці п’ятого порядку:

i_interpolacijna_formyla_nytona42

  1. Скінченні різниці шостого порядку: i_interpolacijna_formyla_nytona43.

Підставляючи отримані значення, значення з таблиці і точку i_interpolacijna_formyla_nytona40, в формулу (4) отримуємо наближене значення функції в заданій точці:

i_interpolacijna_formyla_nytona45

Блок-схема програмної реалізації першої інтерполяційної формули Ньютона:

i_interpolacijna_formyla_nytona25

4 Replies to “Перша інтерполяційна формула Ньютона для рівновіддалених вузлів інтерполяції”

  1. Динчук Василь

    Це що за пастка в циклах? В перших двох умовах “<n-1” а потім бац і “<n” що призводить до помилки виходу за межі масиву наскільки я розумію.

  2. admin Автор

    Доброго вечора Василь. Не зовсім зрозумілим являється Ваш коментар. Якщо мова йде про перших два цикла в блок-схемі, то хочу Вас запевнити, що ніяких «пасток» там не має. По даній блок-схемі було реалізовано delphi-проект Перша інтерполяційна формула Ньютона, який працює безпомилково.

  3. Динчук Василь

    Так, все вірно. Я писав код під C# і трішки заплутався. Дякую за блок-схему і приклад.

  4. admin Автор

    Зрозуміло Василь. Раді, що матеріал був корисним.

Залишити відповідь

Ваша e-mail адреса не оприлюднюватиметься. Обов’язкові поля позначені *

*