常见数学建模方法学习总结及资料汇总

杨利霞

常见数学建模方法学习总结及资料汇总
3 X  h! Q' ?& S6 \数学建模方法
前言
层次分析法（AHP）
( g5 I7 q' \( ?& G) e9 {多属性决策模型
2 c' c% r- o8 ~6 o/ w1 U) l模拟退火模型
种群竞争模型
遗传算法模型
线性规划模型
非线性规划模型
6 Y! i; L+ y8 `+ @9 ?PCA主成分分析
+ Z& P7 ]2 L- {9 S: a聚类分析
多元回归分析
5 l/ V3 I! @/ u, j前言
这是我在进行初级的数学建模时进行的一些原理学习，网上资料显得比较杂乱无章，这里我对一些数学建模方法的文档和视频资料进行了收集并汇总了以下写的好的，在学习时跟着这些资料学习并看看我的汇总就可以达到一个比较好的效果啦~
$ X. N! ~( ]" {

层次分析法（AHP）
应用：决策。（择优选择方案、优劣等级划分…）
$ G% B/ \* Y1 p% B

原理：根据问题的性质，使问题归结为最低层（供决策的方案、措施）到中间层（指标），中间层到最高层（总目标）的相对权值，通过一系列计算，来确定优劣性。
其中，中间层到最高层的相对权值，表示考虑方案的不同指标对最终结果选择的权重。最底层到中间层的相对权值，表示单个指标下，表现多个方案的好坏对比的权重。
# y. Y/ D5 Q* D) I( d9 ]
9 J# S% ^5 i! L. q  t4 f0 r
" c, o, h0 ?- O1 o8 c基本步骤：
8 i4 p" Y+ d$ C9 E; X1、划分层次（层次模型）：分为目标——准则层（指标）——方案层。
2 b, v7 I+ R( L  _" Y2、构造判断矩阵（又称成对比较矩阵）。根据1~9尺度进行相应填写。
3、一致性检验与求权重矩阵 。CR（一致性比率）<0.1则成对比较矩阵通过。一致性检验是检验判断构造的矩阵是否合格。权重矩阵是用来计算最终方案权值的。
9 D+ H/ D" w) C" l+ j! E+ _4、根据权重矩阵求最优方案。


% G! n- K4 t5 T- V, H视频学习：
: @% A7 c7 F% q; G" L6 Q# T9 R先看：
- z- W+ ^" }% b9 r链接: AHQ初步学习
后看: AHQ加深理解
* D6 b3 M2 p& l1 \1 p2 o/ a
' R0 ]9 d2 \2 h( t
多属性决策模型
应用：决策。


原理：类似于一般的加权求总分。主要思想是把量纲不同的数值转化为同一个量纲，然后构建权重矩阵，计算加权分值进行决策。


基本步骤：
" r. d7 i! L* X加权算数平均算子，属性值类型归一化。
1、根据属性类型，划分为效益型、成本型等，对决策矩阵（属性值构成的矩阵）进行归一化。
# _8 i, m: d& X5 v9 A. _, G2、构建成对比较矩阵，得到各个属性指标权重矩阵。（同AHP，也要进行一致性检验）
3、根据权值矩阵与归一化的决策矩阵确定每个方案的最终权值，进行决策。
# l0 p$ J9 O- p" \' @. [9 d
8 e7 Q! i" N4 F' C' V0 U
文档资料：多属性决策模型
7 l& A, X2 W' D- r1 w& R' B% t' V) d$ G视频资料：多属性决策模型
3 D  ?; s2 Y" X" ]

层次分析法和多属性决策模型的区别：层次分析法针对于定性问题定量分析，多属性决策模型则不然。即如果“属性值”并不是一个能用“数值”很好地表示的值的话，采用AHP的1~9尺度是更优方案。如果“属性值”全都是适合的“数值”的话，就用多属性决策模型的方案。

1 r8 i. u3 d  m. @; G
5 i6 M' c: u3 [6 |! r  N5 {, n模拟退火模型
% u* f8 M9 n  }) M$ C5 c$ j6 s+ I应用：无法确定能否在规定时间内找到最优解的NP问题——旅行商问题（TSP，访问每一座城市并回到最初城市的最短回路）；背包问题（在限定的总重量下，如何选择物品使总价值最高）；最大流问题（充分利用管道能力，使总流量最大）；图着色问题（用最少的颜色，对地图进行着色）；…


& r. ^- g0 K9 P3 E" h" w: ]7 w原理：一种寻找最优解的算法，传统的算法可能被局限于局部最优解，模拟退火算法给出一定概率使计算中能够跳出局部最优解继续求解，当求出更优解时可以得到答案。

, l3 K& _1 b* F4 q( b+ C
* L$ R) ~' W  u* Q9 f基本步骤：
' `# W% s$ ^" ^+ t& ?1、初始化。（T温度，S初始解状态，迭代次数）
2、如果计算的最新状态趋近于最优解，则把这个最优解作为下一次计算的起点。
' X% ?3 c& A1 ^" c( E3、如果计算的最新状态趋远于最优解，则根据方案取一定概率，使新的状态作为下一次计算的起点。
7 _9 u+ Q$ ~+ |- v4、重复1、2、3步骤，直至达到平衡状态并满足条件。


3 Y1 y) A2 M1 E7 }, z文档资料：
; q: n5 t7 W* U, `% P模拟退火模型通俗理解
' R2 S9 m3 R: x模拟退火模型加深理解

! B$ N6 i( c) n& k; l. Z* ?
种群竞争模型
应用：种群之间的竞争，优胜劣汰；产品的市场竞争模型…

! {5 T  x, U, a( s7 d: Y; V) t  v
0 x# m7 x1 N1 }  v& S2 w! }& h0 O原理：当两个物种共同生存，相互会起到阻滞作用时建立的模型。种群竞争设置的参数为【初始数量，自然增长率，环境容量，种群对资源的消耗量】。
% x: T' ?1 ]6 s) @1 M
4 K) F$ X% L  l
基本步骤：无。这个模型的资料大多都很简略。

( o0 T# @3 S/ ?; U9 [2 F
' @' u- l0 n, C# }$ v6 ?1 b文档资料：种群竞争模型
1 q  |' }, T  H5 Y' k+ }
0 @7 _' h- ]5 w
遗传算法模型
应用：对于同一个物种，有着个体之间性能的差异，遗传算法针对的是物种内的优胜劣汰。
" c) Q0 u/ R) g2 m) \5 d) x& U- A! }

原理：建立遗传算法模型，通过N代的遗传、交叉、变异、复制，找到问题的最优解。
基本步骤：
7 P& Q. \9 o; q/ Z2 W
% y5 @" n) E: q. b& c
1、根据初始化种群个体来初始化编码和数据。
' e4 W: N. l8 R% y- b2、评估个体的适应度，作为优胜劣汰的评价标准。
3 o. A( x2 h1 E; k4 ?# C: r3、通过选择、交叉、变异、复制（为了保证最优个体能存活，对当代最优个体直接复制到下一代）来决定下一代的种群群体。
4、重复步骤1~3，最好的迭代结果次数是种群群体性能趋于稳定时。


文档资料：
2 v/ _0 _7 T) E+ T0 \( L遗传算法通俗理解
遗传算法加深理解
0 V. P  [5 _/ g$ ?/ [+ T8 [* }: n" c

4 X2 n/ L$ |% q" `线性规划模型
应用：对于有线性限制的问题求最优解。
. M# p2 n0 t9 h+ ?$ j4 s
% h$ l! W8 s( x1 a6 `
3 `, N3 l5 k# u原理：根据限制建立线性关系，根据决策函数求最优解。高中的东西了。如果存在最优解，只能在可行域的边界找到（一般为顶点）。

# P$ c  Y9 H) E+ K
基本步骤：
9 d4 m+ U6 ?. m3 G, u& E1、假设/寻找决策变量。（影响目标结果且有条件限制的变量）
2、建立目标函数。（需要求最优解的东西）
5 X$ m- w. E+ |: \7 D3、寻找约束条件。（变量之间的约束关系）
4、求解。
. W& J) S' U2 s0 `3 [  W
6 D9 `( ^! }% t# |- k( |& s
- ]. \3 `' g, x! J; w3 G文档：线性规划
6 h! X! a3 ]2 L5 a$ i

