Статистическая обработка результатов эксперимента
Министерство образования и Северный (Арктический) федеральный университет имени М.В. Ломоносова Кафедра лесохимических производств | |||||||||||||||||||
Куликова Оксана Николаевна | |||||||||||||||||||
(фамилия, имя, отчество | |||||||||||||||||||
Институт |
ТиПХ |
курс |
II |
группа |
508 |
||||||||||||||
|
|
|||||||||||||||||||
КУРСОВАЯ РАБОТА |
|||||||||||||||||||
По дисциплине «МАТЕМАТИЧЕСКИЕ МЕТОДЫ И МОДЕЛИ» |
|||||||||||||||||||
На темы «Статистическая обработка результатов эксперимента» и «Типовые модели структуры потоков» |
|||||||||||||||||||
|
|
|||||||||||||||||||
Работа допущена к защите |
|||||||||||||||||||
(подпись руководителя) |
(дата) |
||||||||||||||||||
Признать, что работа |
|||||||||||||||||||
выполнена и защищена с оценкой |
|||||||||||||||||||
|
|
|||||||||||||||||||
Руководитель |
С.А. Цаплина | ||||||||||||||||||
(должность) |
(подпись) |
(И. О. Фамилия) | |||||||||||||||||
|
|
|||||||||||||||||||
(дата) |
|||||||||||||||||||
Архангельск |
|||||||||||||||||||
2012 |
|||||||||||||||||||
ЛИСТ ЗАМЕЧАНИЙ
______________________________
СОДЕРЖАНИЕ
ВВЕДЕНИЕ 4
1 АНАЛИЗ РЕЗУЛЬТАТОВ ЭКСПЕРИМЕНТА 5
1.1 Оценка надежности аналитической методики 7
1.2 Дисперсионный анализ результатов 11
1.3 Аппроксимация результатов эксперимента 14
2 ОПИСАНИЕ МНОГОФАКТОРНОЙ СИСТЕМЫ 22
2.1 Расчет линейного уравнения связи (полинома I степени) 23
2.2 Расчет полного квадратного уравнения (полинома II степени) 26
3 РАСЧЕТ ТЕХНОЛОГИЧЕСКОГО АППАРАТА 28
3.1 Определение типа химического реактора 28
3.2 Определение объема химического реактора 31
ВЫВОДЫ 33
СПИСОК ИСПОЛЬЗУЕМОЙ ЛИТЕРАТУРЫ 34
ВВЕДЕНИЕ
Целью курсовой работы является выполнение расчетов, связанных с оценкой факторов, определяющих течение технологического процесса, получение его математического описания, выполнение статистического анализа имеющейся информации, определение параметров технологического процесса и промышленного аппарата с использованием типовых моделей структуры потоков.
Перед студентом впервые за время обучения ставится комбинированная задача, решить ее он должен самостоятельно, используя комплексно теоретические знания, полученные по химии, физике, высшей математике, инженерной графике, вычислительной технике.
Основная часть курсовой работы разбита на 3 раздела и включает 7 расчетных заданий. Выполняются следующие расчеты:
- оценка надежности аналитической методики по данным опыта;
- дисперсионный анализ результатов опытов;
- аппроксимация результатов эксперимента;
- расчет коэффициентов линейного уравнения (полинома I степени);
- расчет коэффициентов полного квадратного уравнения (полином II степени);
- определение типа химического реактора по С – выходной кривой;
- расчет объема химического реактора.
1 АНАЛИЗ РЕЗУЛЬТАТОВ ЭКСПЕРИМЕНТА
Эксперимент – это совокупность действий, направленных на установление взаимосвязи между входными и выходными параметрами системы. Входные параметры (воздействия, оказываемые на систему) иначе называют факторами, а выходные параметры (воздействия, оказываемые системой на окружающую среду) – откликами.
Приняты следующие обозначения: Х – входной параметр; Y – выходной параметр.
Основными этапами эксперимента являются:
- выбор объекта исследования;
- установление параметров (входных и выходных), значения которых должны быть определены в ходе эксперимента;
- выбор и оценка надежности аналитических методик для определения значений параметров;
- проведение дисперсионного анализа с целью установления существенности влияния выбранных факторов;
- проведение основного эксперимента;
- аппроксимация результатов эксперимента с целью получения уравнения связи, которое наилучшим образом будет их описывать;
- установление значений входных параметров, при которых получается оптимальные значения выходных параметров;
- графическая интерпретация полученных результатов с целью более наглядного их представления.
Эксперимент состоит из ряда опытов, в которых могут производиться параллельные определения. В опытах значения хотя бы одного фактора меняется, а параллельные определения осуществляются в строго идентичных условиях.
1.1 Оценка надежности аналитической методики
Аналитическая методика – это строгая последовательность действий, направленных на определение значений какого – либо параметра.
Аналитическая методикам считается надежной, если ее относительная максимальная погрешность не превышает некоторого установленного значения (обычно не более 5 %). При наличии нескольких методик для определения одного и того же параметра желательно выбрать ту, у которой погрешность наименьшая. Погрешность методики оценивается по результатам параллельных определений одного опыта.
Исходные данные для оценки аналитической методики:
X |
X = const | ||||||
Y |
18,0 |
17,7 |
18,3 |
17,1 |
18,9 |
17,3 |
18,7 |
- Определение среднего значения выходного параметра:
,
где m – число параллельных определений;
.
2) Определение выборочной дисперсии, которая характеризует меру отклонения (рассеивания) результатов параллельных определений от их среднего значения:
,
где fu = m – 1 – число степеней свободы выборочной дисперсии; fu = 7 – 1 = 6;
.
3) Определение средней квадратичной погрешности отдельного или единичного результата:
.
4) Проверка результатов на анормальность (на наличие промахов).
Анормальный результат (иначе промах, грубая ошибка) – это резко отклоняющийся результат из серии параллельных определений, полученный в результате грубой ошибки со стороны исследователя.
Обнаружение анормальных результатов проводится 2 способами:
а) с помощью критерия промаха К (грубый способ):
К = 3Su; Yср – К ≤ Y ≤ Yср + К;
К = 3 0,681 = 2,04;
18,00 – 2,04 ≤ Y ≤ 18,00 + 2,04;
15,96 ≤ Y ≤ 20,04
Вывод: анормальных результатов не обнаружено, так как все результаты параллельных определений попадают в данный интервал.
б) с помощью критерия анормальности Н (самый точный способ):
; ;
Нтабл = f(α;m);
; ;
.
Вывод: так как оба расчетных значения критерия анормальности (для минимального и максимального результатов) не превышают табличного, то анормальных результатов среди параллельных определений нет.
5) Определение средней квадратичной погрешности среднего арифметического результата:
;
.
6) Определение табличного значения критерия Стьюдента, который представляет собой нормированную погрешность:
,
где α – уровень значимости, показывающий допустимую долю (или процент) ошибок; в расчетах чаще всего принимают значение α=0,05;
.
7) Определение абсолютной максимальной погрешности опыта:
;
.
8) Определение относительной максимальной погрешности опыта, %:
;
.
Главный вывод: так как относительная максимальная погрешность опыта не превышает 5%, то аналитическую методику можно считать надежной и она может быть использована для определения параметра Y в последующем эксперименте.
9) Установление доверительного интервала, т.е. интервала в котором находится истинное значение параметра Y с вероятностью Р = 1 – α:
;
P=1 – 0,05= 0,95 (95%);
;
10) Установление стабильности параметра Y по коэффициенту вариации, %:
;
.
Если 5%, то параметр Y стабилен, т.е. во времени не изменяется. Процесс, у которого все параметры стабильны, называется стационарным.
Вывод: так как коэффициент вариации не превышает 5%, то параметр Y является стабильным, т.е. не изменяется во времени.
11) Установление необходимого числа параллельных определений для получения результатов с погрешностью, не превышающей 5%:
; ;
Вывод: в каждом опыте требуется производить не менее четырех параллельных определений.
1.2 Дисперсионный анализ результатов опытов
Проведение дисперсионного анализа позволяет выделить факторы, существенно влияющие на систему. Факторы, влияние которых на процесс несущественно, в дальнейшем из рассмотрения исключают.
Для дисперсионного анализа требуется провести несколько опытов (с параллельными определениями), изменяя значения исследуемого фактора.
Исходные данные для дисперсионного анализа результатов опытов:
Опыт |
Параллельные определения | ||||
Y1 |
Y2 |
Y3 |
Y4 |
Y5 | |
1 |
4,8 |
4,7 |
5,0 |
5,3 |
5,2 |
2 |
16,8 |
16,0 |
15,2 |
16,6 |
15,4 |
3 |
20,0 |
19,0 |
21,0 |
19,3 |
20,7 |
- Определение среднего значения параметра в каждом опыте:
,
где m – число параллельных определений в i – том опыте;
;
- Определение выборочной (построчной) дисперсии для каждого опыта – меры отклонения результатов параллельных определений в каждом из опытов от соответствующей им средней величины:
,
где – число степеней свободы дисперсии;
;
;
- Проверка однородности дисперсий и воспроизводимости опытов по критерию Кохрена:
; ;
;
.
Вывод: Gmax < Gтабл, следовательно, дисперсии однородны, а опыты воспроизводимы, т.е. выполнены с заданной степенью точности.
- Определение внутригрупповой дисперсии – средней меры отклонения всей совокупности результатов параллельных определений от соответствующих им средних значений в каждом из опытов:
,
где n – число опытов;
;
Число степеней свободы внутригрупповой дисперсии:
;
.
- Определение среднего значения параметра во всем эксперименте:
;
.
- Определение межгрупповой дисперсии – меры отклонения средних значений параметра в опытах от среднего значения этого параметра во всем эксперименте:
,
где – число степеней свободы межгрупповой дисперсии,
;
.
- Определение критерия Фишера:
; ,
где α – уровень значимости;
; .
Если , то фактор X существенно влияет на систему и, следовательно, должен учитываться при нахождении уравнения связи.
Вывод: так как F > Fтабл, то фактор Х существенно влияет на систему.
1.3 Аппроксимация результатов эксперимента
Аппроксимация проводится с целью установления уравнения связи между входными и выходными параметрами.
Аппроксимацию условно можно разделить на три части:
I часть – построение графика по опытным данным и установление по виду графика возможных видов уравнений связи для их описания;
II часть – определение коэффициентов предполагаемого уравнения связи (или нескольких уравнений);
III часть – оценка надежности полученного уравнения; при наличии нескольких уравнений оценивается надежность каждого уравнения и выбирается более точное.
Для определения коэффициентов уравнений связи можно использовать следующие методы:
Графический – используется в случае линейного уравнения или уравнения, приведенного к линейному виду (линеаризованного) – является самым грубым. Свободный коэффициент определяется как отрезок, отсекаемый прямой от оси ординат, а коэффициент при Х – как тангенс угла наклона прямой к положительной оси абсцисс.
Метод избранных точек – предполагает выбор двух (по числу коэффициентов уравнения) опытных точек, через которые проходит прямая, и решение системы из двух уравнений. Система уравнений составляется путем подстановки выбранных опытных значений X и Y в исходное уравнение.
Метод средних (метод уравновешенной погрешности) – предполагает использование всех опытных данных. Полученная при постановке опытных значений система из n уравнений делится на столько примерно равных частей, сколько коэффициентов. В каждой части уравнения почленно складываются.
Метод наименьших квадратов – основан на условии минимальности суммы квадратов разности опытных значений и соответствующих им расчетных значений.
Исходные данные для аппроксимации результатов эксперимента:
Х |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
|
|
Y |
3,2 |
3,9 |
5,8 |
9,5 |
15,6 |
24,8 |
37,3 |
I часть. Построение графика по опытным данным (рисунок 1).
Рисунок 1 – Зависимость Y=f(X)
Уравнение связи имеет вид .
II часть. Определение коэффициентов уравнения .
Так как уравнение нелинейное, проведем его линеаризацию путем замены переменной X: .
;
В результате получаем данные для определения коэффициентов уравнения:
Х* |
1 |
8 |
27 |
64 |
125 |
216 |
343 |
Y |
3,2 |
3,9 |
5,8 |
9,5 |
15,6 |
24,8 |
37,3 |
Метод средних. Используем все пары значений Y и Х*, полученную систему уравнений делим на 2 части, в каждом уравнении почленно складываем:
Полученное уравнение: .
Графический метод. Строим график зависимости Y = f(X*) (рисунок 2).
Рисунок 2 – Зависимость Y*=f(X)
По графику определяем: 3,2 (отрезок отсекаемый прямой от оси ординат):
;
.
Получаем уравнение: .
Метод избранных точек. Выберем первую и последнюю опытные точки и соответствующие пары значений X* и Y подставим в уравнение :
;
;
;
Получаем уравнение: .
Метод наименьших квадратов. Расчетная система уравнений имеет вид:
Найдём каждую сумму:
Полученные значения подставим в расчетную систему и решим ее:
;
;
Таким образом уравнение имеет вид: .
III часть. Оценка надежности уравнения (используем уравнение полученное самым точным методом).
Способ 1. Данный способ не используем, так как в опытах параллельные определения не проводились.
Способ 2.
- Определение среднего значения параметра в эксперименте: