hssf.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 hssfH
#define hssfH

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

//---------------------------------------------------------------------------
const double Amel=1127.0104803341574386544633680278;
const bool useA0=true; //if true, use A0D+A0C instead of A0C in S-F decomposition as pre-normalizer


struct TSF
{
        //basic characteristics
        int M;            
        int L;            
        int P;            
        double offst;     
        double* F0C;      
        double* F0D;      
  double* logA0C;   
  double* logA0D;   

        double* lp;       

  double F;         
  double Fs;        
  int FScaleMode;   
  int K;            
  double** b;       
  double** h;       
  double* avgb;     
  double* avgh;     

        //other properties

  double rate;      
  double regularity;
  double F0max;     
  double F0min;     
  double F0Cmax;    
  double F0Cmin;    
  double F0Overall; 
  double F0Dmax;    
  double F0Dmin;    
  double* F0;       
  double* logA0;                

        TSF();
        ~TSF();

  //copy function
  void Duplicate(TSF& SF);

  //output routines
        double LogAF(double f);
        double LogAF(double f, int fr);
        double LogAS(int m, int fr);

  //memory handling routines
  void AllocateL(int AnL);
        void ReAllocateL(int newL);
        void AllocateP();
        void AllocateSF();

  //I/O routines
        void SaveSFToFileHandle(FILE* f);
        void SaveToFileHandle(FILE* f);
        void LoadSFFromFileHandle(FILE* f);
        void LoadFromFileHandle(FILE* f);
        void SaveToFile(char* filename);

  //other member functions
  void ShiftFilterByDB(double dB);
};

//--tool procedures----------------------------------------------------------
int Sign(double);

//--general source-filter model routines-------------------------------------
void S_F_b(TSF& SF, atom** Partials);

//--slow-variation SF estimation routines------------------------------------
double P2_DelFtr(double** d, int L, int K, double** x, double F);
double P3_DelSrc(double** d, int L, int M, int K, double** x, double** f, double F);
int SF_SV(double** h, int L, int M, int K, double** a, double** f, double F, double theta, double ep, int maxiter);
double S_F_SV(int M, atom** Partials, double* logA0C, double* lp, int P, int& K, double** h, double* avgh, double** b, double* avgb, double F=0.005, int FScaleMode=0, double theta=0.5, double Fs=44100);
void SF_SV_cf(double* h, double** b, int L, int M, int K, double** a, double** f, double F, double ep, int maxiter);
double S_F_SV_cf(int afres, double* LogAF, double* LogAS, int Fr, int M, atom** Partials, double* A0C, double* lp, int P, int& K, double** h, double** b, double F=0.005, int FScaleMode=0, double Fs=44100);

//--filter-bank SF estimation routines---------------------------------------
int SF_FB(double* hl, int M, int K, double** al, double** fl, double F, int LMode);
double S_F_FB(int M, atom** Partials, double* logA0C, double* lp, int P, int& K, double** h, double* avgh, double** b, double* avgb, double F=0.005, int FScaleMode=0, double Fs=44100);

//--source-filter analysis and synthesis routines----------------------------
void AnalyzeSF_1(THS& HS, TSF& SF, double sps, double offst);
void AnalyzeSF_2(THS& HS, TSF& SF, double*& cyclefrs, double*& cyclefs, double sps, int* cyclefrcount=0, int SFMode=0, double SFF=0.01, int SFFScale=0, double SFtheta=0);
void AnalyzeSF(THS& HS, TSF& SF, double*& cyclefrs, double*& cyclefs, double sps, double offst, int* cyclefrcount=0, int SFMode=0, double SFF=0.01, int SFFScale=0, double SFtheta=0);
void SynthesizeSF(THS* HS, TSF* SF, double sps);


#endif