RNNoise子带能量的理解

void compute_band_energy(float* bandE, const kiss_fft_cpx* X) {
    int i;
    float sum[NB_BANDS] = { 0 };/*NB_BANDS是22个端点，即21个子带*/
    for (i = 0; i < NB_BANDS - 1; i++)
    {
        int j;
        int band_size;  /*对应子带的频点数*/
        band_size = (eband5ms[i + 1] - eband5ms[i]) << FRAME_SIZE_SHIFT;
        for (j = 0; j < band_size; j++) {
            float tmp;
            float frac = (float)j / band_size;  /*理解为窗加权，顶点为端点，覆盖左右两子带*/
            tmp = SQUARE(X[(eband5ms[i] << FRAME_SIZE_SHIFT) + j].r);/*计算子带每个频点的能量，起点为该端点的对应频点  实部*/
            tmp += SQUARE(X[(eband5ms[i] << FRAME_SIZE_SHIFT) + j].i);/*虚部*/
            sum[i] += (1 - frac) * tmp;/*此时sum[i]为上一个循环的sum[i+1](递增窗)，加上右边子带(递减窗)*/
            sum[i + 1] += frac * tmp;/*为下个循环提供左边子带的能量*/
        }
    }
    sum[0] *= 2;/*两头的端点有一半是0，所以单个子带加窗乘2*/
    sum[NB_BANDS - 1] *= 2;
    for (i = 0; i < NB_BANDS; i++)
    {
        bandE[i] = sum[i];
    }
}