LP-torbalı HP filtre çıkışımda bu samimi ne yapıyor?

LP filtreleme ve orijinal filtreden çıkış çıkarma ile bir sinyali HP filtrelemeye çalışıyorum. LP ile filtrelenmiş sinyal elde etmek için, 16.000 örneklik bir sinyal için üçüncü dereceden bir Butterworth filtresi kullanılarak geriye doğru filtrelenmiş olurum. Kesme frekansı 250 Hz'de 0,5 Hz idi. Hiçbir pencere kullanılmadı.

Örneklerin ilk% 85'inde, bu yöntemin çıktısı, aslında HP filtreli sinyal ile tamamen örtüşüyor gibi görünüyor. Bununla birlikte, çıktının sonunda, içtenlikle samimi bir fenomen var.

• Ona ne sebep oluyor?
• Bu yöntemi uygulanabilir yapmak için neyi değiştirebilirim?

Burada iki sinyal kaplanmıştır. Mavi olan HP ile filtrelenir ve kırmızı olan çıkarma ile elde edilir.

İşte örnek başına mutlak fark. Asla tam olarak eşit değildirler; minimum fark 8e-12'dir. Değerlerin% 87'si 0.01'in altındadır.

(Bir LP filtresinin sinden bir HP filtresi almanın daha iyi bir yolu varsa , bunu öğrenmekten memnuniyet duyarım.)$H(z)$

Methinks @ hotpaw2 doğru yolda. Sanırım filtfilt, ters filtre için başlangıç ​​koşullarını tam olarak işlemiyor. Kulağa gerçekten zor bir sorun gibi geliyor - Mathworks'ü yanlış yaptığı için hatalandırabileceğimi sanmıyorum ...

Düzeltme benim için yeterince kolay görünüyor: sinc gidene kadar sinyalinizin sonunu sıfırlayın ve filtreden sonra (şimdi sıfır olmayan!) Dolguyu silin.

Biri için bir ileri geri filtre ve diğeri için yalnızca bir ileri filtre yapıyorsunuz. Böylece iki filtre bir veya her iki kenarda farklı şekilde davranacaktır.

Sinyaliniz güçlü bir şekilde durağansa ve veri penceresi uzunluğunda tam olarak periyodikse, verileri dairesel olarak sarmayı deneyebilir ve filtrelerinizi dürtü yanıt uzunluklarıyla başlatabilirsiniz.