float find_max(const uint16_t *buf, size_t len)
float find_min(const uint16_t *buf, size_t len)
float calculate_average(const uint16_t *buf, size_t len)
float calculate_frequency(const uint16_t *buf, size_t len, float sampling_rate)
void generate_ideal_sine(float *ideal_sine, size_t len, float amplitude, float frequency, float sampling_rate, float phase)
float calculate_rms(const uint16_t *buf, size_t len)
float calculate_snr(const uint16_t *adcBuf, const float *ideal_sine, size_t len)

float calculate_rms(volatile float *buf, size_t len) {
    float sum = 0.0f;
    for (size_t i = 0; i < len; i++) {
        sum += buf[i] * buf[i];
    }
    return sqrt(sum / len);
}

float calculate_snr(volatile uint16_t *adcBuf, volatile float *ideal_sine, size_t len) {
    float signal_rms = calculate_rms(ideal_sine, len);
   
    float noise_sum = 0.0f;
    for (size_t i = 0; i < len; i++) {
        float error = (float)adcBuf[i] - ideal_sine[i];
        noise_sum += error * error;
    }
    float noise_rms = sqrt(noise_sum / len);
    return 20 * log10(signal_rms / noise_rms);
}