CQpub0 Mikami / Mbed 2 deprecated FreqConv

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

Dependencies:   DSProcessingIO mbed

HilbertTransform.hpp@2:0b15cd2b79e3, 2014-12-27 (annotated)

Committer:
CQpub0Mikami
Date:
Sat Dec 27 07:01:43 2014 +0000
Revision:
2:0b15cd2b79e3
Parent:
0:a5a171eda3f8 
3

Who changed what in which revision?

UserRevisionLine numberNew contents of line
CQpub0Mikami 0:a5a171eda3f8 1 //--------------------------------------------------------------
CQpub0Mikami 0:a5a171eda3f8 2 // Hilbert transform filter
CQpub0Mikami 0:a5a171eda3f8 3 // Copyright (c) 2014 MIKAMI, Naoki, 2014/06/23
CQpub0Mikami 0:a5a171eda3f8 4 //--------------------------------------------------------------
CQpub0Mikami 0:a5a171eda3f8 5
CQpub0Mikami 0:a5a171eda3f8 6 #ifndef HILBERT_TRANSFORM_HPP
CQpub0Mikami 0:a5a171eda3f8 7 #define HILBERT_TRANSFORM_HPP
CQpub0Mikami 0:a5a171eda3f8 8
CQpub0Mikami 0:a5a171eda3f8 9 #include "mbed.h"
CQpub0Mikami 0:a5a171eda3f8 10
CQpub0Mikami 0:a5a171eda3f8 11 namespace Mikami
CQpub0Mikami 0:a5a171eda3f8 12 {
CQpub0Mikami 0:a5a171eda3f8 13 template<int order> class Hilbert
CQpub0Mikami 0:a5a171eda3f8 14 {
CQpub0Mikami 0:a5a171eda3f8 15 private:
CQpub0Mikami 0:a5a171eda3f8 16 const float *const hm_; // pointer for filter coefficients
CQpub0Mikami 0:a5a171eda3f8 17 float xn_[order+1]; // buffer for inputs
CQpub0Mikami 0:a5a171eda3f8 18
CQpub0Mikami 0:a5a171eda3f8 19 Hilbert(const Hilbert&);
CQpub0Mikami 0:a5a171eda3f8 20 Hilbert& operator=(const Hilbert&);
CQpub0Mikami 0:a5a171eda3f8 21 public:
CQpub0Mikami 0:a5a171eda3f8 22 Hilbert(const float hk[]) : hm_(hk)
CQpub0Mikami 0:a5a171eda3f8 23 {
CQpub0Mikami 0:a5a171eda3f8 24 for (int k=0; k<=order; k++) xn_[k] = 0.0;
CQpub0Mikami 0:a5a171eda3f8 25 }
CQpub0Mikami 0:a5a171eda3f8 26
CQpub0Mikami 0:a5a171eda3f8 27 // yI: in-phase signal
CQpub0Mikami 0:a5a171eda3f8 28 // yQ: quadrature signal
CQpub0Mikami 0:a5a171eda3f8 29 void Execute(float xin, float& yI, float& yQ)
CQpub0Mikami 0:a5a171eda3f8 30 {
CQpub0Mikami 0:a5a171eda3f8 31 yQ = 0.0;
CQpub0Mikami 0:a5a171eda3f8 32 xn_[0] = xin;
CQpub0Mikami 0:a5a171eda3f8 33
CQpub0Mikami 0:a5a171eda3f8 34 for (int k=0; k<=order/4; k++)
CQpub0Mikami 0:a5a171eda3f8 35 yQ = yQ + hm_[k]*(xn_[2*k] - xn_[order-2*k]);
CQpub0Mikami 0:a5a171eda3f8 36 yI = xn_[order/2]; // in-phase signal
CQpub0Mikami 0:a5a171eda3f8 37
CQpub0Mikami 0:a5a171eda3f8 38 for (int k=order; k>0; k--)
CQpub0Mikami 0:a5a171eda3f8 39 xn_[k] = xn_[k-1]; // move input signals
CQpub0Mikami 0:a5a171eda3f8 40 }
CQpub0Mikami 0:a5a171eda3f8 41 };
CQpub0Mikami 0:a5a171eda3f8 42 }
CQpub0Mikami 0:a5a171eda3f8 43 #endif // HILBERT_TRANSFORM_HPP
CQpub0Mikami 0:a5a171eda3f8 44