ati » from archive
Друзья, а есть тут кто-нибудь, кто помнит математику, относящуюся к DSP? Фильтры и вокруг. Я бы хотел пару дурацких вопросов задать. Pls RT.
Мне собственно надо понять, как отфильтровать всплески сигнала вот такого вида: http://i.imgur.com/UKVX6.png ‎- ati
Пока из всего, что я пробовал из mFilter и signal в R более-менее внятно выглядел НЧ-фильтр Баттерворта, но всё равно не идеально, вот так: http://i.imgur.com/OZkGq.png . Мне бы просто названий методов, приспособленных для этой задачи, дальше бы я сам. ‎- ati
А как должно быть "идеально"? Тупо Fourier transform - подрезать края спектра - reverse Fourier transform не пойдёт? ‎- 9000
навскидку - любой хороший ФНЧ. Если линейность фазы не нужна - БИХ фильтр. У Баттерворта - лучше линейность в полосе пропускания - у Чебышева - круче спад, но пульсации в полосе пропускания ‎- гавс
ой, там по оси y не частота же... короче надо сделать БПФ, посмотреть частоты в "хорошем сигнале" и установить частоту среза повыше раза в 1.5 ‎- гавс
а что кстати за сигнал, можно поподробнее? ‎- гавс
есть еще такая область change detection (н-р http://libgen.info/view.php?id=223779 ) - можно обнаруживать участки со всплесками и фильтровать только их. ‎- гавс
Простите, не ответил @9000 и @nazarovsky -- уехал постик и я не заметил коментов. В итоге выбросы отфильтровал rolling _median_-ом, вручную подобрав ширину окна. Получилось примерно так: http://vozduh.msk.ru/img/weekly_filtered.png , не супергладко, но смысл вроде бы передаёт. ‎- ati
а ширина окна сколько? ‎- гавс
3001 минутных отсчетов. ‎- ati
разные ФНЧ мне не удалось толком настроить, как ни старался. ‎- ati
есть кстати быстрые алгоритмы для скользящей медианы, чтобы не сортировать каждый раз весь массив. ‎- гавс
Я очень надеюсь, что парни, контрибьютящие в http://cran.r-project.org/ именно их и использовали. ‎- ati
а, ну наверное. ‎- гавс
а у чипа есть свой DSP (blackfin, TI OMAP, etc) или все руцями писать? ‎- silpol
@silpol - у меня это уже на сёрвере фильтруется, R-овскими библиотеками. ‎- ati