Знаходження значення функції за допомогою інтерполяційної формули Бесселя

Міністерство освіти і науки України

Вінницький національний технічний університет

Інститут АЕКСУ

Факультет ФЕЛТ

Кафедра Електроніки

КУРСОВА РОБОТА

З дисципліни “Обчислювальна математика”

Алгоритм та програма знаходження значення функції за допомогою інтерполяційної формули Бесселя

2004

ЗАВДАННЯ

на курсову роботу з дисципліни

“Обчислювальна математика”

студентові гр. МП-02 ТвердохлібуА.М. пропонується розробити

алгоритм та програму мовою програмування Турбо Паскаль на знаходження значення функції  за допомогою інтерполяційної формули Бесселя

Основні вхідні дані:

1. Кількість вузлів таблично заданої функції – не більше 100.

2. Похибка обчислень – не більше 0,001.

Основні вихідні дані:

1. Пояснювальна записка до курсової роботи.

2. Виконуваний файл програми.

АНОТАЦІЯ

 

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

ВСТУП

 

Задача знаходження значення функції у міжвузловій точці за допомогою інтерполяційної формули Бесселя має важливе значення при вирішенні як наукових, так і практичних задач, оскільки дає можливість знаходження значення функції у будь-якій точці, в якій це потрібно. В багатьох випадках функція не має аналітичного вигляду, тобто він невідомий, а задана лише декількома точками та значеннями функції в цих точках. Тому для отримання значення функції в проміжних точках застосовуються інтерполяційна формули Гауса (1-а та 2-а), інтерполяційна формула Стірлінга та Бесселя. Останні дві формули є похідними від першої та другої інтерполяційних формул Гауса. Кожна з цих формул має свої переваги та недоліки, що полягають у кількості обчислювальних операцій та в похибці обчислень.

Серед сучасного програмного забезпечення є багато програм чисельного аналізу, до яких можна віднести всесвітньо відомі пакети програм MathCad та MatLab. Вони, як правило, мають зручний інтерфейс та є багатофункціональними. Але їх недоліком є те, що задачі чисельного аналізу певного класу (наприклад, знаходження першої похідної) вирішуються за допомогою лише деякого одного методу. Крім того, вони займають багато дискової пам’яті та вимагають певного часу для того, щоб навчитися ними користуватися. Тому в даній курсовій роботі була поставлена задача розробити програму знаходження значення функції у міжвузловій точці за допомогою інтерполяційної формули Бесселя, яка займала б небагато пам’яті та була б простою у користуванні.

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

ТЕХНІЧНЕ ЗАВДАННЯ

 

1.Основою для проведення роботи є навчальний план кафедри Електроніки ВНТУ.

Замовник - кафедра Електроніки ВНТУ.

Виконавець – студент гр.. МП – 02 Твердохліб А.М.

2. Мета роботи.

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

3. Етапи виконання роботи.

Зміст	Строки  виконання	Чим закінчується етап
1.Отримання і узгодження завдання Розробка ТЗ	2–й тиждень	Технічне завданя
2.Розробка методу рішення	5–й тиждень	Теоретичні відомості про метод, порівняльний аналіз, приклади застосування
3.Алгоритмізація	8–й тиждень	Блок-схеми алгоритмів та їх порівняльний аналіз
4.Розробка і налагодження програми	12–й тиждень	Опис, текст та результати тестування програми
5.Розробка документації на курсову роботу	15–й тиждень	Пояснювальна записка(ПЗ) до курсової роботи
6.Захист курсової роботи	16–й тиждень	ПЗ та виконуваний файл програми

4.1 Кількість вузлів таблично заданої функції – не більше 100.

4.2 Похибка обчислень – не більше 0,001.

4.3 Алгоритм повинен бути оптимізований за часом виконання та розміром.

4.4 Програма повинна бути розроблена за принципами структурного та модульного програмування.

5. Спосіб реалізації результатів.

Робота повинна закінчуватись передачею замовнику пояснювальної записки до курсової роботи та виконуваного файла програми.

6. Техніко економічне обґрунтування.

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

7. Порядок розгляду і приймання роботи:

