让我们生成一个数组:
import numpy as np
data = np.arange(30).reshape(10,3)
data=data*data
array([[ 0, 1, 4],
[ 9, 16, 25],
[ 36, 49, 64],
[ 81, 100, 121],
[144, 169, 196],
[225, 256, 289],
[324, 361, 400],
[441, 484, 529],
[576, 625, 676],
[729, 784, 841]])
Run Code Online (Sandbox Code Playgroud)
然后找到协方差矩阵的特征值:
mn = np.mean(data, axis=0)
data -= mn
C = np.cov(data.T)
evals, evecs = la.eig(C)
idx = np.argsort(evals)[::-1]
evecs = evecs[:,idx]
print evecs
array([[-0.53926461, -0.73656433, 0.40824829],
[-0.5765472 , -0.03044111, -0.81649658],
[-0.61382979, 0.67568211, 0.40824829]])
Run Code Online (Sandbox Code Playgroud)
现在让我们对数据运行matplotlib.mlab.PCA函数:
import matplotlib.mlab as …Run Code Online (Sandbox Code Playgroud) 我有:
>>> import numpy as np
>>> a = np.arange(25).reshape(5, 5)
>>> a
array([[ 0, 1, 2, 3, 4],
[ 5, 6, 7, 8, 9],
[10, 11, 12, 13, 14],
[15, 16, 17, 18, 19],
[20, 21, 22, 23, 24]])
Run Code Online (Sandbox Code Playgroud)
如何获得一个数字和数字之和的数组:
array([[7, 9, 11,], # 1+6=7 , 2+7=9 , 3+8=11
[17, 19, 21,],
[27, 29, 31,]])
Run Code Online (Sandbox Code Playgroud)
我想通过迭代原始数组来做到这一点.就像是:
b[x,y] = [a[x,y]+a[x+1,y] for x in range(0,3) for y in range(1,4)] #ERROR!
Run Code Online (Sandbox Code Playgroud)
虽然这不符合语法原因.有人可以给我正确的语法吗?我不是一名专业程序员,也不熟悉Python.先感谢您.
比方说我有:
numpy.random.seed(20)
a=numpy.random.rand(5000)
b=numpy.random.rand(5000)
Run Code Online (Sandbox Code Playgroud)
我想获得一个位置的索引a[x] > b[x],即所有x的索引
此外,我想得到一个指数(a[x-1] < b[x-1]) && (a[x] > b[x]).
有人可以帮忙吗?我有一种感觉,我必须使用蒙面数组,但我不知道如何.