Вопросы для самопроверки по теме 1.2

1. В чём состоит задача интерполяции функции?

2. Какие критерии согласия обеспечивают совпадение неизвестной функции с интерполирующей?

3. Как называется интерполяция многочленами первой и второй степени?

4. Напишите общие формулы конечных разностей 1-го, 2-го и 3-го порядков.

5. Напишите формулу интерполяционного многочлена Ньютона для пяти узлов.

6. Чему равна третья производная при трёх узлах интерполирования?

1.3. Численное интегрирование

Изучаемые вопросы: Квадратурные формулы Ньютона-Котеса. Квадратурные формулы Гаусса. Задачи оптимизации. Формулы Эйлера и Грегори. Формулы Ромберга. Стандартные программы численного интегрирования. Построение программ с автоматическим выбором шага интегрирования.

Здесь также после изучения материала опорного конспекта и письменных лекций Вам следует решить одну из задач контрольной работы согласно «Методическим указаниям к выполнению контрольной работы.

1.3.1. Приближенное вычисление определенного интеграла

Простейшие формулы для приближённого вычисления определённого интеграла называются квадратурными. В многомерном случае их называют также кубатурными. К простейшим квадратурным формулам относятся формулы прямоугольников, трапеций и формула Симпсона, объединённые общим названием – квадратурные формулы Ньютона-Котеса. Все эти формулы основаны на свойстве аддитивности определённого интеграла, а именно: интеграл по сумме отрезков равен сумме интегралов по этим отрезкам. Поэтому, если нужно вычислить определённый интеграл от некоторой функции вдоль отрезка, то его можно представить в виде суммы интегралов по частичным отрезкам разбиения интервала: , где .

Задача состоит в выборе достаточного числа разбиений отрезка (отрезки, как правило, выбираются одинаковыми), и удачной замене подынтегральной функции. Обычно она заменяется интерполяционным многочленом степени:

, (1)

где – остаточный член интерполяции.

Т. о., на каждом частичном промежутке

,

где– приближённое значение интеграла на частичном промежутке, а– величина ошибки на том же промежутке.

Соответственно, приближённое значение интеграла ,(2)

а ошибка

. (3)

На рис. 1 представлена геометрическая интерпретация определённого интеграла, как площади криволинейной трапеции, ограниченной осью ОХ, графиком функции и прямыми , и интегралана частичном промежутке. (Заштрихованная криволинейная трапеция).

Заметим здесь, что если считать шаг разбиения в методе Симпсона равным целому, без деления пополам, то в расчётах, вместо формулы (2.16) (п.2.4 Учебного пособия), можно использовать следующую:

. (4)

Соответствующие формулы, вместе с оценками погрешностей и примерами вычислений Вы можете найти в Учебном пособии.

Более полное изложение этой темы – в [7], c.86-163.

Вопросы для самопроверки по теме 1.3

1. Напишите формулы прямоугольников, трапеции и Симпсона.

2. Сформулируйте обобщённую теорему о среднем.

1.4. Приближение функций

Из всех вопросов темы 1.4. ­Приближение функций изучается лишь метод наименьших квадратов. Вопросы этой темы не содержатся в контрольной работе, поэтому здесь приводятся только основные теоретические положения.

Метод наименьших квадратов

Пусть известно, что величины исвязаны некоей функциональной зависимостью. Требуется приближенно определить эту функциональную зависимость по экспериментальным данными. Предположим, что в результате измерений получен ряд экспериментальных точек. Мы уже знаем, что через точек всегда можно провести кривую, аналитически выражаемую многочленом- ой степени. Этот многочлен называют интерполяционным. Вообще, замену функции на функцию так, что их значения совпадают в заданных точках

, , (1)

называют интерполяцией.

Однако такое решение проблемы не всегда является удовлетворительным, поскольку из-за случайных ошибок измерения и, возможно, случайной природы самих величин x и y. Т.о., можно записать, что

(2)

где – некоторая случайная ошибка. Поэтому требуется провести кривую так, чтобы она в наименьшей степени зависела от случайных ошибок. Эта задача называется сглаживанием (аппроксимацией) экспериментальной зависимости и часто решается методом наименьших квадратов. Сглаживающую кривую называют аппроксимирующей.

Задача аппроксимации решается следующим образом. В декартовой прямоугольной системе координат наносят точки . По виду расположения этих точек делается предположение о принадлежности искомой функции к определенному классу. Например, линейная , квадратичная и т.п. В общем случае . Неизвестные параметры функции определяются из требования минимума суммы квадратов случайных ошибок, т.е. минимума величины

. (3)

Величина называется также суммарнойневязкой. Необходимым условием минимума функции нескольких переменных является обращение в нуль частных производных невязки:

, . (4)

Решая систему уравнений (4), находят неизвестные параметры и тем самым полностью определяют функцию, которая наилучшим образом (в смысле наименьших квадратов отклонений от исходных точек или наименьшей суммарной невязки) аппроксимирует искомую функцию .

Рассмотрим подробнее линейную зависимость .

Дифференцируя (3), получим следующую систему уравнений

(5)

Из первого уравнения находим , где

, . (6)

Подставляя выражение для во второе уравнение, найдем

, (7)

где

, . (8)

Таким образом,

(9)

есть искомая линейная функция.

Ввиду простоты расчетов аппроксимация линейной зависимости используется довольно часто. Кроме того, многие функции, зависящие от двух параметров, можно линеаризовать путем замены переменных.

Для этого необходимо подобрать такое преобразование исходной зависимости , в результате которого она приобретает линейный вид . Далее решается задача линейной аппроксимации для новой зависимости ивычисленные коэффициенты ипересчитываются в коэффициентыи.

Для ряда часто встречающихся двухпараметрических зависимостей возможные замены переменных (а также, обратные замены для пересчета иви) приведены втабл. 1.

Таблица 1.

Вид зависимости	Замена переменных		Ограничения	Обратная замена переменных
Гиперболическая
Логарифмическая
Показательная
Степенная
Комбинированная

Более полное изложение этой темы – в [7], c.164-200.