Отправлено rst 05 сентября 2006 г. 08:36

Как корректно посчитать 32bit fixed-point IIR в FDAtool?
Ставлю галку Квантование и все галки оптимизации (в том числе - чтоб не переполнялись коэфф. и <=1). Полученный фильтр в float (моделирую на компе) прекрасно пашет, при переходе на s32 - переполнение внутри фильтра. На float-модели замерил диапазон значений внутри фильтра при заданном сигнале - оказалось диапазон ~= +-7, отмасштабировал вх./вых. s32 фильтра на 8 - заработал, но при другом сигнале - другой диапазон!
Линию задержки обнулить не забыл ;)
Так как все-таки так рассчитать IIR32, чтоб нигде в нем не было перехода за 32 бита? (ну или хотя бы чтоб данные в линии задержки не превышали 32)?
TMS320VC5502. Надо ФВЧ - вырезать постоянку (и 100Гц) на ~50дБ, FS=48кГц. Нижняя частота полосы пропускания = не хотелось бы ниже 200Гц, поэтому КИХ получается довольно громоздким - жалко тактов ;)
В общем - как рассчитать IIR так чтоб все числа не превышали некого диапазона?

• Масштабный к-т можно получить, если вычислить sum(|h[n]|). h[n]-ИХ фильтра. — -=ВН=-   (05.09.2006 16:47 193.125.71.140, 426 байт)
  • Вроде в его случае (+) — SM   (05.09.2006 17:35 195.225.131.186, 214 байт)
    • Последний абзец у меня это учитывает:-) — -=ВН=-   (05.09.2006 17:37 193.125.71.140, пустое)
      • Абзец состоит из одного предложения:-) — -=ВН=-   (05.09.2006 17:38 193.125.71.140, пустое)
        • Да я ни грамма не против, что для каждой искать надо, согласно последнему абзецу. Я про то, что из H, перед тем как h искать, надо нули устранить. — SM   (05.09.2006 17:54 195.225.131.186, пустое)
          • Так он автоматом устраняются. Эта структура автоматом образует 2 секции. Первая БИХ с одними полюсами. Вторая КИХ. Так что раздельный поиск все решает:-) — -=ВН=-   (05.09.2006 18:01 193.125.71.140, пустое)
            • Ну если с такой стороны подойти, то конечно. Это я все к тому, чтобы народ не запутался, где какой H брать :) и почему. — SM   (05.09.2006 18:05 195.225.131.186, пустое)
              • Да я не возражаю, наоборот:-) — -=ВН=-   (05.09.2006 18:08 193.125.71.140, пустое)
• Скорее всего (+) — SM   (05.09.2006 10:58 213.141.159.26, 266 байт)
  • Да, действительно (+) — SM   (05.09.2006 11:00 213.141.159.26, 187 байт)
    • Хорошо попробую. Такую структуру выбрал так одна линия задержки хорошо ложится на dual-MAC - одно звено = ~10тактов — rst   (05.09.2006 11:48 195.162.42.34, пустое)
      • Если такты на первом месте - то масштабируйте вход, тем более, если он 16-битный. Исходя из анализа фильтров с убранными нулями. — SM   (05.09.2006 11:57 213.141.159.26, пустое)
        • Попробовал структуру звеньев с двумя линиями задержки - результаты обнадеживают — rst   (06.09.2006 12:14 195.162.42.34, 394 байт)
    • Полностью согласен. (+) — homekvn   (05.09.2006 11:40 212.185.161.237, 1219 байт)
• Какая структура вашего фильтра ? Возможно причина переполнения именно в этом — mag24   (05.09.2006 09:12 62.148.158.178, пустое)
  • На биквадратных звеньях. — rst   (05.09.2006 10:08 195.162.42.34, 241 байт)
    • Блин! конечно же: d0 = x - a1 * d1 * 2 - a2 * d2 * 2 y = d0 + b1 * d1 * 2 + b2 * d2 * 2 — rst   (05.09.2006 10:09 195.162.42.34, пустое)