一份简短又全面的数学建模技能图谱：常用模型&算法总结（一）——主成分分析

浅夏110

本文总结了常用的数学模型方法和它们的主要用途，主要包括数学和统计上的建模方法，关于在数学建模中也挺常用的机器学习算法暂时不作补充，以后有时间就补。至于究竟哪个模型更好，需要用数据来验证，还有求解方法也不唯一，比如指派问题，你可以用线性规划OR动态规划OR整数规划OR图与网络方法来解。
6 m' G: @% M7 L! [% V  K0 G
3 D7 H5 d. i. y6 ]" w( Q总的来说，常用主成分分析来降维，综合评价的方法有因子分析法、层次分析法、....，用蒙特卡罗方法来随机模拟求解；还应掌握数据变换、数据拟合、参数估计、插值等数据处理，线性规划、整数规划、目标规划、动态规划类问题的求解要根据已知信息找出约束条件与目标函数，图论算法也是非常常用的，组合优化算法常用于很难求出最优解的NP问题，还有一些连续离散化的技术 eg.通过插值or拟合or光滑技术【移动平均之类的】可以把离散数据连续化，通过分组【把数据划分成不同的小区间】OR差分就可以把连续数据离散化....。此外的数值分析算法eg方程组求解、矩阵运算、函数积分等算法也经常用得到。至于图像处理方法稍后再补。后面会再补一份用python进行数据分析的资料。
6 C! E. c' x, M
& Y/ f! ?: k" e$ b【1】主成分分析
    目的是希望用较少的变量去解释原来资料中的大部分变异，把相关性很高的变量转化成彼此相互独立或不相关的变量，是一种降维方法。 在描述数据集中的样本时，样本又叫作实例、观测，样本可以由多个属性来描述，这些又可以称为特征、指标、变量、维度，比如描述某企业的员工信息时，数据集就是所有员工信息，每个员工就是一个样本，用来描述样本信息的性别、年龄、工龄、籍贯、工资....就是特征，这些指标可能有某种程度上的相关关系，就会存在信息冗余，就需要特征选择，也就是降维，常用的降维方法有主成分分析pca, SVD奇异值分解，逐步回归。。。
