vibrato.h

Go to the documentation of this file.

/*
    Harmonic sinusoidal modelling and tools

    C++ code package for harmonic sinusoidal modelling and relevant signal processing.
    Centre for Digital Music, Queen Mary, University of London.
    This file copyright 2011 Wen Xue.

    This program is free software; you can redistribute it and/or
    modify it under the terms of the GNU General Public License as
    published by the Free Software Foundation; either version 2 of the
    License, or (at your option) any later version.
*/
#ifndef vibratoH
#define vibratoH


#include <stdio.h>
#include "hs.h"
//---------------------------------------------------------------------------

struct TVo
{
  //basic characteristics
  int M;            
  int L;            
  int P;            
  double h;         
  double* F0C;      
  double* A0C;      
  int afres;        
  double LogAF[4096]; 

  int* peakfr;      
  double* lp;       
  double* Dp;       
  int* Kp;          
  double** fxr;     
  double** LogASp;  
  int K;

  //other properties
  double rate;      
  double regularity;
  double F0max;     
  double F0min;     
  double F0Cmax;    
  double F0Cmin;    
  double F0Overall; 
  double F0Dmax;    
  double F0Dmin;    
  double* F0;       
  double* F0D;      
  double D;         
  double LogAS[100];
  double FXR[100];  
  double FRes[50];  
  double** fres;    

  TVo();
  ~TVo();
  void AllocateL(int AnL);
  void ReAllocateL(int newL);
  void AllocateP();
  void AllocatePK();
  void SaveToFile(char* filename);
  void SaveToFileHandle(FILE* file);
  void LoadFromFileHandle(FILE* file);
        double GetOldP();
};

//--tool procedures----------------------------------------------------------
double QIE(double* y);
double QIE(double* y, double& x);

//--demodulation routines----------------------------------------------------
void DeFM(double* f2, double* f1, double* AA, int npfr, int* peakfr, int Fr, double& fct, int& fcount, double* &frs, double* &f, bool furthersmoothing=false);
void DeFM(double* f2, double* f1, double* AA, int npfr, double* lp, int Fr, double& fct, int& fcount, double* &frs, double* &f, bool furthersmoothing=false);
void DeAM(double* A2, double* A1, int npfr, int* peakfr, int Fr);
void DeAM(double* A2, double* A1, int npfr, double* lp, int Fr);

//--source-filter analysis routine-------------------------------------------
double S_F(int afres, double* LogAF, double* LogAS, int Fr, int M, atom** Partials, double* A0C, double* lp, int P, double F0Overall);
double S_F_SV(int afres, double* LogAF, double* LogAS, int Fr, int M, atom** Partials, double* A0C, double* lp, int P, double F=0.005, int FScaleMode=0, double theta=0.5, double sps=44100);

//--vibrato analysis and synthesis routines----------------------------------
void AnalyzeV(THS& HS, TVo& V, double*& peaky, double*& cyclefrs, double*& cyclefs, double sps, double h, int* cyclefrcount=0,
      int SFMode=0, double SFF=0.01, int SFFScale=0, double SFtheta=0);
void RegularizeV(THS& HS, TVo& V, double sps, double h);
void SynthesizeV(THS* HS, TVo* V, double sps, int UseK=0);
TVo* InterpolateV(double newP, double rate, TVo& V);

//--other functions----------------------------------------------------------
void FindPeaks(int* peakfr, int& npfr, double* F0, int Fr, double periodinframe, double*& peaky);
void FS_QR(int& K, double* FXR, double* data, int Fr, double period, double shift, double* FRES);
void RateAndReg(double& rate, double& regularity, double* data, int frst, int fren, int padrate, double sps, double offst);

#endif