2023年辽宁省数学建模竞赛B题思路详细分析

2023-11-04 13:01

本文主要是介绍2023年辽宁省数学建模竞赛B题思路详细分析,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

摘要略,2023年辽宁省数学建模竞赛B题代码和论文已经完成,代码为全部3问代码,论文包括摘要、问题重述、问题分享、模型假设、符号说明、模型的建立和求解(问题1无监督聚类模型的建立和求解,问题二有监督分类预测模型的建立和求解、问题三迁移预测模型的建立与适用性评估)、模型的评价、参考文献、附录,详情参考见文末

二、 问题重述

2023辽宁省数学建模竞赛B题:“海洋强国”战略部署已成为推动中国现代化建设的重要组成部分,国家对此提出“发展海洋经济,保护海洋生态环境,加快建设海洋强国”的明确要求。 《辽宁省“十四五”海洋经济发展规划》明确未来全省海洋经济的发展战略、发展目标、重大任务、空间部署和保障措施。规划范围包括辽宁省全部海域和大连、丹东、锦州、营口、盘锦和葫芦岛6 个市以及海洋经济发展所依托的相关陆域,规划期限为2021 年至2025 年,展望到2035 年。辽宁省作为中国最北沿海省份,拥有2292.4 公里海岸线(如图1 所示)。

在“海洋强省”建设目标的背景下,完成海洋经济发展规划的目标重在海洋高新技术领域创新。其中关键核心技术之一是攻克水下导航与定位的适配区分类预测技术。

水下航行器在执行水下任务时需要保持自主、无源、高隐蔽性、不受地域和时域限制、高精度的导航与定位。重力辅助导航是满足上述条件的主要方法之一。 在重力辅助导航系统中,影响导航可靠性与精度的关键步骤是选择匹配性高的航行区域,即适配区。适配区的标定与识别技术是最具挑战性的问题之一。选取适配区前需要对研究海域的重力基准图(基础性的是重力异常基准图)进行插值加密处理,基于重力基准图所提供水下航行器航行区域的重力异常变化情况对适配区的选取进行分析。

重力异常(值)的定义为:实际地球内部的物质密度分布不均匀,导致实际观测重力值与理论上的正常重力值总存在偏差,在排除各种干扰因素影响后,仅仅由地球物质密度分布不匀所引起的重力的变化,简称为重力异常。

在重力异常变化显著区域,导航系统可获得高的定位精度;反之,在重力异常变化平坦区域,导航系统会出现定位精度的不敏感。由于不同区域的重力异常特征分布不同,建立可行的适配区分类预测模型,对保障水下航行器的导航精度至关重要。

假设X 为影响区域匹配性的特征属性指标,Y 为刻画区域适配性的输出结果, F 为以X 为输入以Y 为输出的分类预测系统。

基于上述背景分析,请参考附件中的重力异常数据建立数学模型,解决以下问题:

问题一:附件1,给出一组分辨率为1’×1’(相邻两格网点间的距离是1’)的重力异常基准数据A,试通过精细化基准图,合理划分区域,完成各区域的适配性标定(标签Y )。

问题二:根据问题一中各划分区域的适配性标定结果Y ,合理选择区域的特征属性指标(特征X ),试建立有效的区域适配区分类预测模型(系统F )。

问题三:利用附件二中的重力异常基准数据B,试对问题二所建立的系统F进行迁移性预测并讨论该系统F 对新重力异常数据的适用性。

目录:(完整代码和论文见文末)

一、 摘要.... 1

二、 问题重述.... 3

三、 问题分析.... 4

3.1 问题一分析..... 4

3.2 问题二分析..... 4

3.3 问题三分析..... 4

四、 模型假设.... 5

五、 符号说明.... 6

六、 模型的建立和求解.... 6

6.1 问题一模型的建立和求解.... 6

6.1.1数据可视化.... 6

6.1.2 K均值聚类方法:.... 7

6.1.3 聚类求解结果.... 10

6.2 问题二模型的建立和求解.... 11

6.2.1 随机森林算法.... 11

6.2.2 决策树算法.... 12

6.2.3 SVM算法.... 13

6.2.4 使用其他多种分类方法.... 16

6.2.5 不同模型的比较.... 17

6.3 问题三模型的建立和求解.... 19

6.3.1 重力异常数据进行可视化.... 19

6.3.2 迁移预测.... 20

6.3.3 重力异常数据的适用性评估.... 21

七、 模型的评价.... 22

7.1 模型的优缺点.... 22

三、 问题分析

要建立一个数学模型来解决2023辽宁省数学建模竞赛B题,可以采用机器学习或统计分析的方法来完成重力异常数据的区域适配性标定。步骤如下:

3.1 问题一分析

数据准备:

导入附件1中的数据(经度、纬度、重力异常值)。 可能需要对数据进行预处理,如数据清洗,处理缺失值,进行标准化等。 特征工程:

