由附件1“性能数据表”中A1~A14(用装料方式I)整理得下表“sheet2”:
表1 性能数据表
wt%Co/SiO2
Co/SiO2 和 HAP 的质量比
催化剂总量
乙醇浓度
温度
乙醇转化率(%)
C4烯烃选择性(%)
C4烯烃收率(%)
0.5
1
400
1.68
250
4.0
9.6
0.386665445
0.5
1
400
1.68
275
12.1
8.6
1.042826728
0.5
1
400
1.68
300
29.5
10.7
3.159887673
0.5
1
400
1.68
325
43.3
18.9
8.184209991
0.5
1
400
1.68
350
60.5
27.3
16.47764673
0.5
1
400
1.68
400
88.4
41.0
36.27781909
1
0.4925
100
1.68
250
2.5
1.9
0.047111815
1
0.4925
100
1.68
275
5.3
2.6
0.136038297
1
0.4925
100
1.68
300
10.2
3.6
0.366488509
... ...
... ...
(略)
代码1:Matlab程序求参数相关性
clc,clear all,close all
a=xlsread('E:\MATLAB数学建模\数模竞赛\2021届试题\B\附件1.xlsx','sheet2','A2:G75'); figure
plotmatrix(a); covmat=corrcoef(a); figure imagesc(covmat); grid; colorbar;
表2 参数相关性表
参数相关性
wt%Co/SiO2
Co/SiO2 和 HAP 的质量比
催化剂总量
乙醇浓度
温度
乙醇转化率
C4烯烃选择性
乙醇转化率相关性
C4烯烃选择性相关性
wt%Co/SiO2
1
-0.0078
0.1098
0.2266
-0.0147
-0.0255
-0.239
4
4
Co/SiO2 和 HAP 的质量比
-0.0078
1
-0.0413
0.0079
-0.0001
-0.0235
0.128
5
5
催化剂总量
0.1098
-0.0413
1
-0.2017
0.0041
0.3432
0.3634
2
3
乙醇浓度
0.2266
0.0079
-0.2017
1
-0.0546
-0.3225
-0.063
3
6
温度
-0.0147
-0.0001
0.0041
-0.0546
1
0.793
0.6522
1
2
乙醇转化率
-0.0255
-0.0235
0.3432
-0.3225
0.793
1
0.6816
1
C4烯烃选择性
-0.239
0.128
0.3634
-0.063
0.6522
0.6816
1
图1 参数相关性
相关性图形显示:
催化剂总量 对 乙醇转化率 及 C4烯烃选择性 弱相关;
温度 对 乙醇转化率 及C4烯烃选择性 强相关;
乙醇转化率 对 C4烯烃选择性 强相关。
由表2及图1确定研究次序:温度、催化剂总量、乙醇浓度、wt%Co/SiO2、Co/SiO2和HAP 的质量比。
参数:
x1:Co/SiO2
x2:Co/SiO2 和 HAP 的质量比
x3:催化剂总量
X4:乙醇浓度
t:温度
乙醇转化率: y1=f1(tt)= t=t+tt, a1*t^4+b1*t^3+c1*t^2+d1*t+e1(使用万能的幂级数做拟合)
C4 烯烃选择性:y2=f2(tt,y1)= t=t+tt, a2*t^4+b2*t^3+c2*t^2+d2*t+e2+y1^ff
1、温度:A3数据的C4烯烃收率最高,故选A3数据拟合
代码2:Lu脚本代码
!!!using["luopt","math"]; //使用命名空间 目标函数(tt,a1,b1,c1,d1,e1,a2,b2,c2,d2,e2,ff :
i,s,x1,x2,x3,x4,t,y1,y2,py,pc : tyArray,max)= { i=-1, s=0, (++i<max).while{
x1=tyArray(i,0), x2=tyArray(i,1), x3=tyArray(i,2), x4=tyArray(i,3),
t=tyArray(i,4)+tt, py=tyArray(i,5), pc=tyArray(i,6),
y1=a1*t^4+b1*t^3+c1*t^2+d1*t+e1, y2=a2*t^4+b2*t^3+c2*t^2+d2*t+e2+y1^ff,
s=s+[y1-py]^2+[y2-pc]^2 }, sqrt[s/(2*max)] }; main(::tyArray,max)= {
tyArray=matrix{ //存放实验数据 " 1 1 400 0.9 250 9.7 5.5 1 1 400 0.9 275 19.2 8.04 1
1 400 0.9 300 29.3 17.01 1 1 400 0.9 325 37.6 28.72 1 1 400 0.9 350 48.9 36.85
1 1 400 0.9 400 83.7 53.43 1 1 400 0.9 450 86.4 49.9 " }, len[tyArray,0,&max],
Opt1[@目标函数, optmax,1000, optmode,20, optdeep,20, optwaysimdeep, optwayconfra,
optwaylme] //Opt1函数全局优化 };
结果(tt,a1,b1,c1,d1,e1,a2,b2,c2,d2,e2,ff,均方差RMSE):
-261.8675335445256 -3.262424048642732e-007 9.2748570810669e-005
-6.276632898091342e-003 0.446933781199288 15.64595737717238
1.09830167656103e-006 -3.272228670670809e-004
2.28577641473638e-002 -0.9402976319921129 -25.07075625324265
1.234317092139116 0.704111408824131
由以上拟合参数求C4烯烃收率最大值,可求最佳温度。
代码3:Lu脚本代码
!!!using["luopt","math"];; //使用命名空间
init(::tt,a1,b1,c1,d1,e1,a2,b2,c2,d2,e2,ff) = new[real_s,12].SetArray["
-261.8675335445256 -3.262424048642732e-007 9.2748570810669e-005
-6.276632898091342e-003 0.446933781199288 15.64595737717238
1.09830167656103e-006 -3.272228670670809e-004 2.28577641473638e-002
-0.9402976319921129 -25.07075625324265 1.234317092139116 "].in[0 :
&tt,&a1,&b1,&c1,&d1,&e1,&a2,&b2,&c2,&d2,&e2,&ff];; f(x : t, y1, y2 :
tt,a1,b1,c1,d1,e1,a2,b2,c2,d2,e2,ff)= { t=x+tt,
y1=a1*t^4+b1*t^3+c1*t^2+d1*t+e1, y2=a2*t^4+b2*t^3+c2*t^2+d2*t+e2+y1^ff,
-y1*y2/100 }; Opt1[@f, optrange: 250.0, 500.0]; //250~500℃范围内求解,超出范围恐误差加大 g(x :
t, y1, y2 : tt,a1,b1,c1,d1,e1,a2,b2,c2,d2,e2,ff)= { t=x+tt,
y1=a1*t^4+b1*t^3+c1*t^2+d1*t+e1, y2=a2*t^4+b2*t^3+c2*t^2+d2*t+e2+y1^ff, o[y1,"
",y2," ",y1*y2/100,"\r\n"] }; g[426]; //求426℃时乙醇转化率、C4烯烃选择性、C4烯烃收率
结果:
426.3471802046889 -52.65073533983154
93.248565621854752 56.460787404028451 52.648874393061369
故426℃获得最大C4烯烃收率52.65%。
426℃时乙醇转化率为93.25%,C4烯烃选择性为56.46%,C4烯烃收率为52.65%。
表3 A3计算数据
乙醇转化率(%)
C4烯烃选择性(%)
乙醇转化率计算值(%)
C4烯烃选择性计算值(%)
9.7
5.5
9.296473618
5.55136596
19.2
8.0
20.63317954
7.750893896
29.3
17.0
28.01480073
17.65844747
37.6
28.7
37.00061664
28.01122598
48.9
36.9
50.09138418
37.19025191
83.7
53.4
83.29818902
53.38249601
86.4
49.9
86.47081958
49.9053186
图2 A3拟合数据对比
2、催化剂总量:数据量不足,暂不研究
3、乙醇浓度:A7、A8、A12、A9(乙醇浓度0.3、0.9、1.68、2.1)
表4 A7、A8、A12、A9拟合参数
实验
拟合参数(代码2):tt,a1,b1,c1,d1,e1,a2,b2,c2,d2,e2,ff
均方差 RMSE
最大C4烯烃收率
426℃时(代码3)
温度(℃)
收率(%)
乙醇转化率(%)
C4烯烃选择性(%)
C4烯烃收率(%)
A7
-252.8899179528622 1.866666666666673e-007 -5.317552792852978e-005
4.222960774025418e-003 0.3120916811555708 20.56535432567161
-4.586666666666685e-008 1.069646305291495e-005
5.737007359059772e-004 1.677680281411836e-002 5.843953613601379
-4776.690299232196
1.613456224144101e-015
升高
增大
92.92
40.24
37.39
A8
-281.4005327840787 -1.461333333333339e-007 1.788534190327997e-005
1.990680192658441e-003 0.1476147740840605 9.668196322759862
1.209119781600527e-007 -2.107766143733289e-005
1.331148967266221e-003 0.1867554749576292 7.425348864916348
0.3645782561680974
5.617333549722722e-016
461
47.71
62.82
55.92
35.13
A12
-242.4216261300292 -5.333333333333015e-009 7.228338642559277e-006
3.708339912735445e-004 5.947657275415871e-002 0.9248381845068223
-6.188047295963441e-008 1.617363381590294e-005
-1.88202661958118e-004 6.426014442682737e-002 5.716727247155738
-24.54643030158932
7.021666937153402e-016
467.5
36.92
63.00
40.95
25.80
A9
-274.5222725537157 7.253333333333393e-008 2.154728676917969e-006
5.914792247090588e-004 5.000002803548012e-002 2.975978386120882
3.199999999995535e-009 -1.14461149113099e-005
1.69439998784743e-003 0.2195389126636296 8.594734781755825
6.762159126008593e-012
2.27311216819532e-015
486.4
61.53
69.80
43.63
30.45
表5乙醇浓度对C4烯烃收率的影响
wt%Co/SiO2
Co/SiO2 和 HAP 的质量比
催化剂总量
乙醇浓度
温度
乙醇转化率(%)
C4烯烃选择性(%)
C4烯烃收率(%)
1
1
100
0.3
426
92.92
40.24
37.39
1
1
100
0.9
426
62.82
55.92
35.13
1
1
100
1.68
426
63
40.95
25.8
1
1
100
2.1
426
69.8
43.63
30.45
1
1
100
0.3
400
76.0
33.3
25.27911589
1
1
100
0.9
400
56.1
41.4
23.24340866
1
1
100
1.68
400
44.5
36.3
16.16619292
1
1
100
2.1
400
40.8
42.0
17.15112728
发现426℃与400℃时一致:乙醇浓度减小或增大,均可使C4烯烃收率增大,因缺少数据,使C4烯烃收率最大时的乙醇浓度不好拟合(本来要拟合426℃时最佳乙醇浓度)。
4、wt%Co/SiO2:A4、A1、A2、A6(wt%Co/SiO2:0.5、1、2、5)
将A4、A1、A2、A6拟合出426℃数据最佳,但因缺少400℃左右及以上实验数据,恐误差较大,故直接选择350℃数据做拟合。
表6 wt%Co/SiO2对C4烯烃收率的影响
wt%Co/SiO2
Co/SiO2 和 HAP 的质量比
催化剂总量
乙醇浓度
温度
乙醇转化率(%)
C4烯烃选择性(%)
C4烯烃收率(%)
0.5
1
400
1.68
350
60.5
27.25
16.47764673
1
1
400
1.68
350
36.8
47.21
17.37376011
2
1
400
1.68
350
67.9
39.1
26.54080464
5
1
400
1.68
350
55.8
10.65
5.937910401
代码4:Lu脚本代码
!!!using["luopt","math"]; //使用命名空间 目标函数(tt,a1,b1,c1,d1,e1,a2,b2,c2,d2,e2,ff :
i,s,x1,x2,x3,x4,t,y1,y2,py,pc : tyArray,max)= { i=-1, s=0, (++i<max).while{
x1=tyArray(i,0)+tt, x2=tyArray(i,1), x3=tyArray(i,2), x4=tyArray(i,3),
t=tyArray(i,4), py=tyArray(i,5), pc=tyArray(i,6),
y1=a1*x1^4+b1*x1^3+c1*x1^2+d1*x1+e1, y2=a2*x1^4+b2*x1^3+c2*x1^2+d2*x1+e2+y1^ff,
s=s+[y1-py]^2+[y2-pc]^2 }, sqrt[s/(2*max)] }; main(::tyArray,max)= {
tyArray=matrix{ //存放实验数据 " 0.5 1 400 1.68 350 60.5 27.25 16.47764673 1 1 400
1.68 350 36.8 47.21 17.37376011 2 1 400 1.68 350 67.9 39.1 26.54080464 5 1 400
1.68 350 55.8 10.65 5.937910401 " }, len[tyArray,0,&max], Opt1[@目标函数,
optmax,1000, optmode,20, optdeep,20, optwaysimdeep, optwayconfra, optwaylme]
//Opt1函数全局优化 };
结果(tt,a1,b1,c1,d1,e1,a2,b2,c2,d2,e2,ff,均方差RMSE):
-2.84889852050221 5.076959901361098 1.119333563495891
-31.18314598500931 -2.724602384343222 86.10680646706408
-2.915629722656612 -1.403220910994298 12.10486872551318
-3.90874000069633e-015 31.03258015142568 -9062704915.511139
2.87802701599909e-015
由代码3及以上拟合参数求C4烯烃收率最大值,可求最佳wt%Co/SiO2。
求解时修改代码:
Opt1[@f, optrange: 0.001, 10.0]; //0.001~10范围内求解,超出范围恐误差加大
g[2.32]; //求wt%Co/SiO2=2.32时乙醇转化率、C4烯烃选择性、C4烯烃收率
2.321036989494045 -27.19401258771866
79.056540095119573 34.398176048303348 27.194007839616752
故wt%Co/SiO2=2.32时获得最大C4烯烃收率27.19%。
wt%Co/SiO2=2.32时乙醇转化率为79.06%,C4烯烃选择性为34.40%,C4烯烃收率为27.19%。
5、Co/SiO2和HAP 的质量比:数据量不足,暂不研究
以上求解,供研究 “B题 乙醇偶合制备 C4 烯烃” 作参考。