Порядок выполнения работы.

1. Выберите вариант задания из таблицы случайных чисел.

2. Определите среднее арифметическое.

3. Найдите среднюю квадратическую ошибку отдельного измерения.

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

5. Определите среднюю квадратическую ошибку среднего арифметического.

6. Определите характеристики , h и H.

Таблица 1.2

№	Xi		h	H	s0	r0	S2	m3	A3	m4	B4

Контрольные вопросы.

1.Что называется случайной величиной?

2.Виды случайных величин?

3.Дискретная случайная величина? Функция распределения. Плотность распределения? Биномиальное распределение?

4.Распределение Пуассона? Нормальное распределение?

5.Вычисление средних величин?

ЛАБОРАТОРНАЯ РАБОТА №2

ПЛАНИРОВАНИЕ МНОГОФАКТОРНОГО ЭКСПЕРИМЕНТА

Цель работы.

Построение математической модели изучаемого объекта с использованием планирования эксперимента.

Основные положения.

Большое количество экспериментальных задач в химии и химической технологии формулируется как задачи экстремальные: определение оптимальных условий процесса, оптимального состава композиции и т.д. Планирование эксперимента резко повышает точность и уменьшает объем экспериментальных исследований. Поэтому использование метода планирования эксперимента является наиболее эффективным методом получения математических моделей многофакторного процесса. При его реализации можно оценить роль факторов, на которые можно воздействовать (температура, концентрации, давления и др.) при исследовании и оптимизации изучаемого объекта или технологического процесса, получить количественные оценки основных эффектов взаимодействия. Выбор плана определяется постановкой задачи исследования и особенностями объекта. Процесс исследования обычно разбивается на отдельные этапы. Информация, полученная после каждого этапа, определяет дальнейшую стратегию эксперимента. Таким образом, возникает возможность оптимального управления экспериментом. Применение методов планирования значительно повышает эффективность эксперимента [1].

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

Переменные x₁, x₂, …, x_k принято называть факторами. Факторами могут быть какие-либо внешние для объекта исследования воздействия (влажность, температура окружающей среды), или же параметры самого объекта (концентрация, температура, давление, удельная теплоемкость рабочей смеси). Выходные параметры также могут быть разнородными. В зависимости от решаемой задачи выходная величина называется откликом, функцией цели, функцией отклика, параметром оптимизации. Обычно аналитическая связь между входом и выходом (модель объекта) неизвестна, а известны факторы x_i и подлежащие исследованию выходные величины y_i.

Область определения двух факторов x₁, x₂ называется двухфакторным пространством, а эксперимент – двухфакторным экспериментом.

Каждый фактор может принимать определенное количество значений. Эти значения называются уровнями факторов. Например, если факторы могут принимать два значения -1, +1, то число уровней равно 2.

Число опытов в одном эксперименте равно числу различных наборов факторов. Для полного факторного эксперимента необходимое количество опытов N определяется по формуле:

,

где n – количество уровней; k – число факторов.

При проведении двухфакторного эксперимента на двух уровнях результат эксперимента представляет собой квадрат и число опытов равно:

.

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

.

Выбор факторов.

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

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

.

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

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

.

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

.

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

Модель процесса.

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

Модель определяется переменными x_i и постоянными параметрами β_i и в общем случае имеет вид:

.

Модели могут быть линейные относительно x_i:

,

Также они могут быть нелинейного вида:

Построение планов полного факторного эксперимента.

Полным факторным называется такой эксперимент, в котором реализуются все возможные комбинации (наборы) уровней факторов [3]. Если k факторов варьируются на двух уровнях, то число всех возможных наборов – N=2^k. Если k факторов варьируются на трех уровнях, то N=3^k. С увеличением числа факторов k, быстро растет число опытов.

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

Для трехфакторного эксперимента:

В каждой точке может проводиться несколько опытов n_u, которые называются параллельными. Для проведения эксперимента значения факторов разделяются на уровни, задаваемые соответствующими строками.

Основные достоинства планов ПФЭ – простота определения коэффициентов уравнения регрессии, возможность учета произведений взаимодействия факторов без изменения плана основного эксперимента. Преимущества любой матрицы ПФЭ достигаются за счет особого построения плана эксперимента, при котором матрица обладает свойствами ортогональности, нормировки, симметрии и ротатабельности [3].

Свойство ортогональности: сумма построчных произведений элементов любых двух граф равно нулю.

где i, j – номер столбца или номер фактора; i=1,2,…, k (k – общее количество факторов); u – номер набора факторов или номер строки; N – общее число различных наборов или число строк матрицы ПФЭ.

Свойство нормировки: сумма квадратов элементов любой графы равна числу различных опытов – строк N.

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

Свойство ротатабельности: дисперсии предсказанных значений отклика на равных расстояниях от центра плана постоянны и минимальны.

