Автор работы: Пользователь скрыл имя, 02 Ноября 2011 в 16:55, курсовая работа
Моделирование помогает понять и упорядочить результаты эмпирических наблюдений, создать логический каркас научной теории, обнаружить внутренние связи и соотношения между результатами эксперимента. Особую важность задачи моделирования приобретают при рассмотрении сложных объектов, априорные сведения о которых либо отсутствуют, либо незначительны. Многие свойства, которые остаются неучтенными при моделировании, например, из-за того, что они неизвестны, могут в корне изменить картину результатов моделирования, и модель, таким образом, оказывается очень далека от оригинала.
ВВЕДЕНИЕ……………………………………………………………………... 4
1. Осушка газа абсорбционным методом ……..…………..…... 5
2. ОБРАБОТКА Результатов активного эксперимента...........… 6
2.1. Планирование эксперимента ……............................................................ 6
2.2. Определение основных статистических характеристик параллельных опытов……………………………..……………………..…….
8
2.4. Проверка результатов измерений по критерию грубой ошибки…......... 9
2.5. Определение дисперсии воспроизводимости по критерию Кохрена.... 10
3. Построение математической модели абсорбера………… 10
3.1. Определение порядка и расчет коэффициентов модели…….……….. 10
3.2. Проверка модели на адекватность……………………………………..
13
3.2.1. Критерий Фишера.…………………………………………….….. 13
3.2.2. Корреляционная функция остатков………………………………. 15
ЗАКЛЮЧЕНИЕ…………………………………………………………………... 16
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ……………………………. 17
ПРИЛОЖЕНИЕ А. Текст m-файла определение порядка и расчет коэффициентов модели и проверки их на адекватность………………………
18
ПРИЛОЖЕНИЕ Б. Текст m-файла Корреляционная функция остатков……. 19
ПРИЛОЖЕНИЕ В. Критерий грубой ошибки………………………………… 20
ПРИЛОЖЕНИЕ Г. Квантиль распределения Кохрена
Для каждой группы параллельных опытов определяются следующие статистические характеристики:
Максимальное значение - ;
Минимальное значение - ;
Среднее значение - , где число опытов в данной группе (объём выборки);
Дисперсия ;
Среднее квадратичное отклонение .
Результаты вычислений статистических характеристик сведены в таблицу 2.5.
Таблица 2.5. Статистические характеристики
Эксперимент №1 | Эксперимент №2 | Эксперимент №3 | |
47406 | 82678 | 82392 | |
34373 | 58370 | 57371 | |
41445,2 | 72771,5 | 71715,2 | |
26804725 | 70322237 | 47555786 | |
5177,328 | 8385,836 | 6896,07 |
2.3 Проверка результатов измерений по критерию грубой ошибки
Для
оценки выборочных данных по критерию
наличия грубой ошибки (R критерий) для
каждой выборки, полученной в результате
проведения параллельных опытов, вычисляются
величины:
Расчетные значение и сравниваются с . (приложение В). Табличное значение выбирают для уровня значимости α и числа степеней свободы , где - объем выборки.. Если
принимают, что отклонение (или ) определяется случайными явлениями и принадлежит данной генеральной совокупности, в противном случае результат отбрасывается как грубый промах и оценку по R-критерию повторяют для следующего (или ) с пересчитанными значениями основных статистических характеристик.
Эксперимента №1:
9
α=10
1,1513 1,366 =2,26
<
<
Эксперимента №2:
9
α=5
1,1813 1,717 =2,26
< <
Эксперимента №3:
9
α=5
1,577 2,0281 =2,26
< <
Во всех трех выборках Rmax < Rкр и . Принимаем, что отклонение Ymax определяется случайными явлениями и принадлежит данной генеральной совокупности.
2.4 Определение дисперсии воспроизводимости по критерию Кохрена
Дисперсию воспроизводимости ( ) определяют путём сравнения выборочных дисперсий для параллельных опытов. При одинаковом числе опытов во всех выборках для сравнения дисперсий пользуются критерием Кохрена.
Рассчитывается значение
и сравнивается с (приложение Г), выбранным для числа сравниваемых дисперсий и числа степеней свободы , с которым определена каждая дисперсия. Если , сравниваемые дисперсии можно считать однородными, следовательно их можно усреднить. А значит дисперсия воспроизводимости
G=0,6452, Gкр=0,6167
G > Gкр , различие между дисперсиями значительное, поэтому в качестве дисперсии воспроизводимости выбираем меньшую из сравниваемых дисперсий.
S²воспр=26804725
3. Построение математической модели абсорбера
3.1. Определение порядка и расчет коэффициентов модели
Процесс создания модели начинается с выбора типа модели и, как правило, на первом этапе останавливаются на линейном варианте в форме алгебраического многочлена:
где , - неизвестные коэффициенты модели,
G, L - варьируемые входные параметры объёкта.
Поиск
неизвестных коэффициентов
Таблица 3.1. Результаты активного эксперимента
№ | Расход газа
G, m3/ч |
Температура газа T, град. |
Расход абсорбента L, м3/ч |
Концентрация абсорбента X, кг/м3 |
Расход газа Y, м3/ч |
1 | 10000 | 11 | 30 | 47 | 19497 |
2 | 15000 | 11 | 40 | 47 | 21321 |
3 | 20000 | 11 | 50 | 47 | 26611 |
4 | 25000 | 11 | 60 | 47 | 22861 |
5 | 30000 | 11 | 70 | 47 | 26508 |
6 | 10000 | 11 | 30 | 47 | 20782 |
7 | 15000 | 11 | 40 | 47 | 18989 |
8 | 20000 | 11 | 50 | 47 | 19353 |
9 | 25000 | 11 | 60 | 47 | 25396 |
10 | 30000 | 11 | 70 | 47 | 20834 |
11 | 10000 | 11 | 30 | 47 | 26910 |
12 | 15000 | 11 | 40 | 47 | 22512 |
13 | 20000 | 11 | 50 | 47 | 19488 |
14 | 25000 | 11 | 60 | 47 | 24308 |
15 | 30000 | 11 | 70 | 47 | 20448 |
16 | 10000 | 11 | 30 | 47 | 22246 |
17 | 15000 | 11 | 40 | 47 | 19814 |
18 | 20000 | 11 | 50 | 47 | 20238 |
19 | 25000 | 11 | 60 | 47 | 19428 |
20 | 30000 | 11 | 70 | 47 | 26886 |
21 | 10000 | 11 | 30 | 47 | 19326 |
22 | 15000 | 11 | 40 | 47 | 26451 |
23 | 20000 | 11 | 50 | 47 | 22847 |
24 | 25000 | 11 | 60 | 47 | 22163 |
25 | 30000 | 11 | 70 | 47 | 24934 |
Рисунок
2. Поверхность расхода газа Y
при изменении расхода газа G и расхода
абсорбента L.
После вычисления неизвестных коэффициентов рассчитывают остаточную дисперсию:
где - экспериментальное значение выходного параметра для определенных входных сигналов;
- величина выходного параметра, расcчитанного по модели при тех же значениях входных сигналов;
- число коэффициентов в
- число экспериментальных
Таблица
3.2. Поиск неизвестных
Вид модели | Коэффициенты | S2ост | ||||||
b1 | b2 | b3 | b4 | b5 | b6 | |||
b1G+b2L | -0.0039 | 2.0265 | - | - | - | - | 7.9186e+006 | |
b1G+b2L+b3 | -0.0035 | 1.8132 | 2.8059 | - | - | - | 7.6359e+006 | |
b1G^2+b2G+b3L+b4 | 0 | -0.0043 | 2.0505 | 2.6245 | - | - | 7.9314e+006 | |
b1G^3+b2G^2+b3G+b4L+b5 | 0 | 0 | 0.0003 | 0.8159 | 3.2729 | - | 8.1859e+006 | |
b1G^4+b2G^3+b3G^2+b4G+b5L+b6 | 0 | 0 | 0.3826 | 0 | 0 | 0 | 1.8461e+007 |
3.2. Проверка модели на адекватность
3.2.1. Критерий Фишера
Значения коэффициентов в уравнении регрессии, полученные по методу наименьших квадратов, являются оптимальными для выбранной математической модели, однако не всегда корректно останавливаться на этом. Процесс создания модели должен заканчиваться объективной оценкой, насколько точно построенная модель описывает идентифицируемый объект.
Проверка
модели на адекватность производится
путём сравнения суммы
Если ., то с достоверностью в (1 - a)*100% считают модель адекватной объекту, если нет - с той же достоверностью вероятно противоположное утверждение.
В случае, когда адекватность модели не подтвердилась, необходимо вернуться к началу п.2.2 и изменить вид модели. Чаще всего в такой ситуации просто увеличивают порядок модели и весь последующий расчет повторяют. Данные сведены в таблицу 3.3.
Модель адекватна.
Таблица 3.3. Сравнение F
Вид модели | S2ост | F | Fкр | |
b1G+b2L | 7.9186e+006 | 3.3850 | > | 2.9169 |
b1G+b2L+b3 | 7.6359e+006 | 3.5104 | > | 2.9263 |
b1G^2+b2G+b3L+b4 | 7.9314e+006 | 3.3796 | > | 2.9365 |
b1G^3+b2G^2+b3G+b4L+b5 | 8.1859e+006 | 3.2745 | > | 2.9477 |
b1G^4+b2G^2+b3G^2+b4G+b5L+b6 | 1.8461e+007 | 1.4519 | < | 2.9600 |