From f89a1f51993f79c3e95314b649f93f41d2a58bf7 Mon Sep 17 00:00:00 2001
From: bma <baptiste.marechal@femto-st.fr>
Date: Mon, 7 May 2018 14:15:34 +0200
Subject: [PATCH] change highpass frequency corner

---
 fftplot.m | 27 +++++++++++++++++++++++++++
 1 file changed, 27 insertions(+)
 create mode 100755 fftplot.m

diff --git a/fftplot.m b/fftplot.m
new file mode 100755
index 0000000..8f27f9a
--- /dev/null
+++ b/fftplot.m
@@ -0,0 +1,27 @@
+#!/usr/bin/octave-cli --persist
+
+filename = argv(){1};
+col = eval(argv(){2});
+mult = eval(argv(){3});
+data.fs = eval(argv(){4});
+
+[b, a] = butter (3, data.fs/1000, "high");
+
+if length(col) == length(mult)
+	figure
+	hold all
+	grid on
+	cc = 'bkcgmry';
+	for i = [1:length(col)]
+		data.x = filter(b, a, load(filename)(:,col(i)).*mult(i))(floor(end/4):end);
+		N = length(data.x);
+		f = (data.fs/N)*((-(N-1)/2):((N-1)/2));
+		X = abs(fftshift(ifft(data.x)));
+		plot(f, X, cc(mod(i, length(cc))))
+		leg{i} = strcat(filename, ' col', num2str(col(i)));
+		hold on
+	end
+	legend(leg)
+	input("Press to continue...");
+end
+exit
-- 
2.16.4