Вычислительно эффективный КИХ-фильтр с линейной ФЧХ. Вопрос (+) («Телесистемы»: Конференция «Микроконтроллеры и их применение»)

Разработка, производство и продажа радиоэлектронной аппаратуры

| Карта сайта | Пишите нам | В избранное

Требуется программист в Зеленограде
- обработка данных с датчиков; ColdFire; 40 тыс.
e-mail:jobsmp@pochta.ru

Телесистемы | Электроника | Конференция «Микроконтроллеры и их применение»

Вычислительно эффективный КИХ-фильтр с линейной ФЧХ. Вопрос (+)

Отправлено vinogradov aleksei 29 марта 2007 г. 14:19

вот в чем. Пусть у нас есть фильтр со следующей импульсной характеристикой (количество коэффициентов характеристики нечетно):

h(0), h(1), h(2), h(3), h(4).

При этом характеристика симметрична, т.е. h(0)=h(4), h(1) = h(3).

Чувствую, что можно за счет этой симметрии сократить кол-во умножений почти вдвое. Но не могу выдумать вычислительно эффективный алгоритм. Наверняка такой алгоритм существует. Не подскажете, где с ним можно ознакомиться ?

Составить ответ | Вернуться на конференцию

Ответы

Так берите сумму 2-х соответсвующих отсчетов вх. данных и умножайте ее на к-т фильтра. Получите (N-1)/2 суммирований и столько же умножений-накопленний. Плюс еще одно уножение накопление на средний к-т. — -=ВН=- (29.03.2007 14:50:30 193.125.71.140, 371 байт)
- Именно так и делаю. Похоже, я неправильно вопрос задал (+) — vinogradov aleksei (29.03.2007 15:05:31 91.124.250.255, 337 байт)
  - Ну да, ну да, именно так и делаете, второй Ваш вопрос ну просто вытекает из первого, особенно в части числа умножений:-) Данные не двигают. Адреса по кругу бегают- циклические буфера для данных организуют. А с к-тами хитростей нет. — -=ВН=- (29.03.2007 15:19:18 193.125.71.140, пустое)
    - Да, я хвастун. Не делаю. Пока. Меня только хватило на формулу. А вот "конструкцию" буферов только придумываю :) Пока спасибо, буду думать ... — vinogradov aleksei (29.03.2007 15:25:13 91.124.250.255, пустое)
      - Вот сюда посмотри, симметрию добавь по вкусу --> — =AVR= (29.03.2007 15:34:10 80.92.96.19, пустое, ссылка)
        
        Тока осторожно. Там МАС для унсигнед. И если проццедура не обязательно должна ... — -mse- (29.03.2007 15:59:0 83.239.86.98, 1062 байт)
w0*(x0+x4)+w1*(x1+x3)+w2*x2 вместо w0*x0+w1*x1+w2*x2+w3*x3+w4*x4 — GM (29.03.2007 14:47:53 84.252.243.82, пустое)
:)))) Быстрая свёртка.. :))) — Dr_Alex (29.03.2007 14:23:0 213.33.211.66, пустое)
- А если фильтр не очень "длинный" ? Скажем, 31 коэффициент ? — vinogradov aleksei (29.03.2007 14:26:23 91.124.250.255, пустое)
  - А какая разница? — Dr_Alex (29.03.2007 14:27:37 213.33.211.66, пустое)
    - У Айчифера написано, что при количестве отсчетов до 128 быстрая свертка неэффективна ... — vinogradov aleksei (29.03.2007 14:31:30 91.124.250.255, пустое)
      - Не число умножений надо сокращать, а оптимизировать повторное использование коэффициентов, позаботившись предварительно о ПРАВИЛЬНОМ их размещении в регистрах, а не в ОЗУ - при малой длине КИХ это будет весьма эффективно, тем более для такого тормоза как ARM — =AVR=   (29.03.2007 14:37:30 80.92.96.19, пустое)
        
        Для AVR :) — vinogradov aleksei   (29.03.2007 14:38:15 91.124.250.255, пустое)
        
        Те же яйца - load/store architecture — =AVR=   (29.03.2007 14:58:33 80.92.96.19, пустое)
      - Чё правда? :-\ А я думал она где-нить при 4 станет неэффективной.. :-о — Dr_Alex (29.03.2007 14:33:46 213.33.211.66, пустое)