progress in analysis

ccl4/FFTanalisi.m

@@ -24,4 +24,9 @@ f=(0:1:length(x)-1)*1/(length(x)-1)*1/((x(2)-x(1))*3e-5);
 
 Rw=-imag(fft(y)./fft(ysI));
 figure(2)
-plot(f,Rw,f,Rw*0);xlim([0 2000]);ylim([-0.05 2])
+plot(f,Rw,f,Rw*0);xlim([0 2000]);ylim([-0.05 2])
+
+if 0, 
+    AA=[f' Rw]; str=['save fft.dat AA -ascii;'];eval(str);
+end;
+