Question

我正在使用以下代码来识别面孔：

def clasificaImagen(numVS, numEntre, I):

    if numVS > 70: # No puede haber mas de 400 valores singulares
        numVS = 70;

    if numEntre > 76: # Queremos al menos 100 imagenes de prueba
        numEntre = 76;

    # Importamos las imagenes
    personajes = np.array([144,109,200,121,200,77]);
    persuma=np.array([0,144,253,453,574,774,851]);

    # Declaramos nuestras variables
    E = []; #np.mat(4096, numEntre, 6); %zeros(4096, numEntre, 6);
    U = []; #np.mat(4096, numEntre, 6); %zeros(6,1);
    Uk = []; #cell(4096, numVS, 6); %zeros(6, 1);
    aciertos = 0; 
    pruebas = sum(personajes)-(numEntre*6);  
    pguess = 0; # El digito que el programa asigna a cada imagen

    # Calculamos las matrices de entrenamiento para cada digito
    print 'Generando conjuntos de entrenamiento';

    for j in range(0,6):
        # E(j) tiene dimension 76 x numEntre
        E.append(I[:, persuma[j]+1:persuma[j]+numEntre]);
        # U(j) tiene dimension 76 x 76
        u, S, V = np.linalg.svd(E[j]);
        U.append(u);
        # Uk(j) tiene dimension 76 x numVS
        uk = U[j];
        uk = uk[:, 1:numVS];
        Uk.append(uk);

        type(Uk[1])
        print 'Clasificando imagenes';

    for i in range(0,6): 

        rMin = 9999999; 
        P = I[:, persuma[i]+1+numEntre:persuma[i+1]]; 

        for h in range(0,(personajes[i] - numEntre)-1):
            z = P[:, h]; 

            for j in range(0,6):
                ukMult = Uk[j]*Uk[j].T;
                ukMult = ukMult*z;
                ukMult = z - ukMult;
                res = np.linalg.norm(ukMult);
                if res < rMin:
                    pguess = j - 1;
                    rMin = res;

            if pguess == i-1:
                aciertos = aciertos + 1;

    print 'Calculando porcentaje de aciertos';
    return 100*aciertos/pruebas;

问题在于这一行：

ukMult = Uk[j]*Uk[j].T;

发生了什么事？关于代码的评论是西班牙语，但我只是想看看有人能解释矩阵乘法的问题是什么。如有必要，我可以用英文发表评论。

我认为乘法是明确定义的，因为它只是矩阵及其转置的乘积，所以不应该有任何问题。

使用

  np.dot(Uk[j]*Uk[j].T)

也没有用。

Answer 1

调用numpy.dot（）函数时，它需要两个参数，而不是上面显示的参数。尝试这样的事情：

np.dot(Uk[j], Uk[j].T)

文档在这里：http://docs.scipy.org/doc/numpy/reference/generated/numpy.dot.html

他们列出的示例如下所示：

>>> a = [[1, 0], [0, 1]]
>>> b = [[4, 1], [2, 2]]
>>> np.dot(a, b)
array([[4, 1],
       [2, 2]])

另请查看numpy.matrix类：http://docs.scipy.org/doc/numpy-1.10.1/reference/generated/numpy.matrix.html

将Python中循环内的两个矩阵乘以

1 个答案: