别再手动抄数据了！用Python pyvisa库自动读取示波器/万用表数据（附完整代码）-CSDN博客

用Python pyvisa实现实验室仪器自动化：从示波器到万用表的高效数据采集

实验室里最耗时的往往不是创新性工作，而是那些重复性的数据记录任务。想象一下这样的场景：你正在进行电路板测试，需要每隔5分钟记录一次示波器的波形参数和万用表的电压读数，这样的工作可能要持续数小时。传统的手动记录不仅效率低下，还容易引入人为误差。这就是为什么我们需要自动化数据采集——而Python的pyvisa库正是解决这个问题的利器。

1. 搭建自动化测试环境

1.1 硬件准备与连接

在开始编写代码前，我们需要确保硬件连接正确。现代测试仪器通常支持多种接口：

USB连接 ：最简单直接的连接方式，即插即用
LAN接口 ：适合需要远程访问或多设备组网的场景
GPIB接口 ：传统仪器常用的总线接口，稳定性高

以泰克DPO3000系列示波器为例，通过USB连接时，VISA资源字符串通常类似于：

'USB0::0x0699::0x0401::C012345::INSTR'

提示：使用 pyvisa-shell 工具可以快速检测已连接的仪器，避免手动输入资源字符串的麻烦

1.2 软件环境配置

除了安装pyvisa核心库，我们还需要一些辅助工具包：

pip install pyvisa pyvisa-py numpy pandas matplotlib

这里特别推荐安装 pyvisa-py 作为后端，它是一个纯Python实现的VISA库，在跨平台兼容性上表现优异。对于Windows用户，也可以选择NI-VISA或Keysight VISA等商业实现。

配置完成后，用以下代码测试环境是否就绪：

import pyvisa
rm = pyvisa.ResourceManager()
print(rm.list_resources())  # 列出所有可用的仪器资源

2. 仪器控制基础：SCPI命令实战

2.1 理解SCPI命令体系

SCPI(Standard Commands for Programmable Instruments)是控制测试仪器的通用语言。它的命令结构通常遵循以下模式：

:子系统:操作 参数

例如，设置示波器垂直量程的命令可能是：

instrument.write(":CHANNEL1:SCALE 0.1")  # 设置通道1垂直刻度为100mV/div

常见仪器操作命令对照表：

仪器类型	常用SCPI命令	功能描述
示波器	`:AUTOSCALE`	自动调整波形显示
万用表	`:MEAS:VOLT:DC?`	测量直流电压
电源	`:SOURCE:VOLTAGE 3.3`	设置输出电压为3.3V

2.2 高效查询技巧

query() 方法是仪器控制中最常用的方法之一，它结合了写入命令和读取响应的操作。但实际使用中有几个优化技巧：

设置合理超时 ：对于可能耗时的操作，适当增加超时时间

freq = instrument.query(":MEAS:FREQ?", timeout=5000)  # 5秒超时

二进制数据传输 ：波形数据通常较大，二进制格式比ASCII更高效

instrument.write(":WAVEFORM:FORMAT WORD")  # 设置二进制格式
raw_data = instrument.query_binary_values(":WAVEFORM:DATA?", datatype='h')

3. 构建自动化数据采集系统

3.1 定时采集与存储实现

下面是一个完整的定时采集示波器参数并保存到CSV的示例：

import time
import csv
from datetime import datetime

def setup_scope(instr):
    instr.write(":STOP")  # 停止采集
    instr.write(":TIMEBASE:MODE MAIN")
    instr.write(":ACQUIRE:TYPE NORMAL")
    instr.write(":ACQUIRE:COUNT 16")  # 16次平均
    instr.write(":CHANNEL1:COUPLING DC")
    instr.write(":RUN")  # 开始采集

def measure_scope_params(instr):
    params = {}
    params['time'] = datetime.now().isoformat()
    params['freq'] = float(instr.query(":MEASURE:FREQUENCY? CHANNEL1"))
    params['vpp'] = float(instr.query(":MEASURE:VPP? CHANNEL1"))
    params['vavg'] = float(instr.query(":MEASURE:VAVERAGE? CHANNEL1"))
    return params

