CUMCM 2021-B：乙醇偶合制备C4烯烃（多元线性回归分析）

最新推荐文章于 2024-01-19 08:15:00 发布
最新推荐文章于 2024-01-19 08:15:00 发布
阅读量1w 收藏 115
点赞数 5
分类专栏： 北师大 数学模型 文章标签： matlab 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_55818063/article/details/124413334
版权

目录

一、竞赛赛题

二、问题（1）分析

三、问题（2）分析

四、多元线性回归分析（Stata）

1.多元自变量的选定

（1）乙醇转化率

（2）C4烯烃选择性

2.定量变量的描述性统计（Stata）

3.异方差检验

（1）散点图检验 

（2）BP检验（怀特检验的特例）

（3）怀特检验（使用范围广，优先考虑）

4.多重共线性检验（方差膨胀因子：VIF）

5.向后逐步回归

（1）乙醇转化率

（2）C4烯烃选择性

一、竞赛赛题

二、问题（1）分析

CUMCM 2021-B：乙醇偶合制备C4烯烃（1）_Destiny坠明的博客-CSDN博客

三、问题（2）分析

CUMCM 2021-B：乙醇偶合制备C4烯烃（2）_Destiny坠明的博客-CSDN博客

四、多元线性回归分析（Stata）

1.多元自变量的选定

（1）乙醇转化率

        由问题（1）拟合结果可知，乙醇转化率y与温度T呈Logistic阻滞增长关系，故需要将温度T做形式变换后加入到多元线性回归模型中。

        根据问题（1）的阻滞增长模型的拟合结果，r的波动较小，均值为0.0265，其分布散点图如下：

         y0在经过剔除离群值的处理之后，也稳定在一个区间内，均值为0.0386，其分布散点图如下：

        将原式第二项分母中的y0与r用其均值近似代替，并引入β作为补偿系数，可知β≈1.1；再将系数全部移入第一项中，即第一项整体视为多元回归系数β'。

        综上，乙醇转化率（Y1）的多元自变量（X）为：

X1装料方式
X2Co负载量
X3Co/SiO2
X4HAP
X5Co/SiO2 和 HAP装料比
X6乙醇浓度
X7温度的指数形式（阻滞增长）
X8温度的二次形式
X9温度

 

（2）C4烯烃选择性

        由问题（1）拟合结果可知，C4烯烃选择性y与温度T呈二次多项式关系，故需要将温度T做形式变换后加入到多元线性回归模型中。

         综上，C4烯烃选择性（Y2）的多元自变量（X）为：

X1装料方式
X2Co负载量
X3Co/SiO2
X4HAP
X5Co/SiO2 和 HAP装料比
X6乙醇浓度
X7温度的二次形式
X8温度

2.定量变量的描述性统计（Stata）

// 导入数据
import excel "C:\Users\熊锐成\Desktop\XY.xlsx", sheet("XY") firstrow

<XY.xlsx>

装料方式M1M2Co负载量（wt%）Co/SiO2(mg)HAP(mg)Co/SiO2和HAP装料比(mg:mg)乙醇浓度(ml/min)指数型温度二次项温度温度(℃)温度乘乙醇浓度乙醇转化率(%)C4烯烃选择性(%)
A10120020011.68753.7042126625002504202.07 34.05
A10120020011.681461.911347756252754625.85 37.43
A10120020011.682835.574959000030050414.97 46.94
// 定性变量的频数分布，并得到相应字母开头的虚拟变量
tabulate 装料方式 ,gen(M)
// 定量变量的描述性统计
summarize Co负载量wt CoSiO2mg HAPmg 乙醇浓度mlmin 温度 乙醇转化率 C4烯烃选择性

3.异方差检验

（1）散点图检验 

①残差-拟合值散点图（乙醇转化率）

// OLS回归(乙醇转化率)
regress 乙醇转化率 M1 Co负载量wt CoSiO2mg HAPmg CoSiO2和HAP装料比mgmg 乙醇浓度mlmin 指数型温度 二次项温度 温度

// 残差与拟合值的散点图
rvfplot

// 残差与自变量x的散点图
//rvpplot x

        当拟合值处在不同的区间范围内，残差的变化程度差别较大，存在异方差现象。

②残差-拟合值散点图（C4烯烃选择性）

// OLS回归(C4烯烃选择性)
regress C4烯烃选择性 M1 Co负载量wt CoSiO2mg HAPmg CoSiO2和HAP装料比mgmg 乙醇浓度mlmin 二次项温度 温度
// 残差与拟合值的散点图
rvfplot

         当拟合值处在不同的区间范围内，残差的变化程度差别较大，存在异方差现象。

