快捷導(dǎo)航

python人工智能算法之差分進(jìn)化算法的實(shí)現(xiàn)

更新時(shí)間：2023年08月11日 09:30:23 作者：微小冷

DE基于GA，正如進(jìn)化基于遺傳，和遺傳算法相比，差分進(jìn)化引入了差分變異模式，相當(dāng)于開(kāi)辟了一條嶄新的進(jìn)化路徑，下面就來(lái)看看差分優(yōu)化算法是如何實(shí)現(xiàn)的吧

初始化

首先，初始種群隨機(jī)產(chǎn)生解空間(L,R)范圍的向量

(0)表示第0代。

變異

初始化后，進(jìn)行變異操作，其第t tt代變異向量為

vi?(t)=xr1??(t)+rand⋅(xr2??(t)−xr3??(t))

其中r1?,r2?,r3?為隨機(jī)挑選出的與i ii互不相同的序號(hào)，相當(dāng)于隨機(jī)抽選三個(gè)值進(jìn)行基因突變。

交叉

接下來(lái)，開(kāi)始種群繁衍，即種群交叉

其中，r為隨機(jī)數(shù)，cr是預(yù)設(shè)的交叉概率因子；rj?為隨機(jī)生成的某個(gè)整數(shù)，代表第rj?個(gè)分量必為突變分量。

選擇

最后，進(jìn)行下一代的選擇，就是在目標(biāo)個(gè)體和實(shí)驗(yàn)個(gè)體中，選擇最優(yōu)秀的那個(gè)。

xi?(t+1)=argmin(ui?(t),xi?(t))

Python實(shí)現(xiàn)

首先，做一個(gè)DE算法的迭代函數(shù)

import numpy as np
from random import sample, randint, random
def evolve(xs,cr,func):
    # 變異
    vs = []
    N = len(xs[0])
    for _ in range(len(xs)):
        x = xs.pop(0)
        r1, r2, r3 = sample(xs, 3)
        xs.append(x)
        vs.append(r1 + random()*(r2-r3))
    # 交叉
    us = []
    for i in range(len(xs)):
        us.append(vs[i] if random() < cr else xs[i])
        j = randint(0, N-1)
        us[i][j] = vs[i][j]
    # 選擇
    xNext = []
    for x,u in zip(xs, us):
        xNext.append(x if func(x)<func(u) else u)
    return xNext

其中，變異過(guò)程要求隨機(jī)則取三個(gè)不同的個(gè)體，并且這個(gè)個(gè)體不能是當(dāng)前個(gè)體，所以采取的方案是先彈出當(dāng)前個(gè)體，再?gòu)睦锩骐S機(jī)挑選三個(gè)，最后再將當(dāng)前個(gè)體放回去。

交叉過(guò)程通過(guò)兩個(gè)隨機(jī)數(shù)，來(lái)決定下一代的成分。首先是通過(guò)浮點(diǎn)型隨機(jī)數(shù)來(lái)決定是用原有的解還是用新解；然后用一個(gè)整型隨機(jī)數(shù)，來(lái)判定采取變異解的元素序號(hào)。

然后寫出主循環(huán)

# N為個(gè)體數(shù)；nDim為解維度；nIter為迭代次數(shù)
def de(N, nDim, xRange, cr, func, nIter):
    xs = [uniRand(*xRange, nDim) for _ in range(N)]
    for _ in range(nIter):
        xs = evolve(xs, cr, func)
    fs = [func(x) for x in xs]
    xBest = xs[np.argmin(fs)]
    msg = f"當(dāng)前最優(yōu)結(jié)果為{np.min(fs)}，參數(shù)為"
    msg += ", ".join([f"{x:.4f}" for x in xBest])
    print(msg)

測(cè)試

最后寫個(gè)函數(shù)測(cè)試一下

def test(xs):
    _sum = 0.0
    for i in range(len(xs)):
        _sum = _sum + np.cos((xs[i]*i)/5)*(i+1)
    return _sum
???????if __name__ == "__main__":
    de(20, 5, (-5,5), 0.5, test, 100)

最終效果為