Отправлено AntZ 24 июля 2003 г. 06:08

Делаю upsampling с Fs=64000 до Fs=1024000 (шестнадцать раз). Добиваю сампл 15 нулями, получаю каждый шестнадцатый сампл с реальным значением и пятнадцать нулей.

Далее пока делаю FIR на 36 тапов (можно чуть больше, но ненамного - частота высокая). С IIR пока не заморачиваюсь - там то переполнение, то отвратная точность на 16 bit fixed-point. FIR такой:
0-8kHz Passband, 3db Ripple
32-512kHz StopBand 30db Attenuation

Результаты - просто отвратные, фильтр слабый, причем полумерами ничего не добиться. Сглаживает плохо, много импульсных выбросов, особенно на месте "родных" самплов. Глядя на отфильтрованный сигнал хочется уволится с работы и поехать на Чукотку пасти оленей.

Я вот что думаю, а что если сделать интерполяцию в несколько этапов, панпример в два.

1) делаем upsampling c Fs=64000 до Fs=256000 (четыре раза),
фильтр
Passband 0-8Khz
Stopband 32-128Khz
Фильтр не сложный, на и частота умеренная

2) upsampling Fs eще в четыре раза, c 256000 до 1024000
фильтр
Passband 0-8kHz
Stopband 128kHz-512kHz (так как 8-128 уже "чистый" диапазон)
Фильтр тоже должен быть не сложным, да и сигнал уже слегка "подготовленный"

Может быть есть еще какие способы сделать upsampling в 16 раз, при очень ограниченных ресурсах? Кто делал подобные вещи?

• Ответ: немного про полифазные фильтры (-) — TAP   (24.07.2003 12:12, пустое, ссылка)
  • Ответ: И все таки полифазным быстрее :) (+) — TAP   (24.07.2003 17:17, 155 байт)
• А почему не — Fat Robot   (24.07.2003 11:39, 97 байт)
  • А потому, — AntZ   (24.07.2003 12:41, 204 байт)
    • Линейно можно — Fat Robot   (24.07.2003 13:21, 17 байт)
      • Cпектр будет - срашно представить, сплошные изломы — AntZ   (25.07.2003 07:45, пустое)
      • Да и кубически тоже можно — shenon   (24.07.2003 14:51, 164 байт, ссылка)
      • Или CIC фильтр — Fat Robot   (24.07.2003 13:25, 114 байт)
        • После добивания нулями, конечно — Fat Robot   (24.07.2003 13:27, пустое)
• Ответ: Стандартный прием - down/upsampling 2... — PGN   (24.07.2003 10:50, 204 байт)
• Делал на 54xx - 8 кГц - 1.536 МГц ... — проходивший   (24.07.2003 10:31, 581 байт)
  • Вы случайно не модуляцию/дкмодуляцию делали? Narrowband передача на высоких частотах? — AntZ   (24.07.2003 11:02, пустое)
• Ответ: — AntZ   (24.07.2003 06:28, 172 байт)
  • Давным давно я делал полифазу для звука на ADSP-2181и усвоил одну простую истину - на точности 16 бит не получицца... Пришлось смириться и переделывать под 32, но по ресурсам уже не лезло.... вот такая беда со мной случилась :) — Зок Мёдов   (24.07.2003 13:12, пустое)
  • Ответ: делал я с 8 до 400 кГц Получилось нормально Нада было 40 дб и Получил стока (-) — TAP   (24.07.2003 07:37, пустое)
  • Ответ: полифаза не сложнее Надо просто сделать многокаскадный полифазный фильтр (-) — TAP   (24.07.2003 07:34, 212 байт)
    • Времени мало, и памяти тоже. Сделал двухэтапный upsampling - с двумя фильтрами. Первый FIR 24 на Fs=512kHz, второй FIR 16 на 1.024kHz, результаты хорошие, особенно по памяти - все в одном буффере — AntZ   (24.07.2003 08:36, 41 байт)
      • Ответ: Поэкспериментируйте с одинаковыми фильтрами на каждой ступени - сэкономите также на программной памяти... — PGN   (24.07.2003 11:01, 154 байт)
        • В конце концов так и сделал - однинаковый фильтр на всех ступенях — AntZ   (26.07.2003 05:13, пустое)
        • Спасибо, может так и сделаю, пока нашел временную затычку, за два прохода делаю — AntZ   (24.07.2003 11:18, пустое)
      • Ошибся, первый Fs=256, не 512, upsampling два раза в четыре раза. Получилось просто, дешево и сердито. Полифаза, наверное, лучше, но надо дольше возиться — AntZ   (24.07.2003 08:39, пустое)