Files
HendiControlFirmware/Control/brewpi/max31865.py
T
2019-04-15 21:01:17 +00:00

110 lines
1.8 KiB
Python
Executable File

#!/usr/bin/python3
import spidev
from atemperatureSensor import ATemperatureSensor
class Max31865(ATemperatureSensor):
def __init__(self):
super(Max31865, self).__init__('Max31865')
# Open SPI bus
self.spi = spidev.SpiDev()
self.spi.open(0, 0)
self.spi.max_speed_hz = 100000
self.spi.mode = 0b01
self.write_reg(0x00, 0xA3)
self.temp_correction = -0.2
def read_reg(self, addr):
reg = self.spi.xfer([addr, 0xFF])
return reg[1]
def write_reg(self, addr, data):
self.spi.xfer([addr+0x80, data])
def read_digits(self):
self.write_reg(0x00, 0xA3)
msb = self.read_reg(0x01)
lsb = self.read_reg(0x02)
digits = float(256*msb + lsb)/2
return digits
def temperature(self):
return self.__tempRaw()
def __tempRaw(self):
digits = self.read_digits()
R = Max31865.calc_R(430, digits)
T = Max31865.vanDusen_temp(R) + self.temp_correction
# print("R={:.3f} Ohm, T={:.2f} degC".format(R, T))
return T
@staticmethod
def calc_R(Rref, digits):
k = digits / 8192 / 4
R = k*Rref
return R
@staticmethod
def vanDusenLut (R0, Tmin, Tmax, dT):
a = +3.90830e-03
b = -5.77500e-07
c = -4.18301e-12
Rv = []
Tv = []
T = Tmin
while T <= Tmax:
R = R0*(1 + a*T + b*T**2)
if T < 0.0:
R += R0*c*(T - 100)*T**3
Rv.append(R)
Tv.append(T)
T += dT
return Rv, Tv
@staticmethod
def vanDusen_temp(Rmeas):
RLut, TLut = Max31865.vanDusenLut(100, -150, +500, 1.0)
# Find candidate
i = 0
for R in RLut:
if R > Rmeas:
break
i += 1
# Linear interpolation
R0 = RLut[i-1]
R1 = RLut[i]
T0 = TLut[i-1]
T1 = TLut[i]
dR = R1 - R0
dT = T1 - T0
k = (Rmeas-R0)/dR
T = T0 + dT*k
return T
# Main
if __name__ == '__main__':
max31865 = Max31865()
for i in range (0, 8):
print("Reg[0x{}]: 0x{:02X}".format(i, max31865.read_reg(i)))
max31865.temperature()
print("End of program")