3.1Теоретические сведения
Анализ влияния факторов на исходный параметр
При экспериментальных исследованиях часто изучаются объекты, которые в общем случае могут быть представлены в виде «черного ящика», на вход (входы) которого воздействуют величины, называемые факторами, независимыми переменными, или регрессорами. Целью экспериментальных исследований является получение зависимости между входными величинами и выходной величиной, называемой функцией отклика. Функция отклика является, в общем случае функцией многих переменных и о ней имеются самые общие представления (иногда интуитивные). Конечной целью экспериментального исследования является математическая модель, адекватно описывающая поведение объекта.
Входные переменные, или факторы, могут быть разбиты на два класса: контролируемые (измеряемые) и неконтролируемые (не измеряемые). В свою очередь контролируемые переменные могут быть управляемыми и неуправляемыми. Управляемыми являются такие факторы, которые целенаправленно можно изменять в ходе эксперимента. Факторы, для которых изменение не возможно, называются неуправляемыми.
Если все факторы трактуются как качественные, то применяется дисперсионный анализ, если один из факторов качественный, а другие количественные – то корреляционный анализ.
Другой широко распространенной задачей обработки данных является представление их совокупности некоторой функцией у(х). Задача регрессии заключается в получении параметров этой функции такими, чтобы функция приближала бы облако исходных точек (заданных векторами VX и VY) с наименьшей среднеквадратичной погрешностью.
Статистическая обработка данных с помощью системы MathCad
С помощью системы MК можно проводить наиболее распространенные статистические расчеты с данными, представленными векторами их значений. В MК существуют также встроенные статистические функции скалярного аргумента х:
rnd(x) – функция генерации случайных чисел с равномерным распределением;
corr(VX, VY) – коэффициент корреляции двух векторов – VX и VY;
cvar(X, Y) – коэффициент ковариации X и Y. Через функцию erf(x) легко вычисляется дополнительная функция ошибок crfc(x) = 1 - erf(x).
Следующая группа функций относится к вычислению основных статистических параметров одномерного массива данных вектора:
mean(V) – среднее значение;
var(V) – дисперсия (вариация);
side(V) – среднеквадратичная погрешность (квадратный корень из дисперсии);
stdav(V) – стандартное отклонение;
hiet(int, V) – вектор частот попадания данных V в заданные интервалы;
int – служит для построения гистограмм.
На рисунке 3.1 представлен документ, в котором организована генерация вектора X, который содержит 200 случайных чисел, представлено их распределение и вычислены основные статистические параметры массива случайных чисел – вектора X.
При достаточно большом количестве случайных чисел вид гистограммы приближенно говорит о законе их распределения.
Функции распределения плотности вероятности в МК представлены следующим набором:
dbeta(x, s1, s2) – β-распределение (s1, s2 > 0 – параметры формы, 0 < х < 1);
dbino(k, n, р) – биномиальное распределение (возвращает значение вероятности Р(х = k), где n и k целые числа, причем 0 < k < n и 0 < р < 1);
dcauchy(x, l, s) – распределения Коши (l – параметр разложения, s > 0 – параметр масштаба);
dchisq(x, d) – Хи-квадрат-распределение (х, d > 0, причем d – число степеней свободы);
dexp(х, r) – экспоненциальное распределение (r, х > 0);
Рисунок 3.1 – Работа со случайными числами
df(х, d1, d2) – распределение Фишера (d1, d2 > 0 – числа степеней свободы, х>0);
dgama(x, s) – гамма-распределение (s > 0 – параметр формы, х > 0);
dgecm(k, р) – геометрическое распределение (0 < р < 1 – вероятность успеха в отдельном испытании, к – целое неотрицательное число);
dlnor(x, µ, σ) – логнормальное распределение (µ – натуральный логарифм среднего значения, σ > 0 – натуральный логарифм среднеквадратичного отклонения, х > 0);
dlogis(x, l, s) – логистическое распределение (l – параметр разложения, s > 0 – параметр масштаба);
dfibino(k, n, р) – отрицательное биномиальное распределение (n > 0 и k > 0 – целые числа, 0 < р < 1);
dnorm(x, p, a) – нормальное распределение (р – среднее значение, а > 0 – среднеквадратичное отклонение);
dpois(k, x) – распределение Пуассона (к – целое неотрицательное число);
dt(x, d) – распределение Стьюдента (d > 0 – число степеней свободы, х – вещественное число);
duftif(x, а, b) – равномерное распределение (а и b – граничные точки интервала, причем a < b и а < х < b);
dweibuill(x, s) – распределение Вейбулла (s > 0 – параметр формы).
Функции распределения дают вероятность того, что случайная величина будет иметь значения, меньшие или равные определенной величине. Они представлены аналогичным набором функций и отличаются от приведенных тем, что первой буквой имени является буква р, а не d (например, pbeta(x, S1, S2) и т.д.). При экстраполяции используются отрезки прямых, проведенных через две крайние точки.
Следующая группа задает обращения (квантили) функций распределения случайных величин. Имена этих функций отличаются от функций распределения лишь первой буквой. Они начинаются с буквы q – qbeta(р, S1, S2) и т.д. и позволяют по заданной вероятности вычислить такое значение х, при котором вероятность равна или меньше заданного значения р.
Последняя группа статистических функций служит для создания векторов с определенными законами распределения значений их элементов. Их имена так же отличаются от предыдущих двух групп первой буквой. Они начинаются с буквы r. Например, rbeta(n, S1, S2).