Аппроксимация как метод исследования числовых характеристик

Любому специалисту в своей практической деятельности приходится изучать зависимости между различными параметрами исследуемых объектов, процессов и систем. Для вычисления значений сложных функций часто используется вычисление значения отрезка ряда, аппроксимирующего функцию.Аппроксимация также используется для обработки экспериментальных или натурных данных. На практике специалист чаще всего получает зависимости между исследуемыми параметрами экспериментально. В этом случае ставится натурный эксперимент, изменяются значения параметров на входе системы, измеряются значения параметров на выходе системы. Результаты измерений заносятся в таблицу. В результате проведения натурного эксперимента получается, так называемая, табличная функция [11].

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

Интерполяция – это замена исходной функции f(x) функцией φ(x) так, чтобы φ(x) точно проходила через точки исходной функции f(x). Интерполяция еще называется точечной аппроксимацией. Точки исходной функции f(x) называются узлами интерполяции.

Аппроксимация – это замена исходной функции f(x) функцией g(х)так, чтобы отклонение f(x) от g(х)в заданной области было наименьшим. Функция g(х)называется аппроксимирующей. Если исходная функция f(x) задана таблично (дискретным набором точек), то аппроксимация называется дискретной. Если исходная функция f(x) задана аналитически (на отрезке), то аппроксимация называется непрерывной или интегральной.

Аппроксимация или приближение – научный метод, позволяющий использовать одни объекты другими, в каком-то смысле близкими к исходным, но более простыми. Позволяет исследовать числовые характеристики и качественные свойства объекта, сводя задачу к изучению более простых или более удобных объектов (например, таких, характеристики которых легко вычисляются или свойства которых уже известны) [2].

Следует рассматривать два случая:

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

2) аппроксимирующая функция выступает в роли физического закона и с её помощью допускается экстраполировать переменные вне диапазона.

Для построения аппроксимирующей функции в Excel строится линия тренда с выводом коэффициента детерминации R²(величина достоверности аппроксимации). Линия тренда добавляется из контекстного меню кривой, построенной по точкам (функция f(x)) выбором пункта «Добавить линию тренда». В диалоговом окне «Линия тренда», вкладке «Тип» выбирается тип аппроксимирующей функции g(х), и устанавливаются флажки в полях «Показывать уравнение на диаграмме» и «Величина достоверности аппроксимации».

На графике появляется формула функции g(х) и значение достоверности аппроксимации R². Если значение R² близко к 1, то аппроксимирующая функция g(х) выбрана удачно, если не достаточно близко, то необходимо подобрать другую функцию.

Заметим также, что, если использовать в качестве линий тренда полиномы высоких степеней, то можно с точки зрения целевой задачи получить «абсурдный» результат прогноза, несмотря на высокое (близкое к 1) значение величины достоверности аппроксимации R². Иногда удается улучшить ситуацию заданием более точных коэффициентов (5-6 знаков после запятой) при установке для формулы на графике числового формата представления, но на внешних точках такая аппроксимация все равно выдаст недопустимые с точки зрения предметной области результаты.

Пример7.1. Спрогнозируйте заданные таблично статистические данные (рис. 36) на последующие периоды времени.

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

В диалоговом окне «Линия тренда», вкладке «Тип» выберем линейную аппроксимацию, и установим флажки в полях «Показывать уравнение на диаграмме» и «Величина достоверности аппроксимации». Уравнение линии тренда имеет вид у = 67,833*х – 134411. Введем эту формулу в ячейку С2, где в качестве х сошлемся на ячейку, содержащую значение года (А2) (рис. 38).

Добавим еще одну линию тренда, щелкнув опять правой кнопкой мыши по линии графика и выбрав пункт добавления линии тренда.В диалоговом окне на вкладке «Тип» выберем логарифмическую аппроксимацию и установим флажки в полях «Показывать уравнение на диаграмме» и «Величина достоверности аппроксимации». Уравнение второй линии тренда имеет вид у=135801*Ln(х) – 1Е+06. Установим для формулы на графике числовой формат представления для более точного определения свободного члена и введем формулу в ячейку D2.

Рис. 38. Экстраполяцияаппроксимирующей функции

Величина достоверности аппроксимации (R²) для обеих моделей составила 0,9745, она близка к 1, и это говорит о том, что и линейная, и логарифмическая аппроксимирующие функции достаточно точно описывают реальные данные.

Найдем отклонения для обеих аппроксимаций, как модуль разности между статистическими и теоретическими данными, т.е. для ячейки Е2 формула будет иметь вид: = ABS (C2-B2), для ячейки F2– = ABS (D2-B2). Также найдем максимум отклонения в обоих случаях.

Пример 7.2. Пусть мы имеем табличную функцию, которая для нескольких значений x определяет соответствующие значения f(x):

Таблица 25.

Требуется найти значение функции в точках x = 2,5 и x = 4,25.

Решение. Введем исходные данные на лист Excel. Дополнительно построим таблицу, где в качестве значений х будут числа от 0 до 6 с шагом 0,25. По данным первой таблицы построим точечную диаграмму и добавим на график линию тренда. В качестве аппроксимирующей кривой лучше всего выбрать полином 3-й степени, что определяется из внешнего вида графика. Величина достоверности аппроксимации (R²) при этом составляет 0,99, а уравнение имеет вид: y = 0,0885x³– 0,8256x² + 1,7254x– 0,0358 (рис. 39).

Рис. 39. Интерполяцияаппроксимирующей функции

Введем это уравнение в новую таблицу в ячейку Е2, где в качестве х укажем значение из ячейки D2 и скопируем формулу в остальные ячейки. Из таблицы находим значения функции в искомых точках: при x=2,5f(x)=0,5005, при x=4,25 f(x)= -0,8215.