WWW.REFERATCENTRAL.ORG.UA - Я ТУТ НАВЧАЮСЬ

... відкритий, безкоштовний архів рефератів, курсових, дипломних робіт

ГоловнаМатематика, Геометрія, Статистика → Наближення сплайнами третього степеня - Реферат

Наближення сплайнами третього степеня - Реферат


Реферат на тему:
Наближення сплайнами третього степеня
Відомо, що інтерполяція по вузлах, що співпадають з нулями многочлена Чебишева практично не відрізняється від найкращого рівномірного наближення многочленами такого ж порядку. Але цей факт дозволяє проілюструвати суттєве обмеження апроксимації многочленами: якщо функція, яку ми наближуємо, має особливості в деяких точках на інтервалі інтерполяції, то вона погано наближатиметься на всьому інтервалі. Крім того при роботі з многочленами в процесі обчислень швидко накопичуються похибки заокруглень. Вони вже великі для . При використанні многочленів високих степенів їх графіки, як правило, мають осциляції. Цієї загальної залежності можна запобігти, якщо використовувати кусково апроксимуючі функції. При цьому необхідно ставити умови достатньої гладкості спряження графіків многочленів. Під цим розуміють вимогу, щоб в точці з'єднання сусідніх ділянок многочлени, які належать лівій та правій ділянкам і похідні від них до певного порядку співпадали.
Цю задачу добре розв'язують так звані поліноміальні сплайни. Нагадаємо визначення сплайна. На відрізку введемо сітку
. (1)
Функцію називають поліноміальним сплайном степеня m дефекту гладкості k на [a,b], якщо виконуються умови:
1. має на [a,b] неперервні похідні до порядку m-k включно;
2. На кожному відрізку многочлен степеня не вище m. Точки називаються вузлами сплайна .
Простим прикладом сплайна є залишковий член інтерполяції.
Нехай функція u(x) визначена на відрізку і нехай на ньому взято n+1 вузол . Якщо в точках виконуються умови , то називають інтерполяційним сплайном для функції u(x), а вузли - вузлами інтерполяції.
Лінійний інтерполяційний сплайн записують у вигляді
, (2)
,
де
,
а кубічний (дефекту 1) у вигляді
(3)
.
В (2) і (3) вузли сплайна і вузли інтерполяції співпадають.
Перші два доданки кубічного сплайна відповідають лінійному сплайну, а кубічна поправка забезпечує додаткову гладкість. Із зображення (3) випливає властивість інтерполяції незалежно від вибору . Через те, що - лінійна функція, то з (3) одержуємо
, . (4)
Для визначення через значення f(x) у вузлах скористаємося умовою неперервності перших похідних сплайна в точках розбиття.
Обчислення в граничній точці відрізків дає
, (4`)
. (4``)
Умова гладкості приводить до співвідношень
(5)
Це система n-1 лінійних алгебраїчних рівнянь відносно невідомих . Тому задаються ще дві умови, а потім розв'язують методом прогонки відповідну систему лінійних алгебраїчних рівнянь.
Якщо в процесі роботи зі сплайнами значення змінюються, то систему рівнянь ми змушені розв'язувати заново. Тому зручнішими для цього виявляються інтерполянти локального характеру. Перепишемо систему для визначення у вигляді
(6)
де введені такі позначення для різницевих похідних на нерівномірній сітці
Із (6) випливають наближені формули
, (7)
. (8)
Із (7) і (8) випливає, що без розв'язування системи лінійних алгебраїчних рівнянь можна наближено обчислити параметри сплайна. Якщо використовується (7), то одержуємо так званий локальний дискретний кубічний сплайн, а якщо (8) - більш точний.
Слід відзначити, що при цьому для обчислення , , , треба знати зовні відрізка [a,b] значення f-2, f-1, fn+1, fn+2, які визначаються кубічною інтерполяцією з використанням умов
(9)
На користь локальної інтерполяції вказує і та обставина, що значення відомі, як правило, з деякою похибкою.
Іноді зручно використовувати таку форму запису кубічного сплайна
(10)
Тут
При побудові кубічного інтерполяційного сплайна дефекту 1 попутно знаходяться значення похідних першого та другого порядків
де одержується при розв'язуванні системи (6), а визначаються формулами (4), (4`).
Кубічні сплайни мають дуже важливу властивість, яка обумовлює високу ефективність сплайн-інтерполяції. Виявляється, що серед всіх функцій , які інтерполюють функцію, кубічний сплайн з крайовими умовами
(11)
мінімізує функціонал
. (12)
Цей факт дає змогу по іншому визначити кубічний інтерполяційний сплайн: це така функція із класу , яка у вузлах сітки приймає значення і мінімізує функціонал (12).
Якщо то при використанні умов (12) при побудові сплайна погіршується точність наближених формул поблизу границі. Коли та відомі, то слід покласти , . Якщо в кінцевих точках відрізка відома перша похідна, то додаткові співвідношення можна одержати, використовуючи (4`), (4``).
Так поклавши в (4``) i=0 та прирівнявши одержаний вираз до , маємо
(13)
Якщо ж невідоме, то можна діяти таким чином. Будується інтерполяційний поліном 3-го степеня для за точками x0, x1, x2, x3. Значення апроксимується виразом і покладається
(14)
Для знаходження параметрів одержуємо систему лінійних алгебраїчних рівнянь з трьохдіагональною матрицею. Матриця системи симетрична з діагональним переважанням і, отже, з визначником відмінним від нуля.
Розглянемо ще один підхід до побудови інтерполяційних сплайнів. Сплайн для візьмемо у вигляді
(15)
,
де . Перепишемо (15) у вигляді
(16)
,
де , = , + +1,
+ - сплайн 1-го
Loading...

 
 

Цікаве