X-Git-Url: http://git.osdn.jp/view?a=blobdiff_plain;f=win32gl%2Ftest%2Ffftfirm.cpp;fp=win32gl%2Ftest%2Ffftfirm.cpp;h=943b9513f809988ee2ae1c2b319f0411ce8f777b;hb=0d6f30dda998b99119681bf9015e1f34be9a4855;hp=4c03f04fb218db29a697f88829a4a794000aea84;hpb=213008ed37a546099e9daa990bf19885943b876b;p=psychlops%2Fcpp.git

diff --git a/win32gl/test/fftfirm.cpp b/win32gl/test/fftfirm.cpp
index 4c03f04..943b951 100644
--- a/win32gl/test/fftfirm.cpp
+++ b/win32gl/test/fftfirm.cpp
@@ -1,6 +1,45 @@
  #include <psychlops.h>
  using namespace Psychlops;
 
+
+void psychlops_main() {
+	Canvas cnvs(1024, 768, Canvas::window);
+
+	FFT1 fftworkspace;
+	Matrix source_mat, filtered_mat, mat;
+	source_mat.load("voltage_spec_vaio_ie.csv");
+	source_mat.transpose();
+
+	//fftworkspace.set(64,64);
+	//fftworkspace.makeNoise();
+	fftworkspace.set(source_mat);
+	fftworkspace.fft();
+	fftworkspace.getSpectrum(filtered_mat, 1.0);
+
+	int q = filtered_mat.getCols() / 4;
+	Interval col;
+	mat = filtered_mat(1, q*2+3<=col<=q*3);
+	mat.transpose();
+	mat.save("voltage_spec.csv");
+
+	double mean = 0.0;
+	for(int i=0; i<source_mat.getCols(); i++) {
+		mean += source_mat(1, i+1);
+	}
+	mean /= source_mat.getCols();
+	int counter = 0;
+	for(int i=1; i<source_mat.getCols()-3; i++) {
+		if(    source_mat(1, i+1)<mean
+			&& source_mat(1, i+2)<mean
+			&& source_mat(1, i+3)>mean
+			&& source_mat(1, i+4)>mean
+		) { counter++; }
+	}
+	std::cout << counter << std::endl;
+}
+
+
+/*
 const int S = 128;
 
  void psychlops_main() {
@@ -42,3 +81,4 @@ const int S = 128;
 
  }
 
+*/