Теория планирования эксперимента

Реализация задачи оптимизации на основе применения теории планирования эксперимента. Полный факторный эксперимент типа 2n. Дробный факторный эксперимент. Матрица планирования и оценка коэффициентов функции отклика в дробном факторном эксперименте.

Рубрика Экономико-математическое моделирование
Вид контрольная работа
Язык русский
Дата добавления 05.04.2015
Размер файла 69,3 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

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

Размещено на http://www.allbest.ru

Содержание

Задание

Введение

1. Постановка задачи оптимизации

2. Полный факторный эксперимент типа 2n

3. Дробный факторный эксперимент

4. Оценки коэффициентов функции отклика

5. Оценки коэффициентов функции отклика в дробном факторном эксперименте

Расчётная часть

Заключение

Список литературы

Задание

 

 

 

 

 

 

 

 

25

25-1

25-2

 

x0

x1

x2

x3

x4

x5

y

y'

y''

y'''

1

1

-1

-1

-1

-1

-1

36

 

 

 

2

1

1

-1

-1

-1

-1

37

 

 

 

3

1

-1

1

-1

-1

-1

38

 

 

 

4

1

1

1

-1

-1

-1

39

 

 

 

5

1

-1

-1

1

-1

-1

40

 

 

 

6

1

1

-1

1

-1

-1

41

 

 

 

7

1

-1

1

1

-1

-1

42

 

 

 

8

1

1

1

1

-1

-1

43

 

 

 

9

1

-1

-1

-1

1

-1

44

 

 

 

10

1

1

-1

-1

1

-1

43

 

 

 

11

1

-1

1

-1

1

-1

42

 

 

 

12

1

1

1

-1

1

-1

41

 

 

 

13

1

-1

-1

1

1

-1

40

 

 

 

14

1

1

-1

1

1

-1

39

 

 

 

15

1

-1

1

1

1

-1

38

 

 

 

16

1

1

1

1

1

-1

37

 

 

 

17

1

-1

-1

-1

-1

1

38

 

 

 

18

1

1

-1

-1

-1

1

39

 

 

 

19

1

-1

1

-1

-1

1

40

 

 

 

20

1

1

1

-1

-1

1

41

 

 

 

21

1

-1

-1

1

-1

1

42

 

 

 

22

1

1

-1

1

-1

1

43

 

 

 

23

1

-1

1

1

-1

1

44

 

 

 

24

1

1

1

1

-1

1

45

 

 

 

25

1

-1

-1

-1

1

1

43

 

 

 

26

1

1

-1

-1

1

1

41

 

 

 

27

1

-1

1

-1

1

1

39

 

 

 

28

1

1

1

-1

1

1

37

 

 

 

29

1

-1

-1

1

1

1

35

 

 

 

30

1

1

-1

1

1

1

33

 

 

 

31

1

-1

1

1

1

1

31

 

 

 

32

1

1

1

1

1

1

29

 

 

 

Введение

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

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

Планирование экспериментов (ПЭ) охватывает широкий круг вопросов - от учета конкретных особенностей определенных объектов исследования до общих концептуальных проблем. Далее будут рассмотрены общие задачи планирования, изучаемые специальной научной дисциплиной - теорией планирования эксперимента (ТПЭ).

1. Постановка задачи оптимизации

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

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

· выделить управляемые и неуправляемые параметры (факторы) системы и среды, оказывающие существенное влияние на критерий эффективности;

· определить ограничения на значения параметров.

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

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

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

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

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

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

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

Одновременное изменение факторов предполагает их совместимость, что означает осуществимость и безопасность всех их сочетаний. Необходимо также обеспечить независимость изменения каждого фактора, что означает возможность установления любого значения фактора вне связи со значениями других факторов.

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

2. Полный факторный эксперимент типа 2n

На начальных этапах оптимизации для определения градиента применяют неполные полиномы второго порядка или линейные полиномы. Вычисление оценок коэффициентов таких полиномов осуществляется на основе обработки результатов реализации наиболее простых планов, в которых каждый фактор принимает только два значения vi, min или vi, max, расположенные симметрично относительно некоторого нулевого уровня или центра плана по данному фактору. Значения уровней варьирования выбирает исследователь, исходя из возможного диапазона изменения каждого фактора и возможности применения линейной аппроксимации функции отклика в выбранном диапазоне изменений параметра. Без ограничения общности можно считать, что кодированные значения xi принимают значения - 1 и +1 соответственно (или просто - или +). Множество всех точек в k-мерном пространстве, координаты которых являются комбинациями "+" и "-", называется полным факторным планом или планом полного факторного эксперимента типа 2k (ПФЭ). Количество точек в этом плане N =2n.

Для примера возьмем полный факторный эксперимент с тремя независимыми переменными x1, х2 и x3, табл. 1.

Таблица 1

Матрица планирования

Вектор результатов

x0

x1

x2

x3

x1 x2

x1 x3

x2 x3

x1 x2 x3

y

+

-

-

-

+

+

+

-

y1

+

-

-

+

+

-

-

+

y2

+

-

+

-

-

+

-

+

y3

+

-

+

+

-

-

+

-

y4

+

+

-

-

-

-

+

+

y5

+

+

-

+

-

+

-

-

y6

+

+

+

-

+

-

-

-

y7

+

+

+

+

+

+

+

+

y8

Второй, третий и четвертый столбцы таблицы соответствуют собственно плану экспериментов, пятый - восьмой столбцы содержит значения произведений независимых переменных. Фиктивная переменная x0 =1 (первый столбец) введена для единообразия записи расчетных формул коэффициентов полинома. Строки соответствуют опытам, например, первая строка характеризует эксперимент, в котором все независимые переменные находятся на нижнем уровне.

Существует несколько способов построения подобных матриц планирования. В частности можно воспользоваться приемом, характерным для записи последовательности двоичных чисел. В столбце последней переменной x3 знаки меняются поочередно, в столбце предпоследней переменной x2 - чередуются через два элемента, третьей справа переменной x1 - через четыре элемента. Аналогично строится матрица для любого количества переменных, порядок перечисления переменных не играет роли. Столбцы с произведениями переменных вычисляются путем умножения значений элементов в соответствующих столбцах простых переменных.

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

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

симметричности. Сумма всех элементов любого столбца, за исключением первого, равна нулю, например

нормированности. Сумма квадратов элементов любого столбца равна числу опытов, так для i-й переменной

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

Включение в матрицу планирования переменных вида xi2 приведет к появлению единичных столбцов, совпадающих друг с другом и со столбцом x0. Следовательно, нельзя будет определить, за счет чего получено значение x0. Поэтому планы ПФЭ 2n не применимы для построения функции отклика в виде полного полинома второй степени.

3. Дробный факторный эксперимент

С ростом количества факторов n число точек плана в ПФЭ растет по показательной функции 2n. Планы ПФЭ позволяют получить несмещенные оценки градиента функции отклика в центральной точке, но в случае применения линейного полинома оказываются недостаточно эффективными по количеству опытов при большом числе независимых переменных, так как остается слишком много степеней свободы на проверку адекватности модели. Например, при n = 5 на проверку адекватности линейной модели остается 26 степеней. Хотя большое количество опытов и приводит к существенному снижению погрешности в оценке коэффициентов, все же такое число степеней свободы для проверки адекватности является чрезмерным.

Таким образом, в случаях, когда используются только линейные приближения функции отклика, количество опытов следует сократить, используя для планирования так называемые регулярные дробные реплики от ПФЭ, содержащие подходящее число опытов и сохраняющие основные свойства матрицы планирования. Реплика, включающая только половину экспериментов ПФЭ, называется полурепликой, включающая четвертую часть опытов - четвертьрепликой и т. д. Краткое обозначение  указанных дробных реплик 2n - 1, 2n-2 соответственно.

Построение регулярной дробной реплики или проведение дробного  факторного эксперимента (ДФЭ) типа 2n-p предусматривает отбор из множества n факторов n-p основных, для которых строится план ПФЭ. Этот план дополняется р столбцами, которые соответствуют остальным факторам. Каждый из этих столбцов формируется по специальному правилу, а именно, получается как результат поэлементного умножения не менее двух и не более n-p определенных столбцов, соответствующих основным факторам. Иначе говоря, в  дробных репликах p линейных эффектов приравнены к эффектам взаимодействия. Но именно такое построение матрицы планирования и позволяет обеспечить ее симметричность, ортогональность и нормированность.

Таблица 2

Матрица планирования

Вектор результатов

x0

x1

x2

x3

y

+

-

-

+

y1

+

-

+

-

y2

+

+

-

-

y3

+

+

+

+

y4

Правило образования каждого из p столбцов ДФП называют генератором плана. Каждому дополнительному столбцу соответствует свой генератор (для плана типа 2n- p должно быть задано p различных генераторов). Генератор задается как произведение основных факторов, определяющее значение элементов соответствующего дополнительного столбца матрицы планирования. Примером записи генератора для плана 23-1 служит выражение x3 = x1x2, табл. 2. Матрица планирования ДФП типа 2n- p содержит n + 1 столбец и N = 2n- p строк.

4. Оценки коэффициентов функции отклика

С помощью матрицы планирования, описанной в табл. 1, можно вычислить оценки коэффициентов неполного полинома третьей степени

y' = ?0 + ?1x1 + ?2x2 +?х3 +?12x1x2 +?13x1x3 + ?23x2x3 + ?123x1x2х3

или линейной функции

y' = ?0 + ?1x1 + ?2x2 +?3х3.

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

Оценки коэффициентов полинома определяются на основе метода наименьших квадратов и для рассматриваемого типа ПФЭ вычисляются по простым соотношениям.

;

.

(1)

Здесь величина y  соответствует значению отклика  в указанной точке факторного пространства при отсутствии повторных опытов или является оценкой математического ожидания

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

Допустима следующая интерпретация оценок коэффициентов:

?0 соответствует значению функции отклика в центре проводимого эксперимента;

?i равен приращению функции при переходе значения фактора i с нулевого уровня на верхний (это вклад соответствующего фактора в значение функции);

?ij равен нелинейной части приращения функции при одновременном переходе факторов i и j с нулевого уровня на верхний и т.п.

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

.

(2)

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

,

где ?u2 - оценка дисперсии функции отклика в точке u, ru - число повторных опытов в этой точке плана [7, стр. 50]. Дисперсия оценок всех коэффициентов одинакова, поэтому ПФЭ рассмотренного типа являются ротатабельным.

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

5. Оценки коэффициентов функции отклика в дробном факторном эксперименте

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

Запись плана в виде 2n-p не дает полной характеристики регулярной дробной реплики, так как основные эффекты можно приравнять к различным эффектам взаимодействия. Правило смешивания, определяющее коррелированные основные эффекты и эффекты взаимодействия, удобно описывать с помощью определяющего контраста реплики. Определяющий контраст полуреплики получается путем умножения генерирующего соотношения на его же левую часть, а так как для любой кодированной переменной xi2=1, то левая часть формулы определяющего контраста всегда равна единице и обозначается I. В частности, для ДФП типа 23-1 и генераторе x3 = x1x2 имеет место определяющий контраст I = x1x2x3 (генератор умножается на переменную x3, следовательно, x3 x3 = I =  x1 x2 x3).

Чтобы определить, с какими параметрами смешана оценка коэффициента данного фактора, следует умножить обе части определяющего контраста на этот фактор. Учитывая равенство xi2=1, получим порядок смешивания оценок коэффициентов при использовании конкретного плана. В рассматриваемом примере для плана 23-1 и определяющего контраста I = x1x2x3 порядок смешивания факторов следующий:

x1 = x12 x2 x3 = x2 x3x2 = x1 x22 x3 = x1 x3x3 = x1 x2 x32 = x1 x2 .

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

b1*= b1 + b23 ;  b2*= b2 + b13 ;  b3*= b3 + b12 .

Планы типа 2n-р являются ортогональными для моделей с взаимодействиями. Поэтому для вычисления оценок коэффициентов получаются простые формулы, как и для случая ПФЭ

.

Планы дробных реплик строят различным образом, но так, чтобы соблюдались основные свойства матрицы планирования. Например, ДФП 23-1 можно представить одной из двух полуреплик, генераторами которых являются x3 = x1x2 и x3 = - x1x2 соответственно. Определяющие контрасты этих полуреплик:

x32  = I = x1x2x3 и x32  = I = - x1x2x3 .

В этих полурепликах смешивание факторов задается соотношениями:

а) x1 = x2x3 x2 = x1x3 x3 = x1x2 ;

б) x1 = - x2x3 x2 = - x1x3 x3 = - x1x2 .

Коэффициенты линейного полинома в каждой полуреплике:

а) b1* = b1 + b23 ;  b2* = b2 + b13 ;  b3* = b3 + b23 ;

б) b1* = b1 - b23 ;  b2* = b2 - b13 ;  b3* = b3 - b23 .

Реализовав обе полуреплики путем совместной обработки результатов экспериментов можно получить раздельные оценки для линейных эффектов и эффектов взаимодействия (такой вариант плана соответствует ПФЭ).

Разрешающая способность полуреплик (возможность раздельного определения коэффициентов уравнения) зависит от генерирующих соотношений. Так, если для плана 24-1 выбрать генерирующее соотношение x4 = x1x2, то получим реплику с контрастом I =  x1x2x4 и разрешающей способностью x1 = x2x4 и т.д. Здесь линейные эффекты определяются совместно с парными взаимодействиями. Очевидно, что в первую очередь следует пренебречь взаимодействием более высоких порядков из-за их более низкой вероятности существования по сравнению с парными. У полуреплики с контрастом I = x1x2х3x4 или равноценным I = - x1x2х3x4 линейные эффекты будут определяться совместно уже только с тройными взаимодействиями, что повышает точность оценок параметров модели (потенциально величина смещения в оценке коэффициента уменьшается). С ростом количества независимых переменных растет разрешающая способность полуреплик, позволяя оценивать раздельно сначала линейные эффекты, затем парные, тройные  взаимодействия и т. д. Но при этом растет и избыточность экспериментов.

Реплики можно строить высокой степени дробности, сокращая тем самым количество экспериментов. Пусть необходимо изучить влияние пяти переменных и известно, что все эффекты взаимодействия пренебрежимо малы. Для линейного приближения следует определить шесть коэффициентов, что потребует применения плана с количеством точек не менее шести. Ближайшее большее число, соответствующее целой степени 2, равно восьми, это дает возможность получить дробную реплику, эквивалентную ПФЭ 23, т. е. реплику 25 - 2 или четвертьреплику. Для построения четвертьреплики необходимы два генерирующих соотношения. В целях построения такой реплики целесообразно пожертвовать тройным и одним из двойных взаимодействий. Пусть этим двойным взаимодействием будет x1x2. Тогда можно построить четыре различные четвертьреплики, каждая из которых задается двумя генерирующими соотношениями:

а) x4 = x1x2 x5 = х1x2x3 ;

б) x4 = x1x2 x5 = - х1x2x3 ;

в) x4 = - x1x2 x5 = х1x2x3 ;

г) x4 = - x1x2 x5 = - х1x2x3 .

Определяющие контрасты каждой четвертьреплики задаются двумя соотношениями:

а) I = х1x2x4 , I = х1x2x3x5 ;

б) I = х1x2x4 , I = - х1x2x3x5 ;

в) I = - х1x2x4 , I = х1x2x3x5 ;

г) I = - х1x2x4 , I = - х1x2x3x5 .

Из этой совокупности четвертьреплик следует выбрать только одну, например, выберем реплику, задаваемую первой парой генерирующих соотношений. Матрица планирования ДФП получается из матрицы ПФЭ 2n-p для n-p основных факторов добавлением р столбцов, элементы которых вычисляются по соответствующим генерирующим соотношениям, табл. 3

Таблица 3

Матрица планирования

Вектор

 результатов

х0

х1

х2

х3

х4

х5

+

-

-

-

+

-

y1

+

+

-

-

-

+

y2

+

-

+

-

-

+

y3

+

+

+

-

+

-

y4

+

-

-

+

+

+

y5

+

+

-

+

-

-

y6

+

-

+

+

-

-

y7

+

+

+

+

+

+

y8

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

I = х1x2x4 = х1x2x3x5 = x3x4x5 .

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

x1 = x2x4 = x2x3х5 = x1x3x4х5,

x2 = x1x4 = x1x3х5 = x2x3x4х5,

    .    .    .    .    .    .    .

x5 = х1x2x4х5= x1x2х3 = x3x4 .

Оценки коэффициентов линейного полинома задаются соотношениями:

b1* = b1 + b24 + b235 + b1345 ,

b2* = b2 + b14 + b135 + b2345 , и  т. д.

Разрешающая способность выбранной четвертьреплики невысокая - все линейные эффекты определяются совместно с парными взаимодействиями. Этой репликой можно пользоваться для оценки линейных эффектов при условии равенства нулю соответствующих парных взаимодействий. Если такой уверенности нет, то следует применить полуреплику (что требует в два раза большего количества точек плана эксперимента по сравнению с четвертьрепликой) с генерирующим соотношением x5 = х1x2x3x4, пользуясь которым, можно разделить все линейные эффекты и парные взаимодействия.

Построение обобщающего определяющего контраста для реплик более высокой степени дробности производится аналогично четвертьреплике: исходные контрасты сначала перемножаются попарно, получаются контрасты первого уровня; затем контрасты первого уровня снова перемножаются попарно, получаются контрасты второго уровня и так далее, пока не будет исчерпана возможность перемножения. Если получается два и более одинаковых контрастов, то из них оставляется только один. Обобщающий определяющий контраст составляется путем перечисления выражений для всех сформированных контрастов.

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

Расчётная часть

 

x0

x1

x2

x3

x4

x5

x1x2

x1x3

x1x4

x1x5

x2x3

x2x4

x2x5

x3x4

x3x5

x4x5

y

y'

1

1

-1

-1

-1

-1

-1

1

1

1

1

1

1

1

1

1

1

36

49,13

2

1

1

-1

-1

-1

-1

-1

-1

-1

1

1

1

-1

1

-1

-1

37

31,88

3

1

-1

1

-1

-1

-1

-1

1

1

-1

-1

-1

1

1

-1

-1

38

33,13

4

1

1

1

-1

-1

-1

1

-1

-1

1

-1

-1

1

1

-1

-1

39

38,88

5

1

-1

-1

1

-1

-1

1

-1

1

-1

-1

1

-1

-1

1

-1

40

35,63

6

1

1

-1

1

-1

-1

-1

1

-1

1

-1

1

-1

-1

1

-1

41

41,38

7

1

-1

1

1

-1

-1

-1

-1

1

-1

1

-1

1

-1

1

-1

42

42,63

8

1

1

1

1

-1

-1

1

1

-1

1

1

-1

1

-1

1

-1

43

48,38

9

1

-1

-1

-1

1

-1

1

1

-1

-1

1

-1

-1

-1

-1

1

44

40,38

10

1

1

-1

-1

1

-1

-1

-1

1

1

1

-1

-1

-1

-1

1

43

43,63

11

1

-1

1

-1

1

-1

-1

1

-1

-1

-1

1

1

-1

-1

1

42

42,38

12

1

1

1

-1

1

-1

1

-1

1

1

-1

1

1

-1

-1

1

41

45,63

13

1

-1

-1

1

1

-1

1

-1

-1

-1

-1

-1

-1

1

1

1

40

39,88

14

1

1

-1

1

1

-1

-1

1

1

1

-1

-1

-1

1

1

1

39

43,13

15

1

-1

1

1

1

-1

-1

-1

-1

-1

1

1

1

1

1

1

38

41,88

16

1

1

1

1

1

-1

1

1

1

1

1

1

1

1

1

1

37

45,13

17

1

-1

-1

-1

-1

1

1

1

1

1

1

1

1

1

1

1

38

47,88

18

1

1

-1

-1

-1

1

-1

-1

-1

-1

1

1

1

1

1

1

39

44,13

19

1

-1

1

-1

-1

1

-1

1

1

1

-1

-1

-1

1

1

1

40

44,88

20

1

1

1

-1

-1

1

1

-1

-1

-1

-1

-1

-1

1

1

1

41

41,13

21

1

-1

-1

1

-1

1

1

-1

1

1

-1

1

1

-1

-1

1

42

46,38

22

1

1

-1

1

-1

1

-1

1

-1

-1

-1

1

1

-1

-1

1

43

42,63

23

1

-1

1

1

-1

1

-1

-1

1

1

1

-1

-1

-1

-1

1

44

43,38

24

1

1

1

1

-1

1

1

1

-1

-1

1

-1

-1

-1

-1

1

45

39,63

25

1

-1

-1

-1

1

1

1

1

-1

1

1

-1

1

-1

1

-1

43

46,63

26

1

1

-1

-1

1

1

-1

-1

1

-1

1

-1

1

-1

1

-1

41

40,38

27

1

-1

1

-1

1

1

-1

1

-1

1

-1

1

-1

-1

1

-1

39

38,63

28

1

1

1

-1

1

1

1

-1

1

-1

-1

1

-1

-1

1

-1

37

32,38

29

1

-1

-1

1

1

1

1

-1

-1

1

-1

-1

1

1

-1

-1

35

35,13

30

1

1

-1

1

1

1

-1

1

1

-1

-1

-1

1

1

-1

-1

33

28,88

31

1

-1

1

1

1

1

-1

-1

-1

1

1

1

-1

1

-1

-1

31

27,13

32

1

1

1

1

1

1

1

1

1

-1

1

1

-1

1

-1

-1

29

20,88

 

39,38

-0,13

-0,25

-0,50

-1,13

-0,63

0,00

0,00

-0,63

2,38

0,00

-1,25

2,50

-2,50

2,75

3,88

39,38

 

x5=x1x2

25-1

x0

x1

x2

x3

x4

x5

x1x2

x1x3

x1x4

x1x5

x2x3

x2x4

x2x5

x3x4

x3x5

x4x5

y

y''

1

1

-1

-1

-1

-1

-1

1

1

1

1

1

1

1

1

1

1

36

48,50

2

1

1

-1

-1

-1

-1

-1

-1

-1

-1

1

1

1

1

1

1

37

49,50

3

1

-1

1

-1

-1

-1

-1

1

1

1

-1

-1

-1

1

1

1

38

50,50

4

1

1

1

-1

-1

1

1

-1

-1

1

-1

-1

1

1

-1

-1

39

28,50

5

1

-1

-1

1

-1

1

1

-1

1

-1

-1

1

-1

-1

1

-1

40

38,75

6

1

1

-1

1

-1

-1

-1

1

-1

-1

-1

1

1

-1

-1

1

41

44,25

7

1

-1

1

1

-1

-1

-1

-1

1

1

1

-1

-1

-1

-1

1

42

45,25

8

1

1

1

1

-1

1

1

1

-1

1

1

-1

1

-1

1

-1

43

41,75

9

1

-1

-1

-1

1

1

1

1

-1

-1

1

-1

-1

-1

-1

1

44

41,25

10

1

1

-1

-1

1

-1

-1

-1

1

-1

1

-1

1

-1

1

-1

43

44,25

11

1

-1

1

-1

1

-1

-1

1

-1

1

-1

1

-1

-1

1

-1

42

42,75

12

1

1

1

-1

1

1

1

-1

1

1

-1

1

1

-1

-1

1

41

36,75

13

1

-1

-1

1

1

1

1

-1

-1

-1

-1

-1

-1

1

1

1

40

44,50

14

1

1

-1

1

1

-1

-1

1

1

-1

-1

-1

1

1

-1

-1

39

29,00

15

1

-1

1

1

1

-1

-1

-1

-1

1

1

1

-1

1

-1

-1

38

27,50

16

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

37

40,00

17

1

-1

-1

-1

-1

1

1

1

1

1

1

1

1

1

1

1

38

44,00

18

1

1

-1

-1

-1

-1

-1

-1

-1

-1

1

1

1

1

1

1

39

49,50

19

1

-1

1

-1

-1

-1

-1

1

1

1

-1

-1

-1

1

1

1

40

50,50

20

1

1

1

-1

-1

1

1

-1

-1

1

-1

-1

1

1

-1

-1

41

28,50

21

1

-1

-1

1

-1

1

1

-1

1

-1

-1

1

-1

-1

1

-1

42

38,75

22

1

1

-1

1

-1

-1

-1

1

-1

-1

-1

1

1

-1

-1

1

43

44,25

23

1

-1

1

1

-1

-1

-1

-1

1

1

1

-1

-1

-1

-1

1

44

45,25

24

1

1

1

1

-1

1

1

1

-1

1

1

-1

1

-1

1

-1

45

41,75

25

1

-1

-1

-1

1

1

1

1

-1

-1

1

-1

-1

-1

-1

1

43

41,25

26

1

1

-1

-1

1

-1

-1

-1

1

-1

1

-1

1

-1

1

-1

41

44,25

27

1

-1

1

-1

1

-1

-1

1

-1

1

-1

1

-1

-1

1

-1

39

42,75

28

1

1

1

-1

1

1

1

-1

1

1

-1

1

1

-1

-1

1

37

36,75

29

1

-1

-1

1

1

1

1

-1

-1

-1

-1

-1

-1

1

1

1

35

44,50

30

1

1

-1

1

1

-1

-1

1

1

-1

-1

-1

1

1

-1

-1

33

29,00

31

1

-1

1

1

1

-1

-1

-1

-1

1

1

1

-1

1

-1

-1

31

27,50

32

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

29

40,00

39,38

-0,13

-0,25

-0,50

-1,13

-2,25

0,00

0,00

-0,63

4,38

0,00

-1,25

4,50

-2,50

4,63

4,63

 

 

x5=x2x3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x4=x3x1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

x1

x2

x3

x4

x5

x1x2

x1x3

x1x4

x1x5

x2x3

x2x4

x2x5

x3x4

x3x5

x4x5

y

y'''

