49566 (Дослідження методів інтерполяції)

_{МІНІСТЕРСТВО} ОСВІТИ ТА НАУКИ УКРАЇНИ

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

Факультет автоматики і комп’ютерних систем управління

Кафедра автоматики та інформаційно-вимірювальної техніки

Дослідження методів інтерполяції

Пояснювальна записка

до курсової роботи з дисципліни "Моделювання на ЕОМ"

Завдання на курсову роботу

Виконав:

студент групи 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)

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

Рисунок 1. Інтерполяція даних

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

(2)

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

1.2 Інтерполяційна формула Лагранжа

Інтерполяція за Лагранжем вживається в загальному випадку для довільно розташованих вузлів.

Інтерполяційний поліном для методу Лагранжа представлений у вигляді:

, (3)

де всі (j=0,…, n) - поліноми ступеня n, коефіцієнти яких можна знайти з допомогою (n+1) рівняння: .

Для полінома, який шукаємо, отримаємо:

(4)

Формулу (1.8) називають інтерполяційний многочлен Лагранжа.

Треба відзначити дві головні властивості поліномів Лагранжа:

(5) 2)

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

Похибка при інтерполяції за Лагранжем може бути оцінена таким чином:

(6)

де . (7)

1.3 Метод розподілених різниць

Інтерполяційні формули можна отримати, вживаючи визначення поділених різниць.

Так відношення

(8)

називається розподіленою різницею першого порядку, а відношення

(9)

розподіленою різницею другого порядку.

Розподілені різниці порядку n можна отримати з рекурентного співвідношення

. (10)

Для інтерполяції можна використати інтерполяційну формулу Ньютона для нерівно віддалених значень аргументу: [1]

(11).

2. Алгоритми методів

Для інтерполяції методом розподілених різниць спочатку необхідно обчислити за формулою (8) розподілені різниці першого порядку, далі за формулою (9) - розподілені різниці 2 порядку, і так визначати розподілені різниці до порядку n-1, де n - кількість відомих значень функцій при відповідних відомих значеннях аргументу, використовуючи формулу (10). Відповідні обчислення просто реалізувати за допомогою циклу, в я кому значення розподілених різниць вищих порядків будуть обчислюватися за відомими значеннями аргументів та функцій і обчислених на попередньому кроці розподілених різниць нижчих порядків.

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

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

3. Алгоритмізація розв’язання задачі

3.1 Вхідні та вихідні дані

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

Таблиця 1. Тарировочна таблиця для термопари

	20	40	60	80	100	120
	-0,67	-0,25	-0,17	0,61	1,06	1,52

Вихідними даними є значення напруги при температурах .

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

Таблиця 2. Вхідні та вихідні данні

№	Назва змінної	Ідентифікатор	Тип
1	Значення Т,	t [i]	Real
2	Значення U, мВ	u [i]	Real
3	Результат U (75 ), мВ	РР	Real
4	Результат U (93 ), мВ	ZZ	Real

3.2 Структура програми

Программа використовує функції стандартних модулів Crt і Graph. Модуль Crt забезпечує роботу з клавіатурою та екраном в текстовому режимі, а модуль Graph забезпечує роботу з графікою. [2]

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

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

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

Рисунок 2. Логічна схема основної програми

3.3 Розробка плану дослідження роботи програми

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

Перевіримо правильність роботи програми. Запускаємо програму. При запуску програми (ехе-файл чи з середовища Турбо Паскаль) з'являється головне меню програми (рисунок 2). Вибір необхідного пункту меню здійснюється натисканням відповідних цифрових клавіш на клавіатурі, вихід з програми - послідовно клавіш. "Esc" та "Enter".

Рисунок 3. Меню програми

Натискаємо клавішу "Esc", за нею - "Enter" - програма закривається.