이 번역 페이지는 최신 내용을 담고 있지 않습니다. 최신 내용을 영문으로 보려면 여기를 클릭하십시오.
연속 웨이블릿 변환으로부터 특정 주파수나 시간으로 국소화하여 복원
R2026a고베 지진 데이터를 특정 주파수로 국소화하여 근사 성분을 복원합니다. 데이터는 1Hz로 샘플링됩니다. CWT에서 [0.030, 0.070]Hz 주파수 범위에 있는 정보를 추출합니다.
load kobe
Fs = 1;데이터의 CWT를 구합니다.
[wt,f] = cwt(kobe,Fs);
변환된 데이터에 다시 신호 평균을 더하여 지진 데이터를 복원합니다.
xrec = icwt(wt,[],f,[0.030 0.070],SignalMean=mean(kobe));
원본 데이터와 [0.030, 0.070]Hz 주파수 범위에 있는 데이터를 플로팅하고 비교합니다.
t = (0:numel(kobe)-1)/60; tiledlayout(2,1) nexttile plot(t,kobe) grid on title("Original Data") nexttile plot(t,xrec) grid on xlabel("Time (mins)") title("Bandpass Filtered Reconstruction [0.030 0.070] Hz")
![Figure contains 2 axes objects. Axes object 1 with title Original Data contains an object of type line. Axes object 2 with title Bandpass Filtered Reconstruction [0.030 0.070] Hz, xlabel Time (mins) contains an object of type line.](cwtftdemo_01_ko_KR.png)
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 tiledlayout(2,1) nexttile plot(datayear,nino) grid on title("Original Data") nexttile plot(datayear,xrec) grid on title("El Nino Data - Years 2-8")
