Changeset 5033 for branches/PSO/HeuristicLab.Algorithms.ParticleSwarmOptimization/3.3/ParticleUpdater.cs

Timestamp:

12/06/10 14:45:27 (13 years ago)

Author:

epitzer

Message:

Simple but complete PSO implementation (#852)

Location:

branches/PSO/HeuristicLab.Algorithms.ParticleSwarmOptimization/3.3

Files:

• . (modified) (1 prop)
• ParticleUpdater.cs (modified) (1 diff)

branches/PSO/HeuristicLab.Algorithms.ParticleSwarmOptimization/3.3

@@ -3 +3 @@
 obj
 *.user
+*.suo

@@ -3 +3 @@
 obj
 *.user
+*.suo

branches/PSO/HeuristicLab.Algorithms.ParticleSwarmOptimization/3.3/ParticleUpdater.cs

@@ -25 +25 @@
 using HeuristicLab.Operators;
 using HeuristicLab.Parameters;
+using HeuristicLab.Persistence.Default.CompositeSerializers.Storable;
+using HeuristicLab.Common;
 
 namespace HeuristicLab.Algorithms.ParticleSwarmOptimization {
-  public class ParticleUpdater : SingleSuccessorOperator { // ParticleUpdater
+  public class ParticleUpdater : SingleSuccessorOperator {
+
     #region Parameter properties
-
-    public ILookupParameter<RealVector> VelocityParameter {
-      get { return (ILookupParameter<RealVector>)Parameters["Velocity"]; }
+    public LookupParameter<IRandom> RandomParameter {
+      get { return (LookupParameter<IRandom>)Parameters["Random"]; }
     }
-
-    public IParameter RandomParameter {
-      get { return (IParameter)Parameters["Random"]; }
+    public LookupParameter<RealVector> PointParameter {
+      get { return (LookupParameter<RealVector>)Parameters["Point"]; }
     }
-
-    public ILookupParameter<RealVector> CurrentPositionParameter {
-      get { return (ILookupParameter<RealVector>)Parameters["CurrentPosition"]; }
+    public LookupParameter<RealVector> VelocityParameter {
+      get { return (LookupParameter<RealVector>)Parameters["Velocity"]; }
     }
-
-    public ILookupParameter<RealVector> BestLocalParameter {
-      get { return (ILookupParameter<RealVector>)Parameters["BestLocal"]; }
+    public LookupParameter<RealVector> PersonalBestPointParameter {
+      get { return (LookupParameter<RealVector>)Parameters["PersonalBestPoint"]; }
     }
-
-    public ILookupParameter<RealVector> BestGlobalParameter {
-      get { return (ILookupParameter<RealVector>)Parameters["BestGlobal"]; }
+    public LookupParameter<RealVector> BestPointParameter {
+      get { return (LookupParameter<RealVector>)Parameters["BestPoint"]; }
     }
-
-    public ILookupParameter<DoubleMatrix> BoundsParameter {
-      get { return (ILookupParameter<DoubleMatrix>)Parameters["Bounds"]; }
+    public LookupParameter<DoubleMatrix> BoundsParameter {
+      get { return (LookupParameter<DoubleMatrix>)Parameters["Bounds"]; }
+    }
+    public LookupParameter<DoubleMatrix> VelocityBoundsParameter {
+      get { return (LookupParameter<DoubleMatrix>)Parameters["VelocityBounds"]; }
+    }
+    public LookupParameter<DoubleValue> OmegaParameter {
+      get { return (LookupParameter<DoubleValue>)Parameters["Omega"]; }
+    }
+    public LookupParameter<DoubleValue> Phi_PParameter {
+      get { return (LookupParameter<DoubleValue>)Parameters["Phi_P"]; }
+    }
+    public LookupParameter<DoubleValue> Phi_GParameter {
+      get { return (LookupParameter<DoubleValue>)Parameters["Phi_G"]; }
     }
     #endregion
 
+    #region Parameter Values
+    public IRandom Random {
+      get { return RandomParameter.ActualValue; }
+    }
+    public RealVector Point {
+      get { return PointParameter.ActualValue;  }
+      set { PointParameter.ActualValue = value;  }
+    }
+    public RealVector Velocity {
+      get { return VelocityParameter.ActualValue; }
+      set { VelocityParameter.ActualValue = value; }
+    }
+    public RealVector PersonalBestPoint {
+      get { return PersonalBestPointParameter.ActualValue; }
+    }
+    public RealVector BestPoint {
+      get { return BestPointParameter.ActualValue; }
+    }
+    public DoubleMatrix Bounds {
+      get { return BoundsParameter.ActualValue;  }
+    }
+    public DoubleMatrix VelocityBounds {
+      get { return VelocityBoundsParameter.ActualValue; }
+    }
+    public double Omega {
+      get { return OmegaParameter.ActualValue.Value; }
+    }
+    public double Phi_P {
+      get { return Phi_PParameter.ActualValue.Value;  }
+    }
+    public double Phi_G {
+      get { return Phi_GParameter.ActualValue.Value; }
+    }
+    #endregion
+
+    #region Construction & Cloning
+
+    [StorableConstructor]
+    protected ParticleUpdater(bool deserializing) : base(deserializing) { }
+    protected ParticleUpdater(ParticleUpdater original, Cloner cloner)
+      : base(original, cloner) {
+    }
+
     public ParticleUpdater()
       : base() {
-      Parameters.Add(new LookupParameter<IRandom>("Random", "Random number generator (to generate alpha an beta)."));
-      Parameters.Add(new LookupParameter<RealVector>("Velocity", "The velocity vector to update."));
-      Parameters.Add(new LookupParameter<RealVector>("CurrentPosition", "Current position"));
-      Parameters.Add(new LookupParameter<RealVector>("BestLocal", "Best local position"));
-      Parameters.Add(new LookupParameter<RealVector>("BestGlobal", "Best global position"));
-      Parameters.Add(new LookupParameter<DoubleMatrix>("Bounds", "The lower and upper bounds for each dimension of the position vector."));
+      Parameters.Add(new LookupParameter<IRandom>("Random", "Random number generator."));
+      Parameters.Add(new LookupParameter<RealVector>("Point", "Particle's current position"));
+      Parameters.Add(new LookupParameter<RealVector>("Velocity", "Particle's current velocity."));
+      Parameters.Add(new LookupParameter<RealVector>("PersonalBestPoint", "Particle's personal best position"));
+      Parameters.Add(new LookupParameter<RealVector>("BestPoint", "Global best position"));
+      Parameters.Add(new LookupParameter<DoubleMatrix>("Bounds", "The lower and upper bounds for each dimension of the position vector for the current problem."));
+      Parameters.Add(new LookupParameter<DoubleMatrix>("VelocityBounds", "Upper and lower bounds for the particle's velocity vector."));
+      Parameters.Add(new LookupParameter<DoubleValue>("Omega", "The weight for the particle's velocity vector."));
+      Parameters.Add(new LookupParameter<DoubleValue>("Phi_P", "The weight for the particle's personal best position."));
+      Parameters.Add(new LookupParameter<DoubleValue>("Phi_G", "The weight for the global best position."));
     }
 
+    public override IDeepCloneable Clone(Cloner cloner) {
+      return new ParticleUpdater(this, cloner);
+    }
+
+    #endregion
+
     public override IOperation Apply() {
-      double alpha = ((IRandom)RandomParameter.ActualValue).NextDouble();
-      double beta = ((IRandom)RandomParameter.ActualValue).NextDouble();
-      RealVector velocity = (RealVector)VelocityParameter.ActualValue;
+      RealVector velocity = new RealVector(Velocity.Length);
+      RealVector position = new RealVector(Point.Length);
+      double r_p = Random.NextDouble();
+      double r_g = Random.NextDouble();
       for (int i = 0; i < velocity.Length; i++) {
-        velocity[i] = velocity[i] + alpha * (BestLocalParameter.ActualValue[i] - CurrentPositionParameter.ActualValue[i]) + beta * (BestGlobalParameter.ActualValue[i] - CurrentPositionParameter.ActualValue[i]);
+        velocity[i] =
+          Velocity[i] * Omega +
+          (PersonalBestPoint[i] - Point[i]) * Phi_P * r_p +
+          (BestPoint[i] - Point[i]) * Phi_G * r_g;
       }
-      VelocityParameter.ActualValue = velocity;
-      for (int i = 0; i < CurrentPositionParameter.ActualValue.Length; i++) {
-        CurrentPositionParameter.ActualValue[i] = CurrentPositionParameter.ActualValue[i] + VelocityParameter.ActualValue[i];
-        if (CurrentPositionParameter.ActualValue[i] < BoundsParameter.ActualValue[0, 0]) {
-          CurrentPositionParameter.ActualValue[i] = BoundsParameter.ActualValue[0, 0];
-        } else if (CurrentPositionParameter.ActualValue[i] > BoundsParameter.ActualValue[0, 1]) {
-          CurrentPositionParameter.ActualValue[i] = BoundsParameter.ActualValue[0, 1];
-        }
+      BoundsChecker.Apply(velocity, VelocityBounds);
+      for (int i = 0; i < velocity.Length; i++) {
+        position[i] = Point[i] + velocity[i];
       }
+      BoundsChecker.Apply(position, Bounds);
+      Point = position;
+      Velocity = velocity;
+
       return base.Apply();
     }

@@ -3 +3 @@
 obj
 *.user
+*.suo