7.1 Курсова робота приймається комісією в складі двох викладачів за участю керівника роботи.

7.2 Програма перевіряється шляхом тестування на комп’ютері тестових завдань, розроблених замовником і виконавцем роботи.

8. Додаткові відомості.

Дане ТЗ може змінюватись і корегуватись за спільною домовленістю замовника та виконавця.

ТЕОРЕТИЧНІ ВІДОМОСТІ. ІНТЕРПОЛяція ФУНКЦІЙ 1. ПОСТАНОВКА ЗАДАЧІ ІНТЕРПОЛЯЦІЇ

 

Нехай деяка функція у=f(х) задана таблицею (табл.1), тобто при значеннях аргументу х=х₀, х₁, ... , х_n функція f(х) приймає відповідні значення

у₀, у₁,... , у_n.

Таблиця 1

Таблиця експериментальних значень

x	x0	x1	x2	....	xn
y	y0	y1	y2	....	yn

Також нехай необхідно визначити значення у=f(х), (х_i-1<<х_i). Величина х= потрапляє між двома табличними значеннями, тому для обчислення значення функції необхідно запропонувати деякий характер її зміни між відомими експериментальними даними.

Інтерполяцію можна розглядати як процес визначення для даного аргументу х значення функції у=f(х) по її декількох відомих значеннях. При цьому розрізняють інтерполяцію у вузькому смислі, коли х знаходиться між x₀ і x_n, і екстраполювання, коли х знаходиться поза відрізком інтерполяції [x₀ , x_n].

Задача інтерполяції полягає в наступному. На відрізку [а, b] задані n+1 точки х₀, х₁, ... , х_n, що називаються вузлами інтерполяції, і значення деякої функції f(x) у цих точках.

f(x₀) = y₀;

f(x₁) = y₁;(1)

f(x_n) = y_n

Потрібно побудувати функцію Р_n(х) (інтерполюючу функцію), яка б задовольняла таким умовам:

P_n(x₀) = y₀;

P_n(x₁) = y₁;(2)

P_n(x_n) = y_n

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

Геометрично це означає, що потрібно знайти криву y=P_n(х) деякого визначеного типу, що проходить через задану систему точок М_i (х_i,у_i) (i=0,1,2,..,n). Очевидно, можна побудувати множину неперервних функцій, що будуть проходити через задані вузлові точки.[1]

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

У загалі залежність, якої підпорядковується функція, може бути апроксимована багаточленом ступеня n:

Р_n(x) = y = a₀ + a₁ ∙ x + a₂ ∙ x² + ... + a_n ∙ xⁿ. (3)

Таку задачу називають задачею параболічної інтерполяції (або інтерполюванням).

Загалом є багато інтерполяційних формул та методів. До них відносяться такі: інтерполяційні формули Гауса (дві), Стерлінга та Бесселя (які є похідними від формул Гауса), Ньютона (дві) та багато інших.

 

2. ПАРАБОЛІЧНА ІНТЕРПОЛЯЦІЯ

Для визначення коефіцієнтів багаточлена (3) необхідно мати n+1 вузлову точку. Аналітичне визначення коефіцієнтів інтерполяційного багаточлена для n+1 точки зводиться до рішення системи лінійних рівнянь n+1 порядку, кожне з яких являє собою вираз (3), записаний для визначеної вузлової точки

y_i = a₀ + a₁ ∙ x_i + a₂ ∙ x_i² + ... + a_n ∙ x_iⁿ,(4)

де i = 1, 2,. . . n+1.

Даним методом побудови інтерполяційного поліному зручно користуватися, маючи персональний комп’ютер і відповідні програми. Даний метод не є єдиним способом побудови інтерполяційного поліному. Інший підхід, яким часто користуються на практиці, називається методом Лагранжа.[2]

 

3. МЕТОД ЛАГРАНЖА

Нехай при х=х₀, х₁, ... , х_n функція f(х) приймає відповідно значення у₀, у₁,... , у_n. Багаточлен ступеня не вище n, що приймає у вузлових точках задані значення, має вид:

Р_n(х)=у=. (5)

Цей багаточлен (5) називається інтерполяційною формулою Лагранжа і має такі властивості: