Sự chân thành này đang làm gì trong đầu ra bộ lọc LP-Turn-HP của tôi?


10

Tôi đang cố gắng lọc tín hiệu HP bằng cách lọc LP và trừ đầu ra khỏi bộ lọc ban đầu. Để thu được tín hiệu được lọc LP, tôi đã lọc ngược về phía trước bằng cách sử dụng bộ lọc Butterworth bậc ba cho tín hiệu 16.000 mẫu. Tần số cắt là 0,5 Hz @ 250 Hz. Không có cửa sổ đã được sử dụng.

Trong 85% mẫu đầu tiên hoặc hơn, đầu ra của phương pháp này dường như trùng lặp hoàn toàn với tín hiệu được lọc thực sự của HP. Tuy nhiên, ở phần cuối của đầu ra, có một hiện tượng giống như siêu chân thành.

  • Điều gì gây ra nó?
  • Tôi có thể thay đổi gì để làm cho phương pháp này khả thi?

Dưới đây là hai tín hiệu chồng chéo. Màu xanh được lọc bởi HP và màu đỏ thu được bằng phép trừ.

Tín hiệu HP và orig-LP

Đây là sự khác biệt tuyệt đối cho mỗi mẫu. Chúng không bao giờ chính xác như nhau; sự khác biệt tối thiểu là 8e-12. 87% giá trị dưới 0,01.

sự khác biệt giữa hai phương pháp

H(z)


Bạn có thể bao gồm mã MATLAB mà bạn đang sử dụng để tạo các ô không? Dữ liệu không quan trọng; thay vào đó, bạn chỉ có thể mô phỏng nó với các mẫu ngẫu nhiên nếu cần.
Jason R

Câu trả lời:


0

Methinks @ hotpaw2 đang đi đúng hướng. Tôi đoán là bộ lọc không xử lý các điều kiện ban đầu cho bộ lọc ngược hoàn toàn đúng. Nghe có vẻ như là một vấn đề thực sự khó khăn - Tôi không nghĩ mình có thể phạm lỗi với Mathworks vì đã hiểu sai ...

Việc khắc phục có vẻ đủ dễ đối với tôi: zero-pad kết thúc tín hiệu của bạn cho đến khi sự chân thành biến mất và xóa phần đệm (bây giờ khác không!) Sau bộ lọc.


0

Bạn đang thực hiện bộ lọc chuyển tiếp ngược cho một bộ lọc và bộ lọc chỉ chuyển tiếp cho bộ lọc kia. Do đó, hai bộ lọc sẽ hoạt động khác nhau ở một hoặc cả hai cạnh.

Nếu tín hiệu của bạn đứng yên và chính xác theo định kỳ về độ dài của cửa sổ dữ liệu, bạn có thể thử bọc dữ liệu theo vòng tròn và bắt đầu bộ lọc của mình vượt quá độ dài đáp ứng xung của chúng.


Xin lỗi nếu tôi không đủ rõ ràng. Tín hiệu mà tôi đang so sánh với (màu xanh) cũng được lọc ngược ( filtfiltbutter(3,0.5/250,high)). Tín hiệu là từ một ECG mà tôi không thể mong đợi là chính xác định kỳ.
Andreas

Tôi có hiểu lầm gì không? Tôi không thể hiểu upvote vì đoạn đầu tiên rõ ràng là sai.
Andreas
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.