蟒蛇優化多種功能與常見的變量

Question

我試圖儘量減少（全球）使用公共變量3個功能，我想將它們組合成一個功能，並最大限度地減少使用L-BFGS-B（我需要設置變量邊界），但它已經證明很難平衡每個參數與權重，即當一個被最小化時，另一個不會。我還試圖用SLSQP方法，儘量減少其中的一個，同時設置其他作爲約束條件，但制約因素往往被忽略/不符合。 這裏有什麼需要儘量減少，所有的數學都在meritscalculation和meritoflength,meritofROC,,heightorderreturned作爲全局值從計算返回。

def lengthmerit(x0): 
    meritscalculation(x0) 
    print meritoflength 
    return meritoflength 

def ROCmerit(x0): 
    meritscalculation(x0) 
    print meritofROC 
    return meritofROC 

def proximitymerit(x0): 
    meritscalculation(x0) 
    print meritofproximity+heightorder 
    return meritofproximity+heightorder 

我希望所有的人都使用通用X0（含邊界）爲自變量，以儘量減少，是有辦法做到這一點？

Answer 1

這是你想要做什麼？

minimize a * amerit(x) + b * bmerit(x) + c * cmerit(x) 
over a, b, c, x: 
    a + b + c = 1 
    a >= 0.1, b >= 0.1, c >= 0.1 (say) 
    x in xbounds 

如果x是說[x0 x1 .. x9]，建立一個新的變量abcx = [a b c x0 x1 .. x9]， 約束a + b + c = 1與懲罰項加入目標函數， 並儘量減少這樣的：

define fabc(abcx): 
    """ abcx = a, b, c, x 
     -> a * amerit(x) + ... + penalty 100 (a + b + c - 1)^2 
    """ 
    a, b, c, x = abcx[0], abcx[1], abcx[2], abcx[3:] # split 
    fa = a * amerit(x) 
    fb = b * bmerit(x) 
    fc = c * cmerit(x) 
    penalty = 100 * (a + b + c - 1) ** 2 # 100 ? 
    f = fa + fb + fc + penalty 
    print "fabc: %6.2g = %6.2g + %6.2g + %6.2g + %6.2g a b c: %6.2g %6.2g %6.2g" % (
       f, fa, fb, fc, penalty, a, b, c) 
    return f 

和bounds = [[0.1, 0.5]] * 3 + xbounds，即每的a b c in 0.1 .. 0.5左右。
長print S的關係告訴你爲什麼a b c方法0一個 - 可能的amerit() bmerit() cmerit()一個比別人更強壯？ Plot s，而不是print s就容易了。

總結：
1）清楚地制定在紙上的問題，因爲在頂部
2）翻譯成蟒。

Answer 2

這裏是一些縮放的結果和加權

目標函數：

merit_function=wa*meritoflength*1e3+wb*meritofROC+wc*meritofproximity+wd*heightorder*10+1000 * (wa+wb+wc+wd-1) ** 2 

輸入：

abcdex=np.array((0.5, 0.5, 0.1, 0.3, 0.1...)) 

輸出：

fun: array([ 7.79494644]) 

    x: array([ 4.00000000e-01, 2.50000000e-01, 1.00000000e-01, 
    2.50000000e-01...]) 


meritoflength : 0.00465499380753. #target 1e-5, usually start at 0.1 
meritofROC: 23.7317956542   #target ~1, range <33 
Heightorder: 0      #target :strictly 0, range <28 
meritofproximity : 0.0    #target:less than 0.02, range <0.052 

我後數實現運行時，所有的權重往往停留在邊界的最小值，並回到手動調整我開始的縮放問題。

是否有我的優化功能，心不是尋找真正的全球最小的可能性？

這是我如何最小化它：

minimizer_kwargs = {"method": "L-BFGS-B", "bounds": bnds, "tol":1e0 } 

ret = basinhopping(merit_function, abcdex, minimizer_kwargs=minimizer_kwargs, niter=10) 
zoom = ret['x'] 

res = minimize(merit_function, zoom, method = 'L-BFGS-B', bounds=bnds, tol=1e-6)