Фильтрация входного сигнала. Опции

[Xtal1]

Всем доброго времени суток.
Прошу прощения, если не в тот раздел, но более подходящего не нашел.

Проблема заключается в следующем:
Есть поток входных значений, для примера даю выборку из 1000 отсчетов

 sampling.rar ( 9.22 килобайт ) Кол-во скачиваний: 164


Нужно фильтровать такой поток значений, дабы получить "болтанку" на выходе не более 2-3 единиц по амплитуде. Тупым усреднением этого получается достичь лишь при количестве усредняемых значений близкому к 1000.
Можно ли добиться подобного результата с помощью каких-то более серьезных алгоритмов фильтрации и получить при этом время отклика фильтра хотя бы 200-400 сэмплов.

[Guest_TSerg_*]

Это все понятно, но медиана без помех введет нелинейность в сигнал и его оценку, да и ТС сильно беспокоился по поводу 2 дискрет точности выходного сигнала.

P.S.
>однако видно, что шум не симметричный
На крайней картинке надо считать это уже не шумом, а сигналом и поэтому - медиана внесет нелинейность.

P.P.S
Для реккурентного SMA ( скользящее среднее ) вообще выборка только одна
SMA[i] = SMA[i-1] + (X[i] - X[i-n])/n

[AndrewN]

ТС сильно беспокоился по поводу 2 дискрет точности выходного сигнала.

На графике измеренного сигнала видно, что максимальные отклонения равны около 40 в обе стороны, т.е. чуть больше 5 младших бит могут врать, о +/-2 голова может не болеть.

Я поменял окно медианы на 9, даже такое короткое хорошо давит эти +/-40. На мой взгляд (поскольку истинное значение нам никогда не узнать, и приходится довольствоваться оценкой) - задача решена.

P.S. А это была хорошая идея про +/-40. Предположив, что 5 и 3 биты сбоят (40 = 0х28 = 10 1000) почти всегда наоборот, я прибавлял или вычитал 40 из выпадающих значений:

CODE

        if (x[k] < 36629)
            res = x[k] + 40;
        else if (x[k] > 36640)
            res = x[k] - 40;
        else
            res = x[k];

Результат превзошёл все ожидания. Убедитесь сами. Похоже, что эти биты и в правду сбоят!

Сообщение отредактировал AndrewN - Jan 30 2014, 19:18

[thermit]

На графике измеренного сигнала видно, что максимальные отклонения равны около 40 в обе стороны, т.е. чуть больше 5 младших бит могут врать, о +/-2 голова может не болеть.

Я поменял окно медианы на 9, даже такое короткое хорошо давит эти +/-40. На мой взгляд (поскольку истинное значение нам никогда не узнать, и приходится довольствоваться оценкой) - задача решена.

- vale.

Вы придумали собственную задачу и бодро взялись ее решать (не решив, кстати говоря). Между тем с вопросом тс можно ознакомиться в 1 посте.