+ F/ Y) B: ~, \: F, O
4 Z, R9 v0 w$ J主成分分析 (一）： 基本思想与主成分估计方法
( Y3 y! ?& e, C5 A# a& h主成分分析（principal component analysis）是1901年Pearson对非随机变量引 入的，1933年Hotelling将此方法推广到随机向量的情形，主成分分析和聚类分析有很 大的不同，它有严格的数学理论作基础。  主成分分析的主要目的是希望用较少的变量去解释原来资料中的大部分变异，将我 们手中许多相关性很高的变量转化成彼此相互独立或不相关的变量。通常是选出比原始 变量个数少，能解释大部分资料中的变异的几个新变量，即所谓主成分，并用以解释资 料的综合性指标。由此可见，主成分分析实际上是一种降维方法。

4 c3 P/ Z+ D3 d3 ?
, u* t6 f" |4 |

1  基本思想及方法

! Y* U2 l! T5 s+ h

1 确定各个主成分的方法
一个主成分不足以代表原来的 p 个变量，因此需要寻找第二个乃至第三、第四主 成分，第二个主成分不应该再包含第一个主成分的信息，统计上的描述就是让这两个主 成分的协方差为零，几何上就是这两个主成分的方向正交。具体确定各个主成分的方法 如下。
8 M9 Z1 I' Y; l5 i; g2 o


2 注意事项
/ W. P. p. E. \& b 1）主成分分析的结果受量纲的影响，由于各变量的单位可能不一样，如果各自改变量纲，结果会不一样，这是主成分分析的大问题，回归分析是不存在这种情况的， 所以实际中可以先把各变量的数据标准化，然后使用协方差矩阵或相关系数矩阵进行分 析。

4 l1 v) I1 X0 h 2）为使方差达到大的主成分分析，所以不用转轴（由于统计软件常把主成分分 析和因子分析放在一起，后者往往需要转轴，使用时应注意）。

6 l' |/ h1 x; N2 U. c6 s1 O 3）主成分的保留。用相关系数矩阵求主成分时，Kaiser主张将特征值小于1的主成 分予以放弃（这也是SPSS软件的默认值）。

$ Z8 j. d2 k: t 4）在实际研究中，由于主成分的目的是为了降维，减少变量的个数，故一般选取 少量的主成分（不超过5或6个），只要它们能解释变异的70％～80％（称累积贡献率） 就行了。  下面我们直接通过主成分估计（principle estimate）进一步阐述主成分分析的基 本思想和相关概念。

2  主成分估计
主成分估计（principal component estimate）是Massy在1965年提出的，它是回归系数参数的一种线性有偏估计（biased estimate），同其它有偏估计，如岭估计（ridge estimate）等一样，是为了克服小二乘（LS）估计在设计矩阵病态（即存在多重共线性） 时表现出的不稳定性而提出的。  

) o: w/ j0 D0 k4 }; r主成分估计采用的方法是将原来的回归自变量变换到另另一组变量，即主成分，选择其中一部分重要的主成分作为新的自变量（此时丢弃了一部分，影响不大的自变量， 这实际达到了降维的目的），然后用小二乘法对选取主成分后的模型参数进行估计， 后再变换回原来的模型求出参数的估计。
: E* @# _% X% r" W0 N8 h
, {, v1 f6 c+ |  s4 N& q: h
0 d. u# ?% n+ `- y+ n3 q0 G# P" y
. a2 p% o$ Y# m5 I
# G  O, h' r8 W6 J% |单参数主成分估计
4 g& V! o, h/ N2 `& o# E" W  t* y理论上表明：主成分估计在设计阵病态时优于LS估计，但（31）在特征值为1的附 近存在跳跃，会影响计算的稳定性，杨虎在1989年给出的单参数主成分估计解决了这个 问题。

4 E( b/ b: K2 H% q
例3  Hald水泥问题，考察含如下四种化学成分

* [7 [) N( s- X3 D9 Q/ U+ K2 P
4 ?9 g  ]! u& ]2 F. F
" b* P" `: E. c- @, Y9 B* ]
相关系数阵的四个特征值依次为2.2357，1.5761，0.1866，0.0016。最后一个特征值接近于零，前三个特征值之和所占比例（累积贡献率）达到0.999594。于是我们略去 第4个主成分。其它三个保留的特征值对应的三个特征向量分别为

6 F4 F& m! |2 A
0 d. @5 I" N4 m5 Z( U
  V+ W0 Y. H/ d4 D+ e3 y- p0 n对Hald数据直接作线性回归得经验回归方程
% B3 M8 k& H3 M1 U9 K7 U7 Q

, B: o1 j  X' n
' k' g! e+ M, {$ T% u# K再由（31）式计算出主成分估计，即可获得如下主成分回归方程
  G' G4 D, }# X4 u  h  V8 v
8 }4 c$ ^% Q3 a/ W; y
  o4 ^2 I0 w3 {9 g* _1 }7 y
