Разработка открытого урока на тему Интерполяционный полином Лагранжа

Міністерство освіти і науки України
Державний вищий навчальний заклад
„Донецький технікум промислової автоматики”








Методична розробка
відкритого заняття


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

з дисципліни : „Числові методи”

спеціальності 5.080405 „Програмування для електронно-обчислювальної техніки і автоматизованих систем”.






Підготувала викладач
Неліна Н.В.











Донецьк – 2009


План заняття.

Група: 1ПО-07

Спеціальність: 5.080405 „Програмування для електронно-обчислювальної техніки і автоматизованих систем”.

Тема програми: Інтерполяція, екстраполяція та апроксимація функцій

Тема заняття: Інтерполяційний поліном Лагранжа.

Ціль заняття:
Дидактична: сформулювати теоретичні знання про вузли інтерполяції та цнтерполяційні поліноми.
Розвиваюча: розширювати знання зі спеціальності, розвивати логічне мислення, просторове уявлення, пам’ять.
Виховна: виховувати цілеспрямованість, здатність працювати в команді.
Методична: удосконалювати методику проведення лекційних занять з різноманітними формами вивчення нового матеріалу.

Тип заняття: Засвоєння нових знань.

Вид заняття: лекція.

Вид лекції: тематична.

Форма проведення заняття: інтеграційне заняття.

Міжпредметні зв’язки:
Забезпечуючи дисципліни: вища математика.
Забезпечувані дисципліни: „Об’єктно-орієнтоване програмування”, «Математичні методи дослідження операцій»
Метод забезпечення: робоча програма, методична розробка заняття, лекція з теми: „Інтерполяційний поліном Лагранжа», комп’ютерна програма „Тест” з тестами, задачі.
Література:
Обов’язкова:
Бахвалов Н.С. Численные методы- М: Наука, 1973-632с.
Демидович Б.П. Основы вычислительной техники - М: Наука, 1970-664с.
Калиткин Н.Н. Численные методы. М: Наука, 1975-512с..
4. Крылов В.И. Вычислительные методы. М: Наука, 1976-304с
Хід заняття.

1. Організаційний момент.
Перевірка підготовки аудиторії до заняття;
Привітання та перелік студентів;
Перевірка готовності студентів до заняття;
Відповіді на запитання студентів.

2. Підготовка студентів до вивчення нового матеріалу.
Ознайомлення студентів з темою та навчальними цілями заняття;
Мотиваційний момент.

3. Актуалізація опорних знань студентів.
1) Робота з картками (4 студента пишуть відповіді, потім обмінюються листками з відповідями і перевіряють один одного).
2) Опитування з допомогою комп’ютерної програми „Тест” (6 студента).
3) Фронтальне опитування.
Що таке похибка?
Які методи розв’язання лінійних рівнянь Вам відомі?
Чим відрізняються метод хорд та метод дотичних?
Чи можна отримати точний корінь рівняння, якщо використовувати числові методи?
Як називається метод, за допомогою якого можна відновити функцію?
Чи можна, використовуючи метод чисельного дифференцюювання, отримати точне значення похідної?
Як будується ітераційна послідовність?
Чим відрізняється метод Жордана-Гаусса від методу ітерацій?

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

План.
1. Способи задання функції.
2. Математична постановка задачі інтерполювання.
3. Інтерполяційний многочлен Лагранжа.
4. Похибка інтерполяційної формули Лагранжа

Способи задання функції.

