【2022高教社杯数学建模】C题:古代玻璃制品的成分分析与鉴别 赛后总结 39页论文及代码

简介: 2022年高教社杯数学建模竞赛C题的研究成果

1 题目

丝绸之路是古代中西方文化交流的通道,其中玻璃是早期贸易往来的宝贵物证。早期的玻璃在西亚和埃及地区常被制作成珠形饰品传入我国,我国古代玻璃吸收其技术后在本土就地取材制作,因此与外来的玻璃制品外观相似,但化学成分却不相同。
玻璃的主要原料是石英砂,主要化学成分是二氧化硅(SiO2)。由于纯石英砂的熔点较高, 为了降低熔化温度,在炼制时需要添加助熔剂。古代常用的助熔剂有草木灰、天然泡碱、硝石和铅矿石等,并添加石灰石作为稳定剂,石灰石煅烧以后转化为氧化钙(CaO)。添加的助熔剂不同,其主要化学成分也不同。例如,铅钡玻璃在烧制过程中加入铅矿石作为助熔剂,其氧化铅(PbO)、氧化钡(BaO)的含量较高,通常被认为是我国自己发明的玻璃品种,楚文化的玻璃就是以铅钡玻璃为主。钾玻璃是以含钾量高的物质如草木灰作为助熔剂烧制而成的,主要流行于我国岭南以及东南亚和印度等区域。
古代玻璃极易受埋藏环境的影响而风化。在风化过程中,内部元素与环境元素进行大量交换,导致其成分比例发生变化,从而影响对其类别的正确判断。如图 1 的文物标记为表面无风
化,表面能明显看出文物的颜色、纹饰,但不排除局部有较浅的风化;图 2 的文物标记为表面风化,表面大面积灰黄色区域为风化层,是明显风化区域,紫色部分是一般风化表面。在部分风化的文物中,其表面也有未风化的区域。
在这里插入图片描述
在这里插入图片描述

图 1 未风化的蜻蜓眼玻璃珠样品 图 2 风化的玻璃棋子样品

现有一批我国古代玻璃制品的相关数据,考古工作者依据这些文物样品的化学成分和其他检测手段已将其分为高钾玻璃和铅钡玻璃两种类型。附件表单 1 给出了这些文物的分类信息,
附件表单 2 给出了相应的主要成分所占比例(空白处表示未检测到该成分)。这些数据的特点是成分性,即各成分比例的累加和应为 100%,但因检测手段等原因可能导致其成分比例的累加和非 100%的情况。本题中将成分比例累加和介于 85%~105%之间的数据视为有效数据。
请你们团队依据附件中的相关数据进行分析建模,解决以下问题:
问题 1
对这些玻璃文物的表面风化与其玻璃类型、纹饰和颜色的关系进行分析;结合玻璃的类型,分析文物样品表面有无风化化学成分含量的统计规律,并根据风化点检测数据,预测其风化前的化学成分含量。
问题 2
依据附件数据分析高钾玻璃、铅钡玻璃的分类规律;对于每个类别选择合适的化
学成分对其进行亚类划分,给出具体的划分方法及划分结果,并对分类结果的合理性和敏感性进行分析。
问题 3
对附件表单 3 中未知类别玻璃文物的化学成分进行分析,鉴别其所属类型,并对

分类结果的敏感性进行分析。
问题 4
针对不同类别的玻璃文物样品,分析其化学成分之间的关联关系,并比较不同类别之间的化学成分关联关系的差异性。

附件
表单 1 玻璃文物的基本信息
表单 2 已分类玻璃文物的化学成分比例,其中
(1)文物采样点为该编号文物表面某部位的随机采样,其风化属性与附件表单 1 中相应文物一致。
(2)部位 1 和部位 2 是文物造型上不同的两个部位,其成分与含量可能存在差异。
(3)未风化点是风化文物表面未风化区域内的点。
(4)严重风化点取自风化层。
表单 3 未分类玻璃文物的化学成分比例

2 摘要及截图

当下,随着现代考古技术的发展、理化分析水平的显著提升,相关领域对古代玻璃制品的主要化学组成成分的研究有了极大的进步。针对古代玻璃制品分类特征信息、化学组成成分之间的关系及其统计学规律等问题进行了探讨,建立了对应分析、Kmeans++ 聚类、随机森林等模型,解决了古代玻璃制品的分类、鉴别及其化学组成成分间的相关 关系分析等问题。

针对问题一, 依据表单 1、2 给出的数据进行预处理,结合玻璃文物分类特征信息 及对应的频数,建立了对应分析模型,总结了文物的表面风化与其玻璃类型、纹饰和颜 色的相关关系。进而又通过对表单数据的运算处理,得到了文物样品表面有无风化与其 组成成分含量的统计规律。最后,根据风化点的检测数据,本文利用多元线性回归模型 预测得出了相应的风化前的化学成分含量。

