plt.savefig('log_histograma.png', bbox_inches='tight')
def poisson(lambd, j):
return (lambd**j) * np.exp(-lambd) / scipy.misc.factorial(j)
def bose_einstein(lambd, j):
return np.power(lambd, j) / np.power(1 + lambd, 1 + j)
if __name__ == "__main__":
# Ejemplo de la utilizacion de estas funciones:
osci = Osciloscopio('USB0::0x0699::0x0363::C065087::INSTR')
# Adquirimos una pantalla y graficamos los datos:
adquirir_y_graficar(osci)
# Adquirimos una pantalla y guardamos los datos en el archivo
# "medicion0.csv" en la carpeta definida en "path":
path = r"D:\Alumnos\Grupo N\Conteo"
filename = "medicion0.csv"
adquirir_y_guardar(osci, path, filename)
# Adquirimos y guardamos 10 pantallas:
n = 10
adquirir_guardar_multiples(osci, path, n)
# Generamos cuentas y eventos a partir de mediciones ya hechas:
import visa
import matplotlib.pyplot as plt
import numpy as np
import datetime
import time
from instrumentos import Osciloscopio
osci = Osciloscopio('USB::0x0699::0x0363::C102223::INSTR') #Cambia con el osc
#%%
control = True
t0 = time.perf_counter()
TrefRub =
while control == True:
print('start loop')
tiempo1, data1 = osci.get_ventana(1)
tiempo2, data2 = osci.get_ventana(2)
ch1 = np.array([tiempo1, data1])
ch2 = np.array([tiempo2, data2])
#measurement_name = 'Tref'+str(datetime.datetime.now().hour)+'_'+str(datetime.datetime.now().minute)+'_'+str(datetime.datetime.now().second)
measurement_name = str(TrefRub+str(datetime.datetime.now().hour)+'_'+str(datetime.datetime.now().minute)+'_'+str(datetime.datetime.now().second)
print('plotting')
plt.figure(2)
plt.title('Voltaje vs Tiempo', size = 16)
plt.plot(tiempo1,data1,'r*-', markersize=1,label='CH1')
plt.plot(tiempo2,data2,'b*-', markersize=1, label='CH2')
plt.xlabel('Tiempo[ms]', size=12)
plt.ylabel('Voltaje[mV]', size=12)
plt.legend()
"""
# Este archivo tiene los codigos necesarios para la adquisición de los datos.
# import visa
import matplotlib.pyplot as plt
import numpy as np
import datetime
import time
from instrumentos import Osciloscopio
from control import ITC4001
import scipy as sp
from scipy import stats
# itc = ITC4001('USB::0x1313::0x804A::M00404166::INSTR')
osci = Osciloscopio('USB0::0x0699::0x0363::C065093::INSTR')
#%% Inicio y configuracion del ITC4001
import visa
rm = visa.ResourceManager()
tec = rm.open_resource('USB::0x1313::0x804A::M00404166::INSTR')
tec.write('CONF:TEMP')
def temp_set(tec, temp):
tec.write('SOUR2:TEMP {}C'.format(temp))
def temp_sp_read(tec):
temp_sp = float(tec.query('SOUR2:TEMP?'))
print(temp_sp)
return (temp_sp)
# -*- coding: utf-8 -*-
"""
Created on Tue Jun 18 10:18:14 2019
@author: Publico
"""
#%%
import numpy as np
import matplotlib.pyplot as plt
from instrumentos import Osciloscopio
osci = Osciloscopio('USB::0x0699::0x0363::C102220::INSTR')
#%% Inicio y configuracion del ITC4001
import visa
rm = visa.ResourceManager()
tec = rm.open_resource('USB::0x1313::0x804A::M00404166::INSTR')
tec.write('CONF:TEMP')
def temp_set(tec, temp):
tec.write('SOUR2:TEMP {}C'.format(temp))
def temp_sp_read(tec):
temp_sp = float(tec.query('SOUR2:TEMP?'))
print(temp_sp)
return(temp_sp)
def temp_read(tec):
temp = float(tec.query('read?'))
print(temp)
return(temp)