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

[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
};

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

[soolo]

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

У нас нет up-samp перед Lanczos.
Я лишь тестер. Пойду бить программера, чтоб фильтровал после преобразования

Да уш... ©

Сам себе противоречу? Пойду убьюсь =) После полусотни тестов - в самый раз

Непонятно зачем вообще вопрос задают, сами уже всё знают...

Вопрос задавался, чтобы узнать пути решения проблемы. Форумчане подсказали, мы реализовали, протестировали и выбрали оптимальное для нас решение.
Теперь же вы начинаете критиковать это решение. Я уже пожалел, что написал его

Сообщение отредактировал soolo - Nov 21 2011, 08:41

[Fast]

. Я уже пожалел, что написал его..

совершенно напрасно жалеете, результат он через труд сопли и кровь добывается, а не хвалебные речи и сладкие конфеты.
если мне разработчик говорит, что его уже тошнит от шлифования качества и сотой переделки, я отвечаю, если не вырвало - он еще не достиг в этом достаточного профессионального уровня. сорри за готичность.

..ведь каждый второй выходной отсчёт равен входному.

не могли бы пояснить этот момент, почему каждый второй идет на выход неизменным, мне кажется они не могут быть одинаковыми при любом раскладе

У нас нет up-samp перед Lanczos.
Я лишь тестер. Пойду бить программера, чтоб фильтровал после преобразования

упс, чуть не пропустил. Фильтровать надо ДО преобразования частоты