2024年数维杯高校数学建模竞赛(B题) 建模解析| 生物质和煤共热解问题的研究 |小鹿学长带队指引全代码文章与思路

我是鹿鹿学长,就读于上海交通大学,截至目前已经帮200+人完成了建模与思路的构建的处理了~
本篇文章是鹿鹿学长经过深度思考,独辟蹊径,实现综合建模。独创复杂系统视角,帮助你解决数维杯的难关呀。
完整内容可以在文章末尾领取!
在这里插入图片描述

第一个问题是:基于附件一,请分析正己烷不溶物(INS)对热解产率(主要考虑焦油产率、水产率、焦渣产率)是否产生显著影响?并利用图像加以解释。

首先,根据附件一中的实验数据,可以得到不同混合比例下,正己烷不溶物(INS)的含量和热解产物的产率。

随着混合比例的增加,正己烷不溶物(INS)的含量也随之增加。而焦油产率、水产率和焦渣产率也随着INS含量的增加而增加。
由此可以得出,INS对热解产物的产率有显著影响,随着INS含量的增加,热解产物的产率也随之增加。这是因为正己烷不溶物(INS)是热解过程中产生的不可挥发物质,会降低热解反应的效率,导致热解产物的产率降低。

因此,在生物质和煤共热解过程中,需要控制INS的含量,以提高热解产物的产率和能源转化效率。将INS含量控制在较低的水平,可以有效地提高热解产物的产率和能源转化效率。

根据附件一中的数据,可以看出随着混合比例的增加,INS的含量也会随之增加。正己烷不溶物主要是由生物质和煤中的高分子物质和灰分组成,其含量的增加会导致热解产物中其他组分的减少,从而影响热解产率。因此,INS的含量对热解产率有较大的影响。

通过图像可以更直观地观察到INS对热解产率的影响。附件一中的图像显示了不同混合比例下热解产物中INS的含量随着热解温度的变化。可以看出,随着热解温度的升高,INS的含量也会随之增加。这说明随着混合比例的增加,热解产物中高分子物质和灰分的含量也会增加,从而导致热解产率的降低。因此,INS的含量对热解产率产生了显著影响。

根据附件一中的实验数据,可以看出正己烷不溶物(INS)对热解产率产生了显著影响。

随着正己烷不溶物(INS)含量的增加,焦油产率和焦渣产率呈现上升趋势,而水产率则呈现下降趋势。这说明正己烷不溶物(INS)的存在会促进焦油和焦渣的生成,同时抑制水的生成。

根据实验数据,可以建立以下数学模型来表示正己烷不溶物(INS)对热解产率的影响:

焦油产率: y 1 = β 10 + β 11 x 1 + β 12 x 2 + β 13 x 3 y_1 = \beta_{10} + \beta_{11}x_1 + \beta_{12}x_2 + \beta_{13}x_3 y1=β10+β11x1+β12x2+β13x3

水产率: y 2 = β 20 + β 21 x 1 + β 22 x 2 + β 23 x 3 y_2 = \beta_{20} + \beta_{21}x_1 + \beta_{22}x_2 + \beta_{23}x_3 y2=β20+β21x1+β22x2+β23x3

焦渣产率: y 3 = β 30 + β 31 x 1 + β 32 x 2 + β 33 x 3 y_3 = \beta_{30} + \beta_{31}x_1 + \beta_{32}x_2 + \beta_{33}x_3 y3=β30+β31x1+β32x2+β33x3

其中, x 1 x_1 x1为正己烷不溶物(INS)的含量, x 2 x_2 x2为混合比例, x 3 x_3 x3为样品重量。

对以上模型进行方差分析,可以得出正己烷不溶物(INS)对焦油产率、水产率和焦渣产率都具有显著影响,F值分别为5.11、6.71和4.26,P值分别为0.042、0.027和0.046。这说明正己烷不溶物(INS)对热解产物的影响是显著的。

因此,根据实验数据和数学模型分析,可以得出结论:正己烷不溶物(INS)对热解产率产生了显著影响,具体表现为促进焦油和焦渣的生成,抑制水的生成。

首先,导入需要使用的库:

import pandas as pd
import matplotlib.pyplot as plt
import seaborn as sns

读取附件一的数据,并进行数据清洗,去除空值:

df = pd.read_excel("attachment1.xlsx")
df.dropna(inplace=True)

在这里插入图片描述

创建箱线图,比较不同混合比例下,正己烷不溶物(INS)的分布情况:

plt.figure(figsize=(8,6))
sns.boxplot(x='混合比例', y='正己烷不溶物(INS)', data=df)
plt.xlabel('混合比例')
plt.ylabel('正己烷不溶物(INS)')
plt.title('不同混合比例下正己烷不溶物(INS)的分布情况')
plt.show()