В полном факторном эксперименте возможен учет нелинейностей типа произведения факторов (уравнение регрессии отличается от линейного наличием слагаемого x_ix_j) и учет нелинейностей типа квадратов факторов.

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

В общем случае для k факторов полином второй степени имеет вид:

Полный квадратичный полином при k=2 содержит 6 членов и имеет вид:

при k=3 – 11 членов:

Для получения квадратичной зависимости каждый фактор должен фиксироваться как минимум на трех уровнях.

Планы ПФЭ не позволяют найти коэффициенты β_ii при квадратах факторов, так как все графы тождественны графе x₀ (x_i=±1, =1). Кроме того, для этих граф нарушаются условия ортогональности и симметричности. В основе построения планов второго порядка, как и планов первого порядка, лежат принципы ортогональности граф и их симметрия. Условие нормировки может не выполняться.

Для выполнения этих принципов к центру плана второго порядка, ядру – 2^k добавляются дополнительные точки факторного пространства. Точки называются звездными. Положение звездных точек определяется из условия ортогональности всех граф матрицы планирования второго порядка. Кроме значений факторов на уровнях ±1, в плане добавляется точка начала координат x_i=0 (i=1, 2,…,k), и на каждой координате выбираются две звездные точки x_i=α.

Общее число опытов определяется соотношением:

k – количество факторов.

На рисунке 2.1 показано проведение опыта в двухфакторном пространстве для плана ОЦКП. Звездные точки отмечены кружками. Величина α называется плечом звездных точекили звездным плечом. Положение звездных точек зависит от числа варьируемых факторов.

Рисунок 2.1 – расположение звездных точек.

Наличие звездных точек обеспечивает ортогональность граф первых степеней факторов. В ОЦКП каждый фактор фиксируется, в общем случае на пяти уровнях – -α, -1, 0, 1, α.

Для обеспечения ортогональности всех граф матрицы ОЦКП необходимо преобразовать квадраты факторов по формуле:

Значение поправки a определяется следующим образом. Сумма элементов столбца, соответствующего квадратам факторов, должна быть равна нулю:

Откуда

(2.1)

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

Таблица 2.1

Но-мер опы-та

x0

x1

x2

x3

x1∙x2

x1∙x3

x2∙x3

x1∙x2∙x3

y

Точки плана ПФЭ 23 (N0=2k точек)

+1

-1

-1

-1

+1

+1

+1

-1

1-a

1-a

1-a

y1

+1

+1

-1

-1

-1

-1

+1

+1

1-a

1-a

1-a

y2

+1

-1

+1

-1

-1

+1

-1

+1

1-a

1-a

1-a

y3

+1

+1

+1

-1

+1

-1

-1

-1

1-a

1-a

1-a

y4

+1

-1

-1

+1

+1

-1

-1

+1

1-a

1-a

1-a

y5

+1

+1

-1

+1

-1

+1

-1

-1

1-a

1-a

1-a

y6

+1

-1

+1

+1

-1

-1

+1

-1

1-a

1-a

1-a

y7

+1

+1

+1

+1

+1

+1

+1

+1

1-a

1-a

1-a

y8

Звездные точки (2k точек)

+1

-α

α2-a

-a

-a

y9

+1

+ α

α2-a

-a

-a

y10

+1

- α

-a

α2-a

-a

y11

+1

+ α

-a

α2-a

-a

y12

+1

- α

-a

-a

α2-a

y13

+1

+ α

-a

-a

α2-a

y14

Нулевая точка

+1

-a

-a

-a

y15

-

N

-

-

N

Для определения неизвестных a и α нужно сформировать и решить систему из двух уравнений. Одно из них для a было записано раннее – это уравнение (2.1). Другое уравнение получим из условия ортогональности для столбцов x₄ и x₅:

.

После простейших преобразований с учетом того, что - общее число опытов в плане, получаем:

.

Соотношение для a при j=1, 2 или 3 может быть записано как (см. таблицу 2.1):

.

Подставив его в предыдущее уравнение, получаем:

Откуда

.

Тогда

.

И плечо звездных точек

Например, для ОЦКП при числе факторов k=3 параметры плана следующие:

Обобщим результаты в таблицу 2.2.

Таблица 2.2

Число факторов k	Ядро (центр) плана	Общее число опытов N	Звездное плечо α	Значение поправки a
	22		1.0000	2/3=0.6667
	23		1.2154	8/15=0.7303
	24		1.4142	4/5=0.80
	25-1		1.5467	=0.7698
	25		1.5960	=0.8627

Располагая данными таблицы 2.2, можно построить ортогональные центрально-композиционные планы второго порядка.

По результатам опытов формируется полином.

при k=2:

при k=3:

Коэффициенты полинома определяются по формуле:

Полиномы можно представить в виде:

- при k=2;

и при k=3:

Где коэффициент определяется следующим образом при k=2:

при k=3:

Ортогональные центрально композиционные планы при не обеспечивают постоянства дисперсии отклика в различных точках факторного пространства.