1

1

-1

-1

-1

-1

-1

1

1

1

1

1

1

1

1

1

1

36

50,88

2

1

1

-1

-1

-1

1

-1

-1

-1

1

1

1

-1

1

-1

-1

37

46,96

3

1

-1

1

-1

1

-1

-1

1

-1

1

-1

1

-1

-1

1

-1

38

46,71

4

1

1

1

-1

-1

-1

1

-1

-1

-1

-1

-1

-1

1

1

1

39

32,46

5

1

-1

-1

1

-1

-1

1

-1

1

1

-1

1

1

-1

-1

1

40

50,54

6

1

1

-1

1

1

-1

-1

1

1

-1

-1

-1

1

1

-1

-1

41

36,79

7

1

-1

1

1

-1

1

-1

-1

1

-1

1

-1

1

-1

1

-1

42

36,54

8

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

43

40,79

9

1

-1

-1

-1

1

1

1

1

-1

-1

1

-1

-1

-1

-1

1

44

24,21

10

1

1

-1

-1

-1

1

-1

-1

-1

1

1

1

-1

1

-1

-1

43

46,96

11

1

-1

1

-1

1

-1

-1

1

-1

1

-1

1

-1

-1

1

-1

42

46,71

12

1

1

1

-1

-1

-1

1

