Important changes to repositories hosted on mbed.com
Mbed hosted mercurial repositories are deprecated and are due to be permanently deleted in July 2026.
To keep a copy of this software download the repository Zip archive or clone locally using Mercurial.
It is also possible to export all your personal repositories from the account settings page.
Dependencies: BSP_DISCO_F746NG BUTTON_GROUP LCD_DISCO_F746NG TS_DISCO_F746NG UIT_FFT_Real mbed
MyClasses/LPC_Analysis.cpp@2:095b360e0f54, 2015-12-13 (annotated)
- Committer:
- MikamiUitOpen
- Date:
- Sun Dec 13 07:20:14 2015 +0000
- Revision:
- 2:095b360e0f54
- Parent:
- 0:5c237fdcba23
- Child:
- 3:7c26b701f363
3
Who changed what in which revision?
User | Revision | Line number | New contents of line |
---|---|---|---|
MikamiUitOpen | 0:5c237fdcba23 | 1 | //------------------------------------------------------- |
MikamiUitOpen | 0:5c237fdcba23 | 2 | // Class for spectrum analysis using linear prediction |
MikamiUitOpen | 0:5c237fdcba23 | 3 | // |
MikamiUitOpen | 0:5c237fdcba23 | 4 | // 2015/12/08, Copyright (c) 2015 MIKAMI, Naoki |
MikamiUitOpen | 0:5c237fdcba23 | 5 | //------------------------------------------------------- |
MikamiUitOpen | 0:5c237fdcba23 | 6 | |
MikamiUitOpen | 0:5c237fdcba23 | 7 | #include "LPC_Analysis.hpp" |
MikamiUitOpen | 0:5c237fdcba23 | 8 | |
MikamiUitOpen | 0:5c237fdcba23 | 9 | namespace Mikami |
MikamiUitOpen | 0:5c237fdcba23 | 10 | { |
MikamiUitOpen | 0:5c237fdcba23 | 11 | LpcAnalyzer::LpcAnalyzer(int nData, int nFft, int order) |
MikamiUitOpen | 0:5c237fdcba23 | 12 | : AnalyzerBase(nData, nFft), ORDER_(order), |
MikamiUitOpen | 0:5c237fdcba23 | 13 | wHm_(nData-1, nData-1), lp_(nData-1, order), |
MikamiUitOpen | 0:5c237fdcba23 | 14 | xData_(new float[nData]), an_(new float[order]), |
MikamiUitOpen | 0:5c237fdcba23 | 15 | xFft_(new float[nFft]), yFft_(new Complex[nFft/2+1]) {} |
MikamiUitOpen | 0:5c237fdcba23 | 16 | |
MikamiUitOpen | 0:5c237fdcba23 | 17 | LpcAnalyzer::~LpcAnalyzer() |
MikamiUitOpen | 0:5c237fdcba23 | 18 | { |
MikamiUitOpen | 0:5c237fdcba23 | 19 | delete[] xData_; |
MikamiUitOpen | 0:5c237fdcba23 | 20 | delete[] an_; |
MikamiUitOpen | 0:5c237fdcba23 | 21 | delete[] xFft_; |
MikamiUitOpen | 0:5c237fdcba23 | 22 | delete[] yFft_; |
MikamiUitOpen | 0:5c237fdcba23 | 23 | } |
MikamiUitOpen | 0:5c237fdcba23 | 24 | |
MikamiUitOpen | 0:5c237fdcba23 | 25 | void LpcAnalyzer::Analyze(const float xn[], float yn[]) |
MikamiUitOpen | 0:5c237fdcba23 | 26 | { |
MikamiUitOpen | 0:5c237fdcba23 | 27 | wHm_.Execute(xn, xData_); // Windowing |
MikamiUitOpen | 0:5c237fdcba23 | 28 | float em; |
MikamiUitOpen | 0:5c237fdcba23 | 29 | lp_.Execute(xData_, an_, em); |
MikamiUitOpen | 0:5c237fdcba23 | 30 | |
MikamiUitOpen | 0:5c237fdcba23 | 31 | // To spectrum |
MikamiUitOpen | 0:5c237fdcba23 | 32 | xFft_[0] = 1.0f; |
MikamiUitOpen | 0:5c237fdcba23 | 33 | for (int n=0; n<ORDER_; n++) xFft_[n+1] = -an_[n]; |
MikamiUitOpen | 0:5c237fdcba23 | 34 | for (int n=ORDER_+1; n<N_FFT_; n++) xFft_[n] = 0.0f; |
MikamiUitOpen | 0:5c237fdcba23 | 35 | fft_.Execute(xFft_, yFft_); // execute FFT |
MikamiUitOpen | 0:5c237fdcba23 | 36 | |
MikamiUitOpen | 0:5c237fdcba23 | 37 | // Translate to dB |
MikamiUitOpen | 0:5c237fdcba23 | 38 | for (int n=0; n<=N_FFT_/2; n++) |
MikamiUitOpen | 0:5c237fdcba23 | 39 | yn[n] = 10.0f*log10f(em/Norm(yFft_[n])); |
MikamiUitOpen | 0:5c237fdcba23 | 40 | } |
MikamiUitOpen | 0:5c237fdcba23 | 41 | } |
MikamiUitOpen | 2:095b360e0f54 | 42 |