（2）BP检验（怀特检验的特例）

①乙醇转化率

// OLS回归(乙醇转化率)
regress 乙醇转化率 M1 Co负载量wt CoSiO2mg HAPmg CoSiO2和HAP装料比mgmg 乙醇浓度mlmin 指数型温度 二次项温度 温度

// BP检验
estat hettest ,rhs iid

         P值小于0.05，说明在95%的置信水平下拒绝原假设，即我们认为乙醇转化率的多元线性回归模型的扰动项存在异方差。

②C4烯烃选择性

// OLS回归(C4烯烃选择性)
regress C4烯烃选择性 M1 Co负载量wt CoSiO2mg HAPmg CoSiO2和HAP装料比mgmg 乙醇浓度mlmin 二次项温度 温度

// BP检验
estat hettest ,rhs iid

         P值小于0.05，说明在95%的置信水平下拒绝原假设，即我们认为C4烯烃选择性的多元线性回归模型的扰动项存在异方差。

（3）怀特检验（使用范围广，优先考虑）

①乙醇转化率

// OLS回归(乙醇转化率)
regress 乙醇转化率 M1 Co负载量wt CoSiO2mg HAPmg CoSiO2和HAP装料比mgmg 乙醇浓度mlmin 指数型温度 二次项温度 温度

// 怀特检验
estat imtest,white

          P值小于0.05，说明在95%的置信水平下拒绝原假设，即我们认为乙醇转化率的多元线性回归模型的扰动项存在异方差。

②C4烯烃选择性

// OLS回归(C4烯烃选择性)
regress C4烯烃选择性 M1 Co负载量wt CoSiO2mg HAPmg CoSiO2和HAP装料比mgmg 乙醇浓度mlmin 二次项温度 温度

// 怀特检验
estat imtest,white

           P值小于0.05，说明在95%的置信水平下拒绝原假设，即我们认为C4烯烃选择性的多元线性回归模型的扰动项存在异方差。

         由于原始数据的异方差扰动严重，我们使用OLS+稳健的标准误作为指标，以防止异方差影响假设检验的计算。

4.多重共线性检验（方差膨胀因子：VIF）

①乙醇转化率

// OLS回归(乙醇转化率)
regress 乙醇转化率 M1 Co负载量wt CoSiO2mg HAPmg CoSiO2和HAP装料比mgmg 乙醇浓度mlmin 指数型温度 二次项温度 温度

// VIF检验
estat vif

         乙醇转化率的多元线性回归模型的平均VIF为200.93，远大于VIF最小经验值，故其存在多重共线性的问题。

②C4烯烃选择性

// OLS回归(C4烯烃选择性)
regress C4烯烃选择性 M1 Co负载量wt CoSiO2mg HAPmg CoSiO2和HAP装料比mgmg 乙醇浓度mlmin 二次项温度 温度

// VIF检验
estat vif

          C4烯烃选择性的多元线性回归模型的平均VIF为116.68，远大于VIF最小经验值，故其存在多重共线性的问题。

         此处我们将采用逐步回归的方法解决多重共线性的问题，以变量显著性为依据，不断筛选自变量，从而构建出最适的多元回归模型。

5.向后逐步回归

（1）乙醇转化率

// 向后逐步回归（乙醇转化率）
stepwise regress 乙醇转化率 M1 Co负载量wt CoSiO2mg HAPmg CoSiO2和HAP装料比mgmg 乙醇浓度mlmin 指数型温度 二次项温度 温度 ,  r pr(0.05)

①设定显著性水平为0.05，使用OLS+稳健的标准误作为指标，向后逐步回归得到多变量回归结果：

②进一步引入乙醇浓度与温度的交叉项，设定显著性水平为0.05，使用OLS+稳健的标准误作为指标，再次向后逐步回归：

// 添加 温度-乙醇 交叉项
stepwise regress 乙醇转化率 M1 Co负载量wt CoSiO2mg HAPmg CoSiO2和HAP装料比mgmg 乙醇浓度mlmin 指数型温度 二次项温度 温度 温度乘乙醇浓度 ,  r pr(0.05)

         加入乙醇浓度与温度的交叉项之后，拟合优度提高，显著因素增多，模型的泛化能力更强。

（2）C4烯烃选择性

// 向后逐步回归（C4烯烃选择性）
stepwise regress C4烯烃选择性 M1 Co负载量wt CoSiO2mg HAPmg CoSiO2和HAP装料比mgmg 乙醇浓度mlmin 二次项温度 温度 ,  r pr(0.05)

设定显著性水平为0.05，使用OLS+稳健的标准误作为指标，向后逐步回归得到多变量回归结果：

 

