2021长三角数学建模竞赛思路

杨利霞

2021长三角数学建模竞赛思路

3 q% f0 V3 X% I' i2 ^  P4 EA题  Go! Fun游长三角
) u4 l' M4 m1 ^) Z; c- u
- }  R( g/ h8 x+ k
A题如果从旅游业相关企业去看，就不足为奇了，这道题是让我们分析不同类型的游客的旅行规划问题
# D# m& E" p; {2 e0 s3 T
; P  R3 k; x* @; |% B0 C

0 t4 ]" @$ {) e+ F2 h/ d
+ D3 i. K6 z- c9 k) z6 Q7 f
0 v5 |) U; m7 r/ Q: T6 l9 p
6 M; k/ k$ \- g4 V1 W, @
+ N/ r* @! D& |7 Q
6 A- G* p+ T% g! V# \
7 H/ ^5 k% X& O. }3 |2 p4 `  v

7 \% {0 o/ M2 K1 _1 L7 A  X$ a



1 l3 Z* u/ o4 k) m& Y
首先需要爬取不同景点信息，包括经纬度，是否收费，收费的门票多少，周围酒店多少价格，都可以通过地图网页爬取数据，公众号推文有python+selenium案例不会就手动整理。为了能够做第一问，这里可以自行设置一下每个景点中需要花费的时间。
. [$ e5 O$ G8 e% A6 C
9 N2 v# n& H8 r3 Z, E
( V9 J% Y9 d, i8 z! T* B, }6 l家境良好且平时学习刻苦的小李同学，现住于杭州，要到长三角旅行，杭州离上海比较近，行程定为2小时，一天旅游时间设定为8小时，当天需要返回酒店，如果当天没有足够时间去下一个景点则直接返回酒店，8小时包含回酒店所需得时间，那么第一问就可以这么做，长三角景区基本位于上海，那我们就分析上海的景区，一般旅游会综合考虑要玩的景点地理位置及酒店价格选择住宿。假设住宿不变，到了上海采用徒步的形式到达景点，一般人走路速度为1.1-1.5m/s，两点距离根据Haversine公式计算，一天吃喝按100元算。
. I# I4 o1 L4 G5 d0 }6 q% i
* ~+ v5 r$ s& x, f* B, x1 w2 q
. ?4 B& q" D7 v5 ?. n0 k第一问仅考虑旅行时间，构建多目标函数，函数1为旅游景点总数最大，函数2为到景点的路上行程最小，自变量1为酒店编号，自变量2为所有景点序列，可通过randperm实现，在内循环内依次遍历randperm产生的序列，直到满足条件终止，约束条件见上一段落，最后给出一个较好的方案即可。程序中需用矩阵时刻记录参数的变化，满足终止条件就输出结果。

. |6 I1 E8 z3 |
5 ~# N5 h; Q/ b% o( L+ e7 b; E第二问，不考虑时间，只考虑旅行资金，第一问模型中删除旅行时间约束，增加约束条件为旅行资金5000元
+ A) [/ C6 t, Q' R  T

第三问，同时考虑旅行资金和时间再算一遍
& p9 O1 y% ]% Z4 ~: \
- W. J/ D( P* y% l* J
第四问，小李的爷爷也来旅游了，时常怀念曾经的峥嵘岁月，于是跟着小李又来了长三角旅行，总资金变多了，但每天的吃喝消费和住宿费增加了，如果是有门票的景点，费用也增加了，带入第三问模型即可。

; l3 R( C( i, x% X1 g; [- U
. a% T9 g) u1 e* W/ k第五问，后来小李兼职了新媒体，B站粉丝数达52.5万人了，这次到长三角旅行，每天需要抽出1个小时旅行时间来更新视频，上述模型中约束条件8小时改为7小时，假设小李的初始资金为5000，视频收入在第二天到账，视频的播放量+点赞数+投币+收藏，每一万播放量考虑上点赞数+投币+收藏大概在30元左右收入，在本问可以到B站搜一下相关的景点视频，如果会爬虫的同学可以直接用程序搜集，主要抓取一些相关景点视频的播放量及up主的粉丝数，最后取平均值，按粉丝比例算一下小李如果旅游该景点的收益多少，同样的按上文模型求解，本问的区别在于小李的资金是动态变化的。
5 A+ C* [& x/ W) k9 G0 Z$ P7 z
9 l4 t, W- S& \* X# y
' G$ H7 v! ]' V  e
% |3 p5 I% T! e% b
; Q0 s( L) S' P  z
* J  y; d! K, E7 w% k* v" z: JB题  锅炉水冷壁温度曲线


附件1中十条曲线如下
& H$ r# h' m/ q8 N! L4 R, ^
: x$ d  {+ F: l5 U* T
) }3 P, d" i' F5 o# ^/ L: g) `

