WhiteNoiseSignal.h

#ifndef SPIKE_WHITENOISESIGNAL_H
#define SPIKE_WHITENOISESIGNAL_H
 
#include <complex>
#include <random>
#include <string>
#include <vector>
 
#include "Spike/Signal/Signal.h"
#include "Spike/TimeFrame/TimeFrame.h"
 
namespace Spike {
 
class WhiteNoiseSignal : public Signal {
private:
  double alpha;  
  double f_low;  
  double f_high; 
 
  std::random_device rd;                 
  std::mt19937 generator;                
  std::normal_distribution<double> dist; 
 
  std::vector<std::complex<double>>
      frequencies; 
 
public:
  WhiteNoiseSignal(double alpha, double f_low, double f_high,
                   const TimeFrame &time_frame);
 
  WhiteNoiseSignal(const std::string &input_file, const TimeFrame &time_frame);
 
  void calculate_signal();
 
  [[nodiscard]] double signal(double t) const;
 
  void set_alpha(double alpha_new) { alpha = alpha_new; };
 
  [[nodiscard]] const std::vector<double> &get_values() const {
    return signal_values;
  };
 
  [[nodiscard]] const std::vector<std::complex<double>> &
  get_frequencies() const {
    return frequencies;
  };
 
  [[nodiscard]] double get_alpha() const { return alpha; };
 
  [[nodiscard]] double get_variance() const {
    return 2. * (f_high - f_low) * alpha;
  };
 
  void print(std::ostream &out) const override {
    out << "WhiteNoiseSignal(alpha: " << alpha << ", f_low: " << f_low
        << ", f_high: " << f_high << ")";
  }
};
 
} // namespace Spike
 
#endif // SPIKE_WHITENOISESIGNAL_H