可以考虑根据经度和纬度创建额外的特征,如距离某个地理坐标点的距离等。

标定区域适配性标签Y:

为了标定区域的适配性标签Y,需要定义一种方法来将数据点分为不同的区域或类别。这是一个无监督学习问题,可以考虑使用聚类分析(如K均值聚类)来将数据点分成不同的区域。

3.2 问题二分析

分类预测模型建立:

选择合适的机器学习模型,如决策树、随机森林、支持向量机、神经网络等,来预测区域适配性标签Y。 划分数据集为训练集和测试集,用训练集来训练模型,然后用测试集来评估模型的性能。 模型评估:

使用合适的性能指标(如准确度、召回率、F1分数等)来评估模型的性能。 可能需要进行交叉验证以确保模型的稳定性和泛化能力。

3.3 问题三分析

模型应用:

一旦模型建立并评估成功,可以使用该模型来预测附件2中的重力异常数据的区域适配性标签Y。

四、 模型假设

五、 符号说明

以下是问题1、问题2和问题3中需要使用的模型符号及其说明的总结

六、 模型的建立和求解

6.1 问题一模型的建立和求解

6.1.1 数据可视化

使用Python的matplotlib库对2023辽宁省数学建模竞赛B题数据驱动的水下导航适配区分类预测问题一进行数据可视化,然后使用K均值聚类进行数据点分区。

# 设置中文字体
plt.rcParams['font.sans-serif'] = 'SimHei'  # 使用中文字体(这里使用的是“黑体”)
plt.rcParams['axes.unicode_minus'] = False  # 用于正常显示负号# 可视化数据
plt.scatter(data['经度'], data['纬度'], c=data['重力异常值'], cmap='viridis')
plt.xlabel('经度')
plt.ylabel('纬度')
plt.title('重力异常数据分布')
plt.colorbar(label='重力异常值')
plt.savefig('gravity_anomaly_visualization.png', dpi=300)
plt.show()

6.1.2 K均值聚类方法

K均值聚类是一种无监督学习方法,用于将数据点分成K个不同的簇,以使每个数据点属于最近的簇。该方法的目标是最小化数据点与其所属簇中心之间的距离的总和。K均值聚类的算法如下:

选择K个初始簇中心,可以随机选择或根据某种启发式方法选择。

将每个数据点分配到与其最近的簇中心。

更新每个簇的中心,计算每个簇中所有数据点的平均值。

重复步骤2和步骤3,直到簇中心不再发生显著变化或达到预定的迭代次数。

(完整代码见文末)

# 使用K均值聚类
X = data[['经度', '纬度', '重力异常值']]
scaler = StandardScaler()
X_scaled = scaler.fit_transform(X)  # 标准化数据# 设置聚类数(这里设置为3,你可以根据实际情况调整)
kmeans = KMeans(n_clusters=3, random_state=0)# 保存包含区域标签的数据
data.to_excel('附件1_重力异常数据带区域标签.xlsx', index=False)

6.1.3 聚类求解结果

# 可视化不同区域的数据
plt.scatter(region1['经度'], region1['纬度'], label='区域1', c='red', marker='o')
plt.scatter(region2['经度'], region2['纬度'], label='区域2', c='blue', marker='s')
...plt.xlabel('经度')
plt.ylabel('纬度')
plt.title('K均值聚类结果可视化')
plt.legend()
plt.savefig('kmeans_clustering_result.png', dpi=300)
plt.show()

上述代码将可视化数据分布,并将其保存为名"gravity_anomaly_visualization.png"的图像。然后,它使用K均值聚类将数据点分为3个不同的区域,并将带有区域标签的数据保存为一个新的Excel文件。可以根据需要更改聚类数或其他参数以获得更合适的结果 。

1.1 问题二模型的建立和求解

2023辽宁省数学建模B题问题二:根据问题一中各划分区域的适配性标定结果Y ,合理选择区域的特征属性指标(特征X ),试建立有效的区域适配区分类预测模型(系统F )。

为了建立一个有效的区域适配区分类预测模型,可以采用监督学习方法,例如决策树、随机森林、支持向量机等,来预测区域的适配性标定结果Y

完整论文和代码:2023辽宁省数学建模高质量原创论文代码成品分享

这篇关于2023年辽宁省数学建模竞赛B题思路详细分析的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/345306

相关文章

uva 10014 Simple calculations(数学推导)

直接按照题意来推导最后的结果就行了。 开始的时候只做到了第一个推导,第二次没有继续下去。 代码: #include<stdio.h>int main(){int T, n, i;double a, aa, sum, temp, ans;scanf("%d", &T);while(T--){scanf("%d", &n);scanf("%lf", &first);scanf

uva 10025 The ? 1 ? 2 ? ... ? n = k problem(数学)

