Построение математических моделей с использованием только аналитических методов, как правило, трудно реализуемо на практике из-за недостатка необходимой информации о моделируемом процессе или объекте. Поэтому, помимо аналитического подхода используют и экспериментальные методы, позволяющие определить или получить необходимые для построения модели данные либо собственно математическое описание объектов в упрощенном (обобщенном) виде.
К недостаткам экспериментальных методов можно отнести следующие моменты:
1) локальность математической модели, получаемой по эксперименту;
2) большие временные затраты на постановку эксперимента;
3) возможные материальные и финансовые затраты на постановку и проведение эксперимента.
Все экспериментальные методы имеют три основных составляющих, которые необходимо пройти для построения математической модели.
1) Собственно планирование проводимых необходимых экспериментов.
2) Проведение экспериментов и фиксирование получаемых результатов.
3) Обработка полученных данных.
Планирование экспериментов основано на представлении объекта как «черного ящика» (рис. 6.1), т.е. в процессе исследования изучается влияние входной переменной (фактора) x на выходную переменную (реакцию, отклик) y.
Для того, чтобы получить статическую характеристику объекта с k входами: Х(х1,...,хk) и одним выходом yв некоторой области Dx, необходимо задавать поочередно на входе объекта значения входным переменным
Х(1),..., Х(N)Î Dx, выдерживая каждое из них в течение времени t≥T, где Т - оценка времени переходного процесса по всем каналам.
Измеряя на выходе соответствующие установившиеся значения y, получим таблицу значений выходной переменной y1,..., yN. Значения Х(i) и соответствующие им yi задают в табличном виде искомую статическую характеристику.
6.2.1 Основная идея регрессионного анализа
Для исследования взаимосвязи между величинами Х (вход) и y (выход) используются методы корреляционного и регрессионного анализа. Результаты корреляционного анализа позволяют сделать вывод о степени зависимости между переменными, а форма зависимости уточняется методами регрессионного анализа.
На значение величины y оказывают влияние стохастические воздействия разного рода, поэтому форма связи между величинами Х и y определяется линией регрессии, показывающей, как в среднем изменяется величина yпри изменении входной величины Х, т.е. приходится говорить о связи средних значений величиныy c X. Эту связь характеризуют условным математическим ожиданием величины y, вычисляемым при условии, что величина Хприняла определенное значение, а аппроксимирующая функция строитсякак функция регрессии f(X, B) » M[y/X], где B - неизвестные параметры уравнения регрессии.
Задача регрессионного анализа ставится следующим образом. Для каждого i-го опыта имеется набор значений (xi1,...,xik) входных параметров X1 ¸Xk и соответствующее им значение выходного параметра yi. Пример опытных данных приведен в табл. 6.1.
Таблица 6.1 – Пример опытных данных.
№
опыта
Входы
Выходы
Y
X1
X2
---
Xk
x11
x12
---
x1k
y1
x21
x22
---
x2k
y2
-----
-----
-----
---
-----
-----
n
xn1
xn2
---
xnk
yn
Необходимо определить зависимость выхода y от входных факторовxi1,...,xik, которая для случая линейной связибудет иметь следующий вид:
y = b0 + b1x1 + b2x2+ ... + bkxk(1)
Задача сводится к определению оценок коэффициентов уравнения регрессии b0, b1, ... , bk, которые с определенной степенью вероятности будут отражать влияние аргументов xi1,..., xik на y.
Для определения bkиспользуется метод наименьших квадратов.
Таким образом, задача сводится к минимизации:
(2)
где yi - фактическое (экспериментальное) значение выходной переменной,
- рассчитанное значение выходной переменной.
Можно выражение (2) представить с учетом (1):
(3)
То есть, процесс сводится к нахождению таких значений коэффициентов b0, b1, ... , bk, которые давали бы минимальное расхождение между расчетными и измеренными значениями y.
Это условие может быть выполнено, если приравнять к нулю частные производные S по каждому из коэффициентов b0, b1, ... , bk и просуммировать их по индексу i. В результате получается система обычных уравнений (4):
(4)
Решая эту систему, получают коэффициенты модели (1) и, как следствие, явный вид уравнения регрессии для случая линейной связи.
xk
yk
8,8362
13,7195
18,8031
23,7117
28,7673
33,5749
38,5742
43,5962
48,6209
53,6476
Таблица 6.2
Проведено 10 экспериментов, результаты которых приведены в табл. 6.2 С учетом введенных ранее обозначений n=10, k=1. Для случая одного входного параметра искомое уравнение регрессии будет иметь вид y = b0 + b1x1
А для определения коэффициентов регрессии необходимо сформировать систему уравнений (5) из общей системы (4).
(5)
Исходные данные внесем в таблицу Excel и определим расчетные цифры для системы (5): суммы xk, yk, xk2, xk*yk. Подстановка результатов расчета в систему (5) дает уточненную систему уравнений:
10*b0 + 5*b1 = 311,85
55*b0 + 385*b1 =2125,67,
решая которую получаем коэффициенты регрессии b0 = 3,93 b1 =4,962.
Стандартная функция Excel ЛИНЕЙН для определения коэффициентов регрессии дает решение b0 = 3,819 b1 =4,975.