化学反应工程四釜串连数据处理matlab程序Word格式.docx
《化学反应工程四釜串连数据处理matlab程序Word格式.docx》由会员分享,可在线阅读,更多相关《化学反应工程四釜串连数据处理matlab程序Word格式.docx(23页珍藏版)》请在冰豆网上搜索。
385,:
)-0.45)/A;
pt=sum(x(1:
).*ET)/sum(ET)%平均停留时间
sigmaf=sum(x(1:
).^2.*ET)/sum(ET)-pt^2
sig=sigmaf/pt^2n=1/sig
A=288.4449pt=344.33458sigmaf=23259.5671sig=.19617n=5.0975
面积:
A=288.4449平均停留时间:
tm=344.33458方差:
无因次方差:
小于1符合实际流型:
相当于釜数N=5.0975
检验归一性
jifen=int('
(0.534*exp(-((x-352.8)/222.3)^2)+0.4597*exp(-((x-218.4)/127.8)^2)+0.483*exp(-((x-710.6)/999.4)^2)-.45)/288.4449'
0,846.72)
jifen=1.0000003063063450964052571949284
计算Peclet准数:
取其中一个釜进行分析为开放体系:
解方程:
per=solve('
8*x^2+2*x-0.19617=0'
)得两组解:
per=[-.32536529140547271941324376838209](舍去)[.75365291405472719413243768382092e-1]
Peclet准数为:
0.75365
作E(t)-t关系图:
fplot('
[0846.72])
作F(t)-t关系图:
symsx
FT=int((0.534*exp(-((x-352.8)/222.3)^2)+0.4597*exp(-((x-218.4)/127.8)^2)+0.483*exp(-((x-710.6)/999.4)^2)-.45)/288.4449)
functiony=myfun(x)
y=217535077020598272/1057161738408951875*pi^(1/2)*erf(10/2223*x-392/247)+538299452482387968/5285808692044759375*pi^(1/2)*erf(5/639*x-364/213)+126367970891923456/151023105486993125*pi^(1/2)*erf(5/4997*x-187/263)-13194139533312/8457293907271615*x%将积分结果,定义为函数myfun
myfun'
[0846.72])第四个反应釜F(t)-t曲线:
2.将第三个实验数据V(t)-t曲线拟合:
用高斯四阶拟合求曲线表达式:
f=0.7926*exp(-((x-212.2)/126.5)^2)+0.7935*exp(-((x-304.4)/256)^2+0.5084*exp(-((x-127)/72)^2)+0.5974*exp(-((x-867.4)/825.4)^2)
0.7926*exp(-((x-212.2)/126.5)^2)+0.7935*exp(-((x-304.4)/256)^2)+0.5084*exp(-((x-127)/72)^2)+0.5974*exp(-((x-867.4)/825.4)^2)'
0,821.76)
A1=933.00202288276715439975963721132
A=A1-821.76*0.58=456.38122
ET=(y3(1:
374,:
)-0.58)/A;
sig=sigmaf/pt^2
n=1/sig
A=456.38122pt=282.00226sigmaf=21904.7178sig=.27544n=3.6305
A=456.38122平均停留时间:
tm=282.00226方差:
相当于釜数N=3.6305
(0.7926*exp(-((x-212.2)/126.5)^2)+0.7935*exp(-((x-304.4)/256)^2)+0.5084*exp(-((x-127)/72)^2)+0.5974*exp(-((x-867.4)/825.4)^2)-.58)/456.38122'
0,821.76)
jifen=1.0000000063165770808881216392106
8*x^2+2*x-0.4399=0'
)per=0.9873
[0821.76])
同理可得结果如下表:
平均停留时间tm
方差
无因次方差
计算釜数N
Peclet准数
反应釜一
135.24159
11935.77325
0.65257438
1.53239
0.1867639
反应釜二
213.2202
19611.27359
0.431369
2.3182
0.138716
反应釜三
282.00226
21904.7178
0.27544
3.6305
0.987297
反应釜四
344.33458
23259.5671
0.19617
5.0975
0.75365
数据个数494
阀开时间0.5sec
气体温度25℃
气体流量40m3/hr
气体压差10mmH2O
液体温度20℃
液体流量201/hr
操作员号A01-student
实验日期2006\04\24
实验时间16\32\17
Timer,mV1,mV2,mV3,mV4,No.
2.198,0.55,0.53,0.57,0.43,1
4.397,0.55,0.53,0.58,0.44,2
6.595,0.55,0.53,0.58,0.43,3
8.794,0.55,0.53,0.57,0.43,4
10.992,0.55,0.53,0.57,0.43,5
13.190,0.63,0.53,0.57,0.42,6
15.389,2.21,0.52,0.57,0.41,7
17.587,3.12,0.53,0.56,0.42,8
19.785,3.12,0.61,0.56,0.42,9
21.984,3.11,0.71,0.56,0.42,10
24.182,3.08,0.79,0.56,0.42,11
26.381,3.04,0.89,0.57,0.43,12
28.579,2.99,0.97,0.58,0.43,13
30.777,2.96,1.03,0.60,0.43,14
32.976,2.91,1.11,0.62,0.43,15
35.174,2.87,1.17,0.63,0.43,16
37.372,2.83,1.24,0.65,0.43,17
39.571,2.78,1.30,0.67,0.43,18
41.769,2.74,1.34,0.69,0.43,19
43.968,2.69,1.40,0.72,0.43,20
46.166,2.65,1.44,0.75,0.43,21
48.364,2.62,1.49,0.77,0.43,22
50.563,2.58,1.54,0.80,0.44,23
52.761,2.55,1.58,0.83,0.45,24
54.960,2.51,1.62,0.86,0.45,25
57.158,2.48,1.65,0.90,0.46,26
59.356,2.44,1.69,0.93,0.47,27
61.555,2.40,1.72,0.96,0.48,28
63.753,2.37,1.74,1.00,0.48,29
65.951,2.33,1.77,1.03,0.49,30
68.150,2.29,1.80,1.07,0.50,31
70.348,2.27,1.82,1.10,0.51,32
72.547,2.23,1.84,1.13,0.52,33
74.745,2.21,1.85,1.17,0.53,34
76.943,2.17,1.87,1.20,0.54,35
79.142,2.15,1.89,1.23,0.55,36
81.340,2.13,1.90,1.25,0.56,37
83.538,2.10,1.91,1.28,0.57,38
85.737,2.07,1.92,1.32,0.58,39
87.935,2.04,1.93,1.34,0.59,40
90.134,2.02,1.94,1.38,0.61,41
92.332,2.00,1.95,1.40,0.62,42
94.530,1.97,1.95,1.43,0.63,43
96.729,1.94,1.96,1.46,0.64,44
98.927,1.91,1.96,1.48,0.66,45
101.126,1.89,1.96,1.51,0.67,46
103.324,1.87,1.96,1.53,0.68,47
105.522,1.84,1.96,1.55,0.69,48
107.721,1.82,1.96,1.59,0.71,49
109.919,1.79,1.96,1.60,0.72,50
112.117,1.77,1.96,1.63,0.74,51
114.316,1.75,1.96,1.65,0.75,52
116.514,1.72,1.95,1.67,0.76,53
118.713,1.70,1.95,1.69,0.78,54
120.911,1.69,1.95,1.71,0.79,55
123.109,1.67,1.94,1.73,0.80,56
125.308,1.65,1.93,1.74,0.81,57
127.506,1.63,1.92,1.75,0.82,58
129.704,1.60,1.91,1.77,0.83,59
131.903,1.59,1.91,1.79,0.85,60
134.101,1.57,1.90,1.81,0.87,61
136.300,1.56,1.90,1.82,0.88,62
138.498,1.54,1.89,1.84,0.89,63
140.696,1.53,1.88,1.85,0.90,64
142.895,1.51,1.87,1.86,0.92,65
145.093,1.50,1.87,1.88,0.93,66
147.291,1.48,1.85,1.89,0.94,67
149.490,1.46,1.84,1.89,0.95,68
151.688,1.44,1.83,1.90,0.96,69
153.887,1.41,1.81,1.91,0.97,70
156.085,1.40,1.80,1.91,0.97,71
158.283,1.39,1.80,1.92,0.99,72
160.482,1.39,1.79,1.93,1.00,73
162.680,1.38,1.78,1.94,1.02,74
164.879,1.37,1.77,1.94,1.03,75
167.077,1.35,1.76,1.95,1.04,76
169.275,1.34,1.75,1.96,1.05,77
171.474,1.33,1.74,1.96,1.06,78
173.672,1.32,1.73,1.96,1.07,79
175.870,1.30,1.72,1.97,1.08,80
178.069,1.29,1.71,1.97,1.09,81
180.267,1.28,1.70,1.97,1.10,82
182.466,1.27,1.68,1.97,1.10,83
184.664,1.26,1.67,1.97,1.11,84
186.862,1.25,1.66,1.97,1.12,85
189.061,1.23,1.65,1.97,1.12,86
191.259,1.22,1.64,1.97,1.13,87
193.457,1.21,1.62,1.97,1.14,88
195.656,1.21,1.62,1.97,1.15,89
197.854,1.20,1.60,1.97,1.15,90
200.053,1.19,1.59,1.97,1.16,91
202.251,1.18,1.59,1.97,1.17,92
204.449,1.17,1.57,1.96,1.18,93
206.648,1.16,1.56,1.96,1.18,94
208.846,1.14,1.55,1.96,1.19,95
211.045,1.14,1.54,1.96,1.19,96
213.243,1.13,1.52,1.96,1.20,97
215.441,1.12,1.51,1.95,1.20,98
217.640,1.10,1.50,1.94,1.20,99
219.838,1.09,1.49,1.94,1.21,100
222.036,1.08,1.47,1.93,1.21,101
224.235,1.07,1.46,1.93,1.21,102
226.433,1.07,1.45,1.92,1.21,103
228.632,1.06,1.44,1.91,1.21,104
230.830,1.05,1.43,1.90,1.21,105
233.028,1.04,1.42,1.90,1.22,106
235.227,1.03,1.41,1.89,1.22,107
237.425,1.03,1.40,1.88,1.22,108
239.623,1.02,1.39,1.88,1.22,109
241.822,1.01,1.38,1.87,1.23,110
244.020,1.00,1.37,1.86,1.23,111
246.219,0.99,1.36,1.85,1.23,112
248.417,0.99,1.35,1.85,1.23,113
250.615,0.98,1.34,1.84,1.23,114
252.814,0.97,1.33,1.83,1.24,115
255.012,0.96,1.32,1.82,1.24,116
257.211,0.95,1.31,1.82,1.24,117
259.409,0.95,1.30,1.81,1.24,118
261.607,0.94,1.28,1.80,1.24,119
263.806,0.93,1.28,1.79,1.24,120
266.004,0.92,1.27,1.78,1.23,121
268.202,0.92,1.25,1.78,1.23,122
270.401,0.91,1.24,1.77,1.23,123
272.599,0.90,1.23,1.77,1.24,124
274.798,0.89,1.22,1.76,1.24,125
276.996,0.89,1.21,1.75,1.24,126
279.194,0.88,1.20,1.74,1.24,127
281.393,0.88,1.20,1.73,1.24,128
283.591,0.87,1.19,1.73,1.24,129
285.789,0.87,1.18,1.71,1.24,130
287.988,0.87,1.17,1.70,1.24,131
290.186,0.86,1.17,1.69,1.23,132
292.385,0.86,1.16,1.68,1.23,133
294.583,0.86,1.15,1.67,1.22,134
296.781,0.85,1.14,1.67,1.22,135
298.980,0.85,1.13,1.65,1.22,136
301.178,0.84,1.12,1.64,1.22,137
303.377,0.84,1.12,1.63,1.21,138
305.575,0.83,1.11,1.62,1.21,139
307.773,0.83,1.10,1.62,1.21,140
309.972,0.82,1.09,1.60,1.20,141
312.170,0.81,1.08,1.59,1.20,142
314.368,0.81,1.07,1.58,1.20,143
316.567,0.81,1.06,1.57,1.20,144
318.765,0.80,1.05,1.56,1.20,145
320.964,0.80,1.04,1.56,1.19,146
323.162,0.79,1.03,1.55,1.19,147
325.360,0.78,1.02,1.54,1.18,148
327.559,0.78,1.01,1.53,1.18,149
329.757,0.78,1.00,1.52,1.17,150
331.955,0.77,1.00,1.50,1.16,151
334.154,0.77,0.99,1.49,1.16,152
336.352,0.77,0.99,1.48,1.15,153
338.551,0.77,0.98,1.47,1.15,154
340.749,0.76,0.98,1.46,1.14,155
342.947,0.76,0.97,1.45,1.14,156
345.146,0.76,0.96,1.44,1.13,157
347.344,0.76,0.96,1.43,1.13,158
349.543,0.75,0.95,1.42,1.13,159
351.741,0.75,0.95,1.42,1.12,160
353.939,0.75,0.94,1.40,1.12,161
356.138,0.74,0.93,1.40,1.11,162
358.336,0.74,0.92,1.39,1.11,163
360.534,0.74,0.92,1.38,1.11,164
362.733,0.73,0.91,1.37,1.10,165
364.931,0.72,0.90,1.37,1.10,166
367.130,0.72,0.90,1.36,1.10,167
369.328,0.72,0.89,1.34,1.09,168
371.526,0.72,0.89,1.33,1.09,169
373.725,0.72,0.89,1.33,1.08,170
375.923,0.72,0.89,1.32,1.07,171
378.121,0.72,0.88,1.31,1.07,172
380.320,0.71,0.87,1.29,1.06,173
382.518,0.71,0.87,1.28,1.06,174
384.717,0.71,0.86,1.27,1.05,175
386.915,0.71,0.85,1.27,1.04,176
389.113,0.70,0.84,1.26,1.04,177
391.312,0.70,0.84,1.26,1.04,178
393.510,0.69,0.84,1.25,1.04,179
395.709,0.69,0.83,1.24,1.03,180
397.907,0.69,0.83,1.23,1.03,181
400.105,0.69,0.83,1.22,1.02,182
402.304,0.69,0.83,1.22,1.01,183
404.502,0.69,0.82,1.21,1.01,184
406.700,0.69,0.82,1.20,1.00,185
408.899,0.69,0.81,1.19,1.00,186
411.097,0.69,0.81,1.18,0.99,187
413.296,0.69,0.80,1.17,0.99,188
415.494,0.69,0.80,1.16,0.98,189
417.692,0.68,0.79,1.16,0.98,190
419.891,0.68,0.79,1.15,0.97,191
422.089,0.68,0.78,1.14,0.97,192
424.287,0.68,0.78,1.13,0.96,193
426.486,0.67,0.78,1.12,0.96,194
428.684,0.67,0.77,1.12,0.95,195
430.883,0.67,0.77,1.11,0.95,196
433.081,0.67,0.76,1.11,0.94,197
435.279,0.67,0.76,1.10,0.94,198
437.478,0.67,0.76,1.09,0.93,199
439.676,0.67,0.76,1