1.3.5 Сплайн-інтерполяція

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

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

 i=1,2, ... ,m (1.25)

де  – постійні, які визначені вказаними умовами (j= 1,2,3,4).

Перші (2m) умов потребують, щоб сплайни стикалися в заданих точках:

,i=1, 2, ... , m,

,  i=0, 1, ... , m-1. (1.26)

Наступні (2m-2) умов потребують, щоб в місцях дотику сплайнів були рівні перші та другі похідні

 i=1, ... , m-1, (1.27)

 i=1, ... , m-1.


Система алгебраїчних рівнянь має розв’язок, якщо кількість рівнянь дорівнює кількості невідомих. Для цього необхідні ще два рівняння. Як правило, використовують такі додаткові умови:

  (1.28)

Отриманий таким чином сплайн зветься “природним кубічним сплайном”.

В багатьох випадках метод сплайнів є найбільш зручним, тому що це дозволяє отримати аналітичну кусково-поліноміальну функцію. Існують сплайни більш вищих порядків. Вживання цього методу можливо і в інших галузях обчислювальної математики, наприклад, в чисельному інтегруванні і розв’язанні диференціальних рівнянь.

1.4 Уточнена постановка задачі

Нехай на відрізку [а; b] визначено певний клас функцій {Р(х)}, наприклад клас алгебраїчних многочленів, а в точках х0, х1,..., хn цього проміжку задано значення деякої функції y=f(x): y0=f(x0), y1=f(x1),….yn=f(xn). Наближену заміну функції f на відрізку [а; b] однією з функцій Р(х) цього класу так, щоб функція P(х) в точках x0,x1, ..., xn набувала тих самих значень, що й функція f, тобто щоб Р(xi)= уi (і = 0, 1, ..., n), називають інтерполюванням, або інтерполяцією. Точки х0, хi, ..., хп називають вузлами інтерполювання, функцію Р(х) — інтерполюючою функцією, а формулy у=P(х), за допомогою якої обчислюють значення функції f у проміжку [а;b], — інтерполяційною формулою.

З геометричного погляду задача інтерполювання полягає в знаходженні кривої у= Р(х) певного класу, яка проходить через точки площини з координатами (хi, уi)

(i = 0, 1, ....,n) (рис.1.1.1).

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

У деяких випадках доцільніше використовувати інші класи інтерполюючих функцій. Якщо, наприклад, функція f періодична, то функцію Р(х) природно вибирати з класу тригонометричних многочленів, а якщо функція f перетворюється в нескінченність у заданих точках або поблизу них, то функцію Р(х) доцільно вибирати з класу раціональних функцій.

Розглядатимемо лише задачу параболічного інтерполювання, яку сформулюємо так: в n+1 різних точках х0, x1,..., хn задано значення

функції f: y0=f(x0), y1=f(x1),…, yn=f(xn) і треба побудувати многочлен

 (1.29)

 

степеня n, який задовольняв би умови

  (1.30)

Для визначення n+1 коефіцієнтів многочлена (1.29), який задовольняє умови (1.30), запишемо систему (n+1)-го лінійних рівнянь виду:

(1.31)

Ця система має єдиний розв'язок, бо її визначник є визначником Вандермонда, який не дорівнює нулю, бо вузли xi=(i=0,1,…,n) різні. А тому й задача параболічного інтерполювання має єдиний розв'язок, тобто існує єдиний алгебраїчний многочлен виду (1.29), що задовольняє умови (1.30). Многочлен Рn(х), який задовольняє умови (1.29), називають інтерполяційним многочленом, наближену рівність f(x)=Pn(x) –

інтерполяційною формулою, а різницю Rn(f,x)=f(x) – Pn(x) — залишковим членом інтерполяційної формули. Хоч інтерполяційний многочлен, що задовольняє умови (1.29), і єдиний, проте можливі різні форми його запису.

Інтерполяційний многочлен будують тоді, коли:

1) функцію задано таблично для деяких значень аргументу, а треба знайти її значення для значень аргументу, яких у таблиці немає;

2) функцію задано графічно, наприклад за допомогою самописного приладу, а треба знайти її наближений аналітичний вираз;

3) функцію задано аналітичнo, але її вираз досить складний і незручний для виконання різних математичних операцій (диференціювання, інтегрування тощо).



Информация о работе «Дослідження зміни температури термопари за допомогою чисельних методів на ЕОМ»
Раздел: Информатика, программирование
Количество знаков с пробелами: 57641
Количество таблиц: 6
Количество изображений: 20

Похожие работы

Скачать
19877
2
11

... Interpol. Для запуску програми необхідно запустити файл Interpol. exe. Після запуску програми з’являється меню, реалізоване в текстовому режимі, яке дозволяє вибрати необхідний метод інтерполяції. Для вибору методу розподілених різниць необхідно натиснути клавішу "1", для вибору методу Лагранжа - клавішу "2". Після вибору методу для графічного представлення результатів інтерполяції функції в ...

Скачать
266076
11
92

... Методичні вказівки до лабораторної роботи № З «Тепловіддача горизонтальної труби при вільному русі повітря». Тернопіль 2003 У даних методичних вказівках подані теоретичні основи, опис експернментальної установки і практичні рекомендації лля проведення лабораторної роботи і обробки дослідних даних Мета роботи - засвоїти знання з теорії" конвсктивнот теплообміну при ...

Скачать
468112
34
0

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

Скачать
50784
6
1

... підвищенням гнучкості технологічних систем, застосуванням машин з числовим програмним керуванням, впровадженням автоматичних багатофункціональних маніпуляторів-роботів. РОЗДІЛ 3. ЕКОНОМІЧНЕ ОБҐРУНТУВАННЯ ОСНОВНИХ НАПРЯМКІВ РАЦІОНАЛІЗАЦІЇ ТРУДОВИХ ПРОЦЕСІВ НА ВАТ "ЕЛЕКТРОТЕРМОМЕТРIЯ" Провівши аналіз організації трудових процесів на ВАТ "Електротермометрія", слід було б відмітити, що в 2009 ...

0 комментариев


Наверх