rm = pyvisa.ResourceManager()
scope = rm.open_resource('USB0::0x0699::0x0401::C012345::INSTR')
setup_scope(scope)

with open('measurements.csv', 'a', newline='') as f:
    writer = csv.DictWriter(f, fieldnames=['time', 'freq', 'vpp', 'vavg'])
    writer.writeheader()
    
    for _ in range(60):  # 每小时采集60次(每分钟一次)
        data = measure_scope_params(scope)
        writer.writerow(data)
        time.sleep(60)  # 等待1分钟

3.2 多仪器协同工作

在实际测试中，我们经常需要同时控制多种仪器。下面的示例展示了如何同步采集示波器和万用表数据：

def setup_dmm(instr):
    instr.write(":CONFIGURE:VOLTAGE:DC")
    instr.write(":VOLTAGE:DC:RANGE:AUTO ON")
    instr.write(":SAMPLE:COUNT 10")  # 每个读数取10次平均

rm = pyvisa.ResourceManager()
scope = rm.open_resource('USB0::0x0699::0x0401::C012345::INSTR')
dmm = rm.open_resource('USB0::0x0957::0x0618::MY543210::INSTR')

setup_scope(scope)
setup_dmm(dmm)

while True:
    scope.write(":SINGLE")  # 单次触发
    while int(scope.query(":ACQUIRE:COMPLETE?")) == 0:
        time.sleep(0.1)
    
    scope_data = measure_scope_params(scope)
    dmm_voltage = float(dmm.query(":MEASURE:VOLTAGE:DC?"))
    
    print(f"时间: {scope_data['time']}, 频率: {scope_data['freq']:.2f}Hz, "
          f"电压: {dmm_voltage:.3f}V")

4. 高级技巧与故障排除

4.1 波形数据的处理与分析

获取原始波形数据只是第一步，如何有效处理这些数据才是关键。以下是一个完整的波形采集和处理流程：

import numpy as np
import matplotlib.pyplot as plt

def get_waveform(instr, channel=1):
    instr.write(f":WAVEFORM:SOURCE CHANNEL{channel}")
    instr.write(":WAVEFORM:FORMAT WORD")  # 16位有符号整数
    instr.write(":WAVEFORM:BYTEORDER LSBFirst")
    
    # 获取波形参数
    xinc = float(instr.query(":WAVEFORM:XINCREMENT?"))
    xorig = float(instr.query(":WAVEFORM:XORIGIN?"))
    yinc = float(instr.query(":WAVEFORM:YINCREMENT?"))
    yorig = float(instr.query(":WAVEFORM:YORIGIN?"))
    
    # 获取波形数据
    raw_data = instr.query_binary_values(":WAVEFORM:DATA?", 
                                       datatype='h', 
                                       container=np.array)
    
    # 转换为实际电压值
    voltages = (raw_data * yinc) + yorig
    times = np.arange(0, len(voltages)) * xinc + xorig
    
    return times, voltages

times, voltages = get_waveform(scope)
plt.plot(times, voltages)
plt.xlabel('时间 (s)')
plt.ylabel('电压 (V)')
plt.title('示波器波形')
plt.grid(True)
plt.show()

4.2 常见问题解决方案

在实际使用中，你可能会遇到以下典型问题：

连接超时
- 检查仪器IP地址或VISA资源字符串是否正确
- 尝试增加超时时间： instrument.timeout = 10000 (10秒)
数据格式错误
- 明确指定数据格式： query_binary_values(..., datatype='f') (32位浮点)
- 对于ASCII数据，注意去除换行符： response.strip()
仪器无响应
- 检查仪器是否被其他程序占用
- 尝试重置仪器： instrument.write("*RST")
性能优化
- 对于高速采集，禁用自动校准： instrument.write(":CALIBRATION:AUTO OFF")
- 减少查询次数，使用复合命令： instrument.query(":MEASURE:FREQUENCY?;VPP?")