content: /**
   * Copyright 2000-2009 DFKI GmbH.
   * All Rights Reserved.  Use is subject to license terms.
   * <p/>
   * This file is part of MARY TTS.
   * <p/>
   * MARY TTS is free software: you can redistribute it and/or modify
   * it under the terms of the GNU Lesser General Public License as published by
   * the Free Software Foundation, version 3 of the License.
   * <p/>
   * This program is distributed in the hope that it will be useful,
   * but WITHOUT ANY WARRANTY; without even the implied warranty of
   * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   * GNU Lesser General Public License for more details.
   * <p/>
   * You should have received a copy of the GNU Lesser General Public License
   * along with this program.  If not, see <http://www.gnu.org/licenses/>.
   */
  package marytts.signalproc.sinusoidal.test;
  
  import java.io.IOException;
  
  import marytts.signalproc.sinusoidal.Sinusoid;
  
  /**
   * Generates test signals that consist of artificially generated harmonics
   *
   * @author oytun.turk
   */
  public class HarmonicsTester extends SinusoidsTester {
  
      public HarmonicsTester(float freqInHz, int numHarmonics) {
          this(freqInHz, numHarmonics, DEFAULT_DUR, DEFAULT_FS);
      }
  
      public HarmonicsTester(float freqInHz, int numHarmonics, float durationInSeconds) {
          this(freqInHz, numHarmonics, durationInSeconds, DEFAULT_FS);
      }
  
      public HarmonicsTester(float freqInHz, int numHarmonics, float durationInSeconds, int samplingRateInHz) {
          this(freqInHz, numHarmonics, durationInSeconds, samplingRateInHz, null, null);
      }
  
      public HarmonicsTester(float freqInHz, int numHarmonics, float[] amps) {
          this(freqInHz, numHarmonics, amps, null);
      }
  
      public HarmonicsTester(float freqInHz, int numHarmonics, float[] amps, float[] phases) {
          this(freqInHz, numHarmonics, DEFAULT_DUR, DEFAULT_FS, amps, phases);
      }
  
      public HarmonicsTester(float freqInHz, int numHarmonics, float durationInSeconds, float[] amps) {
          this(freqInHz, numHarmonics, durationInSeconds, DEFAULT_FS, amps, null);
      }
  
      public HarmonicsTester(float freqInHz, int numHarmonics, float durationInSeconds, int samplingRateInHz, float[] amps) {
          this(freqInHz, numHarmonics, durationInSeconds, samplingRateInHz, amps, null);
      }
  
      public HarmonicsTester(float freqInHz, int numHarmonics, float durationInSeconds, int samplingRateInHz, float[] amps, float[] phases) {
          if (numHarmonics > 0) {
              Sinusoid[] sins = new Sinusoid[numHarmonics];
  
              float currentAmp, currentPhase;
              for (int i = 0; i < numHarmonics; i++) {
                  if (amps != null && amps.length > i)
                      currentAmp = amps[i];
                  else
                      currentAmp = DEFAULT_AMP;
  
                  if (phases != null && phases.length > i)
                      currentPhase = phases[i];
                  else
                      currentPhase = DEFAULT_PHASE;
  
                  sins[i] = new Sinusoid(currentAmp, freqInHz * (i + 1), currentPhase);
              }
  
              init(sins, durationInSeconds, samplingRateInHz);
          }
      }
  
      public HarmonicsTester(float freqInHz, int numHarmonics, float startTimeInSeconds, float endTimeInSeconds) {
          this(freqInHz, numHarmonics, startTimeInSeconds, endTimeInSeconds, DEFAULT_FS);
      }
  
      public HarmonicsTester(float freqInHz, int numHarmonics, float startTimeInSeconds, float endTimeInSeconds, int samplingRateInHz) {
          this(freqInHz, numHarmonics, startTimeInSeconds, endTimeInSeconds, samplingRateInHz, null);
      }
  
      public HarmonicsTester(float freqInHz, int numHarmonics, float startTimeInSeconds, float endTimeInSeconds, int samplingRateInHz, float[] amps) {
          this(freqInHz, numHarmonics, startTimeInSeconds, endTimeInSeconds, samplingRateInHz, amps, null);
      }
  
      public HarmonicsTester(float freqInHz, int numHarmonics, float startTimeInSeconds, float endTimeInSeconds, int samplingRateInHz, float[] amps, float[] phases) {
          if (numHarmonics > 0) {
              Sinusoid[] sins = new Sinusoid[numHarmonics];
  
              float currentAmp, currentPhase;
              for (int i = 0; i < numHarmonics; i++) {
                  if (amps != null && amps.length > i)
                      currentAmp = amps[i];
                  else
                      currentAmp = DEFAULT_AMP;
  
                  if (phases != null && phases.length > i)
                      currentPhase = phases[i];
                  else
                      currentPhase = DEFAULT_PHASE;
  
                  sins[i] = new Sinusoid(currentAmp, freqInHz * (i + 1), currentPhase);
              }
  
              init(sins, startTimeInSeconds, endTimeInSeconds, samplingRateInHz);
          }
      }
  
      public static void main(String[] args) throws IOException {
          HarmonicsTester s = null;
  
          //Single sinusoid, time-invariant
          float f1 = 115.0f;
          int numHarmonics = 8;
          double defaultAbsMaxVal = 13000.0;
          //s = new HarmonicsTester(f1, numHarmonics);
          float[] amps = new float[numHarmonics];
          for (int i = 0; i < numHarmonics; i++)
              amps[i] = 0.5f / numHarmonics; //(float)Math.pow(2.0, -1.0*i);
          s = new HarmonicsTester(f1, numHarmonics, amps);
          //
  
          if (args.length > 1)
              s.write(args[0], args[1], defaultAbsMaxVal);
          else if (args.length == 1)
              s.write(args[0], defaultAbsMaxVal);
          else {
              String outWavFile = "d:\\h" + String.valueOf(numHarmonics) + "_" + String.valueOf(f1) + "_" + String.valueOf(defaultAbsMaxVal) + ".wav";
              s.write(outWavFile, defaultAbsMaxVal);
          }
  
      }
  }
  
  