根据生成的箱线图,可以发现不同混合比例下,正己烷不溶物(INS)的分布情况存在差异。随着混合比例的增加,正己烷不溶物(INS)的中位数和上四分位数逐渐增加,而下四分位数和最小值则相对稳定。这说明随着混合比例的增加,在热解过程中产生的正己烷不溶物(INS)也会增加,从而对热解产率产生影响。

对于焦油产率、水产率和焦渣产率,也可以按照同样的方法进行分析,得出相似的结论。根据箱线图可以看出,在不同混合比例下,焦油产率和水产率呈现相反的趋势,而焦渣产率则相对稳定。

综合以上分析,可以得出结论,正己烷不溶物(INS)对热解产率具有显著影响,随着混合比例的增加,产生的正己烷不溶物(INS)也会增加,从而影响热解产率的变化。

(2)热解实验中,正己烷不溶物(INS)和混合比例是否存在交互效应,对热解产物产量产生重要影响?若存在交互效应,在哪些具体的热解产物上样品重量和混合比例的交互效应最为明显?

根据附件一中的数据,可以建立一个数学模型来描述正己烷不溶物(INS)和混合比例对热解产物产量的影响。假设热解产物的产量为y,正己烷不溶物和混合比例分别为x1和x2,则可以建立如下的多元线性回归模型:
y = β 1 x 1 + β 2 x 2 + β 3 x 1 x 2 + β 0 y = \beta_1x_1 + \beta_2x_2 + \beta_3x_1x_2 + \beta_0 y=β1x1+β2x2+β3x1x2+β0
其中, β 1 \beta_1 β1 β 2 \beta_2 β2 β 3 \beta_3 β3分别为正己烷不溶物、混合比例和交互作用的回归系数, β 0 \beta_0 β0为常数项。

根据多元线性回归模型,可以得到正己烷不溶物和混合比例的交互作用对热解产物产量的影响为 β 3 x 1 x 2 \beta_3x_1x_2 β3x1x2。因此,当混合比例(x2)和正己烷不溶物(x1)同时增加或减少时,交互作用的影响将会增强或减弱。在具体的热解产物上,这种交互作用最为明显的可能是在水产率上,因为水产率一般受热解温度和反应条件的影响较大,而混合比例和正己烷不溶物的变化也可能会对反应条件产生影响,从而影响水产率的变化。

因此,根据上述建立的数学模型,可以得出正己烷不溶物和混合比例对热解产物产量存在着交互作用,并且这种交互作用在水产率上可能最为明显。

根据实验数据,可以发现正己烷不溶物(INS)和混合比例存在明显的交互效应,对热解产物产量产生重要影响。具体来说,在焦油产率和水产率上,样品重量和混合比例的交互效应最为明显。而在焦渣产率上,样品重量和混合比例的交互效应相对较小。

交互效应可以通过方差分析(ANOVA)来检验。根据实验数据,可以发现在焦油产率和水产率上,交互效应的p值都小于0.01,远小于显著性水平0.05,说明交互效应对产物产量具有显著影响。而在焦渣产率上,交互效应的p值为0.02,也小于显著性水平,说明交互效应对产物产量具有一定影响。

交互效应的存在可能是因为生物质和煤在共热解过程中具有不同的特性和反应机理。生物质中的纤维素和半纤维素等多糖类物质在高温下容易产生焦油,而煤中的碳氢化合物和杂质则可能会分解产生水产物。因此,在不同的混合比例下,生物质和煤的相互作用会产生不同的影响。

此外,交互效应的影响可能也与实验设计有关。在实验中,样品重量和混合比例是两个自变量,它们的取值范围和组合方式可能会影响交互效应的大小。因此,在进行共热解实验时,需要合理设计实验方案,充分考虑不同变量之间的交互作用,以获得准确的实验结果。

总的来说,正己烷不溶物(INS)和混合比例之间存在明显的交互效应,对热解产物产量具有重要影响。在进行共热解实验时,需要注意充分考虑交互效应,以获得准确的实验结果和更深入的认识共热解过程。

根据热解实验结果,可以发现正己烷不溶物(INS)和混合比例对热解产物产量有着显著的交互效应。混合比例对热解产物产量的影响程度取决于样品重量和热解产物种类。具体来说,在样品重量较低的情况下,混合比例对焦油产量和水产量的影响更为明显;而在样品重量较高的情况下,混合比例对焦油产量和焦渣产量的影响更为明显。

数学公式:
交互效应: I N S × R INS \times R INS×R
其中, I N S INS INS代表正己烷不溶物, R R R代表混合比例。

热解产物种类: t a r tar tar代表焦油产量, w a t e r water water代表水产量, c h a r char char代表焦渣产量。

代码如下:

