Вопрос:можно ли устранить данную задержку?

Если фильтруете всегда блок фиксированной длинны, то filtfilt() как раз для этого.

Да,блок всегда фиксированной длины.а если использовать filter?как быть в этом случае?

Например можно добавить в начало сигнала 40 (длинна фильтра/2) отсчетов сигнала с обратным знаком и в обратном порядке во времени. То есть если есть вектор [1 2 3 4 5 .... 100], то его нужно превратить в вектор [-40 -39 .... -1 1 2 3.... 100] и после филтрации начало обрезать. В конец тоже нулей добавить, чтобы не потерять хвост.

Но filtfilt() не хуже совсем, она делает 2 прохода - в одну и другую сторону, чтобы скомпенсировать искажение фазы.

deneb13, делайте с переносом в 0 и обратно, в нуле задавите вашу помеху ФВЧ фильтром, это нормальный подход. FFT здесь не нужно. Блочность не должна влиять на обработку, совпадать она должна с потоковой, это всего лишь манипуляции с данными.

Это несомненно позволит добиться нужного результата, но не совпадет с темой курсовой, так как так или иначе помеху просто отфильтровали, а не оценили и вычли. То что требуется - это извращение еше то, конечно, но чего не сделаешь ради образования

PS: хотя если применить ФНЧ и в 0 снести именно помеху, то на выходе как раз получится чисто помеха, ее обратно на место перенести и можно вычитать.

Сообщение отредактировал Taradov Alexander - Mar 9 2011, 12:15

А что есть какая-то разница выделить сигнал ФНЧ фильтром и вычесть из задержанного сигнала или сразу ФВЧ задавить? Это буквально одно и то же. ИХ ФВЧ фильтра так и получается из ФНЧ прототипа задержкой и вычитанием.

Я понимаю, просто курсовая тогда очень странная получается, больше на ДЗ похожа.

именно это я и пытаюсь сделать!вот только при вычитании, помеха до абсолютного 0 не уходит

Нужно смотреть на задержки фазы, вносимые фильтрами. Ну и плного прямо нуля не выдет, немного все-равно останется.

Я вечерком попробую такое сотворить.

Буду благодарна,если расскажете о своих результатах!

Я решил попробовать для начала просто одну помеху подавить таким методом, но наткнулся на эффект, который сходу объяснить не могу, может коллективный разум поможет.

Код:


Тут происходит перенос помехи в 0, фильтрация ФНЧ и перенос результата назад.

Проблема - если частота помехи низкая (f1 = 20 Гц), то imag(sh) около 0, а real(sh) совпадает по фазе и амплитуде с исходным синусом. Но если частоту увеличивать, то постепенно действительная и мнимая части становятся амплитудой 0.5 и не совпадают по фазе с исходным синусом.

Похоже это из-за фильтрации, но я не могу с ходу понять почему именно.

Угу, из-за фильтрации. Тот фильтр, что у Вас, просто давит зеркальную часть смещенного в ноль синуса. Что там происходит при обратном переносе мне не ведомо (лень формулы писать), но явно не то, что хотелось бы. Полосу пропускание фильтра надо бы увеличить, или делать комплексный фильтр.

Да, действительно, это из-за того, что фильр с действительными коэффициентами.

В общем идею это иллюстрирует, только нужно правильно фильтр посчитать.

Но вообще так избаляться от помех - это изврат.

Скажите, а как можно посчитать комплексный фильтр?может литература есть какая по этому вопросу?дело в том,что если в качестве сигнала использовать шпс, то нельзя сильно увеличивать ПП фильтра-сигнал ведь тоже вырезаться будет.