Передискретизация звука методом интерполяции, подскажите качественный алгоритм Опции

[soolo]

Стоит казалось бы банальная задача: Сигнал с опорной частотой 44.1 кГц преобразовать в сигнал с частотой 48 кГц. Алгоритмы с FFT по некоторым причинам не подходят, поэтому было решено использовать интерполяцию. Но коэффициент 1.088 и звуки вплоть до всеми любимых 20 кГц пока что ставят нас в тупик. Алгоритмы Lanczos, B-Spline, Mittchel, Sinc результат дают, но где-то гармоники зашкаливают, где-то частоты урезаются.
В первую очередь интересуют минимальные гармонические искажения, задача вторая по степени важности - минимальный срез АЧХ. Может быть есть более подходящие для обработки звука алгоритмы? Подскажите, пожалуйста!

[Fast]

для интереса, т.е. чистоты эксперимента, попробуйте вот этот ФНЧ для 2x up-samp перед Lanczos
в нем 256 отводов, обе части симметричны, потребуется 128 умножений

Код

const int BL = 256;
const int16_T B[256] = {
       11,     22,      9,    -18,    -15,     14,     19,    -12,    -22,
        9,     25,     -7,    -27,      5,     29,     -4,    -31,      2,
       33,      0,    -35,     -1,     37,      3,    -40,     -5,     42,
        6,    -44,     -8,     47,     10,    -49,    -12,     52,     14,
      -54,    -16,     57,     18,    -60,    -20,     63,     22,    -66,
      -25,     69,     27,    -72,    -30,     76,     33,    -80,    -36,
       83,     39,    -87,    -43,     91,     46,    -96,    -50,    100,
       54,   -105,    -58,    110,     63,   -115,    -67,    121,     73,
     -126,    -78,    133,     84,   -139,    -90,    146,     97,   -154,
     -104,    162,    112,   -171,   -121,    180,    130,   -191,   -141,
      202,    152,   -214,   -165,    228,    179,   -243,   -195,    261,
      213,   -280,   -233,    302,    257,   -328,   -285,    359,    318,
     -395,   -358,    440,    407,   -496,   -470,    567,    553,   -664,
     -666,    800,    833,  -1008,  -1101,   1367,   1606,  -2135,  -2919,
     4946,  14721,  14721,   4946,  -2919,  -2135,   1606,   1367,  -1101,
    -1008,    833,    800,   -666,   -664,    553,    567,   -470,   -496,
      407,    440,   -358,   -395,    318,    359,   -285,   -328,    257,
      302,   -233,   -280,    213,    261,   -195,   -243,    179,    228,
     -165,   -214,    152,    202,   -141,   -191,    130,    180,   -121,
     -171,    112,    162,   -104,   -154,     97,    146,    -90,   -139,
       84,    133,    -78,   -126,     73,    121,    -67,   -115,     63,
      110,    -58,   -105,     54,    100,    -50,    -96,     46,     91,
      -43,    -87,     39,     83,    -36,    -80,     33,     76,    -30,
      -72,     27,     69,    -25,    -66,     22,     63,    -20,    -60,
       18,     57,    -16,    -54,     14,     52,    -12,    -49,     10,
       47,     -8,    -44,      6,     42,     -5,    -40,      3,     37,
       -1,    -35,      0,     33,      2,    -31,     -4,     29,      5,
      -27,     -7,     25,      9,    -22,    -12,     19,     14,    -15,
      -18,      9,     22,     11
};

\
и глянем разницу в картинках

[petrov]

для интереса, т.е. чистоты эксперимента, попробуйте вот этот ФНЧ для 2x up-samp перед Lanczos
в нем 256 отводов, обе части симметричны, потребуется 128 умножений

Полифазный надо делать с нечётным количеством коэффициентов, одна фаза получилась бы с одним единственным коэффициентом(ведь каждый второй выходной отсчёт равен входному), все вычисления на исходной частоте.