Одному чрезвычано круто распальцованому специалисту в области интерполяции требуется совет, даваемый как бы между прочим, чтобы не сильно уходшить параметры распальцовки, выше означенного спеца. («Телесистемы»: Конференция «Цифровые сигнальные процессоры (DSP) и их применение»)

[an error occurred while processing this directive] Одному чрезвычано круто распальцованому специалисту в области интерполяции требуется совет, даваемый как бы между прочим, чтобы не сильно уходшить параметры распальцовки, выше означенного спеца.

(«Телесистемы»: Конференция «Цифровые сигнальные процессоры (DSP) и их применение»)

Отправлено st256 23 июля 2004 г. 09:10

Специалист это - я, а вопрос в следующем:

У меня на входе частота дискретизации 44.1 кГц (аудио сигнал). Надо из нее получить 96кГц. Вопрос: как это сделать по-приличнее? Выслушаю любые советы.

Пока у нас применяют такой алгоритм
1. Из 44.1 кГц получают при помощи FIR (99-порядок) 88.2 кГц.
2. Из 88.2 кГц получают при помощи кубической интерполяции уже 96кГц.
Погрешность интерполяции в этом случае менее -90дБ.

Но, как было заявлено ранее, я имею некий опыт в такого рода вещах. Мне кажется, что для -90 дБ описываемая схема избыточна по требуемой производительности процессора. Вообщем, кто чем может, поможИте :)

Составить ответ ||| Конференция ||| Архив

Ответы

В тему: может у кого завалялась статья? interpolation in Digital Modem - Part II, IEEE Trans. on Comms. June 1993. — alostap (23.07.2004 15:54, пустое)
У Аналога есть звуковой кодек AD1843. Так так помимо всего прочего ресамплер имеется. Как работает подробно не написано похоже что-то вроде связки сигма - дельта ЦАП сигма-дельта АЦП(Но все в цифре разумеется никакого аналога). Может кто-нибудь слышал о такой технологии? — alostap (23.07.2004 15:42, пустое)
Я вот закачал Матлаб модель для расчета фильтра-компенсатора для AD73322 на основе полифазных фильтров. Прошу ногами сильно не бить — dspman (23.07.2004 13:38, пустое, ссылка)
Ответ: А может есть у кого нибудь оригинальная статья Хогенауэра про CIC? — ux (23.07.2004 13:01, пустое)
- Ответ: Ссылка — PGN (23.07.2004 16:13, 242 байт)
  - Ответ: Присылайте (+) — ux (23.07.2004 17:02, 29 байт)
- нету, вообще Вам ехать или шашечки? — fontp (23.07.2004 14:03, 855 байт)
  - Так интегратор по жизни частенко переполняется , но это никому не мешает — alostap (23.07.2004 15:10, пустое)
    - Ответ: это да, но в данном случае результат абсолютно точный — fontp (23.07.2004 15:28, 144 байт)
  - Ответ: Да знаю я всё это! Просто первоисточник интересен, на него все ссылаются, а его нет в инете. — ux (23.07.2004 14:51, пустое)
    - Ответ: таки шашечки ;) — fontp (23.07.2004 14:54, пустое)
      - А вот мат. выражение для импульсной хар-ки фильтра компенсирующего неравномерность CIC? — ux   (23.07.2004 15:07, пустое)
        
        Дело давнее, — -=ВН=-   (23.07.2004 16:54, 490 байт)
        
        это да, типа Ремеза всегда — fontp   (23.07.2004 17:07, пустое)
        
        Поправка — -=ВН=-   (23.07.2004 17:04, 47 байт)
        
        Ответ: xз, в прoграммках синтеза фильтров, обычно, есть опция компенсировать x/sinx — fontp   (23.07.2004 15:20, 228 байт)
        
        Косрут(ну и выражение) сам считаю, а проги какие? — ux   (23.07.2004 15:25, пустое)
        
        ну разные — fontp   (23.07.2004 15:43, 509 байт)
        
        Ответ: Да и не то это(+) — ux   (23.07.2004 15:31, 296 байт)
        
        Ответ: как когда, когда то, когда не то. Это уже зависит на какую частоту вывести выход CIC — fontp   (23.07.2004 15:49, пустое)
