数学建模社区-数学中国
标题:
主成分分析法（二）：计算步骤
[打印本页]
作者:
浅夏110
时间:
2020-6-8 15:26
标题:
主成分分析法（二）：计算步骤
一：主成分分析简述
2 F8 ^- ^8 D+ C+ u
主成分分析是多元统计分析中重要的降维与分析评价方法；多元分析（multivariate analyses）是多变量的统计分析方法，是数理统计中应用广泛的一个重要分支。由于变量的相关性，不能简单地把每个变量的结果进行汇总，这是多变量统计分析的基本出发点。
2 w5 f) K1 w4 n' J* u: R$ T6 |
* M* x @4 T+ }5 w8 L1 G# U
主成分分析（principal component analysis）是1901年Pearson对非随机变量引 入的，1933年Hotelling将此方法推广到随机向量的情形，主成分分析和聚类分析有很 大的不同，它有严格的数学理论作基础。
. A3 |: u/ J+ V) L3 E3 S* Q
& c4 w. X7 k$ s6 q: ?5 F
主成分分析的主要目的是希望用较少的变量去解释原来资料中的大部分变异，将我 们手中许多相关性很高的变量转化成彼此相互独立或不相关的变量。通常是选出比原始 变量个数少，能解释大部分资料中的变异的几个新变量，即所谓主成分，并用以解释资 料的综合性指标。由此可见，主成分分析实际上是一种降维方法。 主成分分析试图在力保数据信息丢失少的原则下，对多变量的截面数据表进行最佳综合简化，即对高维变量空间进行降维处理。
! u9 r3 D, S n3 F' ~/ O
+ k6 q# w4 Z4 D- L
二：主成分分析法的步骤
) c# k1 ?, K: r8 `4 ]
1）对原始数据进行标准化处理
3 A! |, a0 ~2 k `6 q5 c1 ?! ?" C
% [2 L; O6 f9 Y4 J; O
: C$ P5 G5 j7 G% g
其中
, 即分别为第j个指标的样本均值和标准差。
6 a6 N" Y; {" Q
& r$ Q7 Q. t; m& d* H t
3 E! s/ R6 X9 X3 ]0 Z
4 I4 M! K* k' R4 ~
2）计算相关系数矩阵R
; {6 N" l& X) ]8 ~: E# X9 w( Y
! l7 D& z9 G& i' Z: L) r
! C" v% W e& v
& [- c/ x: m7 ^9 q; ]
3）计算特征值和特征向量
/ J) x. M$ B) Q9 i
特征值要由大到小排序
7 H9 y6 l1 P" P8 ^' q. y
4 E$ p) t* a9 D: }/ N3 _% g' ?- n
Q6 L: n' o2 _2 N* v/ u
, ]9 [/ D7 b* z, ?: N( F+ T8 `3 G
4）选择 p （ p ≤ m ）个主成分，计算综合评价值
1 V, Z& M! S7 Z; |9 B0 x
. T+ z d2 v/ x3 p
, m3 |$ c) ~5 ~# s
8 D3 h7 H) q) {0 R1 b, T$ }8 o& D
h
1 U) F5 u% H9 f3 `/ ^
' f C1 X8 b7 t4 R* ^
matlab实现主成分分析的pcacov( )函数
! d2 k0 |- G0 v8 A
8 Z6 V1 M; ] {# Z, E! ^: i, Q
clc,clear
4 W; n5 U6 D" G: q
load gj.txt %把原始数据保存在纯文本文件gj.txt中
/ y% V/ s% U. Y ?
gj=zscore(gj); %数据标准化
) }8 x2 I8 z& R* g* ?; Z
r=corrcoef(gj); %计算相关系数矩阵
2 r) Y1 e1 q: |4 c
[x,y,z]=pcacov(r);
: t C5 i) d* C& H( }
f=repmat(sign(sum(x)),size(x,1),1);
) u2 k% O: U& k- C7 O0 O
x=x.*f;
/ n2 l/ T$ t+ c$ s
df=gj*x(:,1:4)
# h0 W& O9 r* ]' d& B6 A
tf=df*z(1:4)/100;
1 `- g6 W- d! m2 N5 A0 j# W t9 u
[stf,ind]=sort(tf,'descend')
G% e& L: M5 d6 ]+ h& g
7 C6 O6 a: N u0 E5 u
' i3 @- i; M8 x+ o6 j/ H; y2 m: I
————————————————
% s! A3 O6 W7 M" k2 b4 m
版权声明：本文为CSDN博主「wamg潇潇」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
5 R: U+ [6 ^& f% K/ T$ K0 X( R
原文链接：https://blog.csdn.net/qq_29831163/article/details/88900766
. z V5 `9 n9 [4 y
) C! q' p' L" C* p# Q/ u8 o
7 A2 ^/ Q* ^ u* }/ ^+ a) {
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
: C$ P5 G5 j7 G% g
! C" v% W e& v
Q6 L: n' o2 _2 N* v/ u
, m3 |$ c) ~5 ~# s
h