Методология научных исследований. Планирование экспериментов

Сентябрь 13, 2022

Главная
Математика
Методология научных исследований. Планирование экспериментов

Содержание

2. Планирование эксперимента и его задачи. Виды экспериментов
3. Планирование эксперимента – это процедура выбора числа и условий проведения опытов, необходимых и достаточных для решения
4. КЛАССИФИКАЦИЯ ЭКСПЕРИМЕНТОВ В зависимости от условий эксперименты делятся на несколько видов: 1) промышленный –это эксперимент, поставленный
5. 5) пошаговый – эксперимент, состоящий из отдельных серий опытов. Причем условия проведения каждой следующей серии определяются
6. Планирование оптимальных экспериментов
7. Факторы, параметры оптимизации и модели Выбор факторов, параметров оптимизации и моделей осуществляется с учетом цели исследований
8. Параметр оптимизациии - характеристика цели исследования (выходной параметр), заданная количественно. Параметр оптимизации еще называют (критерием оптимизации,
9. Классификация параметров оптимизации: Параметр оптимизации (отклик )– величина, описывающая результат проведенного эксперимента и зависящая от факторов,
10. Требования к параметрам оптимизации. Параметр оптимизации должен быть: 1 количественным и однозначным, т.е.задаваться одним числом и
11. Факторное пространство и поверхность отклика Каждому фиксированному набору уровней факторов соответствует определенная точка в многомерном пространстве
12. Математические модели Cвязь между уровнями факторов и реакцией (откликом) системы описывают математическими моделями вида: Функцию Фi
13. Требования к модели: адекватность, содержательность, простота и др. Адекватность - способность модели предсказывать результаты эксперимента в
14. Основа МПЭ - представление объекта исследования как кибернетической системы, моделью которой является схема, изображенная на рисунке
15. Поскольку полином* содержитCdk+d коэффициентов, подлежащих определению, план D должен содержать по крайней мере Cdk+d различных экспериментальных
16. Планирование эксперимента для линейного приближения по- поверхности отклика. Вначале выбирают локальную область факторного пространства (область определения).
17. Первый этап планирования эксперимента для получения линейной модели основан на варьировании факторов на двух уровнях: хiв
19. Выписывая комбинации уровней факторов для каждой экспериментальной точки квадрата, изображенного на рис. 4, б, получим план
20. Для оценки свободного члена bо (ожидаемого отклика в центре плана) и определения эффектов взаимодействия b12, b13,
21. Реализация эксперимента 1 Подсчитывают количество исходного сырья и его подготавливают. Сырье должно быть однородным. Если нет
22. Результат (г + 1)-го опыта отбрасывается, если Значения дпя t' для различных α и ρ приведены
23. Реализовав параллельные опыты в точках плана и определив соответствующие значения s2 {у}, следует убедиться в однородности
24. Значения F-критерия для уровня значимости α=0,05 (фрагмент)
25. Проверка значимости отличия откликовYmax и Ymin Вычислив средние значения откликов для каждой точки плана, можно проверить,
26. Обработка результатов эксперимента Метод наименьших квадратов — эффективный и простой способ получения оценок коэффициентов регрессии, если
27. Согласно методу наименьших квадратов находятся такие значения оценок bi величин , которые минимизируют сумму квадратов отклонений
28. Пример. Получение математической модели, описывающей зависимость средних значений площадей поперечных сечений единичных срезов от условий реализации
29. После проведения эксперимента в точках плана и получения всех значений выходного параметра (Y1, … Y8) коэффициенты
30. Матрица планирования эксперимента 23
31. Значения площадей сечения единичного среза , мкм2 * В скобках указаны: размеры зерен dз (мкм), соответствующие
32. план эксперимента. Факторы: зернистость круга dз – Х1; № структуры круга – Х2; угол атаки ω
33. Значения коэффициентов полинома y =144,5 x0+ 104,7x1+ 13,5x2+125 x3+9,7 x1 x2+ 90,5x1 x3+ 11,6x2 x3+ 8,4x1
34. Планирование эксперимента на диаграммах состав—свойство
35. Переменные xi (i = 1, 2, . . ., q) таких систем являются пропорциями (относительным содержанием)
36. Переход к симплексной системе координат
37. Каноническая форма полинома степе- степени n (Шеффе) где Полиномы такого вида (так называемые приведенные полиномы) получаются
38. например, полином второй степени в общем случае имеющий вид в приведенной форме с учетом условия запишется
39. Для оценки коэффициентов аппроксимирующего полинома степени n во всех точках плана, соответствующих узлам {q, n}-решетки, реализуются
40. формулы для оценки коэффициентов приведенных полиномов Для n ≤ 4 они имеют вид: Модель первого порядка:
41. Формулы для расчета дисперсии значений исследуемого свойства : для квадратичной модели для неполной кубической модели
42. Контурные кривые поверхностей отклика, описываемые моделями 2-го порядка
44. Скачать презентацию

Слайд 2

Планирование эксперимента и его задачи.
Виды экспериментов

Слайд 3

Планирование эксперимента –
это процедура выбора числа и условий проведения опытов, необходимых

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

Слайд 4

КЛАССИФИКАЦИЯ ЭКСПЕРИМЕНТОВ
В зависимости от условий эксперименты делятся на несколько

видов:
1) промышленный –это эксперимент, поставленный в условиях предприятия с целью улучшения производства;
2) научно-исследовательский –эксперимент, поставленный в научно-исследовательских лабораториях с целью исследования нового или улучшения существующего процесса, явления;
3) лабораторный - эксперимент, поставленный в научно исследовательских лабораториях с целью изучения хорошо известного, существующего процесса, явления;
4) оптимальный (экстремальный) – эксперимент, поставленный с целью поиска наиболее оптимальных условий его реализации. С математической точки зрения, это эксперимент по поиску экстремумов некоторой функции, отсюда и второе название эксперимента;

