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

Поставить 8 CPU.

Или реализовать фильтры на ПЛИС. А иначе никак, ИМХО, "синхронной" фильтрации по 8 каналам не достичь ...
Хотя, если у вас время между отсчетами достаточное большое ... Получаем, например, по 32 отсчета с каждого канала и для каждого из каналов вызываем функцию фильтрации. Таких вызовов будет 8. Главное, чтобы все они уложились в промежуток между входными отсчетами (с некоторым запасом), тогда при получении 33-го входного отсчета мы сможем выдать выходной отсчет, запомнить 33-й входной отсчет и "пустить его в работу" ...

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

Если вам не нужно 1000 герцовых измерений, то на входе АЦП ставьте аналоговый фильтр, который отрезает выше 250, тогда фильтрация при первой децимации не понадобится.

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

Я так понял, что надо на лету в реалтайме обрабатывать.

А я понял, что надо одновременно брать отсчеты по восьми каналам и одновременно обрабатывать. Короче, пусть ТС объяснится ...

Можно и так. Но есть и другой способ:
Записали в буфер и три раза разными дорожками пустили. Все зависит от ваших целей.

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

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


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

Что-то Вы путаете кислое с пресным. На одном процессоре невозможно параллельно обработать многоканальный поток данных - только последовательно.

Можно организовать быстрый ввод по одному отсчету с каждого канала и до поступления очередных отсчетов реализовать FIR-ы по всем каналам.
Сдвиг по времени каждого отсчета определиться максимально допустимой частотой переключения опроса датчиков.
Время - условное понятие. Если известна частота опроса датчиков, то shift легко вычисляется.

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

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

Если входной сигнал сразу отфильтровать (тем самым, ААА-фильтром) под частоту дискретизации после децимации, то и фильтрации не надо.
Среднее арифметическое - тоже фильтр НЧ.
Простейшая децимация - просто отбросить ненужные выборки, как будто их и не было.

Смотря какой фильтр

так у меня уже сделано - брал 2-й, 4-й и 8-й отсчет
просто не хочу нехороших эффектов с наложением.
мне кажется цифровой фильтр все таки должен быть.