Mercurial > ~darius > hgwebdir.cgi > pyinst
changeset 30:9ce709b7da4b
Add Q&D example code for Anritsu MS2034A and Agilent R5071.
author | Daniel O'Connor <darius@dons.net.au> |
---|---|
date | Wed, 21 Sep 2011 14:58:55 +0930 |
parents | 12b6a3a0a481 |
children | c6c86dcb54ba |
files | agilent_r5071.py anritsu_ms2034a_example.py |
diffstat | 2 files changed, 95 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- /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)
--- /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) +