% P' b' g6 j  Q% {, \& }! f

- d4 t6 H0 s- |, g9 E' K

, \$ {4 @; [, T4 x3 x, r
( ^7 o- J4 z( q: B上图代码


- h+ ^4 V+ M) r! l# z9 {' K  y2 kX=xlsread('附件1.xlsx');
* }" g, r4 B' M  i% \6 m$ xfigure
for i=1:10
    subplot(5,4,2*i-1)
1 \  S$ `) Q- M" V3 A9 m    plot(X(:,1),X(:,i+1))
    title(['管道',num2str(i),'温度曲线'])
    subplot(5,4,2*i)
$ ?* z3 r$ |( w$ m- ^    histogram(X(:,i+1))
    xlabel('温度值')
2 [% O0 q: c2 |- k    ylabel('数据频数')
    title(['管道',num2str(i),'温度曲线统计'])
! h! Y. r2 X4 \  Rend
  W" s8 S# |) d& E4 P2 L* `第一问统计数据特征，统计变量一般有方差、均值和最大最小值，也可以做下histogram统计图说明下温度主要分布情况。


- {5 \( K4 `7 Z( G' Q, m第二问对十个水冷壁管道的工作状态进行评价，在实际生产过程中，温度变化尽可能平稳，水冷壁温度不宜过高，根据这两个条件，可以使用第一问两个指标，方差和最大值。根据这两个条件，确定方差和最大值均是越小越好，那么从数据大小的角度将这两个指标定义为负向指标，数据归一化时则需要注意减最小还是减最大，将数据归一化至[0,1]备用，评价算法可以用因子分析、秩和比、熵权法求得权重后乘以归一化数据求和。
9 x- |0 U- P5 r
6 `) v5 N  A' ]  \3 O/ V
# o  p, Q9 m! u5 R, l/ l第三问附件二中包含了111个操作变量和42个状态变量，涉及高纬度数据首先应当对数据降维，可以去看下公众号发布的主成分分析推文，最后取贡献率到刚好到95%以上的k列降维数据，注意降维数据不是指选出主要指标出来，降维数据表示的是153高维数据的投影数据，直接用这k列数据去分别拟合10个管道的温度数据，常用的回归肯定不行，这里是比较复杂的非线性拟合，可以将k列数据和某个管道温度数据先用1stopt遍历下较符合的公式，通过lsqcurvefit函数进行非线性最小二乘拟合。
9 _2 T7 o0 O/ F7 i

  _6 L2 E  s) W4 }- v" n第四问找出主要操作变量，注意是111个操作变量，求10个温度数据与111个操作数据的相关性，与是个温度数据的相关性取个平均值进行比较，选出最高的介个，算法可采用余弦相似度、皮尔逊等。
( J- _; v. ]! S, P3 C

第五问第四问中的相关性矩阵，取与第十个管道从3172节点以后的温度数据与111个操作数据求相关性，可以选5-10个操作变量备用，拟合出一个关系式，构建一个启发式算法，目标函数1为3172节点以后的数据总下降值，目标函数2为依然高于445的数据个数，目标函数3同样也是自变量为调节的操作变量个数，目标函数4为所有调节变量总下调比例；自变量1为调节的操作变量个数，自变量2设置每个操作变量调整比例为[0-0.2]，同样的对3172节点后的操作变量数据进行整体的调节，向上还是向下调节可通过rand随机。启发式算法可采用模拟退火框架和蒙特卡洛框架，外循环为个体间比较，内循环为蒙特卡洛试验，内循环根据当前的调节操作变量个数随机模拟n次变量调整比例下的函数值，也就是变量2放在了蒙特卡洛试验中，通过非支配排序选择结果最好的试验作为当前个体的变量及函数值，外循环的个体间比较也同样通过非支配排序，迭代结束，输出最优解集的parote图及排名靠前的几个方案。（算法可参照公众号算法推文）
2 d# C! s/ ~7 o2 n0 F) p————————————————
: [/ g5 \, a$ E( s& d' @1 C0 h版权声明：本文为CSDN博主「微信公众号：您好啊数模君」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
, d8 c, u% w1 o1 _& J% @% @原文链接：https://blog.csdn.net/qq_39899679/article/details/117093580

* u; U: s0 q. o

1357404354

da大佬，能给一份，参考论文么，最近正在研究这次的赛题
% m: H; `; F4 y! o

1357404354

mlgb

1357404354 发表于 2022-7-30 14:31

% I% o/ p' ~8 E7 r- w0 s% P2 N- h