🗊Презентация Обработка результатов эксперимента. Матричное исчисление

Обработка результатов эксперимента Матричное исчисление

Регрессионный анализ Это статистический метод проверки гипотезы о пригодности модели и значимости коэффициентов. Он основывается на следующих постулатах: Параметр у – есть случайная величина. Дисперсия у – не зависит от абсолютной величины у. Значения факторов – есть не случайные величины. Наиболее простым классом регрессионной модели является y=b0+b1x1+b2x2+….+ Вектор параметров такой модели находят при условии минимизации ее ошибки, т.е. МНК

Матричные исчисления в многомерной статистике Для нахождения коэффициента регрессии в матричном исчис-лении необходимо решить следующие уравнения: , где X– матрица условий эксперимента. Y – матрица результатов эксперимента, В - матрица столбец коэффи-циентов уравнения x 0…N условный фактор для оценки b0 , т.к. для остальных b1 ,b2 и т.д. есть значения .

Решим матрицу относительно В Умножив обе части равенства на транспонированную матрицу: Получим информационную матрицу Фишера Ее структура приведенная информационная матрица

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

Каждый bi определяется по следующей формуле При ортогональном планировании: Аналогичны рассуждения при проверке дисперсии экспериментов

Ковариация Если ковариация положительна, то с ростом значений одной случайной величины, значения второй имеют тенденцию возрастать, а если знак отрицательный — то убывать. Однако по абсолютному значению ковариации нельзя судить о том, насколько сильно величины взаимосвязаны, так как её масштаб зависит от их дисперсий. Масштаб можно нормировать, поделив значение ковариации на произведение среднеквадратических отклонений, т.е.это коэффициент корреляции, который всегда находится в интервале от −1 до 1.

КРИТЕРИИ ОПТИМАЛЬНОСТИ ПЛАНОВ

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

Первая группа 1.1. Ортогональность. а). диагональность матрицы дисперсии – ковариации. б). направление главных осей эллипсоида рассеяния совпадает с направлением координатных осей пространства. в). оценки параметров коэффициентов независимы друг от друга. 1.2. D – оптимальные планы. а). минимум определителя матрицы дисперсии – ковариации. б). минимум объема эллипсоида рассеяния параметров. в). минимум обобщенной дисперсии всех коэффициентов.  1.3. А – оптимальность. а). минимум следа матрицы дисперсии, ковариации. б). минимум суммы квадратов длин осей эллипсоида рассеяния и минимум длины диагонали прямоугольника, описанного около эллипсоида. в). минимальная средняя дисперсия оценок эксперимента.

    1.4. Е – оптимальные планы (minmax). а). minmax собственного значения матриц дисперсии, ковариации. б). min максимальной оси эллипсоида рассеяния. в). отдельные оценки параметров не обладают слишком большими дисперсией и ковариацией.   1.5. minmax оценки дисперсии коэффициентов. а). minmax диагонали элемента матрицы. б). минимум максимальной проекции оси эллипсоида рассеяния на координатные оси пространства параметров. в). minmax дисперсии оценки коэффициентов.

Вторая группа 2.1. ротатабельность. а). , где  - сферическая координата. в). одинаковая точность предсказаний для точек, равноудаленных от центра плана по любому направлению. 2.2. униформность. а). при в). дисперсия предсказания постоянна в некоторой области вокруг центра эксперимента (интервал варьирования может увеличиваться и уменьшаться без болезненно).  2.3. G – оптимальность. а) в) Minmax – максимального значения дисперсии оценки поверхности отклика.  2.4.Q – оптимальность. а) б) Минимум средних дисперсий оценки поверхности отклика.  Критерии насыщенности плана и критерии композиционности плана.

Критерии насыщенности и композиционности плана. Насыщенность плана обеспечивает минимум числа опытов. Минимум задается числом коэффициентов модели. Приближение к нему служит мерой насыщенности плана. В общем случае КнасN(k+1), где k- число констант модели. Процессы, протекающие в многофазных системах, как правило, не удается описать линейными регрессионными закономерностями. Поэтому необходимо использовать математические модели более высокого порядка. Это и есть ортогональное композиционное планирование эксперимента.

ПЛАНЫ ВТОРОГО ПОРЯДКА

Особенности планов второго порядка

Планы типа 3n Планы с числом степеней свободы, достаточным для получения оценок всех коэффициентов данной модели, можно было бы построить, например, введением трех уровней варьирования для каждой из переменных исследуемой системы. Но такие планы обладали бы большой избыточностью степеней свободы, т. е. избыточностью числа опытов. n = 2, коэффициентов – 6, опытов – 9. Для n = 3, 27 опытов, 10 коэффициентов. Поэтому такие планы полностью никогда не реализуют. Действительно реализуемыми остаются например: Это матрица плана второго порядка, построенного для n=3 на базе матрицы ПФЭ З3. Как видно, из 27 опытов в матрице плана осталось всего 12 и добавлены три опыта в центре, т. е. предусматривается всего 15 опытов