题意是    ?  1  ?  2  ?  ...  ?  n = k 式子中给k,? 处可以填 + 也可以填 - ,问最小满足条件的n。 e.g k = 12  - 1 + 2 + 3 + 4 + 5 + 6 - 7 = 12 with n = 7。 先给证明,令 S(n) = 1 + 2 + 3 + 4 + 5 + .... + n 暴搜n,搜出当 S(n) >=

uva 11044 Searching for Nessy(小学数学)

题意是给出一个n*m的格子,求出里面有多少个不重合的九宫格。 (rows / 3) * (columns / 3) K.o 代码: #include <stdio.h>int main(){int ncase;scanf("%d", &ncase);while (ncase--){int rows, columns;scanf("%d%d", &rows, &col

透彻!驯服大型语言模型(LLMs)的五种方法,及具体方法选择思路

引言 随着时间的发展,大型语言模型不再停留在演示阶段而是逐步面向生产系统的应用,随着人们期望的不断增加,目标也发生了巨大的变化。在短短的几个月的时间里,人们对大模型的认识已经从对其zero-shot能力感到惊讶,转变为考虑改进模型质量、提高模型可用性。 「大语言模型(LLMs)其实就是利用高容量的模型架构(例如Transformer)对海量的、多种多样的数据分布进行建模得到,它包含了大量的先验

【生成模型系列(初级)】嵌入(Embedding)方程——自然语言处理的数学灵魂【通俗理解】

【通俗理解】嵌入(Embedding)方程——自然语言处理的数学灵魂 关键词提炼 #嵌入方程 #自然语言处理 #词向量 #机器学习 #神经网络 #向量空间模型 #Siri #Google翻译 #AlexNet 第一节:嵌入方程的类比与核心概念【尽可能通俗】 嵌入方程可以被看作是自然语言处理中的“翻译机”,它将文本中的单词或短语转换成计算机能够理解的数学形式,即向量。 正如翻译机将一种语言

基于UE5和ROS2的激光雷达+深度RGBD相机小车的仿真指南(五):Blender锥桶建模

前言 本系列教程旨在使用UE5配置一个具备激光雷达+深度摄像机的仿真小车,并使用通过跨平台的方式进行ROS2和UE5仿真的通讯,达到小车自主导航的目的。本教程默认有ROS2导航及其gazebo仿真相关方面基础,Nav2相关的学习教程可以参考本人的其他博客Nav2代价地图实现和原理–Nav2源码解读之CostMap2D(上)-CSDN博客往期教程: 第一期:基于UE5和ROS2的激光雷达+深度RG

CSP 2023 提高级第一轮 CSP-S 2023初试题 完善程序第二题解析 未完

一、题目阅读 (最大值之和)给定整数序列 a0,⋯,an−1,求该序列所有非空连续子序列的最大值之和。上述参数满足 1≤n≤105 和 1≤ai≤108。 一个序列的非空连续子序列可以用两个下标 ll 和 rr(其中0≤l≤r<n0≤l≤r<n)表示,对应的序列为 al,al+1,⋯,ar​。两个非空连续子序列不同,当且仅当下标不同。 例如,当原序列为 [1,2,1,2] 时,要计算子序列 [

数学建模笔记—— 非线性规划

数学建模笔记—— 非线性规划 非线性规划1. 模型原理1.1 非线性规划的标准型1.2 非线性规划求解的Matlab函数 2. 典型例题3. matlab代码求解3.1 例1 一个简单示例3.2 例2 选址问题1. 第一问 线性规划2. 第二问 非线性规划 非线性规划 非线性规划是一种求解目标函数或约束条件中有一个或几个非线性函数的最优化问题的方法。运筹学的一个重要分支。2

三相直流无刷电机(BLDC)控制算法实现:BLDC有感启动算法思路分析

一枚从事路径规划算法、运动控制算法、BLDC/FOC电机控制算法、工控、物联网工程师,爱吃土豆。如有需要技术交流或者需要方案帮助、需求:以下为联系方式—V 方案1:通过霍尔传感器IO中断触发换相 1.1 整体执行思路 霍尔传感器U、V、W三相通过IO+EXIT中断的方式进行霍尔传感器数据的读取。将IO口配置为上升沿+下降沿中断触发的方式。当霍尔传感器信号发生发生信号的变化就会触发中断在中断

每日一题|牛客竞赛|四舍五入|字符串+贪心+模拟

每日一题|四舍五入 四舍五入 心有猛虎,细嗅蔷薇。你好朋友,这里是锅巴的C\C++学习笔记,常言道,不积跬步无以至千里,希望有朝一日我们积累的滴水可以击穿顽石。 四舍五入 题目: 牛牛发明了一种新的四舍五入应用于整数,对个位四舍五入,规则如下 12345->12350 12399->12400 输入描述: 输入一个整数n(0<=n<=109 ) 输出描述: 输出一个整数