Рассчитайте амплитуду сигнала после применения БПФ - PullRequest
Новая
       69

Рассчитайте амплитуду сигнала после применения БПФ

0 голосов
/ 16 мая 2018

Я использую Arduino Due и датчик шума, и я применил библиотеку FFT, чтобы извлечь частоту, и она работает отлично.Но я не знаю, как рассчитать амплитуды и распечатать их на консоли?

Вот код: 

#include "arduinoFFT.h"

#define SAMPLES 32             
#define SAMPLING_FREQUENCY 1000 

arduinoFFT FFT = arduinoFFT();

unsigned int sampling_period_us;
unsigned long microseconds;

double vReal[SAMPLES];
double vImag[SAMPLES];

void setup() {
    Serial.begin(9600);

    sampling_period_us = round(1000000*(1.0/SAMPLING_FREQUENCY));
}

void loop() {

    /*SAMPLING*/
    for(int i=0; i<SAMPLES; i++)
    {
        microseconds = micros();    

        vReal[i] = analogRead(0);
        vImag[i] = 0;

        while(micros() < (microseconds + sampling_period_us)){
        }
    }

    /*FFT*/
    FFT.Windowing(vReal, SAMPLES, FFT_WIN_TYP_HAMMING, FFT_FORWARD);
    FFT.Compute(vReal, vImag, SAMPLES, FFT_FORWARD);
    FFT.ComplexToMagnitude(vReal, vImag, SAMPLES);
    double peak = FFT.MajorPeak(vReal, SAMPLES, SAMPLING_FREQUENCY);


    Serial.println(peak);     

    for(int i=0; i<(SAMPLES/2); i++)
    {


        Serial.print((i * 1.0 * SAMPLING_FREQUENCY) / SAMPLES, 1);
        Serial.print(" ");
        //Serial.println(vReal[i], 1);    //View only this line in serial plotter to visualize the bins
    }

    delay(1000);  
    while(1);       
}

1 Ответ

0 голосов
/ 17 мая 2018

Серийный плоттер не имеет номера образца для построения графика.

Если вы хотите правильно просмотреть амплитуды для каждого окна, вам необходимо использовать следующий код 

for(int i=0; i<(SAMPLES/2); i++)
{
//  Serial.print((i * 1.0 * SAMPLING_FREQUENCY) / SAMPLES, 1);
    for(int i=0; i<20; i++)//Adjust the 20 to make space between each samples.
        Serial.print(0);
    Serial.print(vReal[i]); 
}

А затем отрегулируйте размер последовательного плоттера в соответствии с вашим окном.

Надеюсь, это поможет.

...