В практичній діяльності постійно приходиться стикатися з необхідністю вияву форм зв’язків в процесах і явищах , а також необхідністю їх матиматичного опису.
Розглянемо такі форми зв’язку, для яких деяка величина У, що характеризує процес, залежить від сукупності незв’язаних між собою величин 13EMBED Equation.31415 таким чином, що кожному набору (13EMBED Equation.31415 ) відповідає єдине значення величини Y.
Така однозначна відповідність величини У сукупності незалежних змінних 13EMBED Equation.31415 називається функціональною залежністю, а сама змінна величина У- функцією змінних величин 13EMBED Equation.31415.
Формально можна записати 13EMBED Equation.31415.13EMBED Equation.31415
Якщо величина у є функцією одної незалежної величини х, то цей зв’язок можна представити відношенням
13EMBED Equation.31415
Наприклад. Якщо розглянути площу кола 13EMBED Equation.31415, то площа 13EMBED Equation.31415 є функцією незалежного радіуса
13EMBED Equation.31415
Якщо розглянути об’єм фігури, то він буде функцією трьох вимірів
13EMBED Equation.31415
Із курсу матиматичного аналізу відомі три способи задання функціональних залежностей: - аналітичний, графічний, табличний.
Найбільш зручним способом задання функціональної залежності 13EMBED Equation.31415 є аналітичний, тому що він прямо вказує дії і їх послідовність виконання над незалежною змінною х для отримання відповідного значення Y. Наприклад: зв’язок шляху із часом в рівноприскореному русі можна виразити
13EMBED Equation.31415
При аналітичному способі задання є можливість отримати значення Y для будь-якого фіксованого аргумента х з будь-якою точністю.
Але при цьому необхідно виконувати всю послідовність обчислень і сам спосіб не наглядний.
Ці недоліки відсутні у випадку графічного задання функції. Графіком 13EMBED Equation.31415 називають геометричне місце точок площини хОу, координати яких задовольняють рівняння 13EMBED Equation.31415.
Табличний спосіб задання функції поширений в техніці, фізиці, економіці і він частіше виникає в результаті обробки результатів експерименту або випробувань.
Перевагою табличного способу задання функції являється те, що для кожного значення незалежної змінної, розміщеної в таблиці, можна зразу найти відповідне значення функції.
Математична постановка задачі інтерполювання.
При наявності заданих таблично функцій часто виникає необхідність обчислення значень функції 13EMBED Equation.31415в точках, які відрізняються від значень аргумента , зафіксованих в таблиці. Такі задачі практики формалізуються як математичні задачі інтерполювання .
Нехай на відрізку 13EMBED Equation.31415 задана функція 13EMBED Equation.31415своїм n+1 значеннями
13EMBED Equation.31415
в точках 13EMBED Equation.31415, які називають вузлами інтерполяції. Необхідно знайти аналітичний вираз 13EMBED Equation.31415табульованої функції
y0
y1
y2
y3

yn

x0
x 1
x2
x3

xn

