SF.net SVN: matplotlib: [3970] trunk/py4science/examples/noisy_sine.py

[<jd...@us...>]

Revision: 3970
          http://matplotlib.svn.sourceforge.net/matplotlib/?rev=3970&view=rev
Author:   jdh2358
Date:     2007-10-20 10:56:55 -0700 (Sat, 20 Oct 2007)

Log Message:
-----------
added files workbook

Added Paths:
-----------
    trunk/py4science/examples/noisy_sine.py

Added: trunk/py4science/examples/noisy_sine.py
===================================================================
--- trunk/py4science/examples/noisy_sine.py	                        (rev 0)
+++ trunk/py4science/examples/noisy_sine.py	2007-10-20 17:56:55 UTC (rev 3970)
@@ -0,0 +1,25 @@
+from scipy import arange, sin, pi, randn, zeros
+import pylab as p
+
+a = 2       # 2 volt amplitude
+f = 10      # 10 Hz frequency
+sigma = 0.5 # 0.5 volt standard deviation noise
+
+# create the t and v and store them a 2D array X 
+t = arange(0.0, 2.0, 0.02)                # an evenly sampled time array
+v = a*sin(2*f*pi*t) + sigma*randn(len(t)) # a noisy sine wave
+X = zeros((len(t),2))                     # an empty output array
+X[:,0] = t                                # add t to the first column 
+X[:,1] = v                                # add s to the 2nd column
+p.save('data/noisy_sine.dat', X)            # save the output file as ASCII
+
+# plot the arrays t vs v and label the x-axis, y-axis and title
+# save the output figure as noisy_sine.png
+p.plot(t, v, 'b-')
+p.xlabel('time (s)')
+p.ylabel('volts (V)')
+p.title('A noisy sine wave')
+p.grid()
+p.savefig('noisy_sine.png', dpi=150)
+p.savefig('noisy_sine.eps')
+p.show()


This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.