При выборе формы представления результатов эксперимента в детерминированных системах в первую очередь рассматривается степенная зависимость, как обладающую наибольшей простотой и возможностью варьирования значений коэффициентов и показателей, характеризующих объект. Если степенная форма дает неудовлетворительные результаты, то используется подбор типа эмпирической формулы y=f(x). В стохастических системах для поиска зависимости y=f(x) используется математический аппарат корреляционно-регрессионного анализа.
2. Полученные экспериментальные данные следует представить графически в логарифмической зависимости. Если меж-
где а=Аsinά, в=Аcosά.
6. В тех случаях, когда целью моделирования является определение случайной функциональной зависимости между параметрами объекта, её получают в результате специальной статистической обработки экспериментальных данных.
Допустим, по результатам эксперимента на модели необходимо найти зависимость параметра y от параметра x в виде y = f(x).
|
|
Для этого обработка полученных статистических данных ведётся в следующей последовательности.
1. Определяются статистические средние (оценки математических ожиданий) y и x по формулам:
, ,
где xi и yi – фактические результаты, полученные в ходе эксперимента.
2. Устанавливается функциональная зависимость между x и y.
Для дальнейшего упрощения поиска зависимости на основе корреляционно-регрессионного анализа она представляется в линейном виде
y = a0 + a1x.
Если искомая зависимость нелинейная, то её путём замены переменных необходимо свести к линейному виду.
Например, y = a0 + a1/x, заменить x’ = 1/x => y = a0 + a1x’,
y = a0*e a1x, заменить y = ey1 => y1 = ln a0 + a1x.
3. В полученной линейной зависимости определяется на сколько x и y линейно зависимы друг от друга. Оценку линейной зависимости производят с помощью коэффициента корреляции R
R = , где kx, y – момент корреляции между x и y;
- средние квадратичные отклонения по x и y.
Коэффициент R в линейной зависимости изменяется от -1 до +1. Если R отрицательный, зависимость обратно пропорциональная; если R = 0, зависимость отсутствует.
Оценка коэффициента R на основе статистических данных определяется по формуле
.
Если по абсолютному значению окажется меньше 0,3, то зависимость между x и y слабая (). Если окажется, что , зависимость – сильная. В промежутке средняя зависимость.
Допустим, зависимость между x и y оказалась средняя и выше. В этом случае переходят к следующему пункту.
4. Определяется количественное значение коэффициентов а0 и а1 по методу наименьших квадратов. Для этого определяется сумма квадратов отклонений yi экспериментального от теоретического значения (yi = a0 + a1xi) по формуле
|
|
.
Полученное значение Z является функцией двух коэффициентов а0 и а1. Их необходимо подобрать таким образом, чтобы истинное значение yi мало отличалось от модельного значения yi.= а0+ а1хi . А это значит, найти такие значения а0, а1, при которых Z принимала бы минимальное значение. Минимум функции Z можно найти из системы уравнений с частными производными
; .
Поиск а0 и а1 осуществляется следующим образом. Дифференциальные уравнения записываются на основе полученных статистических данных в виде:
.
Из системы коэффициенты определяются по формулам:
; .
.
5. После того как значения коэффициентов определены, устанавливается их значимость в функциональной зависимости, т.е. определяется насколько их значение влияет на зависимость y=f(x). Значимость коэффициентов находится путём сравнения случайных величин tф с tтеор., где tтеор – теоретическая случайная величина, подчиняющаяся t-распределению и задана таблично для выбранной степени свободы V и уровня доверия ; фактическая величина tф для каждого из коэффициентов определяется по формуле
taj = , j – индекс коэффициента (j = 0,1).
Если окажется, что tф tтеор., то соответствующий коэффициент незначителен и им можно пренебречь.
6. С учетом п. 5 окончательно записывается функциональная зависимость y = f(x).
Пример обработки результатов эксперимента. Допустим, что по результатам планирования эксперимента имеются статистические данные по факторам х1, х2 и отклику y, указанные в табл. 1.
Таблица 1
Номер опыта | Значение факторов | Значение отклика | |
х1 | х2 | y | |
Необходимо установить существует ли линейная зависимость вида .
Для этого выполним следующую статистическую обработку данных.
1. Вычислим коэффициенты корреляции между х1, y и х2, y по формуле
.
Полученные коэффициенты , указывают, что между х2 и y существует сильная прямопропорциональная зависимость, а между х1 и y – зависимость отсутствует.
2. На основании полученных результатов регрессию запишем в виде . Коэффициенты , вычислим по формулам:
, ;
, .
3. Установим значимость коэффициентов b0 и b2.
Для этого по формуле
tbj =
определим, что фактическое значение , .
По таблице t- распределения для и найдем .
4. Так как коэффициенты b0, b2 значимы, окончательно искомая зависимость примет вид .