สิ่งนี้เป็นสิ่งที่ทำในเอาท์พุทตัวกรอง LP-turn-HP ของฉัน

ฉันพยายาม HP กรองสัญญาณโดย LP- กรองมันและลบผลลัพธ์จากตัวกรองเดิม เพื่อให้ได้สัญญาณกรอง LP- ฉันย้อนกลับไปข้างหน้ากรองโดยใช้ตัวกรอง Butterworth ลำดับที่สามสำหรับสัญญาณ 16,000 ตัวอย่าง ความถี่ตัดเป็น 0.5 Hz @ 250 Hz ไม่มีการใช้หน้าต่าง

ในตัวอย่างแรก 85% หรือมากกว่านั้นผลลัพธ์ของวิธีนี้ดูเหมือนจะซ้อนทับกันอย่างสมบูรณ์กับสัญญาณที่กรองโดย HP อย่างไรก็ตามในตอนท้ายของการส่งออกมีปรากฏการณ์ซินที่เหนือชั้น

• อะไรเป็นสาเหตุ
• ฉันจะเปลี่ยนอะไรเพื่อทำให้วิธีนี้เป็นไปได้

นี่คือสัญญาณทั้งสองที่ซ้อนกัน สีน้ำเงินถูกกรองโดย HP และสีแดงที่ได้จากการลบ

นี่คือความแตกต่างที่แน่นอนต่อตัวอย่าง พวกเขาจะไม่เท่ากันอย่างแน่นอน; ความแตกต่างขั้นต่ำคือ 8e-12 87% ของค่าต่ำกว่า 0.01

(หากมีวิธีที่ดีกว่าในการรับตัวกรอง HP จากของตัวกรอง LP ฉันยินดีที่จะเรียนรู้)$H(z)$

Methinks @ hotpaw2 อยู่ในเส้นทางที่ถูกต้อง ฉันเดาว่า filtfilt ไม่ได้จัดการเงื่อนไขเริ่มต้นสำหรับตัวกรองย้อนกลับค่อนข้างถูกต้อง ดูเหมือนว่าเป็นปัญหาที่ยากมาก - ฉันไม่คิดว่าฉันผิดที่ Mathworks จะทำผิด ...

การแก้ไขดูเหมือนจะง่ายพอสำหรับฉัน: zero-pad สิ้นสุดสัญญาณของคุณจนกว่า sinc จะหายไปและลบช่องว่างภายใน (ตอนนี้ไม่ใช่ศูนย์!) หลังจากตัวกรอง

คุณกำลังทำการกรองตัวกรองไปข้างหน้า - ข้างหลังสำหรับตัวกรองหนึ่งและตัวกรองสำหรับส่งต่อ - อย่างเดียว ดังนั้นตัวกรองทั้งสองจะทำงานแตกต่างกันที่หนึ่งหรือทั้งสองขอบ

หากสัญญาณของคุณไม่หยุดนิ่งและเป็นระยะอย่างแน่นอนในความยาวของหน้าต่างข้อมูลคุณสามารถลองห่อข้อมูลเป็นวงกลมและเริ่มตัวกรองของคุณเริ่มต้นด้วยความยาวตอบสนองของแรงกระตุ้น