Main Content

연속 웨이블릿 변환으로부터 특정 주파수나 시간으로 국소화하여 복원

라이브 스크립트 열기

고베 지진 데이터를 특정 주파수로 국소화하여 근사 성분을 복원합니다. CWT에서 [0.030, 0.070]Hz 주파수 범위에 있는 정보를 추출합니다.

load kobe

데이터의 CWT를 구합니다.

[wt,f] = cwt(kobe,1);

변환된 데이터에 다시 신호 평균을 더하여 지진 데이터를 복원합니다.

xrec = icwt(wt,[],f,[0.030 0.070],'SignalMean',mean(kobe));

원본 데이터와 [0.030, 0.070]Hz 주파수 범위에 있는 데이터를 플로팅하고 비교합니다.

subplot(2,1,1)
plot(kobe)
grid on
title('Original Data')
subplot(2,1,2)
plot(xrec)
grid on
title('Bandpass Filtered Reconstruction [0.030 0.070] Hz')

CWT에 주파수 대신 시간 주기를 사용할 수도 있습니다. 엘니뇨 데이터를 불러온 후 시간 주기를 연 단위로 지정하여 데이터의 CWT를 구합니다.

load ninoairdata
[cfs,period] = cwt(nino,years(1/12));

2년부터 8년까지의 역 CWT를 구합니다.

xrec = icwt(cfs,[],period,[years(2) years(8)]);

복원된 데이터의 CWT를 플로팅합니다. 2년부터 8년 사이의 대역을 벗어난 영역에는 에너지가 없는 것을 볼 수 있습니다.

figure
cwt(xrec,years(1/12))

원래 데이터와 2년부터 8년까지의 복원 데이터를 비교합니다.

figure
subplot(2,1,1)
plot(nino)
grid on
title('Original Data')
subplot(2,1,2)
plot(xrec)
grid on
title('El Nino Data - Years 2-8')