public GEPIndividual(int headLen, FunctionSet functionSet, double minVal, double maxVal, int optimizeGenerationCount)
{
this.functionSet = functionSet;
this.optimizeGenerationCount = optimizeGenerationCount;
// 计算尾部长度
int tailLen = headLen * (functionSet.MaxParamCount - 1) + 1;
// 随机生成头部
gene.Add(functionSet.RandomNonterminator());
for (int i = 1; i < headLen; ++i)
{
gene.Add(functionSet.RandomFunction());
}
// 随机生成尾部
for (int i = 0; i < tailLen; ++i)
{
gene.Add(functionSet.RandomTerminator());
}
// 随机生成系数
for (int i = 0; i < tailLen; ++i)
{
coeff.Add(RandomUtil.U(minVal, maxVal));
}
}
public GEPIndividual(List <Function> gene, List <double> coeff, FunctionSet functionSet, int optimizeGenerationCount)
{
this.gene = gene;
this.coeff = coeff;
this.functionSet = functionSet;
this.optimizeGenerationCount = optimizeGenerationCount;
}
static void Main()
{
// 数据点
List <Point> points = new List <Point>
{
// minErr: 0.90
//new Point(-5.25, -2.26),
//new Point(-4.22, 1.19),
//new Point(-3.55, 2.72),
//new Point(-0.73, 2.74),
//new Point(0.45, 0.78),
//new Point(1.07, -1.26),
//new Point(-4.95, -1.12),
//new Point(-2.09, 3.6),
// minErr: 0.90
//new Point(-7.46, -1.42),
//new Point(-5.42, -0.68),
//new Point(-2.92, 0.6),
//new Point(-0.86, 2.98),
//new Point(0.66, 5.62),
//new Point(1.34, 7.48),
//new Point(1.72, 9.6),
//new Point(1.92, 10.86),
// minErr: 0.30
//new Point(-4.7, -1.14),
//new Point(-4.54, -0.35),
//new Point(-4.28, 0.44),
//new Point(-3.86, 1.06),
//new Point(-3.32, 1.2),
//new Point(-2.44, 0.83),
//new Point(-0.92, -0.58),
//new Point(-1.73, 0.27),
// minErr: 1.42
//new Point(-2.1, 4.38),
//new Point(-1.86, 2.7),
//new Point(-1.34, 1.16),
//new Point(0.92, 0.36),
//new Point(2.82, 0.82),
//new Point(4.62, 1.62),
//new Point(5.9, 2.82),
//new Point(6.82, 4),
// minErr: 0.60
//new Point(-2.92, 4.65),
//new Point(-2, 2),
//new Point(-1.02, 0.4),
//new Point(0, -1),
//new Point(1.54, -1.83),
//new Point(4.25, -2.54),
//new Point(6.51, -2.93),
//new Point(8.77, -3.07),
// minErr: 0.39
//new Point(-1.82, -6.62),
//new Point(-0.44, -5.08),
//new Point(2, -3),
//new Point(3.36, -1.56),
//new Point(5.52, 0.66),
//new Point(7.1, 2.44),
//new Point(7.9, 3.28),
//new Point(8.58, 4),
// minErr: 0.48
//new Point(1.9, 0.58),
//new Point(2.38, 1.68),
//new Point(2.78, 2.92),
//new Point(3.42, 4.52),
//new Point(4, 6.16),
//new Point(4.52, 7.34),
//new Point(4.98, 8.76),
//new Point(5.36, 9.76),
// minErr: 0.01
//new Point(-4.09, 0.81),
//new Point(-3.02, -0.12),
//new Point(-1.68, -0.99),
//new Point(-0.5, -0.48),
//new Point(0.43, 0.42),
//new Point(1.5, 1),
//new Point(2.51, 0.59),
//new Point(3.91, -0.69),
// minErr: 0.21
new Point(-3.33, 0.14),
new Point(-2.62, 0.23),
new Point(-1.73, 0.41),
new Point(-0.73, 0.74),
new Point(0.28, 1.23),
new Point(1.15, 2.09),
new Point(1.61, 3.23),
new Point(1.94, 4.19),
};
// 函数集
FunctionSet functionSet = new FunctionSet();
functionSet
.AddFunction(new Var())
.AddFunction(new Const())
.AddFunction(new Add())
.AddFunction(new Sub())
.AddFunction(new Mul())
.AddFunction(new Div())
.AddFunction(new Sin())
.AddFunction(new Exp())
;
// 生成初始种群
List <IIndividual> initPop = new GEPPopulationGenerator
{
PopSize = 200,
HeadLen = 10,
OptimizeGenerationCount = 500,
Data = points,
FunctionSet = functionSet,
ConstMinVal = -5,
ConstMaxVal = 5,
}.Generate();
// 进化
new Evolver1
{
Data = points,
MaxCount = 100,
}.Evolve(initPop);
}