Destiny坠明
关注
  • 5
    点赞
  • 踩
  • 115
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
2021 年高教社杯全国大学生数学建模竞赛 B题 乙醇偶合制备 C4 烯烃 求解思路2
Lu的专栏
11-07 5715
由附件1“性能数据表”中A1~A14（用装料方式I）整理得下表“sheet2”： 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....
2021 年高教社杯全国大学生数学建模竞赛题 B 题 乙醇偶合制备 C4 烯烃
07-06
2021 年高教社杯全国大学生数学建模竞赛题 B 题 乙醇偶合制备 C4 烯烃 C4 烯烃广泛应用于化工产品及医药的生产，乙醇是生产制备 C4 烯烃的原料。 在制备过程中，催化剂组合（即：Co 负载量、Co/SiO2 和 HAP 装料比、乙醇浓度 的组合）与温度对 C4 烯烃的选择性和 C4 烯烃收率将产生影响（名词解释见附录）。 因此通过对催化剂组合设计，探索乙醇催化偶合制备 C4 烯烃的工艺条件具有非常 重要的意义和价值。 某化工实验室针对不同催化剂在不同温度下做了一系列实验，结果如附件 1 和 附件 2 所示。请通过数学建模完成下列问题： (1) 对附件 1 中每种催化剂组合，分别研究乙醇转化率、C4 烯烃的选择性与温 度的关系，并对附件 2 中 350 度时给定的催化剂组合在一次实验不同时间的测试结 果进行分析。 (2) 探讨不同催化剂组合及温度对乙醇转化率以及 C4 烯烃选择性大小的影响。 (3) 如何选择催化剂组合与温度，使得在相同实验条件下 C4 烯烃收率尽可能 高。若使温度低于 350 度，又如何选择催化剂组合与温度，使得 C4 烯烃收率尽可 能高。
CUMCM2021-B.pdf
10-12
2021年高教社杯全国大学生数学建模大赛B题
2021年国赛高教杯数学建模B题乙醇偶合制备C4烯烃解题全过程文档及程序
weixin_43292788的博客
06-09 5834
2021年国赛高教杯数学建模B题乙醇偶合制备C4烯烃解题全过程文档及程序
2021全国大学生数学建模B题 乙醇偶合制备 C4 烯烃
weixin_43292788的博客
09-16 9404
2021高教社杯全国大学生数学建模 B 题 乙醇偶合制备 C4 烯烃 C4 烯烃广泛应用于化工产品及医药的生产，乙醇是生产制备 C4 烯烃的原料。在制备过程中，催化剂组合（即：Co 负载量、Co/SiO2 和 HAP 装料比、乙醇浓度的组合）与温度对 C4 烯烃的选择性和 C4 烯烃收率将产生影响（名词解释见附录）。 因此通过对催化剂组合设计，探索乙醇催化偶合制备 C4 烯烃的工艺条件具有非常重要的意义和价值。 某化工实验室针对不同催化剂在不同温度下做了一系列实验，结果如附件 1 和附件 2 所示。请通过数
CUMCM 2021-B：乙醇偶合制备C4烯烃（2）
Destiny坠明的博客
04-20 7407
目录 一、竞赛赛题 二、问题(1)分析 三、问题(2)分析 探究不同催化剂组合及温度对乙醇转化率以及C4烯烃选择性大小的影响 （1）Co负载量 ①Co负载量——乙醇转化率 ②Co负载量——C4烯烃选择性 （2）Co/SiO2 和 HAP 装料比 ①Co/SiO2 和 HAP 装料比——乙醇转化率 ②Co/SiO2 和 HAP 装料比——C4烯烃选择性 （3）乙醇浓度 ①乙醇浓度——乙醇转化率 ②乙醇浓度——C4烯烃选择性 一、竞赛赛题 二、问题(1)分析 CUMCM
2021年数模国赛B题--乙醇偶合制备 C4 烯烃
m0_63630166的博客
08-11 8987
数模国赛题目
第四题乙醇偶合制备 C4 烯烃
herry57的博客
07-22 2027
b=readmatrix('F\数学建模\培训模拟\研究生组题目四\2021-B\第三题数据.xlsx','Sheet','Sheet3','Range','A1')a=readmatrix('F\数学建模\培训模拟\研究生组题目四\1单因素分析.xlsx','Sheet','C4烯烃选择性','Range','A2')a=readmatrix('F\数学建模\培训模拟\研究生组题目四\1单因素分析.xlsx','Sheet','乙醇转化率','Range','A2')n=length(Y);...
2021国赛数据处理
m0_61651594的博客
09-15 1626
当前，国内外正在大力研究C4烯烃的生产技术。C4 烯烃作为于化工产品及医药的生产的重要原料，对化工产业的发展具有重要意义。在生产过程中，催化剂组合与温度对 C4 烯烃的选择性和 C4 烯烃收率有着重要的影响。我们将对催化剂组合以及温度条件进行分析，探索高效的C4 烯烃制备工艺。 对于问题一，我们将附件一中的所给的数据，运用Matlab软件，绘制出在A1-A14,B1-B7的影响下，乙醇转化率、C4烯烃的选择性与温度的关系图表。通过观察在温度变化下折线的走向趋势，我们发现部分图像存在拐点。这说明当温度较低时
如何写好数学建模竞赛答卷
yillc的专栏
09-04 1531
如何写好数学建模竞赛答卷         一、写好数模答卷的重要性     1.  评定参赛队的成绩好坏、高低，获奖级别， 数模答卷，是唯一依据。     2. 答卷是竞赛活动的成绩结晶的书面形式。     3. 写好答卷的训练，是科技写作的一种基本训练。
2021年数模国赛B-乙醇偶合制备 C4 烯烃问题-二等奖
06-22
2021年数模国赛-乙醇偶合制备 C4 烯烃问题-二等奖，博主版权所有，源码私聊获取，附录代码全。 本文针对乙醇偶合制备 C4 烯烃问题，利用 Newton 插值、多元线性回归，BP 神经网络 与粒子群算法进行研究，旨在解决工业中如何运用乙醇更高效地生产 C4 烯烃的问题。 针对问题一，运用 Newton 插值，在样本量较小，且不知道具体函数表达式的情况下， 可以运用插值函数表示自变量与因变量之间的函数关系，为更好地表现温度与乙醇转换率和 C4 烯烃转换率之间的关系，建立了 Newton 多项式：
cumcm2021B附件5_逆变器参数价格.pdf
10-02
cumcm2021B附件5_逆变器参数价格.pdf
第五讲：2021年国赛B题-乙醇制备C4烯烃优化方案
热门推荐
a4646642的博客
09-21 1万+
2021年国赛B题-乙醇制备C4烯烃原题 第五讲：国赛B题第一小问讲解 2021年国赛B题-乙醇制备C4烯烃数据： 数据一： 数据二： 开始建模 观察此题，就是一个统计推断问题，采取回归思想以及方差分析以及差值算法以及相关分析即可解决。            &nbsp
真题解析 | 2021国赛B题：乙醇偶合制备 C4 烯烃
weixin_44099072的博客
08-18 7853
针对第一小题，需要注意的是，对催化剂组合应该是一一分析的，例如对于催化剂组合编号A1，可以画出来乙醇转化率、C4 烯烃的选择性与温度的三维图进行空间分析，也可以画出来二维图进行波动分析，然后再用些定量的指标进行描述性分析，例如用皮尔逊相关系数分析其相关性。...
A5.2021年全国数学建模竞赛B题-赛题分析与评阅要点（乙醇偶合制备C4烯烃分析）
youcans的博客
09-08 7846
2021年全国数学建模竞赛 2021年 B题（乙醇偶合制备 C4 烯烃）赛题解读分析。 本文首发于 2021年9月8日，内容为 2020 年 B题；9月13日修改，内容为 2021年 B题赛题；9月15日起更新，进行 2021年 B题的分析。 [『Python小白的数学建模课 @ Youcans』](https://blog.csdn.net/youcans/article/details/117295222)带你从数模小白成为国赛达人。
2021年数学建模国赛B题优秀论文（Word）（04烯焼制备分析与试验设计）
hjc的博客
08-17 9583
​ 以下为部分，下载链接：2021年B题论文 ​
2021 年高教社杯全国大学生数学建模竞赛 B题 乙醇偶合制备 C4 烯烃 求解思路1
Lu的专栏
11-03 6826
由附件1“性能数据表”中A1~A14（用装料方式I）整理得下表“sheet2”： 表1 性能数据表 wt%Co/SiO2 Co/SiO2 和 HAP 的质量比 催化剂总量 乙醇浓度 温度 乙醇转化率(%) C4烯烃选择性(%) C4烯烃收率(%) 0.5 1 400
MATLAB练习题：数学建模国赛2021年B题《乙醇偶合制备 C4 烯烃》数据预处理
最新发布
数学建模学习交流的博客
01-19 758
首先观察每一行数据的结构，可以看出每一行要提取的指标由空格隔开，因此我们可以考虑使用split函数在空格处对每一行进行拆分，拆分后每一行变成了多个子字符串；最后，我们会把所有提取出来的数值指标整合到一个8行4列的矩阵M中，每行对应一个催化剂组合的四个指标，从而达成我们的分析目标。现在我们需要从s2的每一行中提取上述介绍的四个指标，并将提取的结果保存到一个8行4列的数值矩阵M中。我们的任务就是从提供的催化剂组合中提取关键的数值指标。​讲解视频：可以在bilibili搜索《乙醇浓度2.1ml/min。
国赛 mathorcup数学建模竞赛中的优化问题
Anoushka的博客
07-20 2096
数学建模优化类题目的解题思路主要包括确定目标和约束、建立模型、选择优化方法、实现算法、评估结果和解释应用
【数学建模】CUMCM-2016A 系泊系统的设计 解题思路整理
05-19
题目描述： 一艘质量为$m$的船停靠在距离海岸线$L$的海域， 用一根长为$L$的铁锚将其系住。现有一条绳子， 一端固定在海岸上某点， 另一端系在船上。 要求绳子始终保持水平， 即绳子与水面的夹角为90°。 现在需要设计一个系泊系统，使得满足以下要求： 1. 绳子的长度为$L$； 2. 系泊系统能够承受风浪力矩$mgl\sin\theta$（$\theta$为船在水平方向上的倾斜角度，$g$为重力加速度，$m$为船的质量，$l$为船到绳子固定点的水平距离）； 3. 系泊系统的总重力不超过$G$，其中包括锚和绳子的重力； 请设计一个系泊系统，确定锚点在海岸线上的位置，以及绳子的长度。 思路： 1. 通过分析，可以得到系泊系统的主要受力情况：重力，浮力，张力和风浪力矩。其中，重力和浮力是与锚点位置无关的，可以先计算出来；风浪力矩是与船的倾斜角度有关的，需要在计算时考虑。 2. 假设锚点位置为$(0,0)$，船的位置为$(x,y)$，则绳子的长度$L$可以表示为$\sqrt{x^2+y^2}$。考虑绳子的张力，由于绳子与水面垂直，所以绳子的张力只有水平方向的分量，即$T\sin\theta$，其中，$\theta$为绳子与水平面的夹角。 3. 计算重力和浮力。船的重力为$mg$，浮力为$\rho Vg$，其中，$\rho$为水的密度，$V$为船的体积。绳子和锚的重力可以通过绳子的长度$L$和锚的重量来计算，即$g(L-x)\rho_s A$，其中，$\rho_s$为锚的密度，$A$为锚的横截面积。 4. 考虑风浪力矩的影响。假设风浪力矩对船的影响可以视为一个力矩，大小为$mgl\sin\theta$，方向垂直于船的水平方向，根据右手定则可以确定其方向。由于绳子的张力与水平方向平行，所以风浪力矩只对船的倾斜角度$\theta$有影响。 5. 在计算时，可以将锚点位置和绳子长度作为变量，通过求解最优化问题，确定最佳的锚点位置和绳子长度，使得系泊系统满足要求。可以考虑使用数值优化方法进行求解。 参考代码： ``` import numpy as np from scipy.optimize import minimize def objective(x): L = np.sqrt(x[0]**2 + x[1]**2) theta = np.arctan(x[1] / x[0]) T = (m*g + rho*V*g + (L-x[0])*rho_s*A) / np.sin(theta) M = T * L * np.sin(theta) return -M def constraint1(x): L = np.sqrt(x[0]**2 + x[1]**2) theta = np.arctan(x[1] / x[0]) T = (m*g + rho*V*g + (L-x[0])*rho_s*A) / np.sin(theta) return G - (m + rho_s*A*(L-x[0]) + rho*V)*g - T*np.cos(theta) def constraint2(x): return x[0] def constraint3(x): return x[1] # 系泊系统参数 m = 1000 # 船的质量 V = 100 # 船的体积 rho = 1000 # 水的密度 rho_s = 7800 # 锚的密度 A = 10 # 锚的横截面积 g = 9.8 # 重力加速度 G = 100000 # 系泊系统的总重力 # 初始值 x0 = [10, 10] # 约束条件 cons = [ {'type': 'ineq', 'fun': constraint1}, {'type': 'ineq', 'fun': constraint2}, {'type': 'ineq', 'fun': constraint3} ] # 求解最优化问题 sol = minimize(objective, x0, method='SLSQP', constraints=cons) # 输出结果 print(sol) ``` 参考资料： [1] 赵凌宇. 数学建模算法与应用[M]. 北京：清华大学出版社，2016. [2] 谷超豪，李涛，刘宝印，等. 数学建模算法应用与分析[M]. 北京：科学出版社，2016.

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值