Add experiment with squaring single frequency signals.

diff --git a/misc/resampling.py b/misc/resampling.py
index 40e77010657b017f68976e4d16b33da115fc0aa6..e3eb49c9ce623eaed67420b8f396d22a6e2298ca 100755 (executable)
--- a/misc/resampling.py
+++ b/misc/resampling.py
@@ -60,6 +60,9 @@ def checkUpsampling(signal):
 signal = array([1,3.5,-7,2.8,51,-3,-23,1000])
 signal2 = append(signal, [1])
 
+evenSignalFreqs = array([0., 0., 0., 0., 1., 0., 0., 0.])
+oddSignalFreqs = array([0., 0., 0., 0.5, 0.5, 0., 0.])
+
 checkUpsampling(signal)
 print
 checkUpsampling(signal2)
@@ -67,4 +70,10 @@ print "\n"
 compareIntegrals(signal)
 print "\n"
 compareIntegrals(signal2)
+print "\nSquaring signal with even number of samples..."
+print "Original signal (in frequency domain): " + str(evenSignalFreqs)
+print "Squared signal (in frequency domain): " + str(fft(size(evenSignalFreqs)*ifft(evenSignalFreqs)**2).round(5))
 
+print "\nSquaring signal with odd number of samples..."
+print "Original signal (in frequency domain): " + str(oddSignalFreqs)
+print "Squared signal (in frequency domain): " + str(fft(size(oddSignalFreqs)*ifft(oddSignalFreqs)**2).round(5))