1

基本的な PSO (粒子群最適化) をコーディングしていますが、粒子インスタンスに属性がないというこのエラーが発生しています__getitem__。すべて問題ないと思いますが、particleクラスにエラーがあるようです。particleクラスを見てみましょう。

from numpy import array
from random import random
from math import sin, sqrt, cos, pi
import matplotlib.pyplot as plt
import pylab

## Settings

c1 = 2
c2 = 2
size = 100
bad_size = 10
dim = 10
max_iterations = 20
Error_limit = 0.00001

def functR(k):
    val = 10*dim
    for i in range(dim):
        val = val + (k[i])**2 - 10*cos(2*pi*k[i])
    return val
#print functR([0]*20)

class particle():
    def __init__(self, pos, fitness,vel, pbestpos, pbestfit):
        self.pos = pos
        self.fitness = fitness
        self.vel = vel
        self.pbestpos = pbestpos
        self.pbestfitness = pbestfit

class swarm():
    def __init__(self, size, bad_size, dim):
        #self.gbest = gbest
        self.size = size
        self.bad_size = bad_size
        self.dim = dim
    def create(self):
        particles = []
        for i in range(size + bad_size):
            p = particle()
            p.pos = array([random() for i in range(dim)])
            p.vel = 0.0
            p.fitness = 0.0
            p.pbestpos = p.pos
            p.pbestfit = p.fitness
            #p = particle(pos, fitness,vel, pbestpos, pbestfit)
            particles.append(p)
        return particles

def optimizer():
    s = swarm(size, bad_size, dim)
    new_swarm = s.create()
    gbest = new_swarm[0]
    gbestfit = functR(gbest)
    i = 0
## The iterative loop
    while i < max_iterations:
        for p in s:
            fitness = functR(p.pos)

            if fitness > p.fitness:
                p.fitness = fitness
                p.pbestpos = p.pos

            if fitness > gbestfit:
                gbest = p
    ## Plotting             
            pylab.xlim([0,10])
            pylab.ylim([0,1.5])
            plt.plot(i,gbest.fitness, "bo")

    ## Velocity and Position update
            vel = p.vel + c1 * random() * (p.pbestpos - p.pos) \
                    + c2 * random() * (gbest.pos - p.pos)
            p.pos = p.pos + vel
        plt.show()      
        i  += 1

    print "gbest fitness   :", gbestfit
    print "Best particle   :", gbest.pos

print optimizer()
4

1 に答える 1

1

ここでは、単一のparticle()インスタンスをリストとして扱っています。

val = val + (k[i])**2 - 10*cos(2*pi*k[i])

kのインスタンスである場合particle()、[i]構文__getitem__はそのインスタンスの呼び出しに変換されます。

ここでそのインスタンスを渡しています：

gbest = new_swarm[0]
gbestfit = functR(gbest)

他の場所では、.pos代わりにパラメーターを渡します。

for p in s:
    fitness = functR(p.pos)

したがって、おそらく次の行でも同じことを行うつもりでしたgbestfit:

gbest = new_swarm[0]
gbestfit = functR(gbest.pos)
于 2015-01-31T09:48:05.757 に答える