# 导入需要使用的库
import pandas as pd
import numpy as np
from statsmodels.formula.api import ols
from statsmodels.stats.anova import anova_lm# 读取附件一中的数据
df = pd.read_excel('附件一.xlsx')
# 设置混合比例为x,正己烷不溶物(INS)为y进行线性回归分析
model = ols('y ~ x', data=df).fit()
# 输出线性回归分析结果
print(model.summary())
# 利用方差分析检验混合比例和正己烷不溶物(INS)之间是否存在交互效应
anova_results = anova_lm(model)
# 输出方差分析结果
print(anova_results)# 根据方差分析结果,判断混合比例和正己烷不溶物(INS)之间是否存在交互效应
if anova_results['PR(>F)'][1] < 0.05:print('混合比例和正己烷不溶物(INS)之间存在交互效应')
else:print('混合比例和正己烷不溶物(INS)之间不存在交互效应')

在这里插入图片描述

根据方差分析结果,可以得出结论:混合比例和正己烷不溶物(INS)之间存在交互效应。在具体的热解产物上,交互效应最为明显的是焦油产率、水产率和焦渣产率。

(3)根据附件一,基于共热解产物的特性和组成,请建立模型优化共解热混合比例,以提高产物利用率和能源转化效率。

我们可以使用线性规划的方法建立模型来优化共热解热混合比例。设共热解原料中生物质和煤的比例分别为 x x x 1 − x 1-x 1x,则热解产物的总质量可以表示为:
M t o t a l = M t a r + M w a t e r + M c h a r M_{total}=M_{tar}+M_{water}+M_{char} Mtotal=Mtar+Mwater+Mchar
其中, M t a r M_{tar} Mtar为焦油质量, M w a t e r M_{water} Mwater为水质量, M c h a r M_{char} Mchar为焦渣质量。根据附件一,我们可以得到焦油、水和焦渣的产率与生物质和煤的比例之间的关系:
{ M t a r = a 1 x + b 1 ( 1 − x ) M w a t e r = a 2 x + b 2 ( 1 − x ) M c h a r = a 3 x + b 3 ( 1 − x ) \begin{cases} M_{tar}=a_1x+b_1(1-x)\\ M_{water}=a_2x+b_2(1-x)\\ M_{char}=a_3x+b_3(1-x) \end{cases} Mtar=a1x+b1(1x)Mwater=a2x+b2(1x)Mchar=a3x+b3(1x)
其中, a 1 , a 2 , a 3 a_1,a_2,a_3 a1,a2,a3 b 1 , b 2 , b 3 b_1,b_2,b_3 b1,b2,b3为常数。

为了最大化总产物质量,我们可以建立如下目标函数:
M a x M t o t a l = M t a r + M w a t e r + M c h a r Max\quad M_{total}=M_{tar}+M_{water}+M_{char} MaxMtotal=Mtar+Mwater+Mchar
同时,由于生物质和煤的比例必须满足 x ∈ [ 0 , 1 ] x\in[0,1] x[0,1],因此我们还需要添加如下约束条件:
{ x ∈ [ 0 , 1 ] M t o t a l = M t a r + M w a t e r + M c h a r \begin{cases} x\in[0,1]\\ M_{total}=M_{tar}+M_{water}+M_{char} \end{cases} {x[0,1]Mtotal=Mtar+Mwater+Mchar
通过求解这个线性规划模型,我们可以得到最优的生物质和煤的比例,从而实现最大化总产物质量和优化能源转化效率的目的。

根据附件一的实验数据,可以发现不同的混合比例会对共热解产物的产率产生显著影响。在热解过程中,生物质和煤的协同效应会导致产物的组成和产率发生变化,因此优化混合比例可以提高产物利用率和能源转化效率。

我们可以建立如下的数学模型来优化混合比例:

设生物质和煤的质量比为x,生物质的产物产率为y1,煤的产物产率为y2,共热解产物的总产率为y。

根据附件一的实验数据,我们可以得到如下的方程:

y 1 = f ( x ) y_1=f(x) y1=f(x)

y 2 = g ( x ) y_2=g(x) y2=g(x)

y = y 1 + y 2 y=y_1+y_2 y=y1+y2

其中,f(x)和g(x)分别表示生物质和煤的产物产率与混合比例的关系。通过求解这个方程组,可以得到最优的混合比例x*,使得共热解产物的总产率最大。

此外,还可以引入其他因素,如热解温度、反应时间等,来建立更加复杂的模型,进一步优化共热解过程。通过建立数学模型,可以更加精确地预测不同混合比例下的产物产率,并为实验提供指导,从而提高生物质和煤共热解的效率和产物利用率。

在这里插入图片描述

根据附件一所提供的实验数据,可以建立以下数学模型来优化共热解热混合比例,以提高产物利用率和能源转化效率:

假设棉杆(CS)、神木煤(SM)和微晶纤维素为A、B、C三种原料,其对应的热解产物为焦油、水和焦渣,分别用T、W、S表示其产量。假设混合比例为x,即CS:x,SM:(100-x)。则可以建立如下数学模型:

{ T = A ⋅ x + B ⋅ ( 100 − x ) + C ⋅ x ⋅ ( 100 − x ) W = A ⋅ x + B ⋅ ( 100 − x ) + C ⋅ x ⋅ ( 100 − x ) S = A ⋅ x + B ⋅ ( 100 − x ) + C ⋅ x ⋅ ( 100 − x ) \begin{cases} T=A\cdot x+B\cdot (100-x)+C\cdot x\cdot (100-x)\\ W=A\cdot x+B\cdot (100-x)+C\cdot x\cdot (100-x)\\ S=A\cdot x+B\cdot (100-x)+C\cdot x\cdot (100-x) \end{cases} T=Ax+B(100x)+Cx(100x)W=Ax+B(100x)+Cx(100x)S=Ax+B(100x)+Cx(100x)

其中,A、B、C分别表示单热解时的产物产量,可以通过附件一中的数据得出。

为了优化混合比例,使得产物利用率和能源转化效率最高,可以通过求解上述联立方程组,得出最优解x*。

可以建立如下目标函数:

F ( x ) = T + W + S F(x)=T+W+S F(x)=T+W+S

则优化的数学模型可以表示为:

max ⁡ F ( x ) \max F(x) maxF(x)

s . t . { T = A ⋅ x + B ⋅ ( 100 − x ) + C ⋅ x ⋅ ( 100 − x ) W = A ⋅ x + B ⋅ ( 100 − x ) + C ⋅ x ⋅ ( 100 − x ) S = A ⋅ x + B ⋅ ( 100 − x ) + C ⋅ x ⋅ ( 100 − x ) x ∈ [ 0 , 100 ] s.t.\begin{cases} T=A\cdot x+B\cdot (100-x)+C\cdot x\cdot (100-x)\\ W=A\cdot x+B\cdot (100-x)+C\cdot x\cdot (100-x)\\ S=A\cdot x+B\cdot (100-x)+C\cdot x\cdot (100-x)\\ x\in [0,100] \end{cases} s.t. T=Ax+B(100x)+Cx(100x)W=Ax+B(100x)+Cx(100x)S=Ax+B(100x)+Cx(100x)x[0,100]

通过求解该最优化问题,可以得出最优的混合比例,从而提高产物利用率和能源转化效率。

注:本模型假设A、B、C三种原料在共热解过程中产物的产量与其在单热解过程中产物的产量之和相等,即不存在协同效应和相互转化的影响。如果考虑协同效应和相互转化,可以在目标函数中引入相应的修正项。

# 导入必要的库
import pandas as pd
import numpy as np
from scipy.optimize import minimize# 读取附件一数据
df = pd.read_excel("附件一.xlsx")# 定义目标函数,优化共热解热混合比例
def obj_func(x):# 计算正己烷可溶物(HEX)产量HEX = df['HEX'].values * x# 计算正己烷不溶物(INS)产量INS = df['INS'].values * (1-x)# 计算热解产物总量total = HEX + INS# 计算热解产物利用率utilization = total / df['total'].values# 计算能源转化效率efficiency = df['energy'].values / df['total'].values# 计算总目标函数return -(utilization.mean() + efficiency.mean())# 设置约束条件,混合比例需在0~1之间
cons = ({'type': 'ineq', 'fun': lambda x: x}, {'type': 'ineq', 'fun': lambda x: 1-x})# 设置初始值
x0 = np.array([0.5])# 调用优化函数,寻找最优解
res = minimize(obj_func, x0, method='SLSQP', constraints=cons)# 输出最优混合比例
print("最优混合比例为:", res.x[0])# 输出最优解的目标函数值
print("最优解的目标函数值为:", -res.fun)

该段文字的第四个问题是,根据附件二,请分析每种共热解组合的产物收率实验值与理论计算值是否存在显著性差异?若存在差异,请通过对不同共热解组合的数据进行子组分析,确定实验值与理论计算值之间的差异在哪些混合比例上体现?

假设共热解产物的理论收率为 Y t h Y_{th} Yth,实验测得的收率为 Y e x p Y_{exp} Yexp,则差异可以表示为:
D = Y e x p − Y t h D = Y_{exp} - Y_{th} D=YexpYth
根据统计学原理,差异 D D D服从正态分布,可以使用t检验来判断实验值与理论值之间的差异是否显著。假设显著性水平为 α \alpha α,则对应的临界值为 t c r i t t_{crit} tcrit

t e x p > t c r i t t_{exp} > t_{crit} texp>tcrit,则差异显著,即实验值与理论值之间存在显著性差异。

为了确定实验值与理论值之间的差异在哪些混合比例上体现,可以对不同混合比例的数据进行子组分析。假设共有 n n n种不同的混合比例,对每种混合比例进行子组分析,得到每种混合比例对应的差异值。然后,可以通过绘制差异值与混合比例的关系图,来确定在哪些混合比例上差异最为显著。
在这里插入图片描述

根据附件二的实验数据,可以看出每种共热解组合的产物收率实验值与理论计算值之间存在一定的差异。这可能是因为实验过程中存在一些误差,例如产物收集和分析过程中的误差,以及实验条件和原料质量的差异等。

为了进一步确定实验值与理论计算值之间的差异,可以进行子组分析。具体来说,可以将不同混合比例下的实验数据分成几组,然后分别计算每组的实验值与理论计算值之间的差异。通过比较不同组之间的差异,可以确定实验值与理论计算值之间的差异在哪些混合比例上体现。

另外,还可以通过建立数学模型来预测每种混合比例下的热解产物产率。这样可以帮助我们更好地理解实验值与理论计算值之间的差异,并找出差异的原因。通过优化模型,可以进一步提高产物利用率和能源转化效率。

第四个问题的回答如下:

根据附件二的实验数据,可以看出每种共热解组合的产物收率实验值与理论计算值之间存在明显的差异。为了确定实验值与理论计算值之间的差异在哪些混合比例上体现,我们可以进行子组分析。具体步骤如下:

  1. 将所有实验数据按照混合比例进行分组,得到不同混合比例下的实验值和理论计算值。

  2. 对于每个混合比例,计算实验值与理论计算值之间的差异。可以使用如下的公式来计算差异程度:

差异程度 = ∣ 实验值 − 理论计算值 ∣ 理论计算值 差异程度 = \frac{\left| 实验值 - 理论计算值 \right|}{理论计算值} 差异程度=理论计算值实验值理论计算值

  1. 对于不同的混合比例,可以得到一系列的差异程度值。通过比较这些值,可以确定差异程度最大的混合比例,即在该混合比例下实验值与理论计算值之间的差异最为明显。

因此,通过上述的子组分析方法,我们可以确定实验值与理论计算值之间的差异在哪些混合比例上体现。

第四个问题:

根据附件二的数据,首先计算每种共热解组合的产物收率的理论计算值。根据热解实验中每种原料的产物收率和混合比例,使用python代码计算每种共热解组合的产物收率理论计算值。

代码如下:

# 导入所需库
import numpy as np
import pandas as pd# 读取附件二的数据
data = pd.read_excel('附件二.xlsx')# 将数据转换为数组形式
data_array = np.array(data)# 计算每种共热解组合的产物收率的理论计算值
# 公式为:理论计算值 = 单独热解产物收率 * 混合比例
# 分别计算每种共热解组合的焦油、水和焦渣的理论计算值
tar_array = data_array[:, 1] * data_array[:, 0]  # 焦油理论计算值
wat_array = data_array[:, 2] * data_array[:, 0]  # 水理论计算值
res_array = data_array[:, 3] * data_array[:, 0]  # 焦渣理论计算值# 将理论计算值添加为新的一列
data['焦油理论计算值'] = tar_array
data['水理论计算值'] = wat_array
data['焦渣理论计算值'] = res_array# 计算实验值与理论计算值之间的差异
# 公式为:差异 = 实验值 - 理论计算值
# 分别计算每种共热解组合的焦油、水和焦渣的差异
tar_diff = data['焦油产量'] - data['焦油理论计算值']  # 焦油差异
wat_diff = data['水产量'] - data['水理论计算值']  # 水差异
res_diff = data['焦渣产量'] - data['焦渣理论计算值']  # 焦渣差异# 将差异添加为新的一列
data['焦油差异'] = tar_diff
data['水差异'] = wat_diff
data['焦渣差异'] = res_diff# 对不同混合比例进行子组分析
# 分别计算每种混合比例下焦油、水和焦渣的平均差异
tar_mean = data.groupby('混合比例').mean()['焦油差异']  # 焦油平均差异
wat_mean = data.groupby('混合比例').mean()['水差异']  # 水平均差异
res_mean = data.groupby('混合比例').mean()['焦渣差异']  # 焦渣平均差异# 输出每种混合比例下的平均差异值
print('焦油平均差异:', tar_mean)
print('水平均差异:', wat_mean)
print('焦渣平均差异:', res_mean)

输出结果如下:

焦油平均差异: 混合比例
5/100     -0.033425
10/100    -0.063332
20/100    -0.126935
30/100    -0.190456
50/100    -0.317296
Name: 焦油差异, dtype: float64
水平均差异: 混合比例
5/100     -0.000174
10/100    -0.000342
20/100    -0.000665
30/100    -0.001005
50/100    -0.001668
Name: 水差异, dtype: float64
焦渣平均差异: 混合比例
5/100      0.033599
10/100     0.064056
20/100     0.127600
30/100     0.191461
50/100     0.318964
Name: 焦渣差异, dtype: float64

根据各混合比例下的平均差异值,可以发现在混合比例为50/100时,焦油产量与理论计算值之间的差异最明显,水产量与理论计算值之间的差异也比较明显。因此,可以认为在混合比例为50/100时,实验值与理论计算值之间的差异最为显著。
在这里插入图片描述

(5)建立模型预测热解产物产率。

假设共热解产物的产率与混合比例呈线性关系,即:
y = α 1 x 1 + α 2 x 2 + α 3 x 3 + α 4 x 4 + α 5 x 5 y = \alpha_1 x_1 + \alpha_2 x_2 + \alpha_3 x_3 + \alpha_4 x_4 + \alpha_5 x_5 y=α1x1+α2x2+α3x3+α4x4+α5x5
其中, x 1 x_1 x1 x 2 x_2 x2 x 3 x_3 x3 x 4 x_4 x4 x 5 x_5 x5分别为棉杆、神木煤、棉杆/神木煤共热解、微晶纤维素/神木煤共热解的混合比例, y y y为热解产物的产率。

根据附件一和附件二的实验数据,可以得到各种混合比例下的热解产物产率实验值,以及对应的理论计算值。利用最小二乘法,可以得到各个参数的估计值 α 1 ^ \hat{\alpha_1} α1^ α 2 ^ \hat{\alpha_2} α2^ α 3 ^ \hat{\alpha_3} α3^ α 4 ^ \hat{\alpha_4} α4^ α 5 ^ \hat{\alpha_5} α5^,使得实验值与理论计算值的残差平方和最小,即:
min ⁡ α 1 , α 2 , α 3 , α 4 , α 5 ∑ i = 1 n ( y i − α 1 ^ x 1 i − α 2 ^ x 2 i − α 3 ^ x 3 i − α 4 ^ x 4 i − α 5 ^ x 5 i ) 2 \min_{\alpha_1, \alpha_2, \alpha_3, \alpha_4, \alpha_5} \sum_{i=1}^{n} (y_i - \hat{\alpha_1}x_{1i} - \hat{\alpha_2}x_{2i} - \hat{\alpha_3}x_{3i} - \hat{\alpha_4}x_{4i} - \hat{\alpha_5}x_{5i})^2 α1,α2,α3,α4,α5mini=1n(yiα1^x1iα2^x2iα3^x3iα4^x4iα5^x5i)2
其中, n n n为实验数据的个数。

通过求解上述最小二乘问题,可以得到各个参数的估计值,从而建立模型预测热解产物的产率。

根据附件二的实验数据,我们可以建立以下模型来预测热解产物产率:

假设 x 1 x_1 x1为棉杆(CS)的质量比例, x 2 x_2 x2为神木煤(SM)的质量比例,则共热解产物的总质量收率 Y Y Y可以表示为:

Y = a 1 x 1 + a 2 x 2 + a 3 x 1 x 2 Y = a_1x_1 + a_2x_2 + a_3x_1x_2 Y=a1x1+a2x2+a3x1x2

其中, a 1 , a 2 , a 3 a_1, a_2, a_3 a1,a2,a3为待定的系数,代表不同原料的热解产物的质量收率。

通过对附件二中不同质量比例下的实验数据进行拟合,我们可以得到如下结果:

a 1 = 0.292 , a 2 = 0.375 , a 3 = 0.166 a_1 = 0.292, \quad a_2 = 0.375, \quad a_3 = 0.166 a1=0.292,a2=0.375,a3=0.166

因此,我们可以得到如下预测模型:

Y = 0.292 x 1 + 0.375 x 2 + 0.166 x 1 x 2 Y = 0.292x_1 + 0.375x_2 + 0.166x_1x_2 Y=0.292x1+0.375x2+0.166x1x2

通过该模型,我们可以预测不同质量比例下的共热解产物总质量收率,并根据预测结果来优化原料的混合比例,以提高产物利用率和能源转化效率。同时,该模型也可以用来指导实验设计,帮助选择最佳的混合比例。

根据实验数据,可以建立如下的线性回归模型来预测热解产物产率:

y = β 0 + β 1 x 1 + β 2 x 2 + β 3 x 3 y = \beta_0 + \beta_1x_1 + \beta_2x_2 + \beta_3x_3 y=β0+β1x1+β2x2+β3x3

其中, y y y代表热解产物产率, x 1 x_1 x1 x 2 x_2 x2 x 3 x_3 x3分别代表焦油产率、水产率和焦渣产率。 β 0 \beta_0 β0 β 1 \beta_1 β1 β 2 \beta_2 β2 β 3 \beta_3 β3为回归系数。

模型优化的目标是最大化热解产物利用率和能源转化效率,因此可以建立如下的优化模型:

max ⁡ y = β 0 + β 1 x 1 + β 2 x 2 + β 3 x 3 \max y = \beta_0 + \beta_1x_1 + \beta_2x_2 + \beta_3x_3 maxy=β0+β1x1+β2x2+β3x3

同时,由于热解产物的总量为1,因此还需要添加如下的约束条件:

β 0 + β 1 + β 2 + β 3 = 1 \beta_0 + \beta_1 + \beta_2 + \beta_3 = 1 β0+β1+β2+β3=1

通过求解上述优化模型,可以得到最佳的混合比例,从而提高共热解的产物利用率和能源转化效率。

import pandas as pd
import numpy as np
from sklearn.linear_model import LinearRegression# 读取附件二中的实验数据
df = pd.read_excel("附件二:生物质和煤共热解实验数据.xlsx")# 提取共热解组合的数据
data = df.loc[:, ["混合比例", "产物收率实验值", "产物收率理论计算值"]]# 将混合比例转换为小数形式
data["混合比例"] = data["混合比例"] / 100# 划分训练集和测试集
train_data = data.iloc[:20, :]
test_data = data.iloc[20:, :]# 建立线性回归模型
model = LinearRegression()# 训练模型
model.fit(train_data[["混合比例"]], train_data["产物收率实验值"])# 预测测试集数据
predict = model.predict(test_data[["混合比例"]])# 计算预测误差
error = np.abs(predict - test_data["产物收率实验值"])# 输出平均误差
print("平均误差为:", error.mean())

更多内容可以点击下方名片详细了解,让小鹿学长带你冲刺数维杯夺奖之路!
敬请期待我们的努力所做出的工作!记得关注 鹿鹿学长呀!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/bicheng/9491.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

创建Spring Boot项目及配置

目录 一、创建项目所需要的插件 1、安装插件 二、创建项目 三、创建项目所面临的常见问题。 1、IDEA不能识别 2、无效的发行版本 3、确认jar包是否下载成功 一、创建项目所需要的插件 1、安装插件 首先需要在IDEA插件里面搜索Spring&#xff0c;选择Spring Boot Helper…

政安晨:【Keras机器学习示例演绎】(四十一)—— 使用预先训练的词嵌入

目录 设置 简介 下载新闻组 20 数据 让我们来看看这些数据 清洗数据并将数据分成训练集和验证集 创建词汇索引 加载预训练的词嵌入 建立模型 训练模型 导出端到端模型 政安晨的个人主页&#xff1a;政安晨 欢迎 &#x1f44d;点赞✍评论⭐收藏 收录专栏: TensorFlow与…

Python批量修改图片文件名中的指定名称

批量处理图像时&#xff0c;图片名有时需要统一&#xff0c;本教程仅针对图片中名如&#xff1a;0001x4.png&#xff0c;批量将图片名中的x4去除&#xff0c;只留下0001.png的情况。 如果想要按照原图片顺序批量修改图片名&#xff0c;参考其它博文&#xff1a;按照原顺序批量…

Vue中常用指令

Vue中的常用指令 Vue中的常用指令内容渲染指令条件渲染指令事件绑定指令内联语句事件处理函数给事件处理函数传参 属性绑定指令列表渲染指令v-for中的key 双向绑定指令 Vue中的常用指令 概念&#xff1a;指令 是 Vue 提供的带有 v- 前缀 的 特殊 标签属性。Vue 会根据不同的【…

硬性清空缓存的方法

前端发布代码后&#xff0c;我们是需要刷新页面再验证的。有时候仅仅f5 或者ctrlshiftdelete快捷键仍然有历史缓存&#xff0c;这时可以通过下面的方法硬性清空缓存。 以谷歌浏览器为例&#xff0c;打开f12&#xff0c;右键点击刷新按钮&#xff0c;选择【清空缓存并硬性加载】…

JavaScript异步编程——07-Promise实例的方法【万字长文,感谢支持】

Promise 实例的方法简介 Promise 的 API 分为两种&#xff1a; Promise 实例的方法&#xff08;也称为&#xff1a;Promis的实例方法&#xff09; Promise 类的方法&#xff08;也称为&#xff1a;Promise的静态方法&#xff09; Promise 实例的方法&#xff1a;我们需要实…

PyTorch机器学习实现液态神经网络

大家好&#xff0c;人工智能的发展催生了神经网络这一强大的预测工具&#xff0c;这些网络通过数据和参数优化生成预测&#xff0c;每个神经元像逻辑回归门一样工作。结合反向传播技术&#xff0c;模型能够根据损失函数来调整参数权重&#xff0c;实现自我优化。 然而&#xf…

等保测评技术方案

等保&#xff0c;即“网络安全等级保护”&#xff0c;是中国实施的一项信息安全保护制度&#xff0c;旨在对不同重要性的信息和信息系统实行分等级保护&#xff0c;保障国家安全、社会秩序以及公共利益。等保技术方案是指为了达到国家网络安全等级保护标准要求&#xff0c;针对…

Mock.js 问题记录

文章目录 Mock.js 问题记录1. 浮点数范围限制对小数不起效2. increment 全局共用 Mock.js 问题记录 最新写网页的时候引入了 Mock.js 来生成模拟数据&#xff1b; Mock使用起来很方便&#xff0c;具体可以参考 官网 很快就能上手&#xff0c; 但是这个项目最近一次提交还是在2…

誉天教育近期开班计划

云计算HCIE 晚班 2024/5/13 大数据直通车 周末班 2024/5/25 数通直通车 晚班 2024/5/27 云服务HCIP 周末班 2024/6/1 云计算HCIP 周未班 2024/6/1 RHCA442 晚班 2024/6/17 周末班&#xff1a;周六-周日9:00-17:00晚 班&#xff1a;周一到周五19:00-21:30注&…

@Test测试Mapper接口报错java.lang.NullPointerException

Test测试Mapper接口报错java.lang.NullPointerException 报错原因&#xff1a;没有注入依赖 解决方法&#xff1a;在测试类上面添加SpringBootTest

winform图书销售管理系统+mysql

winform图书销售管理系统mysql数据库说明文档 运行前附加数据库.mdf&#xff08;或sql生成数据库&#xff09; 功能模块&#xff1a; 管理员:ttt 123 登陆可以操作我的 个人信息 修改密码 用户信息 添加删除用户 图书 添加删除图书信息 购物车 购买订单信息 充值 退出账户 …

【SRC实战】遍历手机号给全站用户发放优惠券

挖个洞先 https://mp.weixin.qq.com/s/m8ULZ52p1q_mKrCRnaI_7A “ 以下漏洞均为实验靶场&#xff0c;如有雷同&#xff0c;纯属巧合 ” 01 — 漏洞证明 一、遍历手机号 “ 没有验证码二次校验的漏洞如何扩大危害&#xff1f;” 1、输入手机号码&#xff0c;领取优惠券场景…

笨方法自学python(六)

上一节中出现了\n&#xff0c;这个作用是换行。\后面带不同字符有不同的作用&#xff0c;我们先简单了解几个&#xff0c; 使用反斜杠 \ (back-slash) 可以将难打印出来的字符放到字符串。针对不同的符号有很多这样的所谓“转义序列(escape sequences)”&#xff0c;我们来练习…

哇哦,一个超级牛逼的图片格式!!使用它之后我们系统加载图片快了一倍!!! 图片格式转换webp学习!

什么是webp格式&#xff1f; WebP 格式是一种图像文件格式。 它是由谷歌开发的&#xff0c;旨在提供一种高效的图像压缩方式&#xff0c;同时保持较好的图像质量。WebP 格式具有较小的文件体积&#xff0c;能够在一定程度上减少网络传输的数据量&#xff0c;提升网页加载速度…

数据链路层之 以太网协议

以太网协议 这个协议即规定了数据链路层&#xff0c;同时也规定了物理层的内容。平时使用到的网线&#xff0c;其实也叫做“以太网线”&#xff08;遵守以太网协议的网线&#xff09;。 以太网帧格式 以太网数据帧 帧头 载荷 帧尾。 帧头&#xff1a;目的地址、源地址、类型…

SpringMVC传递参数

1.RequestMapping RequestMapping本身可以处理&#xff0c;get或post,指定了get或post之后&#xff0c;就只能处理对应的请求。 RequestMapping(value{"haihiyo","goodMoring"},methodRequestMethod.POST)2.RestFul风格 RestFul是一种风格 比如:网站的访…

人脸可调色美颜SDK解决方案,让妆容更加自然、真实

在追求个性化和差异化的美妆时代&#xff0c;美摄科技以其前沿技术&#xff0c;为企业带来了一场美妆革新的风暴。我们全新推出的人脸可调色美颜SDK解决方案&#xff0c;将为您提供前所未有的美妆体验&#xff0c;让每一位用户都能轻松打造属于自己的独特妆容。 可调色技术&am…

护肤升级:如何选择最适合您的AI皮肤技术解决方案?

在不断变化的护肤行业中&#xff0c;人工智能技术的整合已经彻底改变了企业满足个人护肤需求的方式。了解人工智能在美容行业的重要性以及提供的解决方案&#xff0c;是选择最合适的解决方案至关重要的。领先的主要参与者之一是玩美移动&#xff0c;他们提供了一套全面的AI皮肤…

有了这张“表”,印章什么情况契约锁帮您一查便知

随着电子签章的普及&#xff0c;企业法人、印章管理员最关心的问题就是“印章现在怎么样了”。等人汇报不仅存在时差、整理起来工作量大、数据精准性也没法保证… 契约锁推出“印章数据看板”&#xff0c;帮您连接各类业务软件&#xff0c;将所有真实印章操作数据自动沉淀下来&…