яка співпадає у вузлах інтерполяції із значенням заданої функції, тобто
13EMBED Equation.31415
Процес обчислення значень функції в точках 13EMBED Equation.31415, які відрізняються від вузлів інтерполяції, називають інтерполюванням функції 13EMBED Equation.31415.
Якщо аргумент x, для якого визначається наближене значення функції, належить заданому відрізку 13EMBED Equation.31415, то задача обчислення наближеного значення функції називається інтерполюванням у вузькому розумінні. Якщо аргумент знаходиться за межами відрізка інтерполювання 13EMBED Equation.31415 , то задача визначення функції в точці х називається екстраполяцією.
З геометричного погляду задача інтерполювання для функції однієї змінної полягає в знаходжені кривої 13EMBED Equation.31415певного класу, яка проходить через точки з координатами
13EMBED Equation.31415
Мал.1. Геометричне пояснення інтерполяції.
З мал.1. видно, що через дані точки можна провести багато кривих. Таким чином, задача пошуку функції 13EMBED Equation.31415 по кінцевому числу її значень стає невизначеною.
Ця задача стає однозначною, якщо як інтерполюючою функцією 13EMBED Equation.31415для функції 13EMBED Equation.31415, що задана (n+1) своїми значеннями, вибрати многочлен 13EMBED Equation.31415 степені не вище n, такий, що
13EMBED Equation.31415
Многочлен Fn(x), що задовольняє цим умовам, називають інтерполяційним многочленом, а відповідні формули – інтерполяційними формулами. Якщо функція F(x) належить класу степеневих функцій, то інтерполювання називається параболічним. Параболічне інтерполювання найзручніше, оскільки многочлени, які прості за формою і не мають особливих точок, можуть набувати довільних значень, їх легко обчислювати, диференціювати та інтегрувати. При інтерполюванні виникає ряд задач:
- вибір найбільш зручного способу побудови інтерполяційної функції для кожного випадку;
- оцінка похибки при заміні f(x) інтерполюючою функцією F(x) на відрізку [a, b];
- оптимальний вібір вузлів інтерполяції для отримання мінімальної похибки.
Інтерполяційний многочлен Лагранжа
Розглядатимемо задачу параболічного інтерполювання, яку сформулюємо так: в n+1 різних точках задається функція 13EMBED Equation.31415 і необхідно побудувати многочлен
13EMBED Equation.31415 (1)
степені n, що задовольняє умові
13EMBED Equation.31415 (i=0, 1, 2, , n) (2)
При такій постановці вузли інтерполяції x0, x1, x2,....., xn можуть довільно відрізнятися один від одного на відрізу [a,b], тобто вузли інтерполяції нерівновіддалені, або h= xi+1- xi((Const. (h-крок інтерполяції).
Для визначення 13EMBED Equation.31415коефіцієнтів многочлена (1), який задовольняє умові (2), запишемо систему (13EMBED Equation.31415) лінійних рівнянь виду:
13EMBED Equation.31415
Ця система має єдиний розв’язок, бо її визначник є визначником Вандермонда, який не дорівнює нулю, тому що вузли xi (i=0, 1, 2, , n) різні, тобто існує єдиний алгебраїчний многочлен виду (1), що задовольняє умові (2).
Многочлен 13EMBED Equation.31415, який задовольняє умові (2), називають інтерполяційним многочленом, наближену рівність 13EMBED Equation.31415- інтерполяційною формулою, а різницю 13EMBED Equation.31415 - залишковим членом інтерполяційної формули.
Для умови (2) (коли 13EMBED Equation.31415 шукатимемо інтерполяційний многочлен 13EMBED Equation.31415 у такому вигляді:
13EMBED Equation.31415 (3)
де коефіцієнти 13EMBED Equation.31415 (i=0, 1, 2, , n) невідомі. Кожен доданок виразу (3) є многочленом степеня 13EMBED Equation.31415, причому при кожному з коефіцієнтів ai (i=0, 1, 2, , n) множника 13EMBED Equation.31415 немає. Визначимо коефіцієнти ai (i=0, 1, 2, , n), використавши умову (2).
Поклавши в (3) 13EMBED Equation.31415, отримаємо 13EMBED Equation.31415
звідки
13EMBED Equation.31415.
Якщо в (3) покласти 13EMBED Equation.31415 , то
13EMBED Equation.31415
звідки
13EMBED Equation.31415.
Анологічно обчислюємо
13EMBED Equation.31415, (i=0, 1, 2, 3, , n)
Підставивши ці значення коефіцієнтів в (3), отримаємо вираз інтерполяційного многочлена
13EMBED Equation.31415 (4)
Многочлен Ln(x) виду (4) називають інтерполяційним многочленом Лагранжа, а наближену рівність
13EMBED Equation.31415 (5)
інтерполяційною формулою Лагранжа.
Інтерполяційний многочлен Лагранжа можна записити більш стисло. Для цього введемо многочлен 13EMBED Equation.31415 степеня виду
13EMBED Equation.31415 (6)
Продиференціювавши по х цей добуток, отримаємо
13EMBED Equation.31415
Поклавши 13EMBED Equation.31415 (i=0, 1, 2, , n), матимемо
13EMBED Equation.31415 (7)
Підставивши ( 6 ) і ( 7 ) в ( 4 ), знайдемо
13EMBED Equation.31415 (8)
Вирази 13EMBED Equation.31415, що є коефіцієнтами при yi у многочлені Лагранжа, називають коефіцієнтами Лагранжа.
Розглянемо два окремих випадки інтерполяційної формули Лагранжа (5).
Нехай 13EMBED Equation.31415 , тобто значення функції 13EMBED Equation.31415 задано в двох вузлаз 13EMBED Equation.31415 та 13EMBED Equation.31415 і вони дорівнюють 13EMBED Equation.31415 та 13EMBED Equation.31415 відповідно. З формули (4) отримаємо
13EMBED Equation.31415 (9)
Формулу (9) називають формулою лінійного інтерполювання. При лінійному інтерполюванні дуга кривої 13EMBED Equation.31415 на відрізку 13EMBED Equation.31415 замінюється відрізком прямої (9), що лежить між точками 13EMBED Equation.31415 і 13EMBED Equation.31415.
Нехай n=2. Функцію f задано в трьох вузлах xi (i=0, 1, 2) значенями yi. У цьому разі формула (4) має вид
13EMBED Equation.31415 (10)
Формулу (10) називають формулою квадратичного інтерполювання. При квадратичному інтерполюванні дуга кривої y=f(x) на відрзку [x0, x2] замінюється дугою параболи, що проходить через точки (хі, yi) (i=0, 1, 2).
Похибка інтерполяційної формули Лагранжа
Похибка інтерполяційної формули Лагранжа обчислюється за формулою
13EMBED Equation.31415 (11)
Якщо прийняти 13EMBED Equation.31415 то для абсолютної похибки інтерполяційної формули Лагранжа дістанему таку оцінку:
13EMBED Equation.31415 (12)
Приклад. З якою точністю можливо обчислення 13EMBED Equation.31415 за допомогою інтерполяційної формули Лагранжа для функції 13EMBED Equation.31415 при виборі вузлів інтерполяції x0=100; x1=121; x2=144?
Рішення: Знайдемо 13EMBED Equation.31415; 13EMBED Equation.31415; 13EMBED Equation.31415. Звідси 13EMBED Equation.31415=13EMBED Equation.31415 при 100(x(144. Згідно з нерівністю (12) отримаємо
13EMBED Equation.31415(1,6(10-3.

5. Питання для самоконтролю.
1. Що таке вузли інтерполювання?
2. Який вигляд має поліном Лагранжа?
3. Як обчислити похибку?

6. Відповіді на питання студентів.

7. Коментар роботи студентів.

8. Домашнє завдання.
Вивчити поняття та формули.
Розв’язати задачу:
0,68 0,80866 0,715
0,80 1,02964
0,88 1,20966
0,93 1,34087
0,99 1,52368

9. Підсумок заняття.








Завдання для оцінювання знань студентів.
Картки.

4. студента отримують картки з задачами, пишуть відповіді, потім обмінюються листками з відповідями і перевіряють один одного.
1. Обчислити інтеграл: 13 EMBED Equation.3 1415Р=13EMBED Equation.31415
х
0,23
0,12
0,9

у
1,345
2,071
1,982

2. Знайти лінійну функцію:



Питання комп’ютерного тесту.
1 варіант

Какое число называется приближенным?
Какая формула называетмя эмперической?
Каким из предложенных методов можно решить линейное уравнение?
Чем отличается метод итерации от метода Жордана-Гуасса для решения систем уравнений?
Какая формула называется квадратурной?
Метод итерации считается законченым, когда ...
В каком методе необходимо делить отрезок пополам для нахождения корня?
Округлить число 2,4783904 до сотих.
Решите систему доступным вам методом

Методом наименьших квадратов найти коэффициент а для х=0,2;0,3;0,42 и у=0,4;1,5;2,3
Решить интеграл с помощью формулы левых прямоугольников (n=4)

2 варіант

Абсолютная погрешность - это ..
Что вычисляют формулы на рисунке?
Для использования метода хорд надо знать ...
Чем отличается метод итерации от метода Жордана-Гуасса для решения систем уравнений?
Формула на рисунке - это формула ...
Каким из предложенных методов можно решить линейное уравнение?
В методе половинного деления на концах отрезка функция должна принимать значения ...
Округлить число 46, 5680274 до сотих
Методом наименьших квадратов найти коэффициент а для х=0,1;0,3;0,4 и у=1,4;1,6;2,3
Решите систему доступным вам методом

11. Решить интеграл с помощью формулы левых прямоугольников (n=4)


3 варіант

Относительная погрешность - это ..
Из каких этапов состоит построение эмпирических формул?
Метод итерации считается законченым, когда ...
Чем отличается метод итерации от метода Жордана-Гуасса для решения систем уравнений?
Для какого метода эта формула находит погрешность?
Для использования метода касательных надо знать ...
В методе половинного деления до какого этапа надо делить отрезок?
Округлить число 3, 0957 до десятих
Методом наименьших квадратов найти коэффициент а для х=0,1;0,6;0,72 и у=0,9;3,5;2,3
Решите систему доступным вам методом

Решить интеграл с помощью формулы правых прямоугольников (n=4)

Ключ до тестів.

1 варіант

1. Число, которое отличается от точного числа А и заменяет его в расчетах
2. Ту, которую получают в результате исследований
3. Графическим
4. Нахождением итерационной последовательности
5. Приближенного интегрирования
6. Итерационная последовательность сходится
7. Метод половинного деления
8. 2, 48
9. x1=0, x2=1, x3=-1
10. 106
11. 0,012

2 варіант

Абсолютная величина разности между числом и его приближеним
Крэффициенты для метода наименьших квадратов
Знаки производных
Нахождением итерационной последовательности
Трапеций
Графическим
Разных знаков
46, 57
60
x1=1, x2=2, x3=3
77,37

3 варіант

Отношение абсолютной погрешности к модулю точного числа
Восстановление общего вида функции и нахождение ее параметров
Итерационная последовательность сходится
Нахождением итерационной последовательности
Трапецій
Производные первого и второго порядка
Пока не будет достигнута заданная точность
3,1
29
x1=0, x2=-1, x3=-2
0,6592


Критерії оцінювання.

11-9 правильних відповідей – „5”.
8-6 –правильних відповідей – „4”.
5-3 правильні відповіді – „3”.
2 і менше правильних відповідей – оцінка „2”.










13PAGE 15


13PAGE 14215



f(x)

f(x)

y

y1

y0

x0

xn

x2

x1

13EMBED Equation.31415

x



Root EntryEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativexEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation NativeEquation Native