不韋 呂
スペクトログラム このプログラムの説明は，CQ出版社「トランジスタ技術」の2021年10月号から開始された連載記事「STM32マイコンではじめるPC計測」の中にあります．このプログラムといっしょに使うPC側のプログラムについても同誌を参照してください．
Dependencies: Array_Matrix mbed SerialTxRxIntr DSP_ADDA UIT_FFT_Real Window
Diff: MySpectrogram/FFT_Spectrogram.cpp
- Revision:
- 1:d4e3f39ce206
- Parent:
- 0:3bf11d2ab6ad
--- a/MySpectrogram/FFT_Spectrogram.cpp Thu Sep 09 08:55:42 2021 +0000
+++ b/MySpectrogram/FFT_Spectrogram.cpp Wed Dec 08 03:15:17 2021 +0000
@@ -1,23 +1,23 @@
//-------------------------------------------------------
// スペクトログラムで使う FFT 解析用クラス
//
-// 2021/05/24, Copyright (c) 2021 MIKAMI, Naoki
+// 2021/11/17, 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) {}
+ FftSpectropgram::FftSpectropgram(int nFft)
+ : N_FFT_(nFft), fft_(nFft), wHm_(nFft),
+ xData_(nFft), wData_(nFft), yFft_(nFft/2+1)
+ { SwEmphasis(false); } // 最初は高域強調なし
- void FftSpectropgram::Execute(const Array<float> &xn, Array<float> &absFt)
+ 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];
+ // データのコピー，高域強調の有無は SwEmphasis() で切り替え
+ (this->*fp)(xn);
// 直流分を除去
float sum = 0;
@@ -25,9 +25,18 @@
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 の実行
+ 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]);
}
+
+ // データを作業領域にコピーする際に高域強調処理を行う
+ void FftSpectropgram::CopyH(const Array<float> &xn)
+ {
+ // 差分の処理
+ for (int n=1; n<N_FFT_; n++) // n=1 から開始
+ xData_[n] = xn[n] - xn[n-1];
+ xData_[0] = xData_[1]; // n=0 に対応
+ }
}
\ No newline at end of file