From d5e1b113bcab2d22260446b87ef3251f7eeb7943 Mon Sep 17 00:00:00 2001 From: Benjamin Franzke Date: Tue, 8 May 2012 09:22:54 +0200 Subject: Add frequency removal to run_tests (bad result) --- remove_sevenday_frequency.m | 9 +++++++++ 1 file changed, 9 insertions(+) create mode 100644 remove_sevenday_frequency.m (limited to 'remove_sevenday_frequency.m') diff --git a/remove_sevenday_frequency.m b/remove_sevenday_frequency.m new file mode 100644 index 0000000..e0233cc --- /dev/null +++ b/remove_sevenday_frequency.m @@ -0,0 +1,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 + -- cgit