matlab - matlab でのウェーブレット分解と再構成

Question

ウェーブレット ツールボックスを使用することにより、時系列を分解および再構築することが可能になります。

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

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

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

このことから、シリーズから特定の周波数を削除するにはどうすればよいでしょうか? たとえば、系列を周波数成分に分解した後、再構築する前に、日周期などの周波数を削除するにはどうすればよいですか?

24 単位サイクルの信号の例:

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」を最終シリーズからどのように削除しますか?

Answer 1

から取得したベクトルcとそれを変更する必要があります。特定の周波数を削除する場合は、それらの周波数に対応するデータを見つけてゼロにします。次に call を呼び出して、変更された時系列を再構築します。lwavedecwaverec