-1

-1

-1

-1

-1

-1

1

1

1

41

32,46

13

1

-1

-1

1

-1

-1

1

-1

1

1

-1

1

1

-1

-1

1

40

50,54

14

1

1

-1

1

1

-1

-1

1

1

-1

-1

-1

1

1

-1

-1

39

36,79

15

1

-1

1

1

-1

1

-1

-1

1

-1

1

-1

1

-1

1

-1

38

36,54

16

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

37

40,79

17

1

-1

-1

-1

1

1

1

1

-1

1

1

-1

1

-1

1

-1

38

40,96

18

1

1

-1

-1

-1

1

-1

-1

-1

1

1

1

-1

1

-1

-1

39

46,96

19

1

-1

1

-1

1

-1

-1

1

-1

1

-1

1

-1

-1

1

-1

40

46,71

20

1

1

1

-1

-1

-1

1

-1

-1

-1

-1

-1

-1

1

1

1

41

32,46

21

1

-1

-1

1

-1

-1

1

-1

1

1

-1

1

1

-1

-1

1

42

50,54

22

1

1

-1

1

1

-1

-1

1

1

-1

-1

-1

1

1

-1

-1

43

36,79

23

1

-1

1

1

-1

1

-1

-1

1

-1

1

-1

1

-1

1

