日本免费高清视频-国产福利视频导航-黄色在线播放国产-天天操天天操天天操天天操|www.shdianci.com

學無先后，達者為師

網站首頁 編程語言 正文

Python解決非線性規劃中經濟調度問題_python

作者：電力系統與算法之美 ? 更新時間： 2022-06-29 編程語言

1、概述

今天重點講非線性規劃中scipy.optimize.minize庫在非線性規劃中的應用。Scipy 是 Python 算法庫和數學工具包，包括最優化、線性代數、積分、插值、特殊函數、傅里葉變換、信號和圖像處理、常微分方程求解等模塊。?

scipy.optimize 模塊中提供了多個用于非線性規劃問題的方法，適用于不同類型的問題。

brent()：單變量無約束優化問題，混合使用牛頓法/二分法。

fmin()：多變量無約束優化問題，使用單純性法，只需要利用函數值，不需要函數的導數或二階導數。

leatsq()：非線性最小二乘問題，用于求解非線性最小二乘擬合問題。

minimize()：約束優化問題，使用拉格朗日乘子法將約束優化轉化為無約束優化問題。

2、scipy.optimize.minimize參數

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中為參數的部分進行約束限制

3、簡單案例引出

(1)Scipy.optimize實現

import numpy as np
from scipy.optimize import minimize
#目標函數
def fun(args1):
    a,b,c,d=args1
    r=lambda x:(a*x[0]*x[0]+b*x[1]*x[1]+c*x[2]*x[2]+d)
    return r
def con(args2):
    x0min,x1min,x2min=args2
    cons=({'type':'eq','fun':lambda x:-x[0]-x[1]**2+2},
          {'type':'eq','fun':lambda x:x[1]+2*x[2]**2-3},
          {'type':'ineq','fun':lambda x:x[0]**2-x[1]+x[2]**2},
          {'type': 'ineq', 'fun': lambda x: -(x[0]+x[1]**2+x[2]**2-20)},
          {'type':'ineq','fun':lambda x:x[0]-x0min},
          {'type':'ineq','fun':lambda x:x[1]-x1min},
          {'type':'ineq','fun':lambda x:x[2]+x2min})
    return cons
def main():
    args1=(1,2,3,8)
    args2=(0,0,0)
    cons=con(args2)
    x0=np.array((1,2,3))    #初值
    res=minimize(fun(args1),x0,method='SLSQP',constraints=cons)
    print('minf(x):',res.fun)
    print(res.success)
    print('x：',[np.around(i) for i in res.x])
    print('x1:',res.x[0])
    print('x2:',res.x[1])
    print('x3:',res.x[2])
    #另一種表述
    print("optimization problem(res):{}".format(res.x))
    print("Xopt={}".format(res.x))
    print("minf(x)={:.4f}".format(res.fun))
 
if __name__ == "__main__":
    main()

輸出結果

minf(x): 13.878994794551044
True
x： [1.0, 1.0, 1.0]
x1: 0.6743061260520056
x2: 1.1513878035150682
x3: 0.961408393062538
optimization problem(res):[0.67430613 1.1513878 ?0.96140839]
Xopt=[0.67430613 1.1513878 ?0.96140839]
minf(x)=13.8790
?
Process finished with exit code 0

(2)遺傳算法包實現?（—sko.GA&sko.DE）

#from sko.DE import DE   #差分進化法
from sko.GA import GA   #遺傳算法
def obj_func(p):
    x1, x2, x3 = p
    return x1 ** 2 + 2*x2 ** 2 + 3*x3 ** 2+8
 
 
constraint_eq = [
    lambda x: -x[0]-x[1]**2+2,
    lambda x: x[1]+2*x[2]**2-3
]
 
constraint_ueq = [
    lambda x: -x[0]**2+x[1]-x[2]**2,
    lambda x: x[0]+x[1]**2+x[2]**2-20
]
 
#de = DE(func=obj_func, n_dim=3, size_pop=50, max_iter=800, lb=[0, 0, 0],constraint_eq=constraint_eq, constraint_ueq=constraint_ueq)
ga = GA(func=obj_func, n_dim=3, size_pop=50, max_iter=800, lb=[0, 0, 0],
        constraint_eq=constraint_eq, constraint_ueq=constraint_ueq)
 
 
#best_x, best_y = de.run()
best_x, best_y = ga.run()
print('best_x:', best_x, '\n', 'best_y:', best_y)

輸出結果

best_x: [1. 1. 1.]?
?best_y: [14.]
?
Process finished with exit code 0

4、電力系統中應用——經濟調度

(1)案例

(2)Scipy.optimize實現?

import numpy as np
from scipy.optimize import minimize
#目標函數（FG1+FG2+FG3）
def fun(args1):
    a0,a1,a2,b0,b1,b2,c0,c1,c2=args1
    v=lambda x:(a0+a1*x[0]+a2*x[0]*x[0]
                +b0+b1*x[1]+b2*x[1]*x[1]
                +c0+c1*x[2]+c2*x[2]*x[2])
    return v
def con(args2):
    D,x0min,x0max,x1min,x1max,x2min,x2max=args2
    cons=({'type':'eq','fun':lambda x:D-x[0]-x[1]-x[2]},
          {'type':'ineq','fun':lambda x:x[0]-x0min},
          {'type':'ineq','fun':lambda x:-x[0]+x0max},
          {'type':'ineq','fun':lambda x:x[1]-x1min},
          {'type':'ineq','fun':lambda x:-x[1]+x1max},
          {'type':'ineq','fun':lambda x:x[2]-x2min},
          {'type':'ineq','fun':lambda x:-x[2]+x2max})
    return cons
def main():
    args1=(4,0.3,0.0007,3,0.32,0.0004,3.5,0.3,0.00045)
    args2=(700,100,200,120,250,150,300)
    cons=con(args2)
    x0=np.array((150,250,200))    #初值
    res=minimize(fun(args1),x0,method='SLSQP',constraints=cons)
    print('FGi-代價:',res.fun)
    print(res.success)
    print('PGi—解：',[np.around(i) for i in res.x])
    print('PG1:',res.x[0])
    print('PG2:',res.x[1])
    print('PG3:',res.x[2])
 
if __name__ == "__main__":
    main()

輸出結果

FGi-代價: 305.9673913046252
True
PGi—解： [176.0, 250.0, 274.0]
PG1: 176.0874477123534
PG2: 250.0
PG3: 273.9125522876465
?
Process finished with exit code 0

(3)粒子群包實現（pyswarm)

pyswarm是一個支持帶約束的粒子群優化包，sko.PSO中的pso僅支持帶上下限的約束，不支持等式和不等式約束，所以不太喜歡。

from pyswarm import pso
 
def object_func(x):
    return (4+0.3*x[0]+0.0007*x[0]*x[0]+3+0.32*x[1]+0.0004*x[1]*x[1]+3.5+0.3*x[2]+0.00045*x[2]*x[2])
 
#不等式約束
 
def cons1(x):
    return [x[0]+x[1]+x[2]-700]
 
lb = [100, 120, 150]#
ub = [200, 250, 300]
 
xopt, fopt = pso(object_func,lb,ub,ieqcons=[cons1], maxiter=100,swarmsize=1000)
print(xopt)
print(fopt)

輸出結果

Stopping search: Swarm best objective change less than 1e-08
[179.34039956 250. ? ? ? ? 270.65960045]
305.97956393103044
?
Process finished with exit code 0

原文鏈接：https://blog.csdn.net/weixin_46039719/article/details/121880540

欄目分類
最近更新