Python 中的最小化问题，fmin_bfgs 不起作用，但 fmin 可以，“矩阵未对齐”；

Question

我在 python 中有一个函数，它接受一个向量并返回一个实数。我使用 scipy.optimize fmin 和 fmin_bfgs 函数来查找为函数提供近似最小值的参数。然而，当我使用 fmin 时，我得到了一个好的答案（相当慢），但是当我切换到 fmin_bfgs 时，我收到一条错误消息“矩阵未对齐”。这是我的函数：

    def norm(b_):
        b_ = b_.reshape(int(M),1)                   #M already given elsewhere
        Yb = np.dot(Y,b_)                           #Y already given elsewhere

        B = np.zeros((int(M),int(M))) 
        for j in xrange(int(M)):
            B[j][j] = -t[j+1]*np.exp(-t[j+1]*Yb[j]) #The t[j] are already known
        P = np.zeros((int(M),1))
        for j in xrange(int(M)):
            P[j][0] = np.exp(-t[j+1]*Yb[j])

        diff = np.zeros((int(M),1))                 #Functions d(i,b) are known
        for i in xrange(1,int(M)-1):
            diff[i][0] = d(i+1,b_) - d(i,b_)
        diff[0][0] = d(1,b_)
        diff[int(M)-1][0] = -d(int(M)-1,b_)

        term1_ = (1.0/N)*(np.dot((V - np.dot(c,P)).transpose(),W))
        term2_ = np.dot(W,V - np.dot(c,P))          #V,c,P,W already known
        term1_ = np.dot(term1_,term2_)
        term2_ = lambd*np.dot(Yb.transpose(),diff)

        return term1_ + term2_

这是我调用 fmin_bfgs 的方式：

    fmin_bfgs(norm, b_guess,fprime=None,
              args=(),gtol=0.0001,norm=0.00000000001,
              epsilon=1.4901161193847656e-08,maxiter=None,
              full_output=0, disp=1, retall=0, callback=None)

当我调用 fmin 时，它工作正常，只是太慢而没有用（我需要优化几次）。但是当我尝试 fmin_bfgs 时，我收到此错误：

回溯（最近一次调用最后）： 文件“C:\Program Files\Wing IDE 101 4.0\src\debug\tserver_sandbox.py”，第 287 行，模块中 文件“C：\ Python27 \ Lib \ site-packages \ scipy \ optimize \ optimize.py”，第491行，在fmin_bfgs old_fval，old_old_fval中） 文件“C:\Python27\Lib\site-packages\scipy\optimize\linesearch.py​​”，第 239 行，line_search_wolfe2 derphi0、c1、c2、amax） 文件“C:\Python27\Lib\site-packages\scipy\optimize\linesearch.py​​”，第 339 行，在 scalar_search_wolfe2 phi0、derphi0、c1、c2 中） 文件“C:\Python27\Lib\site-packages\scipy\optimize\linesearch.py​​”，第 471 行，位于 _zoom derphi_aj = derphi(a_j) 文件“C:\Python27\Lib\site-packages\scipy\optimize\linesearch.py​​”，第 233 行，在 derphi 中返回 np.dot(gval[0], pk) ValueError：矩阵未对齐

您知道为什么会发生这种情况吗？我为该函数提供的所有矩阵均已正确对齐（并且该函数自 fmin 工作以来一直有效）。非常感谢帮助！

Answer 1

似乎其中一个程序最终处理的数字太大而无法处理。遗憾的是它不能告诉我它做得正确。不过我解决了这个问题，所以没有更多问题。抱歉，如果这浪费了您的时间。