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

[soolo]

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

[Fast]

а зачем здесь к полифазному фильтру еще и лагранж(или фарроу) ?
тогда уж upsampling в 2 раза с помощью ФНЧ КИХ-фильтра, а там уже и Лагранж 6-8 узлов (Lanczos, B-Spline) для понижения 88.2->48

[petrov]

а зачем здесь к полифазному фильтру еще и лагранж(или фарроу) ?
тогда уж upsampling в 2 раза с помощью ФНЧ КИХ-фильтра, а там уже и Лагранж 6-8 узлов (Lanczos, B-Spline) для понижения 88.2->48

Хотите умножайте нули после upsampling, обычно это полифазной структурой КИХ фильтра делают.

Лишнего преобразования частоты тоже хотелось бы избежать.

Ничего лишнего и нету с помощью полифазного фильтра промежуточные отсчёты вычисляются для нормальной работы дробного интерполятора Лагранжа(Фарроу).

[Fast]

Хотите умножайте нули после upsampling, обычно это полифазной структурой КИХ фильтра делают.

зачем же нули умножать в 21м веке..
смысл исключительно расширении полосы, чтобы не "накрыть" полезные составляющие спектра сигнала неравномерной частью АЧХ фильтра Лагранжа
в этом случае Лагранж можно успешно применять как на up- так и на down- sampling на любую дробную частоту
полифазные фильтры же мне кажутся громоздкими, но правда, не знаю , насколько они удобны в FPGA-реализации

[petrov]

зачем же нули умножать в 21м веке..
смысл исключительно расширении полосы, чтобы не "накрыть" полезные составляющие спектра сигнала неравномерной частью АЧХ фильтра Лагранжа
в этом случае Лагранж можно успешно применять как на up- так и на down- sampling на любую дробную частоту
полифазные фильтры же мне кажутся громоздкими, но правда, не знаю , насколько они удобны в FPGA-реализации

Вы понимаете что

тогда уж upsampling в 2 раза с помощью ФНЧ КИХ-фильтра

это и есть полифазный фильтр, если нули не умножать, а поменять местами upsampling и КИХ ФНЧ?.

[Fast]

Вы понимаете что это и есть полифазный фильтр, если нули не умножать, а поменять местами upsampling и КИХ ФНЧ?.

да, при упрощении. Но up-samp в 2 раза и re-samp на произвольную частоту с помощью полифазного фильтра все же разные вещи.

[petrov]

да, при упрощении.

Какое нафиг упрощение? О чём вы? Это есть натуральный полифазный фильтр.

up-samp в 2 раза

И не обязательно в 2, в зависимости от необходимого качества может и больше потребоваться.

re-samp на произвольную частоту с помощью полифазного фильтра все же разные вещи.

Про это речи вообще не было.

[Fast]

Какое нафиг упрощение? О чём вы? Это есть натуральный полифазный фильтр.

такой нафиг, что полифазный - это форма реализации. Полифазным может быть как КИХ, так и БИХ, к примеру, и КИХ с невысчитываемыми нулями не становится автоматом полифазным, также как не становится автоматом лагранж интерполятором фарроу, хотя фарроу одна из форм реализации лагранжа. А становится при определенных условиях реализации.

И не обязательно в 2, в зависимости от необходимого качества может и больше потребоваться.

и каким же образом up-samp на полифазнике в N-раз (+Фарроу) нам поднимет качество больше чем up-samp на КИХ-фильтре в 2 раза (+Фарроу)?

p.s. Могу представить только один вариант, если это Фарроу 3го порядка, с пологой АЧХ, для которого надо стабильно повышать в 4, а лучше в 8 раз сигнал. При Лагранже 8-порядка повышай хоть в 4, 8, 128 - качеству взяться неоткуда..

При насыщенном почти равномерном спектре Фарроу 6 и 10 порядка будут еще хуже. А метод идеальный по определению. Можно ли его выполнить побыстрее для большего числа точек не знаю.

при насыщенном и равномерном спектре до частоты Найквиста ничто не поможет , любая интерполяция будет с искажениями. Потому что цифровать надо было с бОльшей частотой. Никакой преселектор, ни какое АЦП с дельта-сигма-модуляцией не имеет прямоугольной характеристики фильтра НЧ. И если мы видим сильный сигнал вблизи найквиста, сто пудово часть его уже "завернулась" в полосу полезного сигнала при оцифровке