: j& h: @0 Y- v) v' t非线性规划模型
7 Q. ~* k# z" r. y$ [% B( M/ v" |应用：对于有非线性限制的问题求最优解（约束函数或者目标函数含有非线性关系）。

  S0 `8 D, X& r
: y6 L/ s+ F% k& C  J  F6 \) v+ b原理：解决非线性规划模型，这种模型没有通用的解决方案，各个方法都有自己适用的范围。最优解可能在可行域的任意一点找到。通常非线性规划模型可以概括为：

# R6 N0 N+ \8 s% _0 W# _  ]  d7 u
其含义为：优化目标为f(x)，受限于一个等式h(x)与一个不等式g(x)。
( @" A( W& q" _7 v# |* C6 S

基本步骤：
: x' z. Z6 _1 o! v: k" a# q1、同线性规划，找出决策变量；建立目标函数；寻找约束条件。
$ c) f2 @* B( x5 W2 L; V2、求解过程中，选定初始点x，构建点x在趋于最优解的搜索方向，以该搜索方向为基础，设定搜索步长，找到下一个迭代点。
3、以新迭代点为初始点，重复进行2步骤。


- Z$ [! H1 S* J. B: s$ Y文档资料（比较多，稍微看看就好）：非线性规划
6 q/ F1 k4 B; b/ m$ W1 r

PCA主成分分析
应用：降维。

. T+ b; S% w2 t, i1 E
5 Z& m. i7 l/ v/ H1 X原理：在损失较少信息的基础上实现降维。周志华西瓜书上提到的最近重构性和最大可分性，实际上是关于协方差与方差的问题。协方差代表维度之间的相关性，方差可以表示维度可以保留的信息。故PCA降维的优化目标为：在N维中选择R维作为“主成分”，这R维之间两两之间协方差尽量小，方差尽量大。


基本步骤：涉及线性代数、矩阵分析。文档资料中有兴趣自行了解。
# f3 `; F( B2 t2 s6 \7 h- }
" Q. y9 D! L( u) z  ^- H& G
0 M  A+ i! t3 O9 C- P# ^& C; J" [! t文档资料：
PCA理解与基本概念
PCA维度的选择


