Free cookie consent management tool by TermsFeed Policy Generator

source: branches/3040_VectorBasedGP/HeuristicLab.Problems.Instances.DataAnalysis/3.3/Regression/Feynman/Feynman36.cs @ 18058

Last change on this file since 18058 was 17805, checked in by gkronber, 4 years ago

#3075 Use the same noise levels and calculation as in our experiments for the IEEE TeC paper. Reordered instances by name first and noise level second. Removed number of samples from the name.

File size: 3.0 KB
Line 
1using System;
2using System.Collections.Generic;
3using System.Linq;
4using HeuristicLab.Common;
5using HeuristicLab.Random;
6
7namespace HeuristicLab.Problems.Instances.DataAnalysis {
8  public class Feynman36 : FeynmanDescriptor {
9    private readonly int testSamples;
10    private readonly int trainingSamples;
11
12    public Feynman36() : this((int) DateTime.Now.Ticks, 10000, 10000, null) { }
13
14    public Feynman36(int seed) {
15      Seed            = seed;
16      trainingSamples = 10000;
17      testSamples     = 10000;
18      noiseRatio      = null;
19    }
20
21    public Feynman36(int seed, int trainingSamples, int testSamples, double? noiseRatio) {
22      Seed                 = seed;
23      this.trainingSamples = trainingSamples;
24      this.testSamples     = testSamples;
25      this.noiseRatio      = noiseRatio;
26    }
27
28    public override string Name {
29      get {
30        return string.Format("I.34.14 (1+v/c)/sqrt(1-v**2/c**2)*omega_0 | {0}",
31          noiseRatio == null ? "no noise" : string.Format(System.Globalization.CultureInfo.InvariantCulture, "noise={0:g}",noiseRatio));
32      }
33    }
34
35    protected override string TargetVariable { get { return noiseRatio == null ? "omega" : "omega_noise"; } }
36
37    protected override string[] VariableNames {
38      get { return new[] {"c", "v", "omega_0", noiseRatio == null ? "omega" : "omega_noise"}; }
39    }
40
41    protected override string[] AllowedInputVariables { get { return new[] {"c", "v", "omega_0"}; } }
42
43    public int Seed { get; private set; }
44
45    protected override int TrainingPartitionStart { get { return 0; } }
46    protected override int TrainingPartitionEnd { get { return trainingSamples; } }
47    protected override int TestPartitionStart { get { return trainingSamples; } }
48    protected override int TestPartitionEnd { get { return trainingSamples + testSamples; } }
49
50    protected override List<List<double>> GenerateValues() {
51      var rand = new MersenneTwister((uint) Seed);
52
53      var data    = new List<List<double>>();
54      var c       = ValueGenerator.GenerateUniformDistributedValues(rand.Next(), TestPartitionEnd, 3, 10).ToList();
55      var v       = ValueGenerator.GenerateUniformDistributedValues(rand.Next(), TestPartitionEnd, 1, 2).ToList();
56      var omega_0 = ValueGenerator.GenerateUniformDistributedValues(rand.Next(), TestPartitionEnd, 1, 5).ToList();
57
58      var omega = new List<double>();
59
60      data.Add(c);
61      data.Add(v);
62      data.Add(omega_0);
63      data.Add(omega);
64
65      for (var i = 0; i < c.Count; i++) {
66        var res = (1 + v[i] / c[i]) / Math.Sqrt(1 - Math.Pow(v[i], 2) / Math.Pow(c[i], 2)) * omega_0[i];
67        omega.Add(res);
68      }
69
70      if (noiseRatio != null) {
71        var omega_noise = new List<double>();
72        var sigma_noise = (double) Math.Sqrt(noiseRatio.Value) * omega.StandardDeviationPop();
73        omega_noise.AddRange(omega.Select(md => md + NormalDistributedRandom.NextDouble(rand, 0, sigma_noise)));
74        data.Remove(omega);
75        data.Add(omega_noise);
76      }
77
78      return data;
79    }
80  }
81}
Note: See TracBrowser for help on using the repository browser.