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

[an error occurred while processing this directive] Помогите советом... /+/

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

Отправлено Николай 02 сентября 2005 г. 09:31

Необходим цифровой фильтр для выделения сигнала в полосе частот 0 - 55 Гц. Частота дискретизации 6400 Гц. Погрешность коэффициента передачи в полосе пропускания не более 1 %, подавление на частоте 100 Гц - не менее 60 дБ (желательно 80 дБ). Разрядность 16 бит. Фильтр КИХ получается очень длинным и медленным по динамике. Может, использовать несколько каскадов БИХ с прореживанием данных после каждого каскада? Как я понимаю, проблема в том, что ширина переходной полосы очень узкая.
Заранее спасибо.

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

Ответы

Ответ: Многоступенчатый бих с прореживанием - идея хорошая.+ — andy_P (05.09.2005 14:24 80.82.63.185, 431 байт)
- Ответ: Минус один - каждая децимация будет мало-мало алиасить в полосу. Те чем меньше стадий децимации, тем лучше — andy_P (05.09.2005 16:54 80.82.63.185, пустое)
  - В принципе, да. Но алиазинга можно избежать, если сделать децимацию только после первого биквада, а после второго и третьего оставить все как есть. Эти второй и третий биквады весь алиазинг предыдущего шага и "сожрут". Возможны и иные схемы, когда децимация делается не по "полной программе" — homekvn (05.09.2005 17:06 212.185.161.237, пустое)
    - Да вроде не сожрут что в полосу 0..55 Hz попало? — andy_P (05.09.2005 17:42 80.82.63.185, пустое)
      - Так в полосу 0..55 и не попадут, поскольку я предлагаю децимацию сделать только после первого каскада и, скажем, не в 58 раз, как это в принципе можно было бы сделать, доведя частоту дискретизации до 110 Гц, а только в 10-12 раз. Тогда алиазинг в полосе 0..55 будет ничтожно малым, поскольку за первым каскадом будут следовать два каскада без децимации вообще. И вот эти два каскада и "сожрут" алиазинг от первого каскада почти полностью. — homekvn (06.09.2005 11:07 212.185.161.237, пустое)
- Можно и ких. Только все-таки MIPS-ов и памяти он больше скушает (+) — homekvn (05.09.2005 16:10 212.185.161.237, 761 байт)
  - Ответ: дак там полоса 55 hz на частоте дискретизации > 6kHz. Делаем срез первой ступени 100. Имеем ~4 октавы или 24 dB с секции (если децимировать на 2 первый раз). Те 2 секций по моему тупому разумению на первой ступени должно хватить за глаза. Наверное, можно уложиться и в одну секцию, если как-о извратиться — andy_P (05.09.2005 16:42 80.82.63.185, пустое)
    - Да КИХ без прореживания через БПФ на такой-то частоте на калькуляторе можно успеть посчитать... Ловкость рук, так сказать... — -=ВН=- (05.09.2005 17:09 194.190.181.231, пустое)
      - Ни-ни (+) — homekvn   (05.09.2005 17:43 212.185.161.237, 1335 байт)
        
        Ну считайте. — -=ВН=-   (05.09.2005 18:03 194.190.181.231, 561 байт)
        
        Согласен, что упорядочивать после Фурье не нужно. Но на мой взгляд задача не стоит того, чтобы ее через Фурье решать. Кстати, Фурье, на мой взгляд, надо будет 32-разрядное взять. Иначе после FFT и IFFT шумы округления будут достаточно большими. — homekvn   (06.09.2005 10:59 212.185.161.237, пустое)
        
        Для 16-ти разрядного проц. его с бл. пл. запятой можно сделать, 16-ти разрядным. — -=ВН=-   (06.09.2005 12:04 194.190.181.231, пустое)
  - Про память да, памяти КИХ сожрет больше. — -=ВН=- (05.09.2005 16:22 194.190.181.231, 39 байт)
    - Ну, согласен. MIPSы с учетом прореживания, а также того, что при вычислении КИХ достаточно 16-бит, будут такими, что то на то и выйдет. Но задержка... Сам ких имеет неслабую задержку (около 300 семплов на тактовой частоте), да еще и многокаскадность внесет свой вклад. А если еще после этого апсемплить надо?.. — homekvn (05.09.2005 16:45 212.185.161.237, пустое)
      - Какую задержку внесет многокаскадность? Есть фильтр с какой-то пер. ф-ей. И все. Никого не интересует как он реализован. Многкаскадно или немногокаскадно. — -=ВН=-   (05.09.2005 17:05 194.190.181.231, 137 байт)
        
        Не совсем так. Это несложно показать (+) — homekvn   (05.09.2005 17:25 212.185.161.237, 1169 байт)
        
        Ну и посчитайте теперь, какой фильтр, но однокаскадный, будет полностью эквивалентен Вашим 2 каскадам. Будет ли его длина 59 отчетов. — -=ВН=-   (05.09.2005 17:35 194.190.181.231, пустое)
        
        В том-то все и дело, что прежде, чем писать, посчитал. Все так, как написал, и есть. Иначе бы и писать не стал (+) — homekvn   (05.09.2005 17:49 212.185.161.237, 394 байт)
        
        Да что Вы посчитали, если эквив. длина Ваших 2 секций=92 отсчетам. В пересчете на исх. частоту дискр. — -=ВН=-   (05.09.2005 17:53 194.190.181.231, пустое)
        
        Да, Ваша правда. Прошу прощения за бузеш. Если сделать фильтр с таким же подавлением, то получается примерно одинаковая задержка (+-1 семпл). — homekvn   (06.09.2005 10:52 212.185.161.237, пустое)
        
        Прошу прощения. Завтра отвечу. Сейчас ухожу. — homekvn   (05.09.2005 18:05 212.185.161.237, пустое)
  - Сорри (я подумал, что ких, а написано бих). Беру пред. слова назад — homekvn (05.09.2005 16:12 212.185.161.237, пустое)
