From 5ea27d71b3a5d365bebdd56fd74b372886bf37b1 Mon Sep 17 00:00:00 2001 From: Francis Russell Date: Mon, 3 Dec 2012 19:14:44 +0000 Subject: [PATCH] Add experiment with squaring single frequency signals. --- misc/resampling.py | 9 +++++++++ 1 file changed, 9 insertions(+) diff --git a/misc/resampling.py b/misc/resampling.py index 40e7701..e3eb49c 100755 --- 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)) -- 2.47.3