# HG changeset patch
# User Daniel O'Connor <darius@dons.net.au>
# Date 1316582935 -34200
# Node ID 9ce709b7da4bf6305aa6a5d1e184e761c97aeb96
# Parent  12b6a3a0a481f110d677c13b73144e66a56ea503
Add Q&D example code for Anritsu MS2034A and Agilent R5071.

diff -r 12b6a3a0a481 -r 9ce709b7da4b agilent_r5071.py
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/agilent_r5071.py	Wed Sep 21 14:58:55 2011 +0930
@@ -0,0 +1,39 @@
+# Docs are at http://ena.tm.agilent.com/e5071c/manuals/webhelp/eng/
+
+import scpi
+import pylab
+
+# Connect
+inst = scpi.instURL('vxi://203.31.81.47')
+
+# Read ID
+inst.write('*IDN?')
+inst.read()
+# 'Agilent Technologies,E5071C,MY46109815,A.09.54\n'
+
+# Set to 8 byte floats
+inst.write(":FORM:DATA REAL")
+
+# Set to little endian
+inst.write(':FORM:BORD SWAP')
+
+# Grab frequency data
+inst.write(":SENS1:FREQ:DATA?")
+freqs = scpi.bindecode(inst.read(), dtype = numpy.float64)
+
+# Grab trace data
+inst.write(":CALC1:TRACE1:DATA:FDATA?")
+dat = scpi.bindecode(inst.read(), dtype = numpy.float64)
+
+# We only want the real part (no imag for this measurement)
+dat = dat[::2]
+
+# Save for later
+numpy.savez('GS_preamp_20110929', freqs, dat)
+
+# Load..
+f = numpy.load('GS_preamp_20110929.npz')
+freqs = f['arr_0']
+dat = f['arr_1']
+
+pylab.plot(freqs.dat)
diff -r 12b6a3a0a481 -r 9ce709b7da4b anritsu_ms2034a_example.py
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/anritsu_ms2034a_example.py	Wed Sep 21 14:58:55 2011 +0930
@@ -0,0 +1,56 @@
+#!/usr/bin/env python
+
+# Copyright (c) 2011
+#      Daniel O'Connor <darius@dons.net.au>.  All rights reserved.
+#
+# Redistribution and use in source and binary forms, with or without
+# modification, are permitted provided that the following conditions
+# are met:
+# 1. Redistributions of source code must retain the above copyright
+#    notice, this list of conditions and the following disclaimer.
+# 2. Redistributions in binary form must reproduce the above copyright
+#    notice, this list of conditions and the following disclaimer in the
+#    documentation and/or other materials provided with the distribution.
+#
+# THIS SOFTWARE IS PROVIDED BY AUTHOR AND CONTRIBUTORS ``AS IS'' AND
+# ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+# IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+# ARE DISCLAIMED.  IN NO EVENT SHALL AUTHOR OR CONTRIBUTORS BE LIABLE
+# FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+# DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
+# OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
+# HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
+# LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
+# OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
+# SUCH DAMAGE.
+#
+
+import numpy
+from matplotlib import pylab
+import vxi
+import scpi
+
+def test(r):
+    # ID instrument
+    r.write('*IDN?')
+    print "ID is " + r.read(5)
+
+    # XXX: Don't do this, on this device it resets the ethernet settings (duuuh)
+    # Reset to defaults
+    #r.write("*RST")
+
+    # Grab trace data
+    r.write("TRAC1? TRACE1")
+    data = r.read(10)
+    #print "Data - " + dat
+    #ary = ascdecode(data)
+    ary = bindecode(data)
+
+    # Plot data
+    pylab.plot(ary)
+    pylab.show()
+    
+if __name__ == '__main__':
+    r = vxi.VXIDevice('192.168.42.51', device='inst0')
+    test(r)
+