 
        
        - •«Численные методы»
- •Оглавление
- •§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.3 Приближение функций с помощью табличного процессора Excel
В табличном процессоре Excel для интерполяции с помощью формулы Лагранжа можно воспользоваться соответствующей таблицей (таблица 5.1).
А также для аппроксимации Excel имеет удобные встроенные средства. Самый наглядный способ – использование графического представления данных:
- Необходимо ввести таблицу данных и построить по ним точечную диаграмму; 
- Выделив диаграмму, через главное меню выполнить команды Диаграмма  Добавить линию тренда (рис. 5.1) Тренд (тенденция) – термин для обозначения линии графика аппраксимируемой функции. 
- В открывшемся диалоговом окне на закладке Тип необходимо выбрать один из видов функции, предлагаемых для аппроксимации (рис. 5.1). Таким образом перебирая различные виды функций можно подобрать линию тренда наиболее подходящую к заданным значениям 
- В том же диалоговом окне, перейдя к закладке Параметры необходимо установить галочку на флажке показывать уравнение на диаграмме – после выполнения команд на диаграмме появляется уравнение, которое соответствует линии тренда и является искомой функцией (рис. 5.2). 
В закладке Параметры можно также установить галочку на флажке поместить на диаграмму величину достоверности аппроксимации R^2. В статистике данная величина называется коэффициентом детерминированности,она показывает насколько точной получилась построенная функция (рис. 5.2).
 
Рисунок 5.1. Получение линии тренда при аппроксимации функции в Excel
 
Рисунок 5.2. Результат аппроксимации функции в Excel
- Интерполяцию или экстраполяцию можно получить, используя функцию ПРЕДСКАЗ в мастере функций (рис. 5.3). 
 
Рисунок 5.3. Нахождение интерполяционного значения функции в табличном процессоре Excel
Задания для самостоятельного решения
| x | 2,0 | 2,2 | 2,4 | 2,6 | 2,8 | 3,0 | 
| f(x) | 8,69 | 11,44 | 14,70 | 18,53 | 22,98 | 28,10 | 
| x | 1,4 | 1,5 | 1,6 | 1,7 | 1,8 | 1,9 | 
| f(x) | -1,67 | -1,80 | -1,93 | -2,05 | -2,18 | -2,30 | 
Практическая работа №6
Тема: «Составление интерполяционного многочлена Лагранжа»
Цели: освоение применения интерполяционного многочлена Лагранжа для вычисления значения функции по промежуточному значению аргумента;
Задание1. Для функции, заданной таблицей узловых значений, составить формулу интерполяционного многочлена Лагранжа, вычислить с его помощью одно значение функции для промежуточного значения аргумента х*.
Исходные данные для задания 1:
Вариант 1. х* = 0,7
| х | -1 | 0 | 3 | 4 | 
| f(x) | -3 | 5 | 2 | -6 | 
Вариант 2. х* = 2,4
| х | 2 | 3 | 5 | 6 | 
| f(x) | 4 | 1 | 7 | 2 | 
Вариант 3. х* = 0,7
| х | 0 | 2 | 3 | 5 | 
| f(x) | -1 | -4 | 2 | -8 | 
Вариант 4. х* = 7,3
| х | 7 | 9 | 13 | 15 | 
| f(x) | 2 | -2 | 3 | -4 | 
Вариант 5. х* = 3,6
| х | -3 | -1 | 3 | 5 | 
| f(x) | 7 | -1 | 4 | -6 | 
Вариант 6. х* = 1,7
| х | 1 | 2 | 4 | 7 | 
| f(x) | -3 | -7 | 2 | 8 | 
Вариант 7. х* = 2,7
| х | 2 | 4 | 5 | 7 | 
| f(x) | 9 | -3 | 6 | -2 | 
Вариант 8. х* = 0,7
| х | -4 | -2 | 0 | 3 | 
| f(x) | 2 | 8 | 5 | 10 | 
Вариант 9. х* = 2,8
| х | 2 | 4 | 7 | 8 | 
| f(x) | -1 | -6 | 3 | 12 | 
Вариант 10. х* = 0,7
| х | 0 | 1 | 4 | 6 | 
| f(x) | 7 | -1 | 8 | 2 | 
Вариант 11. х* = 2,7
| Х | -1 | -1 | 2 | 4 | 
| f(x) | 4 | 9 | 1 | 6 | 
Вариант 12. х* = 3,4
| Х | -1 | 1,5 | 3 | 5 | 
| f(x) | 4 | -7 | 1 | -8 | 
Вариант 13. х* = – 4,5
| Х | -9 | -7 | - 4 | -1 | 
| f(x) | 3 | -3 | 4 | -9 | 
Вариант 14. х* = 1,3
| Х | -8 | -5 | 0 | 2 | 
| f(x) | 9 | -2 | 4 | 6 | 
Вариант 15. х* = –3,6
| Х | -7 | -5 | - 4 | -1 | 
| f(x) | 4 | - 4 | 5 | 10 | 
Вариант 16. х* = 2,7
| Х | 1 | 4 | 9 | 11 | 
| f(x) | -2 | 9 | 3 | -7 | 
Вариант 17. х* = 8,5
| Х | 7 | 8 | 9 | 13 | 
| f(x) | 6 | -2 | 7 | -10 | 
Вариант 18. х* = 0,7
| Х | -4 | 0 | 2 | 5 | 
| f(x) | 4 | 8 | -2 | -9 | 
Вариант 19. х* = 2,8
| Х | -3 | -1 | 1 | 3 | 
| f(x) | 11 | -1 | 6 | -2 | 
Вариант 20. х* = 0,7
| х | 0 | 3 | 8 | 11 | 
| f(x) | 1 | 5 | - 4 | -8 | 
Задание 2. Вычислить значение заданной функции f(x) для промежуточного значения аргумента х* по интерполяционному многочлену Лагранжа с помощью вычислительной таблицы. Сравните результат интерполирования со значением функции, которое вычислите, используя ее аналитическое выражение.
Задание 3. Для таблично заданной функции получить эмпирическую формулу с помощью математических программных средств (табличный процессор Excel). Сравнить полученную формулу с заданным аналитическим выражением.
Исходные данные для задания 2, 3:
Вариант 1.
  х*
