diff options
| author | Benjamin Franzke <benjaminfranzke@googlemail.com> | 2012-05-08 09:22:54 +0200 |
|---|---|---|
| committer | Benjamin Franzke <benjaminfranzke@googlemail.com> | 2012-05-08 09:24:06 +0200 |
| commit | d5e1b113bcab2d22260446b87ef3251f7eeb7943 (patch) | |
| tree | aa35515ecbbdf9e7fdb0cd829f0f58d7f9187846 /remove_sevenday_frequency.m | |
| parent | 78e6538f5e37a7b1542ae8f4e72f7168bec09c67 (diff) | |
| download | dmc-d5e1b113bcab2d22260446b87ef3251f7eeb7943.tar.gz dmc-d5e1b113bcab2d22260446b87ef3251f7eeb7943.tar.bz2 dmc-d5e1b113bcab2d22260446b87ef3251f7eeb7943.zip | |
Add frequency removal to run_tests (bad result)
Diffstat (limited to 'remove_sevenday_frequency.m')
| -rw-r--r-- | remove_sevenday_frequency.m | 9 |
1 files changed, 9 insertions, 0 deletions
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 + |