1. 程式人生 > >協方差實現(python)

協方差實現(python)

import numpy as np

if __name__ == '__main__':
    a = np.array([[1,2,3],[3,1,1]])
    c_ = np.mean(a, axis =0) #按列計算
    x = a - c_  # 廣播機制
    cov = x.T.dot(x)/(a.shape[0]-1) #協方差
    print(cov)