The experiment using this program is introduced on "Interface" No.3, CQ publishing Co.,Ltd, 2015. 本プログラムを使った実験は,CQ出版社のインターフェース 2015年3月号で紹介しています.

Dependencies:   DSProcessingIO mbed

HilbertTransform.hpp

Committer:
CQpub0Mikami
Date:
2014-12-27
Revision:
2:0b15cd2b79e3
Parent:
0:a5a171eda3f8

File content as of revision 2:0b15cd2b79e3:

//--------------------------------------------------------------
// Hilbert transform filter
// Copyright (c) 2014 MIKAMI, Naoki,  2014/06/23
//--------------------------------------------------------------

#ifndef HILBERT_TRANSFORM_HPP
#define HILBERT_TRANSFORM_HPP

#include "mbed.h"

namespace Mikami
{
    template<int order> class Hilbert
    {
    private:
        const float *const hm_; // pointer for filter coefficients
        float xn_[order+1];     // buffer for inputs

        Hilbert(const Hilbert&);
        Hilbert& operator=(const Hilbert&);
    public:
        Hilbert(const float hk[]) : hm_(hk)
        {
            for (int k=0; k<=order; k++) xn_[k] = 0.0;
        }
        
        // yI: in-phase signal
        // yQ: quadrature signal
        void Execute(float xin, float& yI, float& yQ)
        {
            yQ = 0.0;
            xn_[0] = xin;

            for (int k=0; k<=order/4; k++)
                yQ = yQ + hm_[k]*(xn_[2*k] - xn_[order-2*k]);
            yI = xn_[order/2];              // in-phase signal

            for (int k=order; k>0; k--)
                xn_[k] = xn_[k-1];          // move input signals        
        }
    };
}
#endif  // HILBERT_TRANSFORM_HPP