Российский Химико – Технологический Университет

им.Д.И.Менделеева.

Практическая работа №4 по курсу

«Компьютерное моделирование химических процессов»

«Обработка результатов активного эксперимента».

Группа: Ф-31

Фамилия И.О.: Можевитина Е.Н.

Вариант № 7

Преподаватель: Новикова Д. К.

Москва 2007.

1.Рассмотрим модель, в которой зависимость давления от температуры выражается формулой:

P=exp(A+B/T).

Преобразуем уравнение к виду: lnP=A+B/T или

Y=a₀+a₁х

lnp1

lnp2

lnp3

lnp4

lnp5

lnp6

lnp7

lnp8

lnp9

lnp10

1,813519

1,445058

1,302929

0,664851

0,10589

-0,39802

-0,81385

-1,43666

-1,95327

-2,81308

y(эксп) =

1	1/T1
1	1/T2
1	1/T3
1	1/T4
1	1/T5
1	1/T6
1	1/T7
1	1/T8
1	1/T9
1	1/T10

1	0,05
1	0,030303
1	0,027027
1	0,019231
1	0,015873
1	0,013889
1	0,012658
1	0,011236
1	0,010309
1	0,009091

Ф =

10	0,199617
0,199617	0,005439

Ф^ТФ =

0,374029	-13,7277
-13,7277	687,7036

(Ф^ТФ)^-1 =

-2,24047

101,8054

А =

В =

S ²_R = 0,829584 F^Расч = 1,400314

S ²_Ср = 2,411955

P=exp(-2,24047+101,8054/T).

2. Рассмотрим модель, в которой зависимость давления от температуры выражается формулой:

P=exp(A+B/(С+T)).

Преобразуем уравнение к виду: ТlnP=(AС+В)+АТ-С lnР или

Y=a₀+a₁T-a₂х

Т1lnp1

Т2lnp2

Т3lnp3

Т4lnp4

Т5lnp5

Т6lnp6

Т7lnp7

Т8lnp8

Т9lnp9

Т10lnp10

36,27039

47,68693

48,20838

34,57223

6,671094

-28,6573

-64,2939

-127,863

-189,467

-309,439

y(эксп) =

1	Т1	-lnP1
1	Т2	-lnP2
1	Т3	-lnP3
1	Т4	-lnP4
1	Т5	-lnP5
1	Т6	-lnP6
1	Т7	-lnP7
1	Т8	-lnP8
1	Т9	-lnP9
1	Т10	-lnP10

1	20	-1,81352
1	33	-1,44506
1	37	-1,30293
1	52	-0,66485
1	63	-0,10589
1	72	0,398018
1	79	0,813847
1	89	1,436662
1	97	1,953267
1	110	2,813077

Ф =

10	652	2,082623
652	50386	546,3105
2,082623	546,3105	22,14133

Ф^ТФ =

33,61636	-0,54707	10,33635
-0,54707	0,00893	-0,16888
10,33635	-0,16888	3,239877

(Ф^ТФ)^-1 =

АС+В= -669,966

А= 10,29427

- С= -268,175 отсюда получаем : А= 10,29427, С= 268,175 , В= - 948,43527

S ²_R = 163,0528 F^Расч = 88,76584

S ²_Ср = 14473,51

P=exp(10,29427+- 948,43527/(268,175 +T)).

3. Рассмотрим модель, в которой зависимость давления от температуры выражается формулой:

P=exp(A+B/T+C*T+D*lnT)

Преобразуем уравнение к виду:

lnP=A+B/T+C*T+D*lnT

1,813519

1,445058

1,302929

0,664851

0,10589

-0,39802

-0,81385

-1,43666

-1,95327

-2,81308

lnP1

lnP2

lnP3

lnP4

lnP5

lnP6

lnP7

lnP8

lnP9

lnP10

y(эксп) =

1	1/T1	T1	lnT1
1	1/T2	T2	lnT2
1	1/T3	T3	lnT3
1	1/T4	T4	lnT4
1	1/T5	T5	lnT5
1	1/T6	T6	lnT6
1	1/T7	T7	lnT7
1	1/T8	T8	lnT8
1	1/T9	T9	lnT9
1	1/T10	T10	lnT10

1	0,05	20	2,995732
1	0,030303	33	3,496508
1	0,027027	37	3,610918
1	0,019231	52	3,951244
1	0,015873	63	4,143135
1	0,013889	72	4,276666
1	0,012658	79	4,369448
1	0,011236	89	4,488636
1	0,010309	97	4,574711
1	0,009091	110	4,70048

Ф =

10	0,199617	652	40,60748
0,199617	0,005439	10	0,750119
652	10	50386	2788,78
40,60748	0,750119	2788,78	167,5689

Ф^ТФ =

4027,297	-22510,3	10,0644	-1042,68
-22510,3	127713	-55,0707	5799,801
10,0644	-55,0707	0,026143	-2,6275
-1042,68	5799,801	-2,6275	270,4469

(Ф^ТФ)^-1 =

A= -6,58736

B= 31,56082

C= -0,09186

D= 2,890692

S ²_R = 1,24394E-05 F^Расч = 193896

S ²_Ср = 2,411955

P=exp(-6,58736+31,56082/T+-0,09186*T+2,890692*lnT)