caxapa.ru :: Вот это:

Nikolay_Po (25.04.2017 09:44 - 09:57, просмотров: 145) ответил fk0 на Он не думает. Можно посмотреть картинки и думать как закодить в 16-битах fixed point. Ещё раз рекомендую ссылку. В очередной раз.

Вот это: static float _fstepf_2(filter *fi, float x) { /* /usr/fisher/mipsbin/mkfilter -Bu -Bp -o 4 -a 0.03593750000 0.04531250000 -l */ float *v = fi -> v; v[0] = v[1]; v[1] = v[2]; v[2] = v[3]; v[3] = v[4]; v[4] = v[5]; v[5] = v[6]; v[6] = v[7]; v[7] = v[8]; v[8] = ( 6.9752013701e-07 * x) + ( -0.8573132535 * v[0]) + ( 6.7668969461 * v[1]) + (-23.5932085800 * v[2]) + ( 47.4460745570 * v[3]) + (-60.1843751000 * v[4]) + ( 49.3079205940 * v[5]) + (-25.4811693030 * v[6]) + ( 7.5951586426 * v[7]); return (v[0] + v[8]) - 4 * (v[2] + v[6]) + 6 * v[4]; } С 16-бит фиксированной точкой (с целыми числами - без разницы) не заработает. Слишком высок порядок фильтра и ошибка квантования получит большое усиление, до потери ~~стабильности~~ устойчивости.

Ответить

- Во-первых это нерекурсивный фильтр -- с чего здесь ошибке умножаться? Она больше чем на -60 здесь нигде не умножается. Во-вторых там в каталоге fast целочисленная реализация, вроде. В 16 битах рекурсивный фильтр вообще сложно сделать, IMHO. Я бы fk0(623 знак., 25.04.2017 11:17)
  - Да, вот здесь: Nikolay_Po(188 знак., 25.04.2017 11:53)
    - Да, действительно. Из малоресурсозатратного для микроконтроллеров вспоминаются CIC-фильтры в качестве ФНЧ (уж от скользящего среднего там один шаг). - fk0(26.04.2017 00:48)
      - CIC рассматривал. Не подошли, в сравнение не идут с эллиптическими. В чистом виде АЧХ гребенчатая, требуется каскадирование с другим фильтром. По сложности проще каскад эллиптических поставить. Разве что за исключением случаев, где последние не Nikolay_Po(39 знак., 26.04.2017 09:40)
        
        Как-то утверждение не вяжется с указанием на чрезмерную вычислительную сложность БИХ фильтров. CIC как раз легко считается на 8/16-битном контроллере в целых числах, корректирующий фильтр (КИХ) -- тоже. - fk0(26.04.2017 13:13)
        
        КИХ вообще ни о чём, по сравнению с БИХ. Я с какой стороны ни подступался, КИХ проигрывал БИХ-фильтру при доступности 150..250 тактов CPU на отсчёт. Без SIMD и MAC поддержки CPU КИХ для моей задачи оказались не интересны. Не удалось синтезировать Nikolay_Po(307 знак., 26.04.2017 14:23)
        
        КИХ с симметричной (или антисимметричной) ИХ хорош там, где нужна линейная ФЧХ. USSR(68 знак., 27.04.2017 12:55, )
        
        Сколько отсчётов нужно на ФНЧ с полосой пропускания 0.15 и с ослаблением не хуже 40дБ с частоты 0.25 и выше? Неравномерность в полосе пропускания до 1дБ. - Nikolay_Po(27.04.2017 14:05)
        
        Частота в Найквистах? - USSR(27.04.2017 14:51, )
        
        Если в Найквистах, то MATLAB уверяет, что equiripple FIR 28-го порядка справится. - USSR(27.04.2017 14:57, )
        
        Да, в "Котельниковых". Я считал, у меня лишь 16-й порядок КИХ с трудом вмещается по доступному времени процессора и требует ручной оптимизации. Лучше 8й-10й порядок КИХ, а такого недостаточно для целевой характеристики. Эллиптический БИХ при том Nikolay_Po(116 знак., 27.04.2017 18:02)
        
        У нас в изделии КИХ фильтр на 2000 тапов крутится. Но это FPGA. :) - USSR(27.04.2017 18:16, )
  - Понимаю, о чём вы. Просто в PIC24, к примеру, сложение и умножение за такт делается лишь для 16-битных операндов. Деление за 19 тактов, но оно меняется на сдвиг при делителе 2^N. Nikolay_Po(1130 знак., 25.04.2017 11:46)
    - Деление практически всегда заменяется умножением... Прореживание (decimation) это как раз обычно типовая область применения CIC-фильтров, которые в каких-то условиях имеют меньшую вычислительную сложность, чем другие подходы (вроде снижения fk0(491 знак., 26.04.2017 13:33)
      - Касательно CIC фильтров. Nikolay_Po(1619 знак., 26.04.2017 14:45)

Микроконтроллеры