Как правильно фильтровать? Опции

[inventor]

Задача такая - у меня АЦП может работать только на минимальной частоте 1000 Гц.
мне нужны кратные частоты: 125, 250 и 500.
Получить их можно проведя децимацию, предварительно перед которой
отфильтровав все что выше чстоты найквиста.
т.е. для получения данных на частоте 500,
я запускаю измерения на 1000 Гц,
фильтрую все что выше 250
и беру каждый второй отсчет.
Не стал заморачиваться с фильтрами - смотрю они уже есть в армовсеом коде
(те у кого STM в папке SMSIS)
у меня процессор СС3200 - не нашел примеров именно для него,
пока взял библиотечку у STM - она подходит к техасу (ф-я void arm_fir_decimate_fast_q31(///)).
Есть несколько вопросов, у меня может быть от 1-го канала до 8-ми,
вот как их всех одновременно фильтровать?
Код, который делает АРМ я так понимаю УЖЕ оптимизирован
по скорости, объему и прочему, сокращено количество циклов и всяких MAC операций.
как сделать "нормальную" фильтрацию одновременно всех 8-ми каналов,
и сколько навскидку времени уйдет на вызов функции фильтрации
т.е. переписать саму функцию фильтрации, чтобы использовать только 1 функцию
вместо идкущих подряд 8 и что то другое.
у кого какие мысли?

[ViKo]

Зачем "одновременно"? Обрабатывайте блоки от каждого канала по-очереди. Блок сделайте равным коэффициенту децимации. После каждого вызова функции получите одну точку результата для одного канала. Чисто теоретически, imho.

[inventor]

Зачем "одновременно"? Обрабатывайте блоки от каждого канала по-очереди. Блок сделайте равным коэффициенту децимации. После каждого вызова функции получите одну точку результата для одного канала. Чисто теоретически, imho.

то есть сначала для первого отсчета, через 32 отсчета для второго, через 64 для третьего?
Правильно?
насчте синхронной филттрации.
у меня не звук, где это было бы допустимо
там сейсмоданные и желательно знать когда сигналы к чему относятца.
сейчас у меня обычный FIR фильтр, который дает задержку на 32 отсчета,
т.е. на 32 миллисекунды.
хотелось бы чтобы время, которое я записываю с данными датчиков
соответсвовало первому отсчету в пачке.
а таким образом (не одновременная фильтрация)
время будет соответсвовать только примерно.

[ViKo]

то есть сначала для первого отсчета, через 32 отсчета для второго, через 64 для третьего?
Правильно?

Функция обрабатывает блоки данных. Ваша задача сформировать эти блоки. Если у вас 8 каналов, распихивайте принятые с АЦП данные по 8-ми массивам (блокам). И вызывайте функцию для обработки каждого из блоков. Я бы сделал входной блок размером с к-т децимации, и как только он заполнен, вызывал бы функцию.

Вообще, какая там математика - децимация. Возьмите 2(4, 8) выборки, сложите и поделите на 2 (4, 8).

сейчас у меня обычный FIR фильтр, который дает задержку на 32 отсчета,
т.е. на 32 миллисекунды.
хотелось бы чтобы время, которое я записываю с данными датчиков
соответсвовало первому отсчету в пачке.

Сдвиньте результат на 32 шага назад.

[Tarbal]

Вообще, какая там математика - децимация. Возьмите 2(4, 8) выборки, сложите и поделите на 2 (4, 8).

Децимация действительно легко делается, но предшествующая ей фильтрация требует подсчета свертки.
Да еще и с высокой точностью, что обычная целая математика не прокатит.

Сообщение отредактировал Tarbal - Oct 28 2015, 13:58