고유 벡터를 사용하여 MUSIC을 통한 신호의 기본 주파수 추정


14

문맥:

(면책 조항 : 이것은 통신 문제가 아닙니다).

실제 주기적 신호의 기본 주파수를 추정하려고합니다. 이 신호는 원시 신호를 펄스의 신호와 일치 필터링하여 구성되었습니다. (일치 된 필터). 결과 신호는 다음과 같은 특징이 있습니다.

  • 주기적입니다. (기본은 1 / 기간)이며 이것이 내가 추정하려고하는 것입니다.

  • 시간이 오래 걸리지 않습니다. 구체적으로, 주기적 펄스의 진폭은 진폭이 변할 수있다. (예를 들어, 한 펄스는 낮을 수 있고, 다른 펄스는 높을 수 있고, 다음 펄스는 다시 낮을 수 있으며, 그 매체 다음에 계속 될 수 있습니다).

  • 나는 주파수가 고정되어 있다고 생각합니다 (진폭은 변경하지만 밴드는 변경하지 않음).

  • 고조파 왜곡이 있습니다. 여기서 의미하는 것은 (그리고 내가 틀렸다면 나를 교정하십시오), 신호 내의 개별 펄스는 정현파가 아니라 가우시안, 삼각형, 반-포물선 등과 같은 '펑키'모양이라는 것입니다. .

이 신호의 기본 주파수를 추정하려고합니다.

물론, 원시 신호는 때로는 노이즈 일 뿐이지 만 여전히 경로를 통과하여 필터링됩니다. (나중에 더 자세히 설명).

내가 시도한 것 :

이제 저는 다음과 같은 수많은 기본 주파수 추정기를 알고 있습니다

  1. 자동 상관 법
  2. 음과 그 모든 의존성
  3. FFT 방법.

기타,

  • 음 : 아직 음을 시도하지 않았습니다.

  • FFT 방법 : FFT 방법은 모든 고조파와 기본 요소를 제공하지만 기본이 항상 최고점은 아니기 때문에이 고정되지 않은 비즈니스에서는 특히 까다로울 수 있습니다. 매우 빨리, 당신은 많은 피크 중 어느 것이 기본인지 확인하려고 노력하고 있으며, 이것이 어려운 문제가됩니다.

  • 자기 상관 : 자기 상관 방법은 FFT 방법보다 나은 것으로 보이지만 시간 영역 신호의 진폭 불규칙성에 여전히 민감합니다. 자동 상관 법은 중심 로브와 다음으로 가장 높은 로브 사이의 거리를 측정합니다. 그 거리는 기본에 해당합니다. 그러나 고정되지 않은 경우에는이 보조 로브가 너무 낮아서 일부 임계 값 구성표에서 누락 될 수 있습니다.