Слайд 5

5) пошаговый – эксперимент, состоящий из отдельных серий опытов. Причем условия

проведения каждой следующей серии определяются результатами предыдущих.
6) активный - эксперимент, в ходе которого экспериментатор имеет возможность изменять и/или поддерживать на заданном уровне сколь угодно долго значения параметров, задающих условия проведения эксперимента;
7) пассивный - эксперимент, в ходе которого экспериментатор НЕ имеет такой возможности

Слайд 6

Планирование оптимальных экспериментов

Слайд 7

Факторы, параметры оптимизации и модели
Выбор факторов, параметров оптимизации и

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

Слайд 8

Параметр оптимизациии
- характеристика цели исследования (выходной параметр), заданная количественно.
Параметр

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

Слайд 9

Классификация параметров оптимизации:
Параметр оптимизации (отклик )– величина, описывающая результат проведенного

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

Слайд 10

Требования к параметрам оптимизации.
Параметр оптимизации должен быть:
1 количественным и однозначным,

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

Слайд 11

Факторное пространство
и поверхность отклика
Каждому фиксированному набору уровней факторов соответствует определенная точка

в многомерном пространстве факторов, называемом факторным пространством, а набор значений параметра оптимизации в пространстве– поверхность отклика

Факторное пространство (а) и поверхность отклика (б)

Слайд 12

Математические модели
Cвязь между уровнями факторов и реакцией (откликом) системы

описывают математическими моделями вида:
Функцию Фi называют функцией отклика, а ее геометрический образ — поверхностью отклика.
Цель эксперимента – найти функцию отклика.
Т.к. вид этой функции неизвестен, то получают приближенные уравнения.
Эксперимент необходимо поставить так, чтобы при минимальном количестве опытов, варьируя значения независимых переменных по специальным правилам, построить математическую модель системы и найти оптимальные значения свойств системы.

yi = Фi(х1+х2+…+хl )

Слайд 13

Требования к модели: адекватность, содержательность, простота и др.
Адекватность -

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

Слайд 14

Основа МПЭ - представление объекта исследования как кибернетической системы, моделью которой

является схема, изображенная на рисунке и называемая «черным ящиком».

Основные понятия математического планирования эксперимента (МПЭ)

Слайд 15

Поскольку полином* содержитCdk+d коэффициентов, подлежащих определению, план D должен содержать

по крайней мере Cdk+d различных экспериментальных точек

где xiu — значения, которые принимает i-я переменная в u-м опыте (i- 1,2, ..., к; и u - 1,2,..., N).

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

Реализовав опыты в N точках, получим вектор наблюдений, следующего вида, где Уu — отклик, соответствующий u-той точке плана .

Слайд 16

Планирование эксперимента для линейного приближения по-
поверхности отклика.
Вначале выбирают локальную область

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

Слайд 17

Первый этап планирования эксперимента для получения линейной модели основан на варьировании

факторов на двух уровнях: хiв и xiн , симметрично расположенных относительно нулевого (основного) уровня хi0 (рис. 4, а).или в безразмерном виде - кодах: +1,0,-1
Переход от исходных (натуральных) переменных xi к безразмерным кодированным, осуществляют по формуле

Расположение экспериментальных точек с координатами +1 и -1 при полном факторном эксперименте (ПФЭ) для двух исследуемых факторов приведено на рис. 4, а для трех — на рис. 5. точки

Слайд 18

Слайд 19

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

рис. 4, б, получим план D полного факторного эксперимента(ПФЭ) типа 22 (табл. 1.1) и для куба — план D типа 23 (рис. 5, табл. 1.2).