针对问题二, 首先依据预处理后的数据进行统计分析,总结得出了高钾玻璃、铅钡玻璃的分类规律。之后,建立了 Kmeans++聚类模型,并对其聚类数 K 的取值做了分析,得到两种类型的玻璃样本聚类时,K 值都取 3,即在两个类别的玻璃制品之下,又各自分出三个亚类,进而又分析了三个亚类之中各个特征的差异性,验证了肺类的合理性。

针对问题三, 运用了随机森林模型,随机划分 80%的数据作为训练集,20%的数据 作为测试集,进行求解。最终的鉴别结果分别为:高钾玻璃、铅钡玻璃、铅钡玻璃、铅钡玻璃、铅钡玻璃、高钾玻璃、高钾玻璃、铅钡玻璃。之后,选取了两个最大深度为 5 的决策树,对表单三的八件玻璃制品进行鉴别,求解结果与上述结果相同,验证了模型的准确度与敏感性。

针对问题四, 首先验证了两种玻璃类别其化学成分之间的关联关系近似呈线性相关,进而建立了皮尔逊相关系数分析模型,并通过对不同类别之间的化学成分关联关系进行比较,探究了其差异性。分析得出,铅钡玻璃制品在风化后的化学成分关联关系较强,而高钾玻璃制品在风化后的化学成分关联关系相对减弱。

关键词: 对应分析;多元线性回归;Kmeans++聚类;随机森林;皮尔逊相关系数
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述


3 文章及代码下载方式

方法一:https://www.betterbench.top/#/20/detail
方法二:在Github搜BetterBench-Shop,在此页面中ctrl+F查找“高教社杯”即可找到
在这里插入图片描述
在这里插入图片描述

若下载不成功,私信我,将文章链接发给我

目录
相关文章
|
7月前
|
存储 缓存 NoSQL
使用双 Bitmap实现签到和补签
在高并发场景下,采用双 Bitmap 实现签到与补签方案,具有内存占用低、操作效率高、并发安全等优势。通过 Redis 存储 Bitmap 数据,结合位运算实现签到记录与统计,可高效支撑大规模用户签到需求,适用于社交 APP、会员系统等场景。
|
索引 搜索推荐 缓存
使用Elasticsearch进行高效全文搜索的技术探索
【6月更文挑战第3天】本文探索了使用Elasticsearch进行全文搜索的技术,它基于Lucene,是一款开源、分布式搜索引擎。核心原理在于倒排索引,实现快速查找。全文搜索涉及安装配置、创建索引、索引文档、执行查询及处理结果。为了优化性能,可以调整索引映射、选择合适分析器、利用缓存及优化硬件配置。Elasticsearch因其高效性和灵活性在大数据时代广泛应用。
538 62
|
存储 Java 编译器
JNI内存方面说明以及相关类型手动释放内存
JNI内存方面说明以及相关类型手动释放内存
1625 0
|
API 开发工具 Android开发
简述大疆无人机对接
【2月更文挑战第7天】本文介绍了对接大疆无人机的主要目的,包括实时画面获取、飞行数据监测、操控飞行、媒体管理和业务功能开发等,并列举了多种开发接口如MobileSDK、UXSDK、云开发API等。重点讨论了MobileSDK在Android平台的应用,包括SDK集成步骤、直播推流和获取飞机实时数据的细节。另外,UXSDK用于加速应用开发,提供预设UI组件。上云API则简化了无人机与第三方云平台的集成,支持MQTT、HTTPS和WebSocket协议,适用于行业级无人机。对接流程涉及Pilot2和Dock的配置,以及数据传输和业务功能处理。文章还提及了如何对接多个飞机的方法。
12494 0
简述大疆无人机对接
|
数据可视化 算法 定位技术
Python数据可视化matplotlib和pyecharts参数详解
Python数据可视化matplotlib和pyecharts参数详解
701 0
|
算法
数学建模国赛/美赛常见赛题类型及建模方案(纯干货)
主客观概念主要是在指标定权时来划分的。主观评价与客观评价的区别是,主观评价算法在定权时主要以判断者的主观经验为依据,而客观评价则主要基于测量数据的基本特性来综合定权
1277 0
数学建模国赛/美赛常见赛题类型及建模方案(纯干货)
|
IDE 开发工具
Vscode的远程开发之VScode优势(一)
Vscode的远程开发之VScode优势(一)
354 1
|
数据采集 存储 Web App开发
Python-数据爬取(爬虫)
【7月更文挑战第15天】
1027 3
|
C++
【PTA】​ L1-009 N个数求和​ (C++)
【PTA】​ L1-009 N个数求和​ (C++)
739 0
【PTA】​ L1-009 N个数求和​ (C++)
|
网络协议 Linux iOS开发
《熬夜整理》保姆级系列教程-玩转Wireshark抓包神器教程(2)-Wireshark在Windows系统上安装部署
【2月更文挑战第2天】《熬夜整理》保姆级系列教程-玩转Wireshark抓包神器教程(2)-Wireshark在Windows系统上安装部署
629 4