blob: 26cc35b02f27f392a08fc3d5dec22bc9ec379595 (
plain)
1
2
3
4
5
6
7
8
9
|
function [ret, removed] = remove_sevenday_frequency(quantity)
A = fft(sum(quantity'));
kickout = (abs(A) > 2500) & (abs(A) < 20000);
f = fft(quantity);
removed = diag(kickout) * f;
ret = real(ifft(diag(~kickout) * f));
end
|