그런 다음 MUSIC과 같은 부분 공간 방법을 사용하여 기본을 추정 할 수 있다는 것이 나에게 일어났습니다. 이것을 테스트 한 결과, 신호의 기본에 해당하는 주파수에서 아주 좋은 결과를 얻었습니다. (찾는 신호 수를 2로 설정하면 신호의 공분산 행렬의 가장 높은 고유 벡터 2 개 (고유 값의 최고 값에 해당)를 선택하여 버리고 신호를 생성합니다. 나머지 부분의 노이즈 부분 공간, 가설의 복잡한 정현파를 투영하고 역수를 취하고 좋은 의사 스펙트럼 인 voila를 취하십시오.

질문과 문제 :

  1. 즉, 왜 이것이 더 잘 작동하는지 이해하고 싶습니다.
  2. MUSIC에서는 신호 부분 공간을 버리고 잡음 부분 공간을 사용합니다. 신호 부분 공간의 고유 벡터는 실제로 일종의 '최적 적합'인 것으로 보입니다. 실제로 최적의 일치 필터입니다. 신호 서브 스페이스 고유 벡터를 직접 사용하지 않는 이유는 무엇입니까? (나는 더 이상 음악을 알지 못하지만 왜 잡음 부분 공간을 더 잘 사용합니까?)
  3. 마지막으로, 최종 문제는이 방법이 정상이 아닌 신호에 대해 훨씬 더 강력하게 작동하는 것처럼 보이지만 (문제는 위에서 정의한 바와 같이), 시스템에 노이즈 만있는 경우에도 항상 응답을 얻는다는 것입니다! (위에서 언급 한 원시 사전 일치 필터링 된 신호는 주기적 신호가없는 경우 때때로 백색 잡음 일 수 있음).

이것을 막기 위해 어떤 방법이 존재할 수 있습니까? 나는 고유 값을 살펴 보았고 신호가있는 소음 VS 사례 만있는 경우 부패에 더 많은 '곡률'이 있지만 충분히 강하지 않을 수도 있습니다.

보너스:

  1. 공분산 행렬의 고유 벡터는 언제 다른 것입니까? 그것들이 정현파인지 아닌지를 결정하는 것은 무엇입니까? 왜 구형파가 아닌가? 아니면 다른 모양의 신호를 삽입 하시겠습니까?

Mohammad- 몇 가지 수정 / 설명을 해주시겠습니까? 용어를 고집하는 사람이 될 수 있지만 향후 방문자에게는 중요합니다. 'nice & clean'이외에도 고조파 왜곡을 말할 수 있습니다. 반복적이지 않고 주기적으로 말할 수 있습니까? 고정은 시변 통계 또는 시변 스펙트럼을 나타낼 수 있습니다. 당신은 명확히 할 수 있습니까? 자기 상관 방법은 Yule-Walker 방법의 별칭입니다. '신호 수'라고 말할 때이 실제 정현파 또는 복잡한 지수입니까? 최대 값 고유 값을 사용할 수 있습니까? 순위는 선형 대수에서 다른 의미를 갖습니다. '가장 높은 차이'와 동일 ...
Bryan

1
... (계속) 한가지 중요한 점은 (그리고 명확히 할 때이 대답에 주목할 것입니다) MUSIC 방법은 노이즈 부분 공간 방법이라는 것입니다. 따라서 이상적으로 고유 값이 가장 큰 신호 부분 공간 고유 벡터는 사용되지 않습니다. 또한 신호가 주기적이면 정현파의 합입니다. 주기적인 경우, 불연속 정현파의 합인 푸리에 계열로 정의 할 수 있습니다.
Bryan

@Bryan 돌아 오는 데 시간이 걸리는 것에 대해 죄송합니다 (긴 주말). 실제로 전체 질문을 곧 개선하여 알려 드리겠습니다. 감사합니다!
Spacey

@Bryan 나는 마침내 전체 게시물을 개정하고 제안을 추가했으며 많은 맥락 / 문제를 명확히했습니다. 참조하십시오. 다른 것을 명확히 할 수 있으면 알려주십시오.
Spacey

@Mohammad 고유 벡터의 "강도", 즉 고유 값에 의해 신호가 존재하는지 여부를 식별 할 수 있습니까?
Jim Clay

답변:


8

에프(,에스)=영형V(엑스(),엑스(에스))=영형V(엑스(),엑스(에스))=에프(,에스)
에프(,에스)=에프(에스,0)에스
영형V(엑스(에스),엑스())=이자형나는(에스)엑스μ(엑스)

직관은 신호에서 일부 유한 한 관측치 집합에 대해 추정 된 자기 상관 행렬이 무조건 행렬처럼 행동한다는 것입니다. 상호 관계는 절대 위치가 아니라 시간 차이에만 의존하기 때문입니다. 연산자). 이것에 대한 많은 증거가 있으며 이것은 스케치 직관입니다.

정현파로 대각선 화되는 자기 상관 함수 세트는 고정 프로세스에 해당하는 함수이지만, 많은 다른 프로세스의 자기 상관 함수는 일정한 간격 동안 정현파로 대략 대각선 화됩니다. 이러한 프로세스는 일정 기간 동안 고정 프로세스로 근사 할 수있는 프로세스에 해당합니다. 자세한 내용은 여기에 있습니다 .

일반적인 비 정적 프로세스는 정현파에 의해 대각선 화 될 필요가없는 자기 상관 함수를 가질 수 있습니다.

국부적으로 정지 된 공정은 스펙트럼이 느리게 변화하거나 스펙트럼에 약간의 간격이 급격한 변화가있을 것입니다. 말, 동물 소리, 음악 및 기타 많은 자연 소리가이 설명에 적합합니다. 내가 아는 것처럼 부분 공간 식별 알고리즘이 작동하는 이유는 어떤 형태의 로컬 고 정성 (엄격하지 않은)이 일반적으로 우리가 분석하는 신호 유형을 유지하기 때문입니다.


μ

@MarkS 감사합니다. 1) 이것에 근거하여 공분산 행렬의 고유 벡터가 사인 곡선만큼 프로세스가 정지되어 있다고 말할 수 있습니까? 이것이 일종의 문구 성의 척도 가 될 수 있습니까 ? 2) "... 그리고 순환 행렬은 고유 벡터로 이산 정현파를 가지고 있습니다. (그들은 컨볼 루션 연산자이기 때문에 ...") 이것이 무슨 의미인지, 어떤 연산자인지 확실하지 않습니다. 명확하게 말씀해 주시겠습니까? 3) "자가 상관 함수 집합"이라고 말할 때 공분산 행렬의 행에 대해 이야기하고 있습니까? 다시 감사합니다.
Spacey

@Mohammad Cheers : 1) 그렇습니다. 2) 순환 행렬은 벡터의 모든 순환 순열에서 형성되므로 순환 행렬에 다른 벡터를 곱하면 두 벡터 간의 컨볼 루션이됩니다. 3) 자기 상관 함수 Corr (s, t)는 임의의 프로세스 X에 대한 X (s)와 X (t) 사이의 자기 상관입니다. 나는 연속적이고 이산적인 경우를 동시에 처리하기를 원하기 때문에 함수라고 부릅니다. 샘플 자기 상관 행렬은이 함수에 대한 이산 근사치로 볼 수 있습니다.
Mark S

@ Wiener–Khinchin_theorem을 지적 해 주셔서 감사합니다. 저는 푸리에 분석을 그룹에 대해 처음 배웠으며 신호 처리 클래스에서 공식적으로 소개하지 않았습니다.
Mark S
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.