Realtime spectrum analyzer. Using FFT, linear prediction, or cepstrum smoothing. Version using MEMS microphone and CODEC, named "F746_RealtimeSpectrumAnalyzer_MEMS_Mic" is registered. リアルタイム スペクトル解析器.解析の手法:FFT,線形予測法,ケプストラムによる平滑化の3種類.このプログラムの説明は,CQ出版社のインターフェース誌,2016年4月号に掲載.外付けのマイクまたは他の信号源等を A0 に接続する.線形予測法,ケプストラムは,スペクトル解析の対象を音声信号に想定してパラメータを設定している.MEMS マイクと CODEC を使ったバージョンを "F746_RealtimeSpectrumAnalyzer_MEMS_Mic" として登録.

Dependencies:   BSP_DISCO_F746NG BUTTON_GROUP LCD_DISCO_F746NG TS_DISCO_F746NG UIT_FFT_Real mbed

Committer:
MikamiUitOpen
Date:
Sun Dec 20 13:38:22 2015 +0000
Revision:
7:6598a9b70e5a
Parent:
6:b27ef8d98edc
Child:
8:1f4bc859bc84
8

Who changed what in which revision?

UserRevisionLine numberNew contents of line
MikamiUitOpen 2:095b360e0f54 1 //-----------------------------------------------------------
MikamiUitOpen 2:095b360e0f54 2 // Class for sampling input signal
MikamiUitOpen 2:095b360e0f54 3 //
MikamiUitOpen 7:6598a9b70e5a 4 // 2015/12/20, Copyright (c) 2015 MIKAMI, Naoki
MikamiUitOpen 2:095b360e0f54 5 //-----------------------------------------------------------
MikamiUitOpen 2:095b360e0f54 6
MikamiUitOpen 2:095b360e0f54 7 #include "Sampler.hpp"
MikamiUitOpen 2:095b360e0f54 8
MikamiUitOpen 2:095b360e0f54 9 namespace Mikami
MikamiUitOpen 2:095b360e0f54 10 {
MikamiUitOpen 7:6598a9b70e5a 11 Sampler::Sampler(PinName pin, int fs, int nData)
MikamiUitOpen 2:095b360e0f54 12 : TS_(1000000/fs), N_DATA_(nData), aIn_(pin),
MikamiUitOpen 7:6598a9b70e5a 13 xn_(new int16_t[nData]),
MikamiUitOpen 7:6598a9b70e5a 14 buffer_(new int16_t[nData])
MikamiUitOpen 7:6598a9b70e5a 15 {
MikamiUitOpen 2:095b360e0f54 16 Rd[0] = &Sampler::ReadNorm; // inphase
MikamiUitOpen 2:095b360e0f54 17 Rd[1] = &Sampler::ReadInv; // out-of-phase
MikamiUitOpen 2:095b360e0f54 18 sw_ = 0;
MikamiUitOpen 7:6598a9b70e5a 19 Set(false);
MikamiUitOpen 5:98ec9dd54144 20 for (int n=0; n<nData; n++)
MikamiUitOpen 5:98ec9dd54144 21 buffer_[n] = 32767;
MikamiUitOpen 2:095b360e0f54 22 }
MikamiUitOpen 2:095b360e0f54 23
MikamiUitOpen 7:6598a9b70e5a 24 Sampler::~Sampler()
MikamiUitOpen 7:6598a9b70e5a 25 {
MikamiUitOpen 7:6598a9b70e5a 26 delete[] xn_;
MikamiUitOpen 7:6598a9b70e5a 27 delete[] buffer_;
MikamiUitOpen 7:6598a9b70e5a 28 }
MikamiUitOpen 7:6598a9b70e5a 29
MikamiUitOpen 7:6598a9b70e5a 30 void Sampler::Start()
MikamiUitOpen 2:095b360e0f54 31 {
MikamiUitOpen 2:095b360e0f54 32 count_ = 0;
MikamiUitOpen 7:6598a9b70e5a 33 filled_ = false;
MikamiUitOpen 2:095b360e0f54 34 timer_.attach_us(this, &Sampler::Isr, TS_);
MikamiUitOpen 2:095b360e0f54 35 }
MikamiUitOpen 2:095b360e0f54 36
MikamiUitOpen 2:095b360e0f54 37 void Sampler::Set(bool tf)
MikamiUitOpen 2:095b360e0f54 38 {
MikamiUitOpen 2:095b360e0f54 39 trigger_ = false;
MikamiUitOpen 2:095b360e0f54 40 xnM1_ = 32767;
MikamiUitOpen 2:095b360e0f54 41 filled_ = tf;
MikamiUitOpen 2:095b360e0f54 42 }
MikamiUitOpen 2:095b360e0f54 43
MikamiUitOpen 2:095b360e0f54 44 // Interrupt service routine for Ticker
MikamiUitOpen 2:095b360e0f54 45 void Sampler::Isr()
MikamiUitOpen 2:095b360e0f54 46 {
MikamiUitOpen 2:095b360e0f54 47 if (!trigger_)
MikamiUitOpen 2:095b360e0f54 48 {
MikamiUitOpen 2:095b360e0f54 49 int16_t xn = (this->*Rd[sw_])();
MikamiUitOpen 2:095b360e0f54 50 // Detect rising edge
MikamiUitOpen 2:095b360e0f54 51 if ((xn > (xnM1_+512)) && (xn > 2048))
MikamiUitOpen 2:095b360e0f54 52 trigger_ = true;
MikamiUitOpen 2:095b360e0f54 53 else
MikamiUitOpen 2:095b360e0f54 54 xnM1_ = xn;
MikamiUitOpen 2:095b360e0f54 55 return;
MikamiUitOpen 2:095b360e0f54 56 }
MikamiUitOpen 2:095b360e0f54 57
MikamiUitOpen 2:095b360e0f54 58 buffer_[count_] = (this->*Rd[sw_])(); // Read from A0
MikamiUitOpen 2:095b360e0f54 59
MikamiUitOpen 2:095b360e0f54 60 if (++count_ >= N_DATA_)
MikamiUitOpen 7:6598a9b70e5a 61 {
MikamiUitOpen 7:6598a9b70e5a 62 timer_.detach();
MikamiUitOpen 7:6598a9b70e5a 63 Set(true); // Permits spectrum analysis
MikamiUitOpen 7:6598a9b70e5a 64 memcpy(xn_, buffer_, N_DATA_*sizeof(int16_t));
MikamiUitOpen 7:6598a9b70e5a 65 }
MikamiUitOpen 2:095b360e0f54 66 }
MikamiUitOpen 2:095b360e0f54 67 }