LP 전환 HP 필터 출력에서이 sinc가 수행하는 작업은 무엇입니까?


10

LP를 필터링하고 원래 필터에서 출력을 빼서 신호를 HP 필터링하려고합니다. LP 필터링 된 신호를 얻기 위해 16,000 샘플 신호에 대해 3 차 버터 워스 필터를 사용하여 역방향으로 필터링했습니다. 차단 주파수는 250Hz에서 0.5Hz입니다. 창을 사용하지 않았습니다.

샘플의 처음 85 % 정도에서이 방법의 출력은 실제로 HP 필터링 된 신호와 완전히 겹치는 것 같습니다. 그러나 출력이 끝나면 상당히 비슷한 현상이 발생합니다.

  • 무엇이 원인입니까?
  • 이 방법을 실현시키기 위해 무엇을 바꿀 수 있습니까?

겹쳐진 두 신호는 다음과 같습니다. 파란색은 HP- 필터링되고 빨간색은 빼서 얻습니다.

HP 및 orig-LP 신호

다음은 샘플 당 절대 차이입니다. 그들은 결코 정확히 같지 않습니다. 최소 차이는 8e-12입니다. 값의 87 %가 0.01 미만입니다.

두 방법의 차이점

(LP 필터의 에서 HP 필터를 얻는 더 좋은 방법이 있다면 기꺼이 배우겠습니다.)H()


플롯을 생성하는 데 사용하는 MATLAB 코드를 포함 할 수 있습니까? 데이터는 중요하지 않습니다. 필요한 경우 무작위 샘플로 시뮬레이션 할 수 있습니다.
Jason R

답변:


0

Mehotks @ hotpaw2이 올바른 길에 있습니다. 필트 필트가 역 필터의 초기 조건을 제대로 처리하지 못한다고 생각합니다. 정말 어려운 문제인 것 같습니다. Mathworks가 잘못했다고 잘못 생각할 수는 없습니다 ...

수정은 나에게 쉬운 것처럼 보입니다 : sinc가 사라질 때까지 신호의 끝을 제로 패드로 채우고 필터 뒤의 (현재 0이 아닌!) 패딩을 삭제하십시오.


0

한 필터는 역방향 필터를 수행하고 다른 필터는 전달 전용 필터를 수행합니다. 따라서 두 필터는 한쪽 또는 양쪽 가장자리에서 다르게 작동합니다.

신호가 고정되어 있고 데이터 창 길이에서 정확히주기적인 경우 데이터를 원형으로 감싼 후 임펄스 응답 길이에 따라 필터를 시작할 수 있습니다.


내가 명확하지 않으면 죄송합니다. 내가 비교하고있는 신호 (파란색)도 뒤로 필터링됩니다 ( filtfiltbutter(3,0.5/250,high)). 신호는 정확히주기적인 것으로 예상 할 수없는 ECG에서 나옵니다.
Andreas

내가 뭔가를 오해하고 있습니까? 첫 번째 단락이 분명히 잘못 되었기 때문에 공감대를 이해할 수 없습니다.
Andreas
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.