应用FFT后计算信号的幅度

时间:2018-05-16 17:03:04

标签: arduino fft arduino-due

我正在使用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 :(得分:0)

串行绘图仪没有多个样本来绘制功能。

如果要正确查看每个窗口的幅度,则需要使用以下代码

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]); 
}

然后调整串行绘图仪的大小以适合您的窗口。

我希望这会有所帮助。