
- •«Численные методы»
- •Оглавление
- •§1. Теоретические основы численных методов 10
- •§2. Особенности математических вычислений на эвм. Погрешности вычислений 13
- •§3. Численные методы решения алгебраических и трансцендентных уравнений 25
- •§4. Методы решения систем уравнений 38
- •Введение
- •Из истории вычислительной математики
- •§1. Теоретические основы численных методов
- •§2. Особенности математических вычислений на эвм. Погрешности вычислений
- •Задания для самостоятельного решения
- •Практическая работа №1
- •Примеры выполнения заданий работы
- •Контрольные вопросы.
- •§3. Численные методы решения алгебраических и трансцендентных уравнений
- •3.1 Задача решения алгебраических и трансцендентных уравнений
- •3.2 Локализация корней
- •3.3 Метод деления отрезка пополам (метод бисекции, метод дихотомии)
- •3.4 Метод простой итерации
- •Задания для самостоятельного решения
- •Практическая работа №2
- •Примеры выполнения заданий работы
- •Контрольные вопросы.
- •3.5 Методы Ньютона
- •3.6. Решение уравнений с помощью табличного процессора Excel
- •Задания для самостоятельного решения
- •Практическая работа №3
- •Примеры выполнения заданий работы
- •Контрольные вопросы.
- •§4. Методы решения систем уравнений
- •4.1 Система линейных уравнений
- •4.1.1 Прямые методы решения систем линейных уравнений
- •Задания для самостоятельного решения
- •Практическая работа №4
- •Примеры выполнения заданий работы
- •4.1.2 Вычисление определителей и обратной матрицы
- •Задания для самостоятельного решения
- •Контрольные вопросы.
- •4.1.3 Итерационные методы решения систем линейных уравнений. Метод простой итерации
- •4.2. Решение системы уравнений и вычисление определителя с помощью табличного процессора Excel
- •Задания для самостоятельного решения
- •Практическая работа №5
- •Примеры выполнения заданий работы
- •Контрольные вопросы
- •§5. Методы приближения и аппроксимации функций
- •5.1 Понятия интерполяции и экстраполяции
- •5.2 Интерполяционный многочлен Лагранжа
- •Задания для самостоятельного решения
- •5.3 Приближение функций с помощью табличного процессора Excel
- •Задания для самостоятельного решения
- •Практическая работа №6
- •Примеры выполнения заданий работы
- •Контрольные вопросы
- •5.3 Интерполяционные формулы Ньютона
- •Задания для самостоятельного решения
- •§6. Численное интегрирование
- •6.1 Задача численного интегрирования
- •6.2 Методы прямоугольников и трапеций
- •6.3 Метод Симпсона (метод парабол)
- •Задания для самостоятельного решения
- •Практическая работа №7
- •Примеры выполнения заданий работы
- •Контрольные вопросы
- •6.4 Квадратурная формула Гаусса
- •6.5. Вычисление интеграла с использованием табличного процессора Excel.
- •Задания для самостоятельного решения
- •Практическая работа №8
- •Примеры выполнения заданий работы
- •Контрольные вопросы
- •§7. Решение обыкновенных дифференциальных уравнений
- •7.1. Задача численного решения обыкновенных дифференциальных уравнений
- •7.2. Методы Эйлера
- •7.3 Метод Рунге – Кутта
- •Задания для самостоятельного решения
- •Практическая работа №9
- •Примеры выполнения заданий работы
- •Контрольные вопросы
- •§8. Методы оптимизации
- •8.1 Методы одномерной оптимизации
- •Задания для самостоятельного решения
- •8.2 Методы многомерной оптимизации
- •8.3. Решение задач оптимизации с помощью табличного процессора Excel
- •Задания для самостоятельного решения
- •Практическая работа №10
- •Примеры выполнения заданий работы
- •Контрольные вопросы
- •Литература
5.2 Интерполяционный многочлен Лагранжа
Наиболее распространенными является интерполирование многочленами.
Пусть функция y = f(x) задана таблично:
хi |
x0 |
x1 |
… |
xi |
… |
xn |
yi |
y0 |
y1 |
… |
yi |
… |
yn |
Общий вид интерполяционного многочлена Лагранжа:
(5.1)
Погрешность формулы Лагранжа есть R n (x) =f (x) −L n (x) .
Получаем выражение для погрешности формулы Лагранжа:
,
(5.2)
где Пn+1(x)
= (х– х0)(х–
х1)·
… · (х– хn)
Пример. Используя
интерполяционный многочлен Лагранжа
вычислить в указанных точках
и
значения функции
,
заданной в виде таблицы:
|
|
|
0 |
0,5 |
-0,5 |
1 |
1,5 |
0,7 |
2 |
2,5 |
0,3 |
и оценить погрешность
для заданного значения третьей производной
.
Решение:
Будем использовать многочлен Лагранжа 2-ой степени:
.
Погрешность интерполяции с помощью этого многочлена Лагранжа вычисляется как:
.
Значение функции в точках и , используя многочлен Лагранжа, будет равно:
для
для
.
Абсолютная погрешность интерполяции:
при
равна:
при
равна:
Относительная погрешность интерполяции:
при
равна:
при
равна:
Применение формулы Лагранжа приводит к большому числу однотипных вычислений, поэтому вычислять значения функции, по многочлену Лагранжа, удобнее с использованием вычислительной таблицы (табл. 6.1)
Таблица 6.1
х = х* |
х0 |
х1 |
х2 |
… |
рi |
yi |
yi / pi |
х0 |
х*– х0 |
х0– х1 |
х0– х2 |
… |
р0 |
y0 |
y0 / p0 |
х1 |
х1– х0 |
х*– х1 |
х1– х2 |
… |
р1 |
y1 |
y1 / p1 |
х2 |
х2– х0 |
х2– х1 |
х* –х2 |
… |
р2 |
y2 |
y2 / p2 |
… |
… |
… |
… |
… |
… |
… |
… |
|
|
|
|
|
|
|
S=Σ(yi / pi) |
рi – произведение разностей по строкам: рi = (х*– х0)(х0– х1)· … · (х0– хn)
Многочлен Лагранжа имеет вид:
Ln(x) = S ∙ (х*– х0)( х*– х1)( х* –х2)∙ … · (х* – хn).
Задания для самостоятельного решения
Используя интерполяционный многочлен Лагранжа вычислить в указанных точках и
значения функции , заданной в виде таблицы:
|
|
|
0 |
0,0 |
0,1 |
1 |
0,8 |
0,3 |
2 |
1,6 |
0,3 |
и оценить погрешность
для заданного значения третьей производной
.
Используя интерполяционный многочлен Лагранжа вычислить в указанных точках и
значения функции , заданной в виде таблицы:
|
|
|
0 |
3,5 |
0,2 |
1 |
6,5 |
0,5 |
2 |
9,5 |
-0,9 |
и оценить погрешность
для заданного значения третьей производной
.
Используя интерполяционный многочлен Лагранжа вычислить в указанных точках
и
значения функции , заданной в виде таблицы:
|
|
|
0 |
0,1 |
-7,5 |
1 |
1,5 |
-0,67 |
2 |
2,0 |
7,3 |
и оценить погрешность для заданного значения третьей производной .
4. Используя
вычислительную схему Лагранжа вычислить
в указанных точках
и
значения функции
,
заданной в виде таблицы:
|
|
|
0 |
0,0 |
1,0 |
1 |
0,5 |
-1,7 |
2 |
1,0 |
-0,5 |
5. Используя
вычислительную схему Лагранжа вычислить
в указанных точках
и
значения функции
,
заданной в виде таблицы:
|
|
|
0 |
1,0 |
-0,5 |
1 |
3,5 |
-2,2 |
2 |
6,0 |
-0,6 |