[Pybrainsim-activity] SF.net SVN: pybrainsim:[141] trunk/src/PyBrainSim.py

[<rg...@us...>]

Revision: 141
          http://pybrainsim.svn.sourceforge.net/pybrainsim/?rev=141&view=rev
Author:   rgoj
Date:     2010-06-27 16:41:28 +0000 (Sun, 27 Jun 2010)

Log Message:
-----------
* Changes to dipole source example:
** No noise at level of electrode
** Easy to specify number of evenly spaced electrodes

Modified Paths:
--------------
    trunk/src/PyBrainSim.py

Modified: trunk/src/PyBrainSim.py
===================================================================
--- trunk/src/PyBrainSim.py	2010-06-27 16:33:45 UTC (rev 140)
+++ trunk/src/PyBrainSim.py	2010-06-27 16:41:28 UTC (rev 141)
@@ -290,24 +290,16 @@
         head.setSamplingFrequency(1)
         
         # Adding registration sites
-        head.addRegistrationSite([-numpy.pi/2, numpy.pi/2])
-        head.addRegistrationSite([-numpy.pi/4, numpy.pi/2])
-        head.addRegistrationSite([0, numpy.pi/2])
-        head.addRegistrationSite([numpy.pi/4, numpy.pi/2])
-        head.addRegistrationSite([numpy.pi/2, numpy.pi/2])
+        nElectrodes = 20
+        for i in range(nElectrodes):
+            head.addRegistrationSite([i*numpy.pi/(nElectrodes-1),numpy.pi/2])
         
         # Adding generators
         generator = GeneratorNoisy('Gen', head, position = [ 0.3, 0, 0, 0, 0], mean=1, stddev=0.1)
-        generator = GeneratorNoisy('Gen', head, position = [ 0.3, 0, 0, 0, 0], mean=1, stddev=0.1)
-        generator = GeneratorNoisy('Gen', head, position = [ 0.3, 0, 0, 0, 0], mean=1, stddev=0.1)
 
         # Run the simulation just once (or, equivalently for one second with the sampling rate of 1 Hz)
-        simulatedData = numpy.array(head.runSimulation(5))
+        simulatedData = numpy.array(head.runSimulation(1))
         
-        # Add Gaussian noise with mean 0 and variance 1
-        electrodeNoise = numpy.random.normal(0, 0.0000001, simulatedData.shape);
-        simulatedData = simulatedData + electrodeNoise
-        
         print(simulatedData)
         
     


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