МІНІСТЕРСТВО ОСВІТИ ТА НАУКИ УКРАЇНИ
Вінницький національний технічний університет
Факультет автоматики і комп’ютерних систем управління
Кафедра автоматики та інформаційно-вимірювальної техніки
Дослідження методів інтерполяції
Пояснювальна записка
до курсової роботи з дисципліни "Моделювання на ЕОМ"
Завдання на курсову роботу
Виконав:
студент групи 2АВ-06 О. Г Кирста
Перевірив:
к.т.н., доц. В.В. Кабачій
Вінниця, 2009
Анотація
Кирста О.Г. Дослідження методів інтерполяції. Вінниця, ВНТУ. - 2006.
Українська мова, стор., іл.
В курсовій роботі на прикладі інтерполяції напруги вольтметра розглянуто методи інтерполяції: ітераційний та метод розподілених різниць. Методи алгоритмізовано і реалізовано на ЕОМ в середовищі мови програмування Turbo Pascal 7.0. Наводяться блок-схема і лістинг програми.
Зміст
Вступ
1. Короткі теоретичні відомості
1.1 Задача інтерполяції
1.2 Інтерполяційна формула Лагранжа
1.3 Метод розподілених різниць
2. Алгоритми методів
3. Алгоритмізація розв’язання задачі
3.1 Вхідні та вихідні дані
3.2 Структура програми
3.3 Розробка плану дослідження роботи програми
4. Розв’язання задачі в пакеті прикладних програм
5. Аналіз результатів моделювання
6. Інструкція користувачеві
Висновки
Література
Додатки
В наш час, коли надзвичайно швидкими темпами розвивається наука і техніка, людина освоює все нові і нові галузі, все більше проникає як в надра землі так і за її межі, з’являється багато нових і досить складних задач, рішення яких потребує нових методів і нових підходів. Зокрема надзвичайно велика кількість задач електроніки, електротехніки, механіки, кібернетики та ряду інших галузей науки вимагають від вчених інженерів вирішення досить складних математичних задач, які вимагають певного аналізу та нестандартного підходу до вирішення.
З’являються задачі, які не можна розв’язати за допомогою класичної математики і отримати точний розв’язок, і взагалі досить часто про отримання точного розв’язку не доводиться говорити, оскільки отримати його при існуючих умовах просто неможливо. Тож ставляться задачі отримати приблизні розв’язки, але якомога близькі до точних. Тому в таких задачах використовуються різні наближені методи рішення тієї чи іншої задачі.
Чисельні методи направлені на вирішення задач, що виникають на практиці. Вирішення задачі чисельними методами зводиться до арифметичних та логічних операцій над числами, що вимагає застосування обчислювальної техніки. Умови і рішення задач найчастіше є наближеними, тобто мають похибки, причинами яких є невідповідність побудованої математичної моделі реальному об’єкту, похибки методів обчислень, похибки округлення тощо.
Обчислювальна техніка сьогодення є потужним засобом для фактичного виконання обчислень. Завдяки цьому в багатьох випадках стало можливим відмовитись від наближеного трактування прикладних задач і перейти до їх вирішення в точній постановці. Розумне використання сучасної обчислювальної техніки не можливе без вмілого застосування чисельних методів.
Мета інтерполяції - побудова функції
, яка приймає в окремих точках (вузли інтерполяції) значення, (1)що збігається з раніше заданими значеннями в цих точках невідомої функції
. Геометрично це означає, що потрібно знайти криву певного типу, яка проходить через систему точок (рисунок 1).Рисунок 1. Інтерполяція даних
В загальних випадках ця задача має нескінчену множину розв’язків чи зовсім не має розв’язку, але вона стає однозначною, якщо замість довільної функції
шукати поліном ступеня не вище , який задовольняє умову (1), тобто (2)Інтерполяційну формулу
, як правило, використовують для наближеного обчислення значень даної функції для . Така операція зветься інтерполяцією. Треба відзначити, що розрізняють інтерполяцію в вузькому розумінні, коли , та екстраполяцію, коли знаходиться за межами інтервалу , тобто чи . [1]Інтерполяція за Лагранжем вживається в загальному випадку для довільно розташованих вузлів.
Інтерполяційний поліном для методу Лагранжа представлений у вигляді:
, (3)де всі
(j=0,…, n) - поліноми ступеня n, коефіцієнти яких можна знайти з допомогою (n+1) рівняння: .Для полінома, який шукаємо, отримаємо:
(4)Формулу (1.8) називають інтерполяційний многочлен Лагранжа.
Треба відзначити дві головні властивості поліномів Лагранжа:
(5) 2)якщо
лінійно залежить від , то слушний принцип суперпозиції: інтерполяційний поліном суми декількох функцій дорівнює сумі інтерполяційних поліномів доданків.Похибка при інтерполяції за Лагранжем може бути оцінена таким чином:
(6)де
. (7)Інтерполяційні формули можна отримати, вживаючи визначення поділених різниць.
Так відношення
(8)називається розподіленою різницею першого порядку, а відношення
(9)розподіленою різницею другого порядку.
Розподілені різниці порядку n можна отримати з рекурентного співвідношення
. (10)Для інтерполяції можна використати інтерполяційну формулу Ньютона для нерівно віддалених значень аргументу: [1]
(11).Для інтерполяції методом розподілених різниць спочатку необхідно обчислити за формулою (8) розподілені різниці першого порядку, далі за формулою (9) - розподілені різниці 2 порядку, і так визначати розподілені різниці до порядку n-1, де n - кількість відомих значень функцій при відповідних відомих значеннях аргументу, використовуючи формулу (10). Відповідні обчислення просто реалізувати за допомогою циклу, в я кому значення розподілених різниць вищих порядків будуть обчислюватися за відомими значеннями аргументів та функцій і обчислених на попередньому кроці розподілених різниць нижчих порядків.
Після цього за формулою (11) здійснюється інтерполяція в необхідних точках.
При інтерполяції методом Лагранжа використовується формула (4). Спочатку за формулою (5) визначаються коефіцієнти
, які далі множаться на відповідні значення функції, і шуканий інтерполяційний поліном визначається як сума добутків коефіцієнтів на відповідні значення функції.Для розв’язання поставленої задачі потрібні певні вхідні данні, на основі яких будуть проводитись обчислення. Вхідними даними для інтерполяції є значення напруги при певних значеннях температури. Вони подаються у вигляді тарировочної таблиці 1.
Таблиця 1. Тарировочна таблиця для термопари
20 | 40 | 60 | 80 | 100 | 120 | |
-0,67 | -0,25 | -0,17 | 0,61 | 1,06 | 1,52 |
Вихідними даними є значення напруги при температурах
.