-1

44

36,54

24

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

45

40,79

25

1

-1

-1

-1

1

1

1

1

-1

-1

1

-1

-1

-1

-1

1

43

24,21

26

1

1

-1

-1

-1

1

-1

-1

-1

1

1

1

-1

1

-1

-1

41

46,96

27

1

-1

1

-1

1

-1

-1

1

-1

1

-1

1

-1

-1

1

-1

39

46,71

28

1

1

1

-1

-1

-1

1

-1

-1

-1

-1

-1

-1

1

1

1

37

32,46

29

1

-1

-1

1

-1

-1

1

-1

1

1

-1

1

1

-1

-1

1

35

50,54

30

1

1

-1

1

1

-1

-1

1

1

-1

-1

-1

1

1

-1

-1

33

36,79

31

1

-1

1

1

-1

1

-1

-1

1

-1

1

-1

1

-1

1

-1

31

36,54

32

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

29

40,79

 

39,38

-0,13

-0,25

-0,50

-2,25

-2,25

0,00

0,00

1,75

4,63

0,00

2,25

4,13

2,13

4,38

-2,38

 

 

Условия

25

25-1

25-2

y

y'

y''

y'''

36

45,63

48,50

50,88

37

29,88

49,50

46,96

38

31,13

50,50

46,71

39

36,63

28,50

32,46

40

33,63

38,75

50,54

41

39,13

44,25

36,79

42

40,38

45,25

36,54

43

45,88

41,75

40,79

44

38,38

41,25

24,21

43

41,38

44,25

46,96

42

40,13

42,75

46,71

41

43,13

36,75

32,46

40

37,63

44,50

50,54

39

40,63

29,00

36,79

38

39,38

27,50

36,54

37

42,38

40,00

40,79

38

27,63

44,00

40,96

39

37,13

49,50

46,96

40

37,88

50,50

46,71

41

43,38

28,50

32,46

42

39,38

38,75

50,54

43

44,88

44,25

36,79

44

45,63

45,25

36,54

45

51,13

41,75

40,79

43

39,63

41,25

24,21

41

42,63

44,25

46,96

39

40,88

42,75

46,71

37

43,88

36,75

32,46

35

37,38

44,50

50,54

33

40,38

29,00

36,79

31

38,63

27,50

36,54

29

41,63

40,00

40,79

Заключение

Теория ПЭ охватывает практически все встречающиеся на практике варианты исследования объектов. В дальнейшем будут рассмотрены следующие типовые задачи экспериментального исследования:

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

· приближенное аналитическое описание функциональной связи показателей качества с параметрами системы по результатам проведенного эксперимента. Традиционные методики проведения экспериментов из-за зависимости компонентов восстанавливаемого аналитического описания не позволяют определить раздельное влияние каждого фактора на результирующий показатель, т. е. эти методики обеспечивают получение аналитических зависимостей, пригодных лишь для решения интерполяционных задач. В отличие от них ТПЭ дает возможность оценить вклад каждого параметра в значение показателя, т.е. приближенно восстановить закон функционирования объекта по экспериментальным данным. Полученное аналитическое описание объекта можно использовать для предварительного исследования вариантов построения системы или в интересах построения модели старшей системы, включающей данный объект на правах элемента;

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

Кроме указанных, существуют и другие задачи, решаемые с помощью ТПЭ, например:

· испытания образцов техники. Планирование должно позволить оценить степень соответствия показателей качества образцов заданным требованиям при минимальном объеме испытаний;

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

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

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

Список литературы

1. Асатурян В.И. Теория планирования эксперимента. - М.: Радио и связь, 1983.

2. Налимов В.В. Теория эксперимента. - М.: Наука, 1971.

3. Таблицы планов эксперимента для факторных и полиномиальных моделей / Под ред. В.В. Налимова. - М.: Металлургия, 1982.

Размещено на Allbest.ru


Подобные документы

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

    лекция [220,4 K], добавлен 13.11.2009

  • Сущность и особенности планирования эксперимента, кодирование исходных факторов. Составление плана эксперимента для определения зависимости концентрации меди от расхода шихты, содержания кислорода в дутье. Выбор математической модели объекта исследования.

    курсовая работа [1,8 M], добавлен 11.12.2012

  • Общие сведения о планировании эксперимента. Анализ методики составления планов эксперимента для моделей первого и второго порядков. Положения о планировании второго порядка. Ортогональные и рототабельные центральные композиционные планы второго порядка.

    реферат [242,7 K], добавлен 22.06.2011

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

    курсовая работа [42,1 K], добавлен 03.11.2014

  • Составление матрицы плана факторного эксперимента и разработка матрицы его базисных функций. Написание алгебраического полинома плана и корреляционный анализ результатов эксперимента. Функция ошибки и среднеквадратичное отклонение регрессионной модели.

    контрольная работа [698,2 K], добавлен 13.06.2014

  • Формирование иерархии при решении проблемы "выбор фрезы". Третий этап окончательного определения. Глобальные приоритеты выбора. Полный факторный эксперимент. Определение однородности дисперсий. Расчетные значения критериев. Неполная квадратичная модель.

    курсовая работа [1,6 M], добавлен 12.09.2014

  • Основные параметры сетевой модели системы планирования и управления. Правила построения сетевых графиков. Характеристики элементов сетевой модели. Метод пересмотра планов. Численная реализация задачи сетевого планирования. Метод графической оценки.

    реферат [154,4 K], добавлен 19.03.2015

  • Задачи сетевого планирования и управления. Виды операций: составные, параллельные, зависимые и независимые. Полный и независимый резерв времени для критических операций. Приведение модели к каноническому виду. Решение задач двойственным симплекс-методом.

    курсовая работа [2,9 M], добавлен 20.05.2014

  • Факторный анализ. Задачи факторного анализа. Методы факторного анализа. Детерминированный факторный анализ. Модели детерминированного факторного анализа. Способы оценки влияния факторов детерминированном факторном анализе. Стохастический анализ.

    курсовая работа [150,0 K], добавлен 03.05.2007

  • Получение функции отклика показателя качества Y2 и формирование выборки объемом 15 и более 60. Зависимость выбранного Y от одного из факторов Х. Дисперсионный анализ и планирование эксперимента. Проведение корреляционного и регрессионного анализа.

    курсовая работа [827,2 K], добавлен 19.06.2012

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.