Ответ: И я вякну... — PGN (23.07.2004 12:49, 856 байт)
- Дело в том, что я не помню, на что и что именно я пытался ответить. — st256 (23.07.2004 15:05, 84 байт)
Опоздал... — -=ВН=- (23.07.2004 12:05, 302 байт)
Делал похожее = — alostap (23.07.2004 11:09, 327 байт)
- Добивали нулями? — st256 (23.07.2004 11:47, пустое)
  - Ну не совсем добивал. Просто у меня было 8 векторов коэффициентов - каждый из которых представлял децимированный вариант ИХ КИХ .. И далее линию задержки с входными отсчетами я умножал скалярно на соотв. вектор. — alostap (23.07.2004 14:48, пустое)
    - Конечно эквивалетно добиванию нулями. — alostap (23.07.2004 14:49, пустое)
А если полифазный фильтр с передискретизацией 320/147? 44.1*320/147 аккурат 96 дает. — dspman (23.07.2004 10:48, пустое)
- Хоть и считать там немного на каждую выборку, но коэффициентов-то будет бооольшое количество... — SM (23.07.2004 12:03, пустое)
  - Это точно. Но если иметь быстрый способ вычиления на лету sin(x)/x, то можно память на мипсы поменять. Или часть коэффициентов положть равной нулю ( может тогда вообще без фильтра обойдемся:)) — dspman (23.07.2004 12:12, пустое)
    - И кроме sin(x)/x окна еще хорошего такого, чтоб на -90 дб :)) — SM (23.07.2004 12:18, пустое)
- От 15 миллионов самплов в секунду DSP может копыта отбросить — AntZ (23.07.2004 11:47, пустое)
  - Так не всеже 320 выходных сэмплов на один входной нужно рассчитывать! Ведь 146 из 147 дальше не нужны! — alostap (23.07.2004 15:02, пустое)
    - Оставшихся тоже хватит чтобы поиметь Кама-Сутру по полной программе — AntZ (23.07.2004 15:27, пустое)
  - Ответ: — st256 (23.07.2004 11:48, пустое, ссылка)
- А вот это идея! Спасибо. Подумаю обязательно. — st256 (23.07.2004 11:38, 1 байт)
  - Реализуется двумя FIR фильтрами - один на входной частоте работает, второй на выходной, вроде не много ресурсов должен взять, Выбором оконной функции для импульсных характеристик фильтра можно регулировать уровень боковых лепестков АЧХ. — dspman (23.07.2004 11:49, пустое)
    - А Вы реально что-то подобное делали? Я делал ресамплинг 9600->8000 через полифазу, но там случай намного проще - upsampling 5/downsampling 6. Именно поэтому сомнения и закрались — AntZ (23.07.2004 13:06, пустое)
      - Ага, и такой тоже делал. Я покапаюсь в матлабовских модельках, может чего найду для общественности. — dspman (23.07.2004 13:24, пустое)
    - А как по-английски "полифазный фильтр"? — st256 (23.07.2004 11:52, пустое)
      - Polyphase. — dspman   (23.07.2004 11:59, пустое)
        
        Спасибо! — st256   (23.07.2004 12:08, пустое)
        
        есть ещё чуднАя программка MultiRateFilterDesign — fontp   (23.07.2004 13:15, 343 байт)
        
        Есть такая чудная книга "Digital Signal Processing" Монолакис, Проакис. Там это все хорошо описано. — dspman   (23.07.2004 12:16, пустое)
        
        Ответ: А где взять можно? — ux   (23.07.2004 13:11, пустое)
        
        Сам бы с радостью где-нибудь взял. Я работал с бумажной копией. — dspman   (23.07.2004 13:17, пустое)
Способ №7: Посмотреть описание специализированных микросхем ресемплинга - могут выясниться интересные мелочи. — АПРК (23.07.2004 10:24, пустое)
- :) так я Вам такую микросхему и описал :)))) — st256 (23.07.2004 11:39, пустое)
В аналогичном случае применял почти то-же самое (+) — SM (23.07.2004 09:43, 225 байт)
- Мне видится так:(+) — ARK (23.07.2004 10:15, 639 байт)
  - Ответ: (+) — SM (23.07.2004 10:19, 145 байт)
- Ответ: — st256 (23.07.2004 09:56, 281 байт)
  - Ответ: (+) — SM (23.07.2004 10:07, 269 байт)
    - Вдогонку (+) — SM (23.07.2004 10:21, 197 байт)

Перейти к списку ответов ||| Конференция ||| Архив ||| Главная страница ||| Содержание ||| Без кадра

E-mail: info@telesys.ru

Ответы

Отправка ответа