안녕하세요.

fft를 IDL로 처음 시도하고 있는데요, 관련해서 질문이 있습니다.

old IDL user group에서 글을 읽다가

홍성재님이 작성하신 글(2009-02-03 23:03:12, Hit : 456, Vote : 0)과 답변을 참고로 제가 fft를 다시 해보았는데요, 질문이 있습니다.

===========================================================================================

X = FINDGEN(96)*!dtor
T1 = 24*SIN(X*24)+12 ; 24시간 주기, 24의 진폭을 가진 sine 함수
T2 = 12*SIN(X*48)+12 ; 12시간 주기, 12의 진폭을 가진 sine 함수
T3 = 6*SIN(X*96)+12 ; 6시간 주기, 6의 진폭을 가진 sine 함수
T = T1+T2+T3 ; T라는 T1, T2, T3의 합으로 이루어진 새 함수를 만들었습니다.

DN = 96 ; time series의 갯수
lag = FINDGEN(DN/2)+1
a_cor = A_CORRELATE(T, lag) ; 1~48 lag의 autocorrelation을 구합니다.

F = fft(T)
F=F^2
power = ABS(F) ; 파워를 구했습니다.

DEVICE, decomposed=0
;WINDOW, 2, xsize=1200, ysize=600

PLOT, lag, power, $ ; 1~48 time step에 대한 파워스펙트럼을 plot
YTITLE=’Power Spectrum of T’, BACKGROUND=255, COLOR=0, $
XTITLE=’Cycles / Hour’, /YLOG

END

====================================================

이상우 (2009-02-05 19:56:15)

FFT를 T에 직접 적용하고 파워스펙트럼을 구해야 원하시는 결과가 나옵니다. 보니까 a_correlate를 적용한 결과에 FFT가 적용되어 있던데, 이러면 T자체에 대한 결과는 나오지 않겠죠. 그리고 파워스펙트럼은 FFT결과를 제곱해야 나오는 것으로 압니다.

======================================================================================

여기서 질문입니다.

결과 그림이 x축(cycle/hour)에서 6,12,24 부분에 피크가 존재해야 하고 24부분에 가장 큰 피크가 존재해야 하는데 그렇지 않은 거 같습니다. 그리고 살짝 오른쪽으로 쉬프트가 된거 같습니다. 그리고 앞 부분에 가장 큰 피크 값은 무엇을 의미하는지 궁금합니다.

이상우님께서 홍성재님 질문에 대한 답변대로 FFT를 T에 대해서 했고, FFT 결과를 제곱했습니다.

• 이 게시글은 willms에 의해 10 years, 2 months 전에 수정됐습니다.