Diff: SDR_Library/QuadOscIir.hpp

Revision:

2:4bec6b2be809

Parent:

1:30d9fb51dec1

--- a/SDR_Library/QuadOscIir.hpp	Mon Sep 23 07:32:10 2019 +0000
+++ b/SDR_Library/QuadOscIir.hpp	Sat Aug 29 11:19:03 2020 +0000
@@ -1,7 +1,7 @@
 //--------------------------------------------------------------
-//  IIR フィルタを利用する二相発振器のクラス
+//  IIR フィルタを利用する直交信号発生器のクラス
 //
-//  2019/04/11, Copyright (c) 2019 MIKAMI, Naoki
+//  2020/07/30, Copyright (c) 2020 MIKAMI, Naoki
 //--------------------------------------------------------------
 
 #include "mbed.h"
@@ -15,59 +15,44 @@
     {
     public:
         // コンストラクタ
-        //      f0：     周波数，単位：Hz
-        //      ts：     標本化間隔，単位：μs
-        //      amp：    振幅
+        //      f0  周波数，単位：Hz
+        //      ts  標本化間隔，単位：μs
+        //      amp 振幅
         QuadOscIir(float f0, float ts, float amp = 1)
             : TS_(ts*1.0E-6f), A0_(amp) { Set(f0); }
 
-        // sin/cos 発生，実行時間の実測値：91 ns
+        // sin/cos 発生
         void Generate(float &sinX, float &cosX)
         {
-            RecursivePart();
-            sinX = b1S_*un2_;
-            cosX = A0_*un1_ + b1C_*un2_;
-        }
-    
-        // sin 発生，実行時間の実測値：69 ns
-        float GenerateSin()
-        {
-            RecursivePart();
-            return b1S_*un2_;
-        }
+            float vn = a1_*vn1_ - vn2_;
+            sinX = b1S_*vn1_;
+            cosX = A0_*vn + b1C_*vn1_;
 
-        // cis 発生，実行時間の実測値：80 ns
-        float GenerateCos()
-        {
-            RecursivePart();
-            return A0_*un1_ + b1C_*un2_;
+            vn2_ = vn1_;    // v[n-2] ← v[n-1]
+            vn1_ = vn;      // v[n-1] ← v[n]
         }
 
         // 出力周波数の設定，周波数を変更し再スタートする場合も使える
         //      f0：周波数，単位：Hz
         void Set(float f0)
         {
-            static const float pi2 = 6.283185f;
-            a1_ = 2.0f*cosf(pi2*f0*TS_);
-            b1S_ = A0_*sinf(pi2*f0*TS_);
-            b1C_ = -A0_*cosf(pi2*f0*TS_);
-            un1_ = 1.0f;
-            un2_ = 0.0f;
+            static const float PI2 = 6.283185f;
+            a1_ = 2.0f*cosf(PI2*f0*TS_);
+            b1S_ = A0_*sinf(PI2*f0*TS_);
+            b1C_ = -A0_*cosf(PI2*f0*TS_);
+            vn1_ = 1.0f;
+            vn2_ = 0.0f;
         }
 
     private:
         const float TS_;        // 標本化間隔
         const float A0_;        // 振幅
         float a1_, b1S_, b1C_;
-        float un1_, un2_;
+        float vn1_, vn2_;
 
-        // u[n] = a1*u[n-1] - u[n-2] の計算
-        void RecursivePart()
-        {
-            float un = a1_*un1_ - un2_;
-            un2_ = un1_;    // u[n-2] ← u[n-1]
-            un1_ = un;      // u[n-1] ← u[n]
-        }
+        // コピー･コンストラクタ，代入演算子の禁止のため
+        QuadOscIir(const QuadOscIir&);
+        QuadOscIir& operator=(const QuadOscIir&);
     };
 }
-#endif  // QUAD_PHASE_OSC_IIR_HPP
+#endif  // QUAD_PHASE_OSC_IIR_HPP
\ No newline at end of file