- •Институт холода и биотехнологий
- •Основы научных исследований, организации и планирования эксперимента Учебно-методическое пособие
- •Санкт-Петербург
- •Введение
- •1. Методология математического моделирования технологических машин и оборудования пищевых производств
- •1.1. Основные понятия планирования эксперимента (активный эксперимент)
- •Полный факторный эксперимент для построения модели зависимости качества измельчения от конструктивных параметров рабочих органов
- •Планирование эксперимента для построения модели зависимости качества измельчения от конструктивных параметров рабочих органов
- •1.2. Особенности обработки данных при пассивном эксперименте
- •Системы нормальных уравнений для различных форм связи
- •Замена переменных для линеаризации моделей
- •Данные пассивного эксперимента о выходе продукции
- •Корреляционная матрица
- •Результаты математико-статистической обработки
- •2. Решение задачи оптимизации технологических машин и оборудования пищевых производств методами линейного программирования
- •2.1. Алгоритм решения простых задач условной оптимизации
- •2.2. Решение задачи линейной оптимизации средствами Excel и Mathcad
- •Условия задачи линейного программирования
- •3. Варианты домашних заданий
- •Экспериментальные данные по испытаниям трубчатого теплообменника
- •Список литературы
- •Содержание
- •Институт холода и биотехнологий
- •Учебно-методическое пособие
Полный факторный эксперимент для построения модели зависимости качества измельчения от конструктивных параметров рабочих органов
Номера опытов
|
Номер пункта опыта
|
Значения нормированных факторов |
Средние значения функций отклика Yi |
|||||
X1 |
X2 |
X1X2 |
Y1 |
Y2 |
Y3 |
Y4 |
||
1, 5, 9 |
1 |
1 |
1 |
1 |
1560 |
1600 |
1580 |
1560 |
2, 6, 10 |
2 |
–1 |
1 |
–1 |
1465 |
1495 |
1480 |
1495 |
3, 7, 11 |
3 |
1 |
–1 |
–1 |
1470 |
1488 |
1452 |
1452 |
4, 8, 12 |
4 |
–1 |
–1 |
1 |
1390 |
1410 |
1400 |
1390 |
Для получения статистически достоверной математической модели при анализе экспериментальных данных проверялась однородность дисперсий выборок функций отклика Y1,Y2,Y3 и Y4. Для этой цели вычислялось значение критерия Фишера по формуле
(1.1)
Для данных табл. 1.4 это значение оказалось равным 4,0, что меньше табличного значения, равного 9,28, для доверительной вероятности 95 %.
Такое соотношение между расчетным и табличным значениями критерия Фишера свидетельствует о воспроизводимости эксперимента, что позволяет вычислить коэффициенты модели по формулам
а0 = ; аj = , (1.2)
где Xj – столбцы элементов матрицы планирования; – среднеариф-метическое значение для каждой из функций отклика по трем параллельным опытам; М – число разных экспериментов (М = 4).
Вычисления по приведенным формулам для выполненного эксперимента с помощью пакета прикладных программ Eхсеl дали следующие оценки коэффициентов уравнений регрессии:
для y1: a0 = 1471,25; a1 = 41,25; a2 = 43,75;
для y2: a0 = 1498,25; a1 = –49,25; a2 = 45,75.
для y3: a0 = 1478; a1 = –52; a2 = 38;
для y4: a0 = 1474,25; a1 = –53,25; a2 = 31,75.
В соответствии с проведенными вычислениями оценки уравнения регрессии записываются в следующем виде:
y1 = 1471,25 – 41,25 1 + 43,75 2; |
y2 = 1498,25 – 49,25 1 + 45,75 2; |
y3 = 1478 – 52 1 + 38 2; |
y4 = 1474,25 – 53,25 1 + 31,75 2.
|
Уравнения регрессии, построенные по этим формулам, только в том случае соответствуют реальному процессу, когда каждый их член вносит вклад, значимо отличающийся от случайных колебаний функций отклика. Это условие выполняется, если абсолютная величина коэффициента больше его доверительного интервала, определяемого при помощи критерия Стьюдента со степенью свободы
f = M (k – 1), ai 2 a,
так как
где k – число параллельных опытов.
Вычисления интервалов достоверности коэффициентов регрессии
свидетельствуют о том, что в уравнениях регрессии значимы лишь те члены, для которых выполняются условия:
для – 0,3637;
для – 26,3373;
для y3 – 16,3373;
для y4 – 30,3373.
Последним шагом, предшествующим использованию полученных экспериментально соотношений, является проверка их адекватности. Эта проверка позволяет судить о том, не отброшены ли в процессе обработки результатов величины, существенные для достоверного воспроизведения полученной зависимостью реального процесса, и правильно ли выбрана искомая математическая модель. Она заключается в вычислении расчетного значения критерия Фишера и его сравнении с табличным. Расчетное значение определяется соотношением
при
где ур – рассчитанное по уравнениям значение функций отклика у1 и у2; – дисперсия среднего значения функции отклика; k – число коэффициентов в уравнении.
Пользуясь данными табл. 3.8–3.9 [2] и результатами произведенных ранее вычислений для расчетного критерия Фишера, находим Fp = 14,2222. Для доверительной вероятности 95 % и соответствующих степеней свободы числителя и знаменателя табличное значение критерия Фишера оказалось равным 12,6744.
Соотношение Fp > F свидетельствует о том, что выбранная на первом этапе линейная математическая модель неадекватна реальной зависимости эксплуатационных характеристик рабочих органов измельчительной машины от варьируемых в эксперименте факторов. Для нахождения зависимости, адекватной реальной, достроим план эксперимента до центрального композиционного ротатабельного, добавляя так называемые «звездные точки» для значений нормированных факторов Хi = 2 и проводя дополнительные эксперименты в центре плана. Соответствующая матрица планирования и некоторые вспомогательные величины приведены в табл. 1.5.
Таблица 1.5