SquareWave.hh

 
//
// Copyright (C) 2011-2022 Yuki YOKOKURA
// This program is free software;
// you can redistribute it and/or modify it under the terms of the FreeBSD License.
// For details, see the License.txt file.
 
#ifndef SQUAREWAVE
#define SQUAREWAVE
 
#include "Matrix.hh"
 
namespace ARCS {    // ARCS名前空間
double SquareWave(const double freq, const double phase, const double time);    
double SquareWave(const double freq, const double phase, const double time, const double starttime);    
 
template <size_t M>
void SquareWave(const double freq, const double phase, const double time, Matrix<1,M>& output){
    double r = sin(2.0*M_PI*freq*time + phase);
    if(0 < r){
        output =  Matrix<1,M>::ones();
    }else{
        output = -Matrix<1,M>::ones();
    }
}
 
template <size_t M>
Matrix<1,M> SquareWaveVec(const double freq, const double phase, const double time){
    Matrix<1,M> ret;
    SquareWave(freq, phase, time, ret);
    return ret;
}
 
}
 
#endif