Слайд 20

Для оценки свободного члена bо (ожидаемого отклика в центре плана) и

определения эффектов взаимодействия b12, b13, ..., b123, ... план эксперимента D расширяют до так называемой матрицы планирования X добавлением соответствующего «фиктивной переменной» единичного столбца X0 и столбцов произведений Х1Х2, Х1Х3, Х2,Х3, Х1Х2Х3….. (см., например, табл. 1.3).

Слайд 21

Реализация эксперимента
1 Подсчитывают количество исходного сырья и его подготавливают. Сырье

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

среднее арифметическое значений параметров оптимизации в v-ой точке плана;

— число степеней свободы, равное количеству опытов минус единица.

Слайд 22

Результат (г + 1)-го опыта отбрасывается, если
Значения дпя t' для различных

α и ρ приведены в специальных таблицах

Фрагмент таблицы

Слайд 23

Реализовав параллельные опыты в точках плана и определив соответствующие значения

s2 {у}, следует убедиться в однородности дисперсий, т. е. в том, что среди всех дисперсий нет таких, которые бы значительно превышали все остальные.
Требование однородности дисперсий является одним из требований регрессионного анализа, и поэтому проверке однородности дисперсий следует уделять должное внимание. Для сравнения двух дисперсий применяется F-критерий (критерий Фишера), представляющий отношение большей дисперсии к меньшей. Согласно критерию Фишера гипотеза об однородности двух дисперсий s21 {у} и s22 {у} (s21 {у} > s22 {у} ) определенных с f1 и f2 степенями свободы соответственно, не отвергается, если расчетное отношение Fp = s21 {у} / s22{у} не превышает табличного значения Fτ для числа степеней свободы f1 и f2 и для уровня значимости α (обычно α= 0,05). Таблица F-критерия приводятся в справочной литературе, а для уровня значимости α= 0,05 приведена ниже.

Слайд 24

Значения F-критерия для уровня значимости α=0,05 (фрагмент)

Слайд 25

Проверка значимости
отличия откликовYmax и Ymin
Вычислив средние значения

откликов для каждой точки плана, можно проверить, значимо ли отличаются друг от друга Ymax и Ymin.
Значимость различия двух средних можно проверить с помощью t-критерия (критерия Стьюдента) по формуле:

Средние значимо не отличаются, если экспериментальное значение t -критерия не превосходит табличного tT (см.таблицу) для f = r max + r min степеней свободы и обычно 5%-ного уровня значимости, т. е. вероятность того, что при f = r max + r min степенях свободы значение величины t будет больше по абсолютной величине, чем tT, уровня значимости равного 0,05. Если значение t меньше табличного, то с вероятностью Р ~ 1 — а = 0,95 можно считать, что разницы между результатами двух опытов нет.

Слайд 26

Обработка результатов эксперимента
Метод наименьших квадратов — эффективный и простой способ

получения оценок коэффициентов регрессии,
если независимые переменные х достаточно точно поддерживаются на определенных уровнях, а наблюдаемые значения отклика –у1, у2, ..., уN по данным N опытов плана D представляют собой независимые и нормально распределенные случайные величины.
Полученный в результате опытов ограниченный статистический материал дает возможность определить лишь оценки b0, bl .. ..., bк' теоретических коэффициентов регрессии b0, b1 ..., bк в уравнении (*), справедливых для некоторой гипотетической совокупности, состоящей из всех мыслимых опытов. Тогда уравнение регрессии,полученное на основании N опытов, запишется следующим образом:

где у — значение выхода, предсказанное
уравнением для ряда входных условий.

Слайд 27

Согласно методу наименьших квадратов находятся такие значения оценок bi величин ,

которые минимизируют сумму квадратов отклонений (невязок) ɛи опытных точек от величин, предсказанных регрессионным уравнением (**), т. е. минимизирующие функцию

Слайд 28

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

единичных срезов от условий реализации операции шлифования (интенсивности съема Q=vдt, характеризуемой углом атаки ; структуры шлифовального круга - № структуры; зернистости круга -№ зернистисти).
Предположим, что для описания зависимости выходного параметра Y от трех независимых факторов X1, X2, X3 достаточно неполного полинома третьей степени:
y =b0 x0+ b1 x1+ b2 x2+b3 x3+b12 x1 x2+ b13 x1 x3+ b23 x2 x3+ b123 x1 x2 x3;
где b0 – свободный член, а b0, b1, b2… b123 -коэффициенты регрессионного уравнения;
x1, x2, x3 – псевдокоординаты факторных переменных X1, X2, X3 соответственно, принимающие значения -1 или +1, под которыми закодировано максимальное и минимальное значение фактора

