remove_sevenday_frequency.m


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(1-kickout) * f));
end