使用PyMeasure读取锁相放大器时出现异常结果的问题排查

问题描述 #

通过PyMeasure的SR830驱动，经GPIB连接电脑与锁相放大器，配置内部参考源、设置频率后尝试读取幅值（magnitude）和相位（phase）参数，代码如下：

from pymeasure.instruments.srs import SR830
from time import sleep

LI = SR830("GPIB address")

LI.reference_source="Internal"
LI.frequency=668.4
# LI.time_constant=4/LI.frequency
sleep(3*LI.time_constant)
LI.auto_gain()
    
sleep(10)
# 丢弃初始读取数据，避免锁相放大器的数据异常输出
throwaway=LI.magnitude
throwaway2=LI.magnitude

def lock_in_test(repeats=1):
    itr=0
    while itr<repeats:
        res=LI.magnitude
        sleep(3)
        phase=LI.phase
        print(f"mag={res},phase={phase}")
        sleep(1)
        itr+=1

lock_in_test(10)

设备屏幕显示预期结果：相位量级为0.1，幅值约7e-6，但实际输出如下：

mag=9.0,phase=7.67835e-06
mag=7.67788e-06,phase=7.67788e-06
mag=10.37,phase=7.70722e-06
mag=10.37,phase=7.61545e-06
mag=10.37,phase=7.73097e-06
mag=10.37,phase=7.68859e-06
mag=10.37,phase=7.67648e-06
mag=10.37,phase=7.69558e-06
mag=10.37,phase=7.65413e-06
mag=10.37,phase=7.68114e-06

存在两个核心问题：

• 幅值从某一时刻开始重复输出错误值
• 调用LI.phase时读取到的是幅值数据，注释该语句后幅值读取恢复正常
  已尝试降低GPIB总线速率，问题仍存在。

原因分析 #

这类问题通常源于PyMeasure SR830驱动的属性缓存机制，或连续读取单个参数时的GPIB指令时序冲突：

1. PyMeasure对仪器属性的读取默认会缓存结果，若驱动未正确处理连续读取的刷新逻辑，会导致后续读取到旧数据或错误的寄存器值
2. 分开读取magnitude和phase时，两次GPIB指令的间隔可能触发仪器内部的数据输出队列异常，导致返回数据串位

解决方法 #

1. 一次性读取多参数（推荐） #

SR830支持通过read_measurement方法一次性读取所有测量参数，避免分开读取的时序冲突。修改测试函数如下：

def lock_in_test(repeats=1):
    itr=0
    while itr<repeats:
        # 一次性读取幅值(X/Y的模)和相位，返回格式为 [X, Y, R, theta]
        X, Y, mag, phase = LI.read_measurement()
        print(f"mag={mag},phase={phase}")
        sleep(4)  # 合并原有的sleep(3)+sleep(1)
        itr+=1

2. 手动发送GPIB指令绕过驱动缓存 #

直接使用底层GPIB指令读取参数，避免驱动的缓存干扰：

def lock_in_test(repeats=1):
    itr=0
    while itr<repeats:
        # 发送读取幅值指令（R参数对应指令'OUTP? 3'）
        LI.write("OUTP? 3")
        mag = float(LI.read())
        sleep(0.1)  # 短间隔确保指令执行完成
        # 发送读取相位指令（theta参数对应指令'OUTP? 4'）
        LI.write("OUTP? 4")
        phase = float(LI.read())
        print(f"mag={mag},phase={phase}")
        sleep(4)
        itr+=1

3. 禁用属性缓存 #

若坚持使用属性读取，可手动刷新缓存或禁用驱动的缓存机制：

# 在读取前清除缓存
def lock_in_test(repeats=1):
    itr=0
    while itr<repeats:
        # 强制刷新幅值缓存
        LI.magnitude.refresh()
        res=LI.magnitude
        sleep(3)
        # 强制刷新相位缓存
        LI.phase.refresh()
        phase=LI.phase
        print(f"mag={res},phase={phase}")
        sleep(1)
        itr+=1

4. 优化仪器配置与读取时序 #

• 确保time_constant配置合理，若注释了手动设置，需确认仪器自动设置的时间常数是否匹配测量需求，避免数据未稳定就读取
• 减少不必要的sleep时间，改用等待仪器就绪的逻辑（如检查仪器状态寄存器）

内容的提问来源于stack exchange，提问作者Nitzan R