两个方程的区别在于后者具有更小的均方误差，因而更稳定。此外前者所有系数都无法 通过显著性检验。  计算的MATLAB程序如下：
& J# L4 l5 c& I$ Z1 h% N, v6 l9 |7 lclc,clear
! u+ ^2 \1 Q7 j3 w4 Uload sn.txt  %把原始的x1,x2,x3,x4,y的数据保存在纯文本文件sn.txt中
[m,n]=size(sn);num=3; %num为选取的主成分的个数
mu=mean(sn);sigma=std(sn);
) M; J2 W  Y% ?/ rsnb=zscore(sn); %数据标准化
b=snb(:,1:end-1); %x1,x2,x3,x4的数据赋给b
r=cov(b);  %标准化数据的协方差阵就是相关系数阵
[x,y,z]=pcacov(r);
f=repmat(sign(sum(x)),size(x,1),1);
x=x.*f;
%以下是普通的小二乘法回归
' d- W9 P- l6 J) q& jr=[ones(m,1),b]\snb(:,end);  %标准化数据的回归方程系数
6 b0 g2 E& l+ f6 W- l3 ^; B, x+ C' {* Hbzh=mu./sigma;
8 }# X8 M+ G! f/ Wch10=mu(end)-bzh(1:end-1)*r(2:end)*sigma(end)  %原始数据的常数项
4 t! K& ~8 w0 N( P, D$ afr=r(2:end);fr=fr';
ch1=fr./sigma(1:end-1)*sigma(end) %原始数据的x1,x2等等系数
%以下是主成分回归
pval=b*x(:,1:num);
rp=[ones(m,1),pval]\snb(:,end);  %主成分数据的回归方程系数
0 a( z4 f1 P% Q: ^7 A1 Y- bbeta=x(:,1:num)*rp(2:num+1);           %标准化数据的回归方程系数
ch20=mu(end)-bzh(1:end-1)*beta*sigma(end)  %原始数据的常数项
& J) N9 k; ^$ C3 R% h5 cfr=beta';
ch2=fr./sigma(1:end-1)*sigma(end) %原始数据的x1,x2等等系数
check1=sqrt(sum((sn(:,1:end-1)*ch1'+ch10-sn(:,end)).^2)/(m-n))
check2=sqrt(sum((sn(:,1:end-1)*ch2'+ch20-sn(:,end)).^2)/(m-num-1))
8 P) G- J4 B- k6 U. g; J9 G


4 @1 z/ Z/ T$ |: {. F  R
0 k- Z4 N& r$ O$ v8 u7 I
$ ^" |: K; I  X& u0 |0 _主成分分析法（二）：计算步骤
主成分分析法的步骤
1）对原始数据进行标准化处理


) c/ n* i8 e2 d4 o+ c/ O           其中 , 即分别为第j个指标的样本均值和标准差。



0 e8 h( G% \0 V6 T. B2）计算相关系数矩阵R
, a  ]1 W) a: k$ S. |2 [. ?

3）计算特征值和特征向量
   特征值要由大到小排序
" d; `4 `  ~" n5 g: r. t& {


& c9 j+ B" R/ s" e% ~, _4）选择 p （ p ≤ m ）个主成分，计算综合评价值




h

& U  j4 x6 U3 Q0 \matlab实现主成分分析的pcacov( )函数
& A/ L0 I: q2 I3 M* P9 \clc,clear
load gj.txt   %把原始数据保存在纯文本文件gj.txt中
gj=zscore(gj); %数据标准化
' s) F6 i7 \5 C+ `, A5 T; ir=corrcoef(gj);  %计算相关系数矩阵
/ r7 N& \. u# d& Z[x,y,z]=pcacov(r);
$ j$ D: i! ~  d7 u8 P) C0 |" gf=repmat(sign(sum(x)),size(x,1),1);
x=x.*f;
df=gj*x(:,1:4)
% V0 B  z6 Q7 W5 R; Wtf=df*z(1:4)/100;
[stf,ind]=sort(tf,'descend')

4 g/ C/ [' K3 z% O+ s

————————————————
* M5 v2 Q) M& k. e5 ^- U+ [& l; q9 y版权声明：本文为CSDN博主「wamg潇潇」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
6 E1 o$ {3 O9 r% J* u原文链接：https://blog.csdn.net/qq_29831163/article/details/88900766
, Z  l8 x5 G+ j% P. \; Z3 [: G

; h4 s7 c& Q( K
4 `9 ]& e- Y7 S+ D" g1 a; r