不韋 呂
スペクトログラム このプログラムの説明は，CQ出版社「トランジスタ技術」の2021年10月号から開始された連載記事「STM32マイコンではじめるPC計測」の中にあります．このプログラムといっしょに使うPC側のプログラムについても同誌を参照してください．
Dependencies: Array_Matrix mbed SerialTxRxIntr DSP_ADDA UIT_FFT_Real Window
Diff: MySpectrogram/FFT_Spectrogram.cpp
- Revision:
- 0:3bf11d2ab6ad
- Child:
- 1:d4e3f39ce206
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/MySpectrogram/FFT_Spectrogram.cpp Thu Sep 09 08:55:42 2021 +0000
@@ -0,0 +1,33 @@
+//-------------------------------------------------------
+// スペクトログラムで使う 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]);
+ }
+}
\ No newline at end of file