不韋 呂 / Mbed 2 deprecated TrG_Spectrogram

スペクトログラム このプログラムの説明は，CQ出版社「トランジスタ技術」の2021年10月号から開始された連載記事「STM32マイコンではじめるPC計測」の中にあります．このプログラムといっしょに使うPC側のプログラムについても同誌を参照してください．

Dependencies:   Array_Matrix mbed SerialTxRxIntr DSP_ADDA UIT_FFT_Real Window

MySpectrogram/FFT_Spectrogram.cpp

Committer:
MikamiUitOpen
Date:
2021-09-09
Revision:
0:3bf11d2ab6ad
Child:
1:d4e3f39ce206

File content as of revision 0:3bf11d2ab6ad:

//-------------------------------------------------------
//  スペクトログラムで使う FFT 解析用クラス
//
//  2021/05/24, Copyright (c) 2021 MIKAMI, Naoki
//-------------------------------------------------------

#include "FFT_Spectrogram.hpp"

namespace Mikami
{
    FftSpectropgram::FftSpectropgram(int nData, int nFft)
        : N_DATA_(nData), N_FFT_(nFft),
          fft_(nFft), wHm_(nFft, nData-1), b1_(1.0f),
          xData_(nFft), wData_(nFft), yFft_(nFft/2+1) {}

    void FftSpectropgram::Execute(const Array<float> &xn, Array<float> &absFt)
    {
        // 高域強調
        for (int n=0; n<N_DATA_-1; n++)
            xData_[n] = xn[n+1] - b1_*xn[n];

        // 直流分を除去
        float sum = 0;
        for (int n=0; n<N_FFT_; n++) sum = sum + xData_[n];
        float ave = sum/N_FFT_;
        for (int n=0; n<N_FFT_; n++) xData_[n] = xData_[n] - ave;

        wData_ = wHm_.Execute(xData_);	// 窓掛け
        fft_.Execute(wData_, yFft_);	// FFT の実行
        for (int n=0; n<=N_FFT_/2; n++) // 絶対値に変換
            absFt[n] = 100.0f*abs(yFft_[n]);
    }
}