Построение статической модели абсорбера

Автор работы: Пользователь скрыл имя, 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
ПРИЛОЖЕНИЕ Г. Квантиль распределения Кохрена

Содержимое работы - 1 файл

курсач моделирование1.doc

— 771.50 Кб (Скачать файл)

   2.2. Определение основных статистических характеристик параллельных опытов

   Для каждой группы параллельных опытов определяются следующие статистические характеристики:

   Максимальное  значение - ;

   Минимальное значение - ;

   Среднее значение - , где число опытов в данной группе (объём выборки);

   Дисперсия ;

   Среднее квадратичное отклонение .

   Результаты  вычислений статистических характеристик сведены в таблицу 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кр , различие между дисперсиями значительное, поэтому в качестве дисперсии воспроизводимости выбираем меньшую из сравниваемых дисперсий.

воспр=26804725 

   3. Построение математической модели абсорбера

   3.1. Определение порядка и расчет коэффициентов модели

   Процесс создания  модели  начинается  с выбора типа модели и, как правило,  на первом этапе останавливаются  на линейном варианте в форме алгебраического  многочлена:

                                          y=Gb1+Lb2

   где  , - неизвестные коэффициенты модели,

   G, L - варьируемые входные параметры объёкта.

   Поиск неизвестных коэффициентов осуществляют с помощью регрессии по методу наименьших квадратов на основе данных активного эксперимента (таблица 3.1; рисунок 2). Данные поиска сведены в таблице 3.2.

   Таблица 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. Критерий Фишера

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

   Проверка  модели на адекватность производится путём сравнения суммы квадратов  отклонений экспериментальных данных от рассчитанных по модели (остаточная дисперсия ) с показателем точности проводимых измерений (дисперсия   воспроизводимости ).  Подобная процедура известна в теории  регрессионного анализа под названием "критерий Фишера" (F -критерий). С помощью этого критерия производят сравнение двух дисперсий, рассчитывая отношение большей дисперсии к меньшей и сравнивая полученный результат с табличным значением. , выбранным для уровня значимости и чисел степеней свободы и для и соответственно ( , ). 

                                                      ,

   Если  ., то с достоверностью в (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

Информация о работе Построение статической модели абсорбера