Достаточно трех биквадов для подавления 80 дБ — homekvn (02.09.2005 16:09 212.185.161.237, 980 байт)
- Ответ: Может, конечно, я и ошибся где-то, но у меня фильтр с коэффициентами, которые Вы привели, режет уже 35 Гц на -3 дБ. /-/ — Николай (05.09.2005 16:56 80.73.170.86, пустое)
  - Послал на Ваш е-мейл программу рассчета и маленький комментарий. Будут вопросы - пишите. — homekvn (05.09.2005 18:01 212.185.161.237, пустое)
- Ответ: Спасибо. — Николай (02.09.2005 18:08 80.73.170.86, пустое)
  - Маленькое уточнение — homekvn (02.09.2005 20:10 84.146.6.148, 973 байт)
    - Ответ: Использую BLACKFIN. Библиотечной функции с 24 или 32 битами я не нашел. — Николай (05.09.2005 10:14 80.73.170.86, пустое)
      - Вот нашел (+) — homekvn (05.09.2005 12:04 212.185.161.237, 401 байт)
        
        Ответ: Нашел. Очень полезная вещь, спасибо. — Николай (05.09.2005 14:45 80.73.170.86, пустое)
      - Надо где-то поискать. Мой коллега, который на Блекфине пишет, говорит, что видел эту 32-разрядную функцию. Поищите в инете. В крайнем случае, придется написать на ассемблере самому. — homekvn (05.09.2005 11:41 212.185.161.237, пустое)
Да динамики будут прим. одинаковы, что у КИХ, что у БИХ. — -=ВН=- (02.09.2005 10:19 194.190.181.231, 335 байт)
- Ответ: Многоступенчатый БИХ, наверно, вы имели ввиду (в первом предложении)? Или все-таки КИХ? /-/ — Николай (02.09.2005 10:35 80.73.170.86, пустое)
  - КИХ имел ввиду. — -=ВН=- (02.09.2005 11:09 194.190.181.231, 681 байт)
Ответ: Если владеете немецким(+) — Vadim Kudryavtsev (02.09.2005 10:15 210.94.41.89, 489 байт)
- Ответ: Нелинейность ФЧХ, в принципе, не помешает /-/ — Николай (02.09.2005 10:42 80.73.170.86, пустое)
  - подскажите, где она вообще может помешать? — 2b|!2b?.. (02.09.2005 16:59 212.46.19.12, пустое)
    - Ответ: — homekvn (02.09.2005 17:14 212.185.161.237, 503 байт)
- Пследняя фраза особенно хороша, типа был приказ лепить КИХ только и исключительно с линейной фазой:-) — -=ВН=- (02.09.2005 10:26 194.190.181.231, пустое)
а делают вообще — 2b|!2b?.. (02.09.2005 10:07 212.46.19.12, 99 байт)
- Ответ: да делают. Результаты нормальные - прореживается :-) — andy_P (05.09.2005 14:26 80.82.63.185, пустое)
- Ответ: А я еще не пробовал /-/ — Николай (02.09.2005 10:09 80.73.170.86, 1 байт)

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

E-mail: info@telesys.ru

Ответы

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