一、背景：

現在項目上有一個用python 實現非線性規劃的需求。非線性規劃可以簡單分兩種，目標函數為凸函數 or 非凸函數。

凸函數的 非線性規劃，比如fun=x^2+y^2+x*y，有很多常用的python庫來完成，網上也有很多資料，比如CVXPY

非凸函數的 非線性規劃（求極值），從處理方法來說，可以嘗試以下幾種：

1.純數學方法，求導求極值；

2.使用神經網絡，深度學習來處理，可參考反向傳播算法中鏈式求導的過程；

3.尋找一些python庫來做，本文介紹scipy.optimize.minimize的使用方法

二、庫方法介紹

官方文檔：https://docs.scipy.org/doc/scipy/reference/generated/scipy.optimize.minimize.html

來看下改方法的入參

scipy.optimize.minimize(fun, x0, args=(), method=None, jac=None, hess=None, hessp=None, bounds=None, constraints=(), tol=None, callback=None, options=None)

解釋：

fun: 求最小值的目標函數

x0:變量的初始猜測值，如果有多個變量，需要給每個變量一個初始猜測值。minimize是局部最優的解法，所以

args:常數值,后面demo會講解，fun中沒有數字，都以變量的形式表示，對于常數項，需要在這里給值

method:求極值的方法，官方文檔給了很多種。一般使用默認。每種方法我理解是計算誤差，反向傳播的方式不同而已，這塊有很大理論研究空間

constraints:約束條件，針對fun中為參數的部分進行約束限制

三、demo

1.計算 1/x+x 的最小值

# coding=utf-8from scipy.optimize import minimizeimport numpy as np #demo 1#計算 1/x+x 的最小值 def fun(args):  a=args  v=lambda x:a/x[0] +x[0]  return v  if __name__ == "__main__":  args = (1) #a  x0 = np.asarray((2)) # 初始猜測值  res = minimize(fun(args), x0, method='SLSQP')  print(res.fun)  print(res.success)  print(res.x)

執行結果:函數的最小值為2點多，可以看出minimize求的局部最優

2.計算 (2+x1)/(1+x2) - 3*x1+4*x3 的最小值 x1,x2,x3的范圍都在0.1到0.9 之間

# coding=utf-8from scipy.optimize import minimizeimport numpy as np # demo 2#計算 (2+x1)/(1+x2) - 3*x1+4*x3 的最小值 x1,x2,x3的范圍都在0.1到0.9 之間def fun(args): a,b,c,d=args v=lambda x: (a+x[0])/(b+x[1]) -c*x[0]+d*x[2] return vdef con(args): # 約束條件 分為eq 和ineq #eq表示 函數結果等于0 ； ineq 表示 表達式大于等于0  x1min, x1max, x2min, x2max,x3min,x3max = args cons = ({'type': 'ineq', 'fun': lambda x: x[0] - x1min},/    {'type': 'ineq', 'fun': lambda x: -x[0] + x1max},/    {'type': 'ineq', 'fun': lambda x: x[1] - x2min},/    {'type': 'ineq', 'fun': lambda x: -x[1] + x2max},/   {'type': 'ineq', 'fun': lambda x: x[2] - x3min},/    {'type': 'ineq', 'fun': lambda x: -x[2] + x3max}) return cons if __name__ == "__main__": #定義常量值 args = (2,1,3,4) #a,b,c,d #設置參數范圍/約束條件 args1 = (0.1,0.9,0.1, 0.9,0.1,0.9) #x1min, x1max, x2min, x2max cons = con(args1) #設置初始猜測值  x0 = np.asarray((0.5,0.5,0.5))  res = minimize(fun(args), x0, method='SLSQP',constraints=cons) print(res.fun) print(res.success) print(res.x)