Как эффективно вычислить матрицу ядра Гаусса в numpy? ⇐ Python

[Anonymous]

Код: Выделить всё

def GaussianMatrix(X,sigma):
row,col=X.shape
GassMatrix=np.zeros(shape=(row,row))
X=np.asarray(X)
i=0
for v_i in X:
j=0
for v_j in X:
GassMatrix[i,j]=Gaussian(v_i.T,v_j.T,sigma)
j+=1
i+=1
return GassMatrix
def Gaussian(x,z,sigma):
return np.exp((-(np.linalg.norm(x-z)**2))/(2*sigma**2))

Это мой текущий путь. Есть ли способ использовать матричную операцию для этого? X — точки данных.

Подробнее здесь: https://stackoverflow.com/questions/297 ... y-in-numpy