Слайд 29

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

параметра (Y1, … Y8) коэффициенты b0, b1 …b123 рассчитываются как среднее значение выходных параметров, взятых со знаком столбца сочетания факторов, соответствующих вычисляемому коэффициенту. Так, например, коэффициент b12 определяется по формуле:

Для проверки адекватности полученной матмодели ставятся
эксперименты в дополнительных одной-трех точках плана (на рис. 6 это точки 9 и 10), называемых проверочными. Полученные в этих точках экспериментальные значения выходного параметра сравниваются по известному алгоритму с рассчитанными для этих же условий по регрессионному уравнению и делается заключение о пригодности данной матмодели для дальнейшего исследования.

Слайд 30

Матрица планирования эксперимента 23

Слайд 31

Значения площадей сечения единичного среза
, мкм2
* В скобках указаны: размеры

зерен dз (мкм), соответствующие № зернистости и объемные доли зерна Кз в круге, соответствующие номеру № его структуры

Слайд 32

план эксперимента.
Факторы: зернистость круга dз – Х1; № структуры круга –

Х2; угол атаки ω –Х3 и составим матрицу плана полнофакторного эксперимента 23

Слайд 33

Значения коэффициентов полинома
y =144,5 x0+ 104,7x1+ 13,5x2+125 x3+9,7 x1 x2+ 90,5x1

x3+ 11,6x2 x3+ 8,4x1 x2 x3;

полином

Далее осуществляется проверка адекватности полинома

Слайд 34

Планирование эксперимента
на диаграммах состав—свойство

Слайд 35

Переменные xi (i = 1, 2, . . ., q) таких

систем являются пропорциями (относительным содержанием) i-x компонентов смеси и удовлетворяют условию
***

Геометрическое место точек, удовлетворяющих условию нормированности суммы переменных (***), представляет собой (q — 1)- мерный правильный симплекс (треугольник для q = 3, тетраэдр для q = 4 и т. д.). Каждой точке такого симплекса соответствует смесь определенного состава.

Слайд 36

Переход к симплексной системе координат

Слайд 37

Каноническая форма полинома степе-
степени n (Шеффе)
где
Полиномы такого вида (так называемые

приведенные полиномы) получаются из обычных полиномов * соответствующей степени для q переменных введением соотношения **

Слайд 38

например, полином второй степени в общем случае имеющий вид
в приведенной форме

с учетом условия

запишется

Для оценки коэффициентов приведенного полинома были предложены планы Точками таких планов являются узлы {q, n}-симплексных реше- решеток

q = 3: а — линейная {3,1}; б — квадратичная {3,2};
в — неполно-кубическая; г — кубическая 3,3); 9 — четвертой степени {3,4};
е - восьмой степени {3,8}; q = 4: ж - квадратичная {4,2}; з - кубическая {,4,3}

Слайд 39

Для оценки коэффициентов аппроксимирующего полинома степени n во всех точках плана,

соответствующих узлам {q, n}-решетки, реализуются опыты и определяются отклики системы у.

Обозначение откликов в точках симплексных решеток

Обозначение откликов в матрице планирования для {3, 3}-решетки

Слайд 40

формулы для оценки коэффициентов приведенных полиномов
Для n ≤ 4 они

имеют вид: Модель первого порядка:

для трехкомпонентной смеси:

для q-компонентной смеси

где

Модель второго порядка:

где

для q -компонентной смеси

где

Для трехкопонентной смеси

Слайд 41

Формулы для расчета дисперсии значений исследуемого свойства :
для квадратичной модели
для

неполной кубической модели

Слайд 42

Методология научных исследований. Планирование экспериментов

Содержание

Планирование эксперимента и его задачи.Виды экспериментов

Планирование эксперимента –это процедура выбора числа и условий проведения опытов, необходимых

КЛАССИФИКАЦИЯ ЭКСПЕРИМЕНТОВ В зависимости от условий эксперименты делятся на несколько

5) пошаговый – эксперимент, состоящий из отдельных серий опытов. Причем условия

Планирование оптимальных экспериментов

Факторы, параметры оптимизации и модели Выбор факторов, параметров оптимизации и

Параметр оптимизациии - характеристика цели исследования (выходной параметр), заданная количественно. Параметр

Классификация параметров оптимизации:Параметр оптимизации (отклик )– величина, описывающая результат проведенного

Требования к параметрам оптимизации.Параметр оптимизации должен быть: 1 количественным и однозначным,

Факторное пространствои поверхность откликаКаждому фиксированному набору уровней факторов соответствует определенная точка

Математические модели Cвязь между уровнями факторов и реакцией (откликом) системы

Требования к модели: адекватность, содержательность, простота и др. Адекватность -