= 3,8       
 
  
| х | 1,3 | 2,1 | 3,7 | 4,5 | 6,1 | 7,7 | 8,5 | 
| f(x) | 1,7777 | 4,5634 | 13,8436 | 20,3952 | 37,3387 | 59,4051 | 72,3593 | 
Вариант 2.
  х*
= 3,5       
 
  
| х | 1,2 | 1,9 | 3,3 | 4,7 | 5,4 | 6,8 | 7,5 | 
| f(x) | 0,3486 | 1,0537 | 1,7844 | 2,2103 | 2,3712 | 2,6322 | 2,7411 | 
Вариант 3.
  х*
= 0,5       
 
  
| х | -3,2 | -0,8 | 0,4 | 2,8 | 4,0 | 6,4 | 7,6 | 
| f(x) | -1,9449 | -0,6126 | 0,3097 | 1,8068 | 2,0913 | 1,4673 | 0,6797 | 
Вариант 4. х* = 4,8
| х | 1,3 | 2,1 | 3,7 | 4,5 | 6,1 | 7,7 | 8,5 | 
| f(x) | 1,7777 | 4,5634 | 13,8436 | 20,3952 | 37,3387 | 59,4051 | 72,3593 | 
Вариант 5. х* = 4,1
| х | 1,3 | 2,1 | 3,7 | 4,5 | 6,1 | 7,7 | 8,5 | 
| f(x) | 1,7777 | 4,5634 | 13,8436 | 20,3952 | 37,3387 | 59,4051 | 72,3593 | 
Вариант 6. х* = 3,9
| х | 1,2 | 1,9 | 3,3 | 4,7 | 5,4 | 6,8 | 7,5 | 
| f(x) | 0,3486 | 1,0537 | 1,7844 | 2,2103 | 2,3712 | 2,6322 | 2,7411 | 
Вариант 7. х* = 3,3
| х | -3,2 | -0,8 | 0,4 | 2,8 | 4,0 | 6,4 | 7,6 | 
| f(x) | -1,9449 | -0,6126 | 0,3097 | 1,8068 | 2,0913 | 1,4673 | 0,6797 | 
Вариант 8. х* = 4,0
| х | 1,3 | 2,1 | 3,7 | 4,5 | 6,1 | 7,7 | 8,5 | 
| f(x) | 1,7777 | 4,5634 | 13,8436 | 20,3952 | 37,3387 | 59,4051 | 72,3593 | 
Вариант 9. х* = 2,9
| х | 1,3 | 2,1 | 3,7 | 4,5 | 6,1 | 7,7 | 8,5 | 
| f(x) | 1,7777 | 4,5634 | 13,8436 | 20,3952 | 37,3387 | 59,4051 | 72,3593 | 
Вариант 10. х* = 5,3
| х | 1,2 | 1,9 | 3,3 | 4,7 | 5,4 | 6,8 | 7,5 | 
| f(x) | 0,3486 | 1,0537 | 1,7844 | 2,2103 | 2,3712 | 2,6322 | 2,7411 | 
Вариант 11. х* = 4,1
| х | -3,2 | -0,8 | 0,4 | 2,8 | 4,0 | 6,4 | 7,6 | 
| f(x) | -1,9449 | -0,6126 | 0,3097 | 1,8068 | 2,0913 | 1,4673 | 0,6797 | 
Вариант 12. х* = 6,6
| х | 1,2 | 1,9 | 3,3 | 4,7 | 5,4 | 6,8 | 7,5 | 
| f(x) | 0,3486 | 1,0537 | 1,7844 | 2,2103 | 2,3712 | 2,6322 | 2,7411 | 
Вариант 13. х* = 4,4
| х | 1,3 | 2,1 | 3,7 | 4,5 | 6,1 | 7,7 | 8,5 | 
| f(x) | 1,7777 | 4,5634 | 13,8436 | 20,3952 | 37,3387 | 59,4051 | 72,3593 | 
Вариант 14. х* = 5,2
| х | -3,2 | -0,8 | 0,4 | 2,8 | 4,0 | 6,4 | 7,6 | 
| f(x) | -1,9449 | -0,6126 | 0,3097 | 1,8068 | 2,0913 | 1,4673 | 0,6797 | 
Вариант 15. х* = 3,7
| х | 1,2 | 1,9 | 3,3 | 4,7 | 5,4 | 6,8 | 7,5 | 
| f(x) | 0,3486 | 1,0537 | 1,7844 | 2,2103 | 2,3712 | 2,6322 | 2,7411 | 
Вариант 16. х* = 0,4
| х | 1,3 | 2,1 | 3,7 | 4,5 | 6,1 | 7,7 | 8,5 | 
| f(x) | 1,7777 | 4,5634 | 13,8436 | 20,3952 | 37,3387 | 59,4051 | 72,3593 | 
Вариант 17. х* = 7,5
| х | -3,2 | -0,8 | 0,4 | 2,8 | 4,0 | 6,4 | 7,6 | 
| f(x) | -1,9449 | -0,6126 | 0,3097 | 1,8068 | 2,0913 | 1,4673 | 0,6797 | 
Вариант 18. х* = 2,5
| х | 1,2 | 1,9 | 3,3 | 4,7 | 5,4 | 6,8 | 7,5 | 
| f(x) | 0,3486 | 1,0537 | 1,7844 | 2,2103 | 2,3712 | 2,6322 | 2,7411 | 
Вариант 19. х* = 6,8
| х | 1,3 | 2,1 | 3,7 | 4,5 | 6,1 | 7,7 | 8,5 | 
| f(x) | 1,7777 | 4,5634 | 13,8436 | 20,3952 | 37,3387 | 59,4051 | 72,3593 | 
Вариант 20. х* = – 1,7
| х | -3,2 | -0,8 | 0,4 | 2,8 | 4,0 | 6,4 | 7,6 | 
| f(x) | -1,9449 | -0,6126 | 0,3097 | 1,8068 | 2,0913 | 1,4673 | 0,6797 | 
