Ejemplo n.º 1
0
Signal SignalSpace::inner_product(const Signal& a, const Signal& b,
                                  double t_s, double t_e,
                                  double dt, int cnt)
{
    assert(t_e > t_s);
    double t = 0.0, sum = 0.0;
    int i = 0, len = 0;
    Signal c(1 / dt, t_e - t_s);
    len = c.get_len();

    for (t = t_s, i = 0; t < t_e; t += dt, ++i) {
        if (cnt && i % cnt == 0)
            sum = 0.0;
        sum += a.at_t(t) * b.at_t(t) * dt;
        if (i < len)
            c[i] = sum;
    }
    return c;
}
Ejemplo n.º 2
0
Signal Signal::operator+(const Signal& obj) const
{
  unsigned length = 0;
  double sfreq = 0.0, dura = 0.0;

  assert(type != SIG_COS && obj.type != SIG_COS);

  if (type == SIG_SAMP && obj.type == SIG_SAMP)
    assert(samp_freq == obj.samp_freq && len == obj.len);

  if (type == SIG_SAMP && obj.type == SIG_SAMP) {
    length = len;
    sfreq = samp_freq;
    dura = duration;
  } else if (type == SIG_SAMP) {
    length = len;
    sfreq = samp_freq;
    dura = duration;
  } else {
    length = obj.len;
    sfreq = obj.samp_freq;
    dura = obj.duration;
  }

  Signal c(sfreq, dura);

  for (unsigned i = 0; i < length; ++i) {
    if (type == SIG_SAMP && obj.type == SIG_SAMP)
      c.sample[i] = sample[i] + obj.sample[i];
    else if (type == SIG_SAMP)
      c.sample[i] = sample[i] + obj.at_t(i / obj.samp_freq);
    else
      c.sample[i] = at_t(i / samp_freq) + obj.sample[i];
  }
  return c;
}