Python 如何使用ADC在Raspbery Pi中获得尽可能高的采样率?

Python 如何使用ADC在Raspbery Pi中获得尽可能高的采样率?,python,raspberry-pi,pyqtgraph,adc,Python,Raspberry Pi,Pyqtgraph,Adc,我在一个使用Raspberry Pi 3 B的项目中工作,我通过ADC MPC3008从红外传感器(Sharp GP2Y0A21YK0F)获取数据,并使用PyQtgraph库实时显示数据 ADC的数据表显示,在5.0V电压下,采样率为200khz。然而,我每秒只得到大约30个样本 是否可以使用树莓pi实现200khz 如果是,我应该学习或实施什么才能获得它 如果没有,我应该做什么才能获得可能的最高采样率,我如何才能找出什么是最高的采样率 这是我的密码: # -*- coding: utf-8 -

我在一个使用Raspberry Pi 3 B的项目中工作,我通过ADC MPC3008从红外传感器(Sharp GP2Y0A21YK0F)获取数据,并使用PyQtgraph库实时显示数据

ADC的数据表显示,在5.0V电压下,采样率为200khz。然而,我每秒只得到大约30个样本

是否可以使用树莓pi实现200khz

如果是,我应该学习或实施什么才能获得它

如果没有,我应该做什么才能获得可能的最高采样率,我如何才能找出什么是最高的采样率

这是我的密码:

# -*- coding: utf-8 -*-

import time
import Adafruit_GPIO.SPI as SPI
import Adafruit_MCP3008
from collections import deque
import serial
import pyqtgraph as pg
from pyqtgraph.Qt import QtCore, QtGui
import numpy as np

SPI_PORT   = 0
SPI_DEVICE = 0
mcp = Adafruit_MCP3008.MCP3008(spi=SPI.SpiDev(SPI_PORT, SPI_DEVICE))

win = pg.GraphicsWindow()
win.setWindowTitle('pyqtgraph example: Scrolling Plots')

nsamples=600 #tamanho das matrizes para os dados
tx_aq = 0 #velocidade da aquisição
intervalo_sp = 0.5 #intervalo para secao de poincare

# 1) Simplest approach -- update data in the array such that plot appears to scroll
#    In these examples, the array size is fixed.
p1 = win.addPlot()
p1.setRange(yRange=[0,35])

p2 = win.addPlot()
p2.setRange(yRange=[-100,100])

p3 = win.addPlot()
p3.setRange(yRange=[-100,100])
p3.setRange(xRange=[-0,35])

#p3.plot(np.random.normal(size=100), pen=(200,200,200), symbolBrush=(255,0,0), symbolPen='w')
'''
p3.setDownsampling(mode='peak')
p3.setClipToView(True)
p3.setRange(xRange=[-100, 0])
p3.setLimits(xMax=0)
'''

data1= np.zeros((nsamples,2),float) #ARMAZENAR POSICAO
vec_0=deque()
vec_1=deque()
vec_2=deque()
ptr1 = 0

data2= np.zeros((nsamples,2),float) #ARMAZENAR VELOCIDADE
diff=np.zeros((2,2),float)
diff_v=deque()

data3= np.zeros((nsamples,2),float)
data3_sp=np.zeros((1,2),float)

ptr3=0

curve1 = p1.plot(data1)
curve2 = p2.plot(data2)
curve3 = p3.plot(data3)

#Coeficientes da calibração do IR
c1=-7.246
c2=44.17
c3=-95.88
c4=85.28

tlast=time.clock()
tlast_sp=time.clock()
#print tlast