( u7 t# q( F8 E  c! z, o聚类分析
4 \4 p$ S( g% D0 l$ t应用：无监督分类。
5 t2 ?. |8 L  z; w5 y
" _, U0 o* U8 c  \2 r
! B8 @. Z& |& e9 m% A1 _8 M原理：对没有标示性的数据，按照算法，进行簇划分。


基本步骤：
聚类算法分很多类别。最常使用是基于距离和基于密度的聚类方案，具体方案选择结合实际。
( O. E  S$ A" a2 G- Z  ^7 v在基于距离聚类中比较典型的是k-means聚类，通过多次迭代更新中心点，找到多个中心点（簇的中心）与多个簇，实现：簇内点紧凑，簇间分离程度高。需要设置K（簇数）的值，效率高。
基于密度的聚类有DBSCAN和MeanShift（均值漂移）等，主要思路是根据数据点的紧凑程度，将聚合程度较高的数据形成簇。
基于密度的聚类相比于基于距离的聚类，可以找到数据点分布呈不规则形状的簇，效率要差一些，相比之下，基于距离的聚类，同一个簇中各个数据点都不会太偏离中心点，效率高。
; ~% E: }" l5 d* |5 p

鸢尾花分类的聚类代码，用来辅助学习与理解。
1
from sklearn import datasets
from sklearn import model_selection
0 _0 E6 L7 G, y5 K7 `2 b4 lfrom sklearn.cluster import KMeans
+ H, j" H- @0 _1 N9 H: Iimport matplotlib.pyplot as plt

$ b1 u  K/ p7 P. {0 U
# F  w8 D! F5 k4 c# step1.加载数据
( }& z" V  k4 c6 V4 Riris = datasets.load_iris()

- ~9 \* j+ s" [- ]$ M
# 2.取特征空间中的4个维度
7 Q- N# O8 J* o( T% X7 g1 DX = iris.data[:,:4]  
1 j) L) ?) U' \; @) S* C
" ~0 r& ?1 S! _* g
+ [0 E0 T: D  O9 M( Y9 n! r/ K# 3.搭建模型，构造KMeans聚类器
estimator = KMeans(n_clusters=3)

1 y: F# C! q) f) }
#开始聚类训练   
estimator.fit(X)
, \* w2 O8 g' n. i/ L  B

# 获取聚类标签
label_pred = estimator.labels_


) a- @* S+ ^4 j2 W# 绘制数据分布图（以花萼长度和宽度为展示依据）
# plt.scatter(X[:, 0], X[:, 1], c="red", marker='o', label='see')
# print(X[:,0])
! p# Q0 C( L' ?# plt.xlabel('calyx length')  
# plt.ylabel('calyx width')  
# plt.legend(loc=2)  
4 {- c* ?8 a2 F7 F8 ]+ K# plt.show()  
2 F6 e/ \7 H) V0 {/ S1
" v- `0 q' c  q/ x1 B2
3
( Q  `- M* c! Y4
5
6
7
8
7 Q9 A" C4 c1 K( [9
10
11
. ~* O- e8 y  e12
13
" D. x( u- K) \9 W/ i3 @# _14
15
16
17
- n/ N9 X7 {8 U/ K  P9 |6 {7 I) _18
6 h* m2 |# M% A. ]" [19
20
  r: Z$ O, a  S# B& |. c21
22
8 [) u) b4 i8 N; b23
, E/ ^! k( @: ^, G9 I* g2 j0 o- j24
25
( j# D6 \& G  P7 @" s, x26
: Q8 O. T; ?; q% r3 I27
% v; J# x+ Y  c) P以下为聚类结果图，针对两个维度进行了显示。

( p# ]5 T' d; h- y' [" N
; W. L% j8 E) Q) ^) \
% Q8 k+ z7 Z) j! W9 E: S7 h
! B" p9 g  ]& y5 w9 f/ L多元回归分析
, l  V7 L; |4 \, a# V. _应用：主要用于描述关系或进行预测。

+ @+ H/ g# W% v" n8 E8 i2 D' b
7 z, k/ O$ J8 s5 m- A+ i8 Y原理：建立因变量和多个自变量之间的线性或者非线性关系。常见的回归分析是线性回归，在线性回归中因变量y也可以不作为结果，而作为另一个需要得到的结果g(x)的因变量，即：
* ^0 k8 W7 i9 V+ k( t5 ]. N# \，这样的的线性回归称为广义线性回归。
基本步骤：主要分为线性回归和非线性回归。就不总结了大伙自己看看吧。
1 x& F* j3 s  e8 A3 Y

文档资料：挺多的不找了，大部分都是线性回归，非线性回归资料少。


- P# J% b1 `0 ~# A6 W* Y; x下面两个图分别是线性回归和非线性回归的图。
3 h3 o7 t3 X' m: p
% R* l% g# P0 o( X' D

. Q; T/ }; h! g- |) r7 ]) s" x
7 i3 y3 d( c6 V; E6 i8 y! @————————————————
版权声明：本文为CSDN博主「Rulcy」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
$ m3 j+ A1 G" f2 m. Z3 [/ b/ q原文链接：https://blog.csdn.net/weixin_44371912/article/details/106934687
8 L0 J$ W& q2 W

, Y  M2 e+ L" ]0 M. g4 `