不韋 呂 / Mbed 2 deprecated TrG_FFT_Analyzer

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

Dependencies:   Array_Matrix mbed SerialTxRxIntr DSP_ADDA UIT_FFT_Real Window

MyFFT_Analyzer/FFT_Analyzer.cpp@1:d9dbfbe95c8d, 2021-11-03 (annotated)

Committer:
MikamiUitOpen
Date:
Wed Nov 03 01:29:16 2021 +0000
Revision:
1:d9dbfbe95c8d
Parent:
0:e5fc70976c00 
2

Who changed what in which revision?

UserRevisionLine numberNew contents of line
MikamiUitOpen 0:e5fc70976c00 1 //-------------------------------------------------------
MikamiUitOpen 0:e5fc70976c00 2 // FFT を使ってスペクトル解析を行うクラス
MikamiUitOpen 0:e5fc70976c00 3 //
MikamiUitOpen 0:e5fc70976c00 4 // 2021/05/22, Copyright (c) 2021 MIKAMI, Naoki
MikamiUitOpen 0:e5fc70976c00 5 //-------------------------------------------------------
MikamiUitOpen 0:e5fc70976c00 6
MikamiUitOpen 0:e5fc70976c00 7 #include "FFT_Analyzer.hpp"
MikamiUitOpen 0:e5fc70976c00 8
MikamiUitOpen 0:e5fc70976c00 9 namespace Mikami
MikamiUitOpen 0:e5fc70976c00 10 {
MikamiUitOpen 0:e5fc70976c00 11 FftAnalyzer::FftAnalyzer(int nFft)
MikamiUitOpen 0:e5fc70976c00 12 : N_FFT_(nFft), fft_(nFft), wBlk_(nFft),
MikamiUitOpen 0:e5fc70976c00 13 xData_(nFft), wData_(nFft), yFft_(nFft/2+1) {}
MikamiUitOpen 0:e5fc70976c00 14
MikamiUitOpen 0:e5fc70976c00 15 void FftAnalyzer::Execute(const Array<float> &xn, Array<float> &absFt)
MikamiUitOpen 0:e5fc70976c00 16 {
MikamiUitOpen 0:e5fc70976c00 17 xData_ = xn; // データを作業領域にコピー
MikamiUitOpen 0:e5fc70976c00 18
MikamiUitOpen 0:e5fc70976c00 19 // 直流分を除去
MikamiUitOpen 0:e5fc70976c00 20 float sum = 0;
MikamiUitOpen 0:e5fc70976c00 21 for (int n=0; n<N_FFT_; n++) sum = sum + xData_[n];
MikamiUitOpen 0:e5fc70976c00 22 float ave = sum/N_FFT_;
MikamiUitOpen 0:e5fc70976c00 23 for (int n=0; n<N_FFT_; n++) xData_[n] = xData_[n] - ave;
MikamiUitOpen 0:e5fc70976c00 24
MikamiUitOpen 0:e5fc70976c00 25 wData_ = wBlk_.Execute(xData_); // 窓掛け
MikamiUitOpen 0:e5fc70976c00 26 fft_.Execute(wData_, yFft_); // FFT の実行
MikamiUitOpen 0:e5fc70976c00 27 for (int n=0; n<=N_FFT_/2; n++) // 絶対値に変換
MikamiUitOpen 0:e5fc70976c00 28 absFt[n] = 100.0f*abs(yFft_[n]);
MikamiUitOpen 0:e5fc70976c00 29 }
MikamiUitOpen 0:e5fc70976c00 30 }