...
Pensaba usar la ecuación de Steinhart-Hart, pero efectivamente es para más de 30ºC, entonces tomaré más puntos, y luego con Matlab, o Excel creo que también se puede, aproximar a polinomios o exponenciales.
Otra opción, que al menos a mí me gusta mas. Supongo que al termiostor lo vas a disponer formando un divisor de tensión y leer esa tensión con un microcontrolador. Como ese valor de tensión te introduce otra alinealidad (es una función racional de la resistencia) , hacé directamente una tabla Temp-Tensión en lugar de Temp-Resistencia y la corrección hacerla solamente con polinomios evaluados por el método de Horner, que es más rápido y ocupa mucha menos memoria que usar logaritmos y pow() .
Por otro lado si el rango es amplio vas a tener el problema de la resolución, porque solamente en la franja alrededor de la temperatura central vas a tener mediciones óptimas, para temperaturas altas y bajas vas a tener todo amontonado alrededor de Vcc y gnd.
Una manera de zafar de eso es usar salidas open-collector e ir conmutando resistencias. Como la respuesta del divisor de tensión a NTC resulta bastante lineal dentro de +/-10° alrededor de la frecuencia central, solamente tenés que programar cambios cada 20°. O menos si por ejemplo aproximás por una cúbica.
Saludos.