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

Revision:
4:99d4d5ea06a2
Parent:
3:7c26b701f363
--- a/MyClasses/CepstrumAnalysis.cpp	Mon Dec 14 13:44:42 2015 +0000
+++ b/MyClasses/CepstrumAnalysis.cpp	Tue Dec 15 02:39:31 2015 +0000
@@ -1,7 +1,7 @@
 //-------------------------------------------------------
 //  Class for spectrum analysis using cepstrum
 //
-//  2015/12/14, Copyright (c) 2015 MIKAMI, Naoki
+//  2015/12/15, Copyright (c) 2015 MIKAMI, Naoki
 //-------------------------------------------------------
 
 #include "CepstrumAnalysis.hpp"
@@ -10,9 +10,8 @@
 {
     CepstrumAnalyzer::CepstrumAnalyzer(int nData, int nFft, int nLifter)
             : AnalyzerBase(nData, nFft, nFft), N_LIFTER_(nLifter),
-              xFft_(new float[nFft]), yFft_(new Complex[nFft/2+1]),
-              lifter_(new float[nLifter]), db_(new Complex[nFft/2+1]),
-              cep_(new float[nFft]),
+              yFft_(new Complex[nFft/2+1]), lifter_(new float[nLifter]),
+              db_(new Complex[nFft/2+1]), cep_(new float[nFft]),
               cepLft_(new float[nFft]), cepFt_(new Complex[nFft/2+1])
     {
         // Generate lifter
@@ -23,7 +22,6 @@
 
     CepstrumAnalyzer::~CepstrumAnalyzer()
     {
-        delete[] xFft_;
         delete[] yFft_;
         delete[] lifter_;
         delete[] db_;
@@ -34,8 +32,7 @@
 
     void CepstrumAnalyzer::Analyze(const float xn[], float yn[])
     {
-        wHm_.Execute(xn, xFft_);    // Windowing and zero-padding
-        fft_.Execute(xFft_, yFft_); // Execute FFT
+        fft_.Execute(xn, yFft_);    // Execute FFT
 
         // Translate to dB and to complex
         for (int n=0; n<=N_FFT_/2; n++)
@@ -56,4 +53,3 @@
             yn[n] = cepFt_[n].real();
     }
 }
-