Що це sinc робить у моєму фільтрному виводі HP на LP?

Я намагаюся HP-фільтрувати сигнал, LP-фільтруючи його і віднімаючи вихідний вихідний фільтр. Для отримання сигналу, відфільтрованого LP, я фільтрую назад вперед, використовуючи фільтр Баттерворта третього порядку для сигналу зразка 16000. Частота відсікання становила 0,5 Гц при 250 Гц. Жодне вікно не використано.

У перших 85% або більше зразків висновок цього методу, здається, повністю перекривається з фактично відфільтрованим HP сигналом. Однак наприкінці випуску є надзвичайно гріхоподібне явище.

• Що це викликає?
• Що я можу змінити, щоб зробити цей метод здійсненним?

Ось два сигнали накладені. Синій - фільтрується НР, а червоний - відніманням.

Ось абсолютна різниця на зразок. Вони ніколи не бувають абсолютно рівними; мінімальна різниця - 8е-12. 87% значень нижче 0,01.

(Якщо є кращий спосіб отримати фільтр HP з фільтра LP , я з радістю дізнаюся його.)$H(z)$

Methinks @ hotpaw2 на вірному шляху. Я думаю, що filtfilt не справляється з початковими умовами зворотного фільтра. Це звучить як справді важка проблема - я не думаю, що я можу винуватий Mathworks за те, що він помилився ...

Виправлення здається мені досить простим: нульовий проміжок кінця сигналу, поки синхрон не піде, і видаліть (тепер ненулевий!) Прокладку після фільтра.

Ви робите фільтр "вперед-вперед" для одного, а другий - лише для вперед. Таким чином, два фільтри будуть вести себе по-різному в одному або обох краях.

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