caxapa.ru :: Простейший а-ля медианный фильтр. Инициализируешь 32-х битную целую переменную "Фильтр" первым измеренным значением АЦП, умноженным на 2^10, к примеру. Каждый период отсчётов, берёшь значение "Фильтра" и делишь его на 2^10. Сравниваешь полученное значение и значение свежего отсчёта АЦП. Если код АЦП больше делённого значения фильтра, то инкрементируешь "Фильтр". Если мешьше - декрементируешь. Если равно - ничего не делаешь. Выход фильтра

Простейший а-ля медианный фильтр. Инициализируешь 32-х битную целую переменную "Фильтр" первым измеренным значением АЦП, умноженным на 2^10, к примеру. Каждый период отсчётов, берёшь значение "Фильтра" и делишь его на 2^10. Сравниваешь полученное значение и значение свежего отсчёта АЦП. Если код АЦП больше делённого значения фильтра, то инкрементируешь "Фильтр". Если мешьше - декрементируешь. Если равно - ничего не делаешь. Выход фильтра - величина "Фильтр" поделённая

на 2^10 - уже получена для сравнения.

Таким образом, на выходе получишь практически медианное значение за период ~1024 отсчёта, с точностью +/-1/1024. Показатель степени выбираешь по результатам тестирования - размениваешь скорость реакции на стабильность замеров.

P.S. Пожалуй, с 10-й степенью я загнул. Скачок с 2^14 до нуля займёт 2^14*2^10 отсчётов, 9 часов где-то. Или я что-то путаю. Сам использовал такой фильтр для компенсации дрейфа порогового значения при ловле импульсов небольшой амплитуды. Время установления было минуты при 19кГц выборки.

Ответить

- Спасибо! Cкpипaч(401 знак., 24.05.2023 22:58 - 29.05.2023 08:53)
  - Когда используете знаковые целые (int32_t, например), не используйте сдвиги. Просто пишите *256 или /256. Иначе нарвётесь на знак, рано или поздно. Сдвиг со знаком без специальных мер не работает. А компилятор достаточно умён, чтобы заменить 256 на сдвиги с должным учётом знака. Nikolay_Po(293 знак., 24.05.2023 23:09 - 23:16)
    - Починили всего-то в c++20 - RxTx(25.05.2023 22:16)
      - А по-подробнее? Сдвиги знак учитывать стали? Оно (ссылка)? Nikolay_Po(1 знак., 25.05.2023 22:27, ссылка)
        
        ISO/IEC JTC1 SC22 WG21 N 4860 [expr.shift] 7.6.7 Shift operators RxTx(59 знак., 25.05.2023 23:16, ссылка)
        
        Спасибо! - Nikolay_Po(26.05.2023 10:51)
    - Собственно 2^4 это константа, которую можно выбрать другой. С тем чтобы она была выше возможной помехи. - Cкpипaч(24.05.2023 23:15)
      - Вижу, понимание есть. Дерзайте! Всё получится. - Nikolay_Po(24.05.2023 23:17)
- Ссылочку на источник бы. - RxTx(26.05.2023 00:28)
  - Да где-то в какой-то конференции, типа stack overflow, в ответах. Нашёл, когда захотел получать медианное значение по большому количеству отсчётов, а ни памяти, ни процессорного времени для вычисления "в лоб" не хватило бы. - Nikolay_Po(26.05.2023 10:57)
    - Посмотри историю браузера будь ласка :) Мне интересно каким образом это может быть медианой. - RxTx(27.05.2023 12:37)
- Спасибо! Красивое! Dingo(354 знак., 25.05.2023 05:27)
- Думать лень про написаное.. но разве это не то же самое y = y*k + x*(1-k)... для целочисленного коэффициента 2^n всё решается сдвигами и суммированием. - POV(24.05.2023 22:48)
  - Нет, не то же самое. В вашем примере, при большом отклонении, к значению аккумулятора прибавится большее значение. Вклад помехи зависит от её величины. В моём алгоритме (нашёл его, когда искал быстрый медианный фильтр) возмущение аккумулятора помехой не больше единицы, вне зависимости от величины помехи. - Nikolay_Po(24.05.2023 23:02)
- Брр. Вроде бы врубился. Посмотри, правильно ли я понял идею? Cкpипaч(204 знак., 24.05.2023 21:57)
  - Идеально! Только порядок придётся снизить, иначе слишком долго будет отрабатывать большие изменения. Нужно на практике проверять. - Nikolay_Po(24.05.2023 23:04)
- Фраза "инкрементируешь фильтр" - опечатка? Можно попросить вас формулу написать? (сейчас как-раз пытаюсь "изобрести БИХ в целых числах", похоже вы описываете именно его, но как-то не сходится). - Cкpипaч(24.05.2023 20:52)
  - Назови переменную Filter. И инкрементируй/декрементируй её по условию. Формулу как писать, с условием, вроде помню, но на движке Сахары этого не изобразить. - Nikolay_Po(24.05.2023 20:59)
    - Переменную правильно назвать аккумулятор с декрементом и инкрементом. - Visitor(24.05.2023 21:23)
      - Как хош называй. Да, аккумулятор - подходит. - Nikolay_Po(24.05.2023 21:25)
- А классика, то ли Титце и Шенк, то ти Хоровитц и Хилл так статанализ для какой то установки делали, я по мотивам разбор однотональных частот процем сделал, инкремент, если в зону временную попала, декремент, если нет. И штук 8 таких фильтров в ПИКе, работало. - Visitor(24.05.2023 20:07)
  - АОН-ДТМФ? - Alex68(24.05.2023 20:15)
    - Нет, там 2 тона, просто не обработать, это для проводных линий пейджинга было. Типа внутриканальная адресация. - Visitor(24.05.2023 20:19)
      - Bell 201\212 или чего-то там такое? - Alex68(24.05.2023 20:21)
        
        Да, сообщения им шли;-) Потом и это сделали Белл 202. - Visitor(24.05.2023 20:32)

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