Arduinoos

Part 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11

So far so good, we got our magnitude spectrum. We may now analyze it just like any other spectrum in order to get meaningful information. In our case, we may want to identify the main frequencies from the spectrum. Firstly we will run a peak picking algorithm in order to locate major peaks, and then keep the most intense one. Those who want to perform spectral comparison will have to keep all significant peaks.

There are lots of peak picking algorithm! But the one we need has to be small and fast. This one is pretty trivial but it works great.

Secondly, because of the poor peak shapes (lack of large vectors of data) we need to interpolate the peak apex in order to get an accurate frequency. I am using a non iterative quadratic interpolation which gives good results too.

Both routines are merged in one function

double majorPeak(double *vD, uint8_t samples) {
	double maxY = 0;
	uint8_t IndexOfMaxY = 0;
	for (uint8_t i = 1; i < (samples - 1); i++) {
		if ((vD[i-1] < vD[i]) && (vD[i] > vD[i+1])) {
			if (vD[i] > maxY) {
				maxY = vD[i];
				IndexOfMaxY = i;
			}
		}
	}
	double delta = 0.5 * ((vD[IndexOfMaxY-1] - vD[IndexOfMaxY+1]) / (vD[IndexOfMaxY-1] - (2.0 * vD[IndexOfMaxY]) + vD[IndexOfMaxY+1]));
	double interpolatedX = (IndexOfMaxY + delta) / samplingDuration;
	return(interpolatedX);
}

Links:

• Quadratic Interpolation of Spectral Peaks

Next post on same subject