主成分分析（pca）

import numpy as np

x = np.array([2.5,0.5,2.2,1.9,3.1,2.3,2,1,1.5,1.1])
y = np.array([2.4,0.7,2.9,2.2,3.0,2.7,1.6,1.1,1.6,0.9])
dx = np.zeros(10)
dy = np.zeros(10)
C = np.zeros([2,2])

sx = np.sum(x)/10
sy = np.sum(y)/10

def one(dn,n,sn):
    for i in range (0,10):
        dn[i] = n[i] - sn
    return dn

dx = one(dx,x,sx)
dy = one(dy,y,sy)


C[0][0] = np.sum(dx*dx)/9
C[0][1] = np.sum(dx*dy)/9
C[1][0] = np.sum(dy*dx)/9
C[1][1] = np.sum(dy*dy)/9

eig = np.linalg.eig(C)

eigvalues = eig[0]
eigvectors = eig[1]

final = np.zeros(10)
data = np.zeros([10,2])
for i in range(0,10):
    data[i][0] = dx[i]
for i in range(0,10):
    data[i][1] = dy[i]

tzxl = np.zeros(2)
tzxl[0] = eigvectors[0][1]
tzxl[1] = eigvectors[1][1]
tzxlt = tzxl.transpose()
final = np.dot(data,tzxlt)

print(final)

[-0.82797019  1.77758033 -0.99219749 -0.27421042 -1.67580142 -0.9129491
  0.09910944  1.14457216  0.43804614  1.22382056]