#!?/usr/local/bin/python3.6
#?coding=utf-8

#?作者：huwenhao
#?Github主頁：?https://github.com/huwenhao1127/
import?numpy?as?np
import?os

np.set_printoptions（threshold=np.inf）
os.environ[‘TF_CPP_MIN_LOG_LEVEL‘]?=?‘2‘


#?PCA（principal?component?analysis）
#?data:?raw?data
#?K:?new?dimensions
def?pca（data?k）:
????X?=?np.asmatrix（data）
????RawDim?NumData?=?data.shape[0]?data.shape[1]

????#?Make?X?a?mean-normalized?m*n?data?matrix
????for?i?in?range（RawDim）:
????????X[i]?=?X[i]?-?X[i].mean（）

????#?Compute?symmetric?matrix?of?X
????C?=?X*X