AM中波放送用SDR.CICフィルタのみを使用.CQ出版社「トランジスタ技術」誌,2021年4月号に掲載

Dependencies:   mbed

Committer:
MikamiUitOpen
Date:
Wed Dec 01 00:37:44 2021 +0000
Revision:
3:878a48f15e89
Parent:
2:4bec6b2be809
4

Who changed what in which revision?

UserRevisionLine numberNew contents of line
MikamiUitOpen 1:30d9fb51dec1 1 //-------------------------------------------------------------
MikamiUitOpen 1:30d9fb51dec1 2 // ダウンサンプリングで使う3段の CIC フィルタ用クラス
MikamiUitOpen 1:30d9fb51dec1 3 //
MikamiUitOpen 2:4bec6b2be809 4 // 2020/08/02, Copyright (c) 2020 MIKAMI, Naoki
MikamiUitOpen 1:30d9fb51dec1 5 //-------------------------------------------------------------
MikamiUitOpen 1:30d9fb51dec1 6
MikamiUitOpen 1:30d9fb51dec1 7 #include "mbed.h"
MikamiUitOpen 1:30d9fb51dec1 8
MikamiUitOpen 1:30d9fb51dec1 9 #ifndef CIC3_FILTER_CLASS_HPP
MikamiUitOpen 1:30d9fb51dec1 10 #define CIC3_FILTER_CLASS_HPP
MikamiUitOpen 1:30d9fb51dec1 11
MikamiUitOpen 1:30d9fb51dec1 12 namespace Mikami
MikamiUitOpen 1:30d9fb51dec1 13 {
MikamiUitOpen 1:30d9fb51dec1 14 class Cic3Stage
MikamiUitOpen 1:30d9fb51dec1 15 {
MikamiUitOpen 1:30d9fb51dec1 16 public:
MikamiUitOpen 1:30d9fb51dec1 17 // コンストラクタ
MikamiUitOpen 2:4bec6b2be809 18 // rate ダウンサンプリングの率,1/10 にダウンサンプリングする場合は 10
MikamiUitOpen 2:4bec6b2be809 19 // amp 直交信号発生器の出力の振幅
MikamiUitOpen 1:30d9fb51dec1 20 Cic3Stage(int rate, float amp)
MikamiUitOpen 1:30d9fb51dec1 21 : G0_(1.0f/(amp*rate*rate*rate)),
MikamiUitOpen 1:30d9fb51dec1 22 vn1_(0), vn2_(0), vn3_(0), vn3M1_(0), yn1M1_(0), yn2M1_(0) {}
MikamiUitOpen 1:30d9fb51dec1 23
MikamiUitOpen 2:4bec6b2be809 24 // 積分器
MikamiUitOpen 2:4bec6b2be809 25 void Integrate(int16_t xn)
MikamiUitOpen 1:30d9fb51dec1 26 {
MikamiUitOpen 2:4bec6b2be809 27 vn1_ += xn; // 積分器1段目
MikamiUitOpen 2:4bec6b2be809 28 vn2_ += vn1_; // 積分器2段目
MikamiUitOpen 2:4bec6b2be809 29 vn3_ += vn2_; // 積分器3段目
MikamiUitOpen 1:30d9fb51dec1 30 }
MikamiUitOpen 1:30d9fb51dec1 31
MikamiUitOpen 2:4bec6b2be809 32 // くし形フィルタ
MikamiUitOpen 2:4bec6b2be809 33 float CombFilter()
MikamiUitOpen 1:30d9fb51dec1 34 {
MikamiUitOpen 2:4bec6b2be809 35 int32_t yn1 = vn3_ - vn3M1_; // くし形フィルタ1段目
MikamiUitOpen 2:4bec6b2be809 36 int32_t yn2 = yn1 - yn1M1_; // くし形フィルタ2段目
MikamiUitOpen 2:4bec6b2be809 37 int32_t yn3 = yn2 - yn2M1_; // くし形フィルタ3段目
MikamiUitOpen 1:30d9fb51dec1 38
MikamiUitOpen 2:4bec6b2be809 39 vn3M1_ = vn3_; // 現在の値を保存,くし形フィルタ1段目
MikamiUitOpen 2:4bec6b2be809 40 yn1M1_ = yn1; // 現在の値を保存,くし形フィルタ2段目
MikamiUitOpen 2:4bec6b2be809 41 yn2M1_ = yn2; // 現在の値を保存,くし形フィルタ3段目
MikamiUitOpen 1:30d9fb51dec1 42
MikamiUitOpen 1:30d9fb51dec1 43 return G0_*yn3;
MikamiUitOpen 1:30d9fb51dec1 44 }
MikamiUitOpen 1:30d9fb51dec1 45
MikamiUitOpen 1:30d9fb51dec1 46 private:
MikamiUitOpen 1:30d9fb51dec1 47 const float G0_;
MikamiUitOpen 2:4bec6b2be809 48 int32_t vn1_, vn2_, vn3_; // 積分器で使う変数
MikamiUitOpen 2:4bec6b2be809 49 int32_t vn3M1_, yn1M1_, yn2M1_; // くし形フィルタで使う変数
MikamiUitOpen 1:30d9fb51dec1 50
MikamiUitOpen 2:4bec6b2be809 51 // コピー・コンストラクタ,代入演算子の禁止のため
MikamiUitOpen 1:30d9fb51dec1 52 Cic3Stage(const Cic3Stage&);
MikamiUitOpen 1:30d9fb51dec1 53 Cic3Stage& operator=(const Cic3Stage&);
MikamiUitOpen 1:30d9fb51dec1 54 };
MikamiUitOpen 1:30d9fb51dec1 55 }
MikamiUitOpen 2:4bec6b2be809 56 #endif // CIC3_FILTER_CLASS_HPP