def getdata():
    global vec_0, vec_1, vec_2, tlast
    timenow=time.clock()

    if timenow-tlast>=tx_aq:
        #name=input("HUGO")

        tlast=timenow

        t0=float(time.clock())
        str_0 =mcp.read_adc(0)
        t1=float(time.clock()) 
        str_1 =mcp.read_adc(0)
        t2=float(time.clock())
        str_2 =mcp.read_adc(0)

        d0x=(float(str_0))*(3.3/1023)
        d0= c1*d0x**3+c2*d0x**2+c3*d0x+c4
        vec_0=(t0, d0)

        d1x=(float(str_1))*(3.3/1023)
        d1= c1*d1x**3+c2*d1x**2+c3*d1x+c4
        vec_1=(t1, d1)

        d2x=(float(str_2))*(3.3/1023)
        d2= c1*d2x**3+c2*d2x**2+c3*d2x+c4
        vec_2=(t2, d2)

        functions()

def diferenciar():
    global data2


    diff=(data1[-1,1]-data1[-3,1])/(data1[-1,0]-data1[-3,0])

    data2[:-1] = data2[1:]
    data2[-1,1] = diff
    data2[-1,0] = data1[-2,0]


def organizar():
    global data1, data3

    data1[:-1] = data1[1:]
    vec_x1=np.array(vec_1)
    data1[-1]=vec_x1

def EF(): #ESPACO DE FASE
    global data3, ptr3

    data3[:-1] = data3[1:]
    data3[-1,0]=data1[-1,1]
    data3[-1,1]=data2[-1,1]

def SP():
    global timenow_sp, tlast_sp

    timenow_sp=time.clock()

    if timenow_sp-tlast_sp>=intervalo_sp:

        tlast_sp=timenow_sp

        data3_sp[0,0]=data3[-2,0]
        data3_sp[0,1]=data3[-2,1]
        p3.plot(data3_sp, pen=None, symbol='o', symbolPen=None, symbolSize=4, symbolBrush=('r'))
        #print data3_sp

def plotar():
    global ptr1
    curve1.setData(data1)    
    ptr1 += 1
    curve2.setData(data2)
    #curve2.setPos(ptr1, 0)

    #p3.plot(data3)

def functions():

    diferenciar()
    organizar()
    EF()
    SP()
    plotar()

def update1():
    global data1, curve1, ptr1

    getdata()


# update all plots
def update():
    update1()

timer = pg.QtCore.QTimer()
timer.timeout.connect(update)
timer.start(50)



## Start Qt event loop unless running in interactive mode or using pyside.
if __name__ == '__main__':
    import sys
    if (sys.flags.interactive != 1) or not hasattr(QtCore, 'PYQT_VERSION'):
        QtGui.QApplication.instance().exec_()
我正试图找到一个解决办法,但到目前为止我失败了


你们能帮我一下吗,或者至少告诉我在哪里可以找到这方面的信息?

这种采样率在像Raspberry Pi这样的通用计算机上是无法实现的,尤其是在
MCP3008
上。原因是MCP系列ADC在
~2.7Mhz
SPI时钟在
5V
时达到峰值

为了以
200KHz
的速率进行读取,您需要一个专用板

但是,您可以尝试使用高达96 kHz的PCM1803A

96kHz采样很容易通过I2S ADC实现。我有96kHz,24位 立体声输入工作使用一个简单的I2S编解码器上的突破板。 可能会有更高的采样率,但我使用的编解码器 (PCM1803A)在96kHz时最大

下面还讨论了这一点:

仅使用SPI ADC,Pi将无法达到150ksps。不 即使只有一个频道。我想我听说过最好的是50 Ksps 频率上会有一定的抖动 取样

2通道*150ksps=300ksps

对于开销,假设每个样本大约32位,您将看到 9.6mbps的原始数据

只有Pi和ADC是不行的

您需要一个外部微控制器/adc将数据发送到Pi 通过USB或以太网

以及

基本问题是:

  • Raspberry Pi不是为高速数据采集而设计的
  • MCP系列ADC的最高频率为2.7Mhz,SPI时钟为5V
  • RPi的SPI延迟
Pi上的SPI接口无法准确地 以精确的间隔从ADC读取100000个样本


你的回答很有启发性。非常感谢你!