Коэффициенты модели где A, B, C, D – некоторые коэффициенты, взятые из таблицы. Ошибки оценок коэффициентов модели могут быть вычислены через соответствующие дисперсии Геометрический смысл плана

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

Параметры ортогональных композиционных панов

Расчет коэффициентов моделей

Влияние композицоинности плана Если в центре плана поставить не один эксперимент, то может быть определено такое число точек, что постоянство информации будет иметь место не только в равноотстоящих от центра плана точках, но и на большей части отрезка радиус-вектора 0    1. Информационные контуры для центральных комплексных планов второго порядка:

Пример использования симметричного ортогонального композиционного плана второго порядка Изучали механические свойства одного из алюминиевых деформированных сплавов в зависимости от содержания лития, температуры и времени старения. В качестве отклика выбрали предел прочности. По результатам измерений, проведенных ранее, было установлено, что дисперсия 4, число степеней свободы 10. Был реализован ПФЭ 23, по результатам кото- рого рассчитаны коэффициенты для линейной модели:

Дисперсия в определении коэффициентов: Дисперсия в определении коэффициентов:   Дисперсия неадекватности: , где N- количество опытов. K’-количество коэффициентов.

Проведем дополнение модели В таком случае можно дополнить данную модель звездными точками и опытами в центре плана. Для того, чтобы сделать матрицу планирования ортогональной, мы должны ввести столбец новой переменной для нахождения коэффициентов регрессии при квадратичных планах. Для 22=0,67; для 23=0,73.

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

По результатам эксперимента, используя коэффициенты рассчитываем коэффициенты регрессии: По результатам эксперимента, используя коэффициенты рассчитываем коэффициенты регрессии: По критерию Стьюдента: статистически значимыми будут Уравнение регрессии: Критерий Фишера расчетный гипотеза об адекватности модели не отвергается.

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

Планы Хартли и Вестлейка Ортогональные и ротатабельные центральные композиционные планы второго порядка строятся на основе матриц ПФЭ 2n (для n <5), ДФЭ 2n-1 (для n =5) и ДФЭ 2п-2 для (n >8), обеспечивающих взаимно независимое (несмешанное) определение всех линейных и парных коэффициентов. Планы Хартли для n = 2, 3, 4, 5 могут быть построены на базе полуреплик соответствующих размерностей; для n = 6 в качестве линейной основы плана может использоваться 1/4-реплика ДФЭ 26-2. Очевидно, что разрешающая способность планов Хартли и Вестлейка должна быть ниже, чем для планов ортогональных или ротатабельных, однако в ряде практических задач этой разрешающей способности может оказаться достаточно для получения удовлетворительной модели исследуемой системы. Планы Хартли можно сделать почти ортогональными, если в них, как и в ортогональных центральных композиционных планах, сделать замену переменных согласно уравнению:

звездное плечо определяется по формуле: звездное плечо определяется по формуле: N0 – число точек в исходном ПФЭ. r – число повторений в звездных точках и в центре плана. N = N0+(2n+1)r – общее число опытов композиционного плана. Планы Вестлейка дают достаточно хороший результат, если получается удовлетворительная модель системы. Но все коэффициенты оказываются взаимно коррелируемыми. Линейная основа плана Хартли для п = 4 Это полуреплика 24-1, кото­рую можно использовать для построения композиционного плана

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

Параметры D-оптимальных планов

Пример расчета для плана второго порядка Для алюминиевого сплава, легированного цин­ком (10,0-10,5%), магнием (2,0-2,5%), медью (2,5-3,0%), марганцем (0,4%), титаном (0,05%) и цирконием (0,13%), прове­ряли оптимальность (с точки зрения прочности и пластичности при статическом растяжении) режима термообработки, применяемого обычно для алюминиевых высокопрочных сплавов системы Аl—Zn—Mg—Сu: закалка с 465° С в воде и старение при 140° С, 16 часов. Расчет коэффициентов регрессии проводили по общему уравнению В=(XТ X) -1 XТ Y

Матрица исходных данных

Результаты На поле диаграммы выделены области, соответствующие принятому режиму старения (Л) и режимам, при которых достигаются максимальные значения прочности. Как видно, при режимах старения, соответствующих областям А' и А", сплав имеет не только более высокую, чем при старении 140°, 16 час, прочность, но и более высокую пластичность.

Выводы Таким образом, проверяемый режим термообработки не является оптимальным для исследуемого сплава. На основании полученных результатов в качестве оптимального был рекомендован режим старения 120° С, 20 час. Для случаев применения сплава, когда требования к его пластичности не являются определяющими, был рекомендован режим старения 130° С, 4 час, обеспечивающий достаточно высокий уровень прочности и пластичности при минимальной продолжительности процесса старения.