2015-01-13 6 views
0

FFT에서 얻은 데이터에 곡선을 맞추고 싶습니다. 이 작업을 수행하는 동안 FFT가 binned 데이터를 제공한다는 것을 기억했기 때문에이를 커브 피팅으로 다르게 처리해야하는지 궁금했습니다.fft : 비닝 데이터 적용

빈이 구조에 비해 좁은 경우 데이터를 다르게 처리해야 할 필요는 없다고 생각하지만 사실 그렇지 않습니다.

binned 데이터에 맞는 올바른 방법은 bin 및 fit 값의 차이가 아니라 각 bin에 대해 bin 영역과 fitted curve 아래 영역을 최소화하여 각 bin의 에너지가 최소화되도록하는 것입니다. 커브로 표시되는 빈의 범위에있는 에너지와 일치합니다.

내 질문은 : 이것에 대해 정확하게 생각하고 있습니까? 그렇지 않다면 어떻게해야합니까?

또한이 주제에 대한 정보를 둘러 보았을 때 예를 들어 "최대 로그 가능성"이 발생했지만 내 상황에 어떻게 적용되는지 이해하기에 충분한 정보를 찾지 못했습니다.

추신 :이 질문에 맞는 사이트인지 알려지지 않은 경우 더 좋은 곳이 있으면 알려 주시기 바랍니다.

+0

안녕하세요!지금까지 시도한 것을 보여주십시오. 특히 fft 스펙트럼을 생성하는 코드 예제와이를 적용하고자하는 커브의 예제를 보는 것이 도움이됩니다. [이] (http://stackoverflow.com/help/mcve)를 참조하십시오. – gg349

+0

질문은 여전히 ​​이론적 인 단계에 있기 때문에 아마 http://dsp.stackexchange.com에서 더 나을 것입니다. –

+0

예제가 좋지 않을 것이라고 생각합니다. 도움이됩니다. 이론적 인 질문이기 때문에 : "나는 이것을 어떻게 얻을 수 있는가?"가 아니라 "어떤 방법으로 정답을 주느냐". 의 DSP, 물론, 내가 어떻게 그걸 간과 했나요? – Oebele

답변

0

옌스이 최대 로그 우도 추정을 사용하여 수행 할 수 있습니다. 이것은 측정 된 데이터를 산출 할 가능성이 가장 큰 매개 변수 집합을 찾는 방법입니다.이 기법은 통계에서 유래했습니다.

필자는 이것을 binned 데이터에 적용하는 방법에 대한 이해할 수있는 소스를 찾았습니다. 슬프게도 여기에 수식을 입력 할 수 없으므로 전체 소스에 대한 설명을 참조하십시오. 슬라이드 4는 this slide show입니다.

편집 :
잡음이 신호의 경우이 방법은 아주 잘 작동하지 않았다. 조금 더 견고한 방법은 최소 제곱합으로 영역 사이의 차이가 최소화되는 방법입니다.

이 방법을 방어 할만한 문헌을 찾지 못했지만 최대 로그 우도 추정에서 발생하는 것과 유사하며 노이즈가없는 테스트 케이스의 경우 매우 유사한 결과를 산출합니다.

0

unwindowed FFT의 경우 bin 간의 올바른 보간법은 Sinc (sin (x)/x) 또는주기적인 Sinc (Dirichlet) 보간 커널을 사용하는 것입니다. 대역 제한 신호의 샘플을 FFT하기 때문에 연속 스펙트럼을 재구성합니다.

+0

나는 보간을 찾고 있지는 않지만 실제 피크 모양의 적합을 찾고있다. – Oebele

+0

실제 스펙트럼 피크 모양은 FFT로 입력 된 데이터에서 사용할 수있는 정보가 주어지면 Sinc 재구성에 의해 생성됩니다. – hotpaw2

+0

나는 여기서 무엇을 의미하는지 표현하는 데 어려움이있다. 그러나 나는 그것을 시도해 보자 : 그것은 사실 일 수 없다. 예를 들어 가우시안 레이저 펄스를 측정하면, 스펙트럼은 sinc가 아닌 gaussian이되어야하기 때문이다. 따라서 sinc 행동은 DFT 또는 이산화의 '인공물'이어야합니다. DFT의 효과와 일치하기 때문에 sinc을 피팅하는 것은 물리적 관점에서 이해가되지 않습니다. 물리적 인 관점에서 예상되는 스펙트럼 모양을 설명하는 함수에 적합하길 원합니다. – Oebele

0

스펙트럼을 (FFT에서) 보간하는 매우 간단하고 효과적인 방법은 제로 패딩을 사용하는 것입니다. 그것은 FFT 이전에 windowing을 사용하거나 사용하지 않은 채로 작동합니다.

  1. (N을 길이 N의 사용자의 입력 벡터를 가지고 M은
  2. 길이의 FFT를 수행 0으로 원래의 N 값을 넘어 모든 값을 설정 정수
  3. 인 길이 M의 *의 N,로 확장 * M)을 OUPUT 쓰레기통

  • 계산 크기 당신이 얻을 것은 보간 스펙트럼이다.

    안부,

  • +0

    이 함수는 sinc 함수로 스펙트럼을 보간하지만 추가 정보는 제공하지 않습니다. 내 목표는 피크 보간의 정확한 위치를 찾는 것입니다. 보간법은이를 위해 도움이되지 않습니다. – Oebele

    +0

    제공 할 추가 정보가 없습니다. Sinc 재구성은 FFT가 주어진 샘플에서 사용할 수있는 모든 정보를 사용합니다 (밴드 제한 신호에서 나온 것으로 가정). – hotpaw2

    +0

    @Ale 추가 백색 가우시안 노이즈로 사인파의 주파수를 낮추고 싶다면 피사 렌코 (Pisarenko)와 아칸소 (AR)와 같은 방법을 찾아 보는 것이 좋습니다. 그것이 다른 것이라면 문제를 다시 설명하십시오. – Jens