А какие алгоритмы ресемплинга (8000->44100, 44100->8000) сейчас самые кошерные? По цена/качество?

Я в этом профан, но попробую пальцем в небо ткнуть: линейная, квадратичная, кубическая и прочая интерполяция не достаточно крута?

АПВС?

http://en.wikipedia.org/wiki/Multi-rate_digital_signal_processing

8000

эмм, ты где это откопал?

И раз уж пошла такая пьянка, разве с помощью обратного преобразования фурье можно сделать ресемплирование?

эмм, ты где это откопал?

Ну например в IP телефонии бывает применяется даже частота дискретизации 4000 Гц. Качество конечно дерьмо, несмотря на все заявления, что для голоса этого достаточно.

«Через FFT» - это частый путь в обработке сигналов, сам тащусь, иногда так неожиданно. Например через FFT делают качественную реверберацию, соответствующую заданному помещению по снятому в этом помещении записи хлопка/щелчка.

Ну например записи телефонного канала в этом хранятся. У них же лимит до 3400 гц, поэтому семплировать их на 8000 - самое подходящее их экономических соображений, особенно когда их нужно хранить сотни нефти.

да, пардон, не подумал, что-то я совсем повернулся на музыке.

Например через FFT делают качественную реверберацию, соответствующую заданному помещению по снятому в этом помещении записи хлопка/щелчка.

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

Ну и между делом, я прошёл по приведённой выше ссылке на Вики, увидел, что при ресемплировании используют Фурье, но не для того, чтобы произвести само ресемплирование, а чтобы отфильтровать лишние частоты, которые могу создать нежелательные шумы после изменения частоты дискретизации.

А в итоге всё равно используют интерполяцию.

Ну нормальные простые фильтры легко делаются без FFT, с ним даже хуже получается. Ну например ФНЧ реализуется примитивным алгоритмом с потреблением нескольких байт памяти. Непонятно даже, нафига там FFT, если нужно сделать простой ФНЧ. Про интерполяцию до меня не дошло... как она алгоритмически реализуется, особенно для случаев, когда 2 разные частоты дескретизации не соотносятся в целое число раз. Подъём до каких-нибудь 7 мгц и опускание обратно меня чё-то мало устраивает, нафига столько памяти жрать...

libsamplerate. имеет несколько настроек для изменения параметра цена/качество. скорость зависит от железа. на процессорах без аппаратного FP можно не пытаться.

как она алгоритмически реализуется, особенно для случаев, когда 2 разные частоты дескретизации не соотносятся в целое число раз

для получения коэффициента интерполяции k/n (k,n - целые) делается последовательно интерполяция в k раз + децимация в n раз

Про интерполяцию до меня не дошло... как она алгоритмически реализуется, особенно для случаев, когда 2 разные частоты дескретизации не соотносятся в целое число раз

Ну в случае линейной интерполяции по двум точкам (t1, x1) и (t2, x2), мы можем полагаем, что они связаны линейной функцией x = f(t).

Тогда, чтобы вычислить значение нового (x', t') где t1 <= t' <= t2 мы задаём функцию:

x' = (x2 - x1) / (t2 - t1) * (t' - t1) + x1

для получения коэффициента интерполяции k/n (k,n - целые) делается последовательно интерполяция в k раз + децимация в n раз

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

Ну это приводит к тому, что в какой-то момент ты работаешь с данными в мегагерцах. Нах надо.

Какой-то шум добавляется при каждом проходе в любую сторону. Ошибки округления...

любая обработка в какой-то степени корежит исходную информацию. передискретизация не исключение (по сути это набор фильтров). но именно _шум_ не добавляется.

Можно подробности про реверберацию?

ты уперся в конкретную реализацию. на самом деле все много проще и делается парой простеньких фильтров (CIC или полифазный), буфером FIFO и PLL в случае аппаратной реализации.

это можно делать в целочисленной арифметике

FPU?

ну оно по-разному называется. где-то FPU, где-то VFP, и т.д.

Мне честно пофиг в какой, главное - реализацию в студию и чтобы сильно не пожирало ЦП (-; Или статью о реализации.

http://www.informit.com/articles/article.aspx?p=361985&seqNum=4

здесь вся теория CIC фильтров http://ieeexplore.ieee.org/xpl/login.jsp?tp=&arnumber=1163535&url=htt...

Посмотри, как это сделано в «foobar 2000» и «Album Player (Aplayer)». Это два наиболее ценимых за качество звука аудиоплеера для PC. В сопроводительной документации наверняка этот вопрос раскрыт.

обработка в какой-то степени корежит исходную информацию.

но именно _шум_ не добавляется.

А что же тогда такое «шум», если не искажение исходного сигнала?

это можно делать в целочисленной арифметике

Прошу прощения, но мне не совсем понятно, а как целочисленная арифметика позволяет избежать ошибок округления?

http://ru.wikipedia.org/wiki/Шум

а в данном случае речь идет об искажении АЧХ главным образом. хотя, в некоторых реализациях появятся ошибки округления, они похожи на шум

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

Понятно, спасибо.

http://audiophilesoft.ru/publ/my/foo_resamplers/11-1-0-34

спасибо