- •Методические указания по выполнению лабораторных и расчетно-графических работ
- •1 Лабораторная работа №1. Методы интерполяции.
- •1.1 Расчеты в точке
- •1.2 Расчеты в точке с большей кривизной
- •1.3 Проверка проведенных вычислений по компьютерной программе
- •1.4 Исследование влияния числа узлов сетки на точность интерполяции
- •2 Лабораторная работа №2. Численные методы дифференцирования функции одной переменной. Формулы Стирлинга
- •3 Лабораторная работа №3. Численные методы решения нелинейных уравнений. Методы дихотомии, хорд, касательных
- •3.1 Задача отделения корня.
- •3.2 Нахождение корня методом дихотомии
- •3.3 Нахождение корня методом хорд
- •3.4 Нахождение корня методом касательных (Ньютона)
- •4 Лабораторная работа №4. Численные методы решения систем нелинейных уравнений. Метод Ньютона.
- •5 Лабораторная работа №5. Численные методы поиска экстремума функций одной переменной. Методы дихотомии, «золотого сечения», Ньютона.
- •5.1 Поиск минимума методом дихотомии
- •5.2 Поиск минимума методом «золотого сечения»
- •5.3 Поиск минимума методом Ньютона
- •Выводы: При поиске минимума функции на отрезке
- •6 Лабораторная работа №6. Численное решение обыкновенных дифференциальных уравнений методами Эйлера, модифицированного метода Эйлера с пересчетом, Рунге-Кутты.
- •6.5 Решение с помощью программного комплекса «чмриз».
- •7 Лабораторная работа №7. Вычисление определенных интегралов методом Монте-Карло.
Методические указания по выполнению лабораторных и расчетно-графических работ
1 Лабораторная работа №1. Методы интерполяции.
Задание. Протабулировать произвольно заданную нелинейную функцию наотрезкес постоянным шагом,где n – число точек (узлов) таблицы n = 4…10. По результатам табулирования функции построить таблицу
|
|
|
… |
|
|
|
|
… |
|
Используя линейную и квадратичную интерполяцию, вычислить значения функции в двух произвольных, не лежащих в узлах таблицы, точках и,в которых функция имеет разную кривизну. Вычислить в процентах относительную погрешность интерполяции по отношению к точному значению функции в точках и,и сделать соответствующие выводы.
Решение.
Пусть функция имеет вид – т.е. многочлен четвертого порядка. (Если выбранная функция – многочлен, то порядок его должен быть задан выше двух, иначе при использовании квадратичной интерполяции получим абсолютно точный результат). Протабулируем функцию на отрезкес постоянным шагом,где n=4 – число точек (узлов) таблицы. По результатам табулирования функции построим таблицу 1.1 и график функции (рисунок 1.1)
Таблица 1.1
|
|
|
1,25 |
|
|
|
|
2 |
|
Рисунок 1.1 График функции
Выберем две точки и ,в которых функция имеет разную кривизну.
1.1 Расчеты в точке
Вычислим точное значение функции в этой точке: .
Вычислим значение функции, используя формулы линейной интерполяции:
, т.е.
Точка находится между узламии.В этом случае
Относительная погрешность линейной интерполяции в процентах равна
.
Вычислим теперь значение функции в точке , используя формулыквадратичной интерполяции. Найдем приближенное значение функции с помощью формулы квадратичной интерполяции
.
Для нахождения коэффициентов составим систему уравнений (2.14) (смотри теоретический материал)
Так как точка ближе к узлу,чемк узлу , то выбираем для расчета следующие три узла: ,,.Соответственно ,,.
Тогда система запишется в виде:
Решая эту систему, находим . Искомое значение функции .
Относительная погрешность квадратичной интерполяции в процентах равна
.
1.2 Расчеты в точке с большей кривизной
Вычислим точное значение функции в этой точке: .
Вычислим значение функции, используя формулы линейной интерполяции:
, т.е.
Точка находится между узламии.В этом случае
.
Относительная погрешность линейной интерполяции в процентах равна
Вычислим теперь значение функции в точке , используя формулыквадратичной интерполяции. Найдем приближенное значение функции с помощью формулы квадратичной интерполяции
.
Для нахождения коэффициентов составим систему уравнений (2.14) (смотри теоретический материал)
Так как точка ближе к узлу,чемк узлу , то выбираем для расчета следующие три узла: ,,.Соответственно ,,.
Тогда система запишется в виде:
Решая эту систему, находим . Искомое значение функции .
Относительная погрешность квадратичной интерполяции в процентах равна
.
1.3 Проверка проведенных вычислений по компьютерной программе
1.4 Исследование влияния числа узлов сетки на точность интерполяции
Проведенные выше расчеты показали, что при выбранном числе узлов сетки порешность интерполяции слишком большая. Для повышения точности интерполяции следует увеличить число узлов сетки, или, что, то же самое, увеличить шаг при табулировании заданной функциина отрезке, (). Выберем число узлов n = 7 так, чтобы шаг был равным . По результатам табулирования функциипостроим таблицу 1.2.
Таблица 1.2
|
|
1,0156 |
|
1,0156 |
1,25 |
2,2656 |
|
|
|
0,5 |
|
1,5 |
2 |
2,5 |
|
Проведем для тех же точек и аналогичныерасчеты по компьютерной программе.
Выводы:
1) Применение квадратичной интерполяции для данной функции – многочлена четвертой степени, дает более точные результаты, чем линейная интерполяция.
2) При увеличении кривизны кривой погрешность интерполяции увеличивается.
3) Для повышения точности интерполяции (если погрешности интерполяции не устраивают расчетчика) следует увеличить число узлов сетки, особенно в областях с большой кривизной кривой .