fmin_cg: erhält einen Valueerror wegen Typenkonflikt

Question

den folgenden Code Ausführen gibt einen Valueerror:

def f(theta): 
    theta = theta.reshape(2, 1) 
    return np.linalg.norm(x*theta -y)**2 

def fprime(theta): 
    theta = theta.reshape(2,1) 
    return (x.T*(x*theta - y)) 

x = np.matrix('1, 2; 3, 4; 5, 6') 
y = np.matrix('4; 2; 1') 
thetainit = np.matrix('0; 0') 

scipy.optimize.fmin_cg(f, np.ravel(thetainit), fprime=fprime) 

Aber die Dimensionen sind in Ordnung und das x0 Argument für fmin_cg abgeflacht ist mit np.ravel wie in der Dokumentation beschrieben. Hier die Fehlermeldung:

1169  gnorm = vecnorm(gfk, ord=norm) 
    1170  while (gnorm > gtol) and (k < maxiter): 
-> 1171   deltak = numpy.dot(gfk, gfk) 
    1172 
    1173   try: 

Ich schätze jede Hilfe.

Answer 1

Dieser Code funktioniert:

def f(theta): 
    theta = theta.reshape(2, 1) 
    return np.linalg.norm(x*theta -y)**2 
# return np.linalg.norm(x*theta -y)**2 

def fprime(theta): 
    theta = theta.reshape(2,1) 
    g = np.array(x.T*(x*theta - y)).flatten() 
    return g 

x = np.matrix('1, 2; 3, 4; 5, 6') 
y = np.matrix('4; 2; 1') 
thetainit = np.matrix('0; 0') 

scipy.optimize.fmin_cg(f, np.ravel(thetainit), fprime=fprime) 

Der entscheidende Punkt hinzuzufügen war 'abflachen()'!