Sınırlama ile optimizasyon

0

Soru

Kısıtlı bir optimizasyon problemini çözmek istiyorum.

maks {ln (c1) + ln (c2)}

s. t. 4 (c1) + 6 (c2) ≥ 40

Bu kod yazdım :

import numpy as np
from scipy import optimize

def main():
    """
    solving a regular constrained optimization problem    
    max ln(cons[0]) + ln(cons[1]) 
    st. prices[0]*cons[0] + prices[1]*cons[1] <= I         
    """    
   
    prices = np.array([4.0, 6.0])
    I = 40.0
   
    util = lambda cons: np.dot( np.log(cons))  #define utility function
    budget = lambda cons: I - np.dot(prices, cons)   #define the budget constraint
    
    initval = 40.0*np.ones(2)    #set the initial guess for the algorithm
    
    res = optimize.minimize(lambda x: -util(x), initval, method='slsqp', 
                            constraints={'type':'ineq', 'fun':budget}, 
                            tol=1e-9)
    assert res['success'] == True
           
    
    print(res)

Ne yazık ki, kodum herhangi bir çözüm yazdırmıyor. Bunu anlamam için bana yardım EDEBİLİR MİSİN?

mathematical-optimization python scipy
2021-11-17 08:52:35
1

En iyi cevabı

1

Kodunuz o zamandan beri bir TypeError verir np.dot iki argüman bekler, tanımınıza bakın utils işlev. Bu nedenle, kullanın

# is the same as np.dot(np.ones(2), np.log(cons))
utils = lambda cons: np.sum(np.log(cons))

yerine.

2021-11-17 21:09:35

Diğer dillerde

Bu sayfa diğer dillerde

Русский
..................................................................................................................
Italiano
..................................................................................................................
Polski
..................................................................................................................
Română
..................................................................................................................
한국어
..................................................................................................................
हिन्दी
..................................................................................................................
Français
..................................................................................................................
Česk
..................................................................................................................
Português
..................................................................................................................
ไทย
..................................................................................................................
中文
..................................................................................................................
Español
..................................................................................................................
Slovenský
..................................................................................................................