matlab에서 웨이브 렛 분해 및 재구성

Question

웨이블릿 도구 상자를 사용하면 시계열을 분해하고 재구성 할 수 있습니다.

load leleccum; 
s = leleccum(1:3920); 

% decomposition 
[c,l] = wavedec(s,3,'db5'); 

% Reconstruction 
a0 = waverec(c,l,'db5'); 

이렇게하면 시리즈에서 특정 주파수를 어떻게 제거 할 수 있습니까? 예를 들어 시리즈를 주파수 구성 요소로 분해 한 후에 재구성하기 전에 일별주기와 같은 빈도를 제거하려면 어떻게해야합니까?

t = 1:365; 
raw = 20+(10-2).*rand(1,length(t)); 
signal_1 = 10*sin(2*pi*t/24); 

y = raw + signal_1; % example data series 

% decomposition 
[c,l] = wavedec(y,3,'db5'); 

어떻게 마지막 시리즈에서 24 즉 'signal_1'의 주기성에 신호를 제거하는 것입니다 : 24 단위주기 신호의

예?

Answer 1

wavedec에서 얻는 c 및 l 개의 벡터를 수정해야합니다. 특정 주파수를 제거하려면 해당 주파수에 해당하는 데이터를 찾아 제로 아웃하십시오. 그런 다음 waverec으로 전화하여 수정 된 시계열을 재구성하십시오.