Free cookie consent management tool by TermsFeed Policy Generator

source: branches/3116_GAM_Interactions/HeuristicLab.Problems.Instances.DataAnalysis/3.3/Regression/Feynman/Feynman55.cs @ 17932

Last change on this file since 17932 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.2 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 Feynman55 : FeynmanDescriptor {
9    private readonly int testSamples;
10    private readonly int trainingSamples;
11
12    public Feynman55() : this((int) DateTime.Now.Ticks, 10000, 10000, null) { }
13
14    public Feynman55(int seed) {
15      Seed            = seed;
16      trainingSamples = 10000;
17      testSamples     = 10000;
18      noiseRatio      = null;
19    }
20
21    public Feynman55(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("II.6.11 1/(4*pi*epsilon)*p_d*cos(theta)/r**2 | {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 ? "Volt" : "Volt_noise"; } }
36
37    protected override string[] VariableNames {
38      get { return new[] {"epsilon", "p_d", "theta", "r", noiseRatio == null ? "Volt" : "Volt_noise"}; }
39    }
40
41    protected override string[] AllowedInputVariables { get { return new[] {"epsilon", "p_d", "theta", "r"}; } }
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 epsilon = ValueGenerator.GenerateUniformDistributedValues(rand.Next(), TestPartitionEnd, 1, 3).ToList();
55      var p_d     = ValueGenerator.GenerateUniformDistributedValues(rand.Next(), TestPartitionEnd, 1, 3).ToList();
56      var theta   = ValueGenerator.GenerateUniformDistributedValues(rand.Next(), TestPartitionEnd, 1, 3).ToList();
57      var r       = ValueGenerator.GenerateUniformDistributedValues(rand.Next(), TestPartitionEnd, 1, 3).ToList();
58
59      var Volt = new List<double>();
60
61      data.Add(epsilon);
62      data.Add(p_d);
63      data.Add(theta);
64      data.Add(r);
65      data.Add(Volt);
66
67      for (var i = 0; i < epsilon.Count; i++) {
68        var res = 1.0 / (4 * Math.PI * epsilon[i]) * p_d[i] * Math.Cos(theta[i]) / Math.Pow(r[i], 2);
69        Volt.Add(res);
70      }
71
72      if (noiseRatio != null) {
73        var Volt_noise  = new List<double>();
74        var sigma_noise = (double) Math.Sqrt(noiseRatio.Value) * Volt.StandardDeviationPop();
75        Volt_noise.AddRange(Volt.Select(md => md + NormalDistributedRandom.NextDouble(rand, 0, sigma_noise)));
76        data.Remove(Volt);
77        data.Add(Volt_noise);
78      }
79
80      return data;
81    }
82  }
83}
Note: See TracBrowser for help on using the repository browser.