thermistor.sourceforge.netNTC Calculation

1 Abstract The project offers support for NTC thermistor calculations. The Steinhart-Hart equation is a mathematical model for these thermistors that seems to fit for a wide range of temperatures with high precision. Software to calculate the characteristic Steinhart-Hart coefficients based on temperature-resistance tables for given thermistors as well as functions allowing conversion of temperature values to resistance and vice versa is provided. 2 Description A model for the resistivity of a semiconductor as a function of the temperature was found by Steinhart and Hart 1968 ( [1] ). The Steinhart-Hart law describes the absolute temperature T (in Kelvins) as a function of the NTC thermistor’s resistivity (in Ω) according to the formula Steinhart-Hart polynom 1 / T = a 0 + a 1 · ln r + a 3 · (ln r) 3 The figure below shows the typical graph of an NTC thermistors characteristic, giving the reciprocal temperature (in 1 / K) over the natural logarithm of the resistance (in Ω). The constants a 0 , a 1 and a 3 , also called Steinhart-Hart coefficients, vary depending on the type of thermistor. To support developer when creating temperature measurement applications, thermistor manufacturer often supply these constants for their products. They also publicate tables where resistivity of thermistor products for a wider range of temperature values are listed. This project provides software to calculate temperature value for a given resistance of an NTC thermistor with given Steinhart-Hart coefficients, calculate resistance value for a given temperature for an NTC thermistor with given Steinhart-Hart coefficients and evaluate Steinhart-Hart coefficients for an NTC thermistor descibed by a temperature-resistance table. Apart from the standard Steinhart-Hart equation other forms have been found. For application with lower CPU power a simplified form of the Steinhart-Hart equation can be used. Simplified Steinhart-Hart polynom 1 ⁄ T = a 0 + a 1 · ln r On the other hand a quadratic term can be inserted into the formula to increase accuracy giving the extended Steinhart-Hart equation Extended Steinhart-Hart polynom 1 / T = a 0 + a 1 · ln r + a 2 · (ln r) 2 + a 3 · (ln r) 3 An introduction to thermistors and the Steinhart-Hart polynom can be found at Wikipedia [ 2 ]. 3 Software Software is provided by this project for the calculations given in section algorithms. The classes / modules can be used for conversion from resistance to temperature, conversion from temperature to resistance, calculation of Steinhart-Hart coefficients for an NTC thermistor whose characteristics is given as T-R table. The latter can be done for standard, simplified or extended Steinhart-Hart polynom. Download the archives for C or Java here . 3.1 Java-Classes For Java the following classes are available Class Description NtcException Exception used in this thermistor framework NtcTable Represents an T-R table for an NTC thermistor. Provides methods to get name and description of the NTC thermistor as well as methods to add entries and to lookup the resistance for a given temperature value. NtcThermistorModel Represents an NTC thermistor modelled using the Steinhart-Hart polynom 1 / T = a 0 + a 1 · ln r + a 3 · (ln r) 3 . The model will be created using a NtcTable object. It provider methods to read the themistor name, get the NtcTable object, get the Steinhart-Hart coefficients as well as methods for conversion from temperature to resistance and vice versa. NtcThermistorExtendedModel Represents an NTC thermistor modelled using the extended Steinhart-Hart polynom 1 / T = a 0 + a 1 · ln r + a 2 · (ln r) 2 + a 3 · (ln r) 3 . Derived from NtcThermistorModel thus providing the same methods. NtcThermistorSimplifiedModel Represents an NTC thermistor modelled using the simplified Steinhart-Hart polynom 1 ⁄ T = a 0 + a 1 · ln r . Derived from NtcThermistorModel, provides the same methods. Note : Generics are used in the project so Java 5 or higher is required according to that. A small test class NtcTest is part of the package together with an example text file containing an T-R table of a theoretical NTC (simu.txt). To make developers life easy the whole application is deployed as an (executable) jar file containing all sources, classes and a manifest file with main class NtcTest. After download the test application can be invoked by java -jar thermistor.0.1.jar The output produced can be seen here . 3.2 C Modules Three C modules are available Module Description coeff.c Program used to calculate Steinhart-Hart coefficients from an NTC table (given as text file). ttor.c T emperature to r esistance calculation. rtot.c R esistance to t emperature calculation. For usage under Windows they have been compiled (using Cygwin’s C compiler). For other OS please compile them with your system C-compiler. An example text file containing an T-R table of a theoretical NTC (simu.txt) is also provided for testing purposes. Example outout of the "coeff" program is given here . 4 Algorithms 4.1 Conversion from resistance to temperature Based on a measured resistance value of an NTC thermistor the extended Steinhart-Hart equations allows a simple calculation of the temperature 1 / T = a 0 + a 1 · ln r + a 2 · (ln r) 2 + a 3 · (ln r) 3 Here r is the resistance in Ω and T the absolute temperature in K (K = Kelvin). With absolute zero T abs = -273.15 ° C , the formula for the temperature t in °C finally leads to t = 273,15 ° C + [ a 0 + a 1 · ln r + a 2 · (ln r) 2 + a 3 · (ln r) 3 ] -1 By setting appropriate coefficients to zero, the calculation according to the simplified or standard Steinhart-Hart equation can be done as well. 4.2 Conversion from temperature to resistance To calculate resistance values in Ω from temperature in °C the root r of the extended Steinhart-Hart equation must be found (for the standard polynom simply set a 2 =0 in the formulas below) 1 / T = a 0 + a 1 · ln r + a 2 · (ln r) 2 + a 3 · (ln r) 3 with y = ln r , we get 1 / T = a 0 + a 1 · y + a 2 · y 2 + a 3 · y 3 We introduce the substitutions T = t + T abs b = a 2 / a 3 c = a 1 / a 3 d = (a 0 - 1 ⁄ T) / a 3 p = c - ⅓· b 2 q = 2 / 27 · b 3 - 1 / 3 · b · c + d u = [ -q/2 + ( q 2 /4 + p 3 /27 ) ½ ] ⅓ v = [ -q/2 - ( q 2 /4 + p 3 /27 ) ½ ] ⅓ getting r = e u + v - b/3 for the resistance r (in Ω) of a given temperature t (in ° C). 4.3 Calculation of Steinhart-Hart coefficients Sometimes calculation of coefficients is done using special temperature values. Inserting four value pairs in the range of interest into the extended Steinhart-Hart poylonm leads to a system of linear algebraic equations (Three value pairs for the standard Steinhart-Hart polynom). Temperatures typically used are for example 0° C, 15° C, 25° C and 70° C. By solving this system the values for a 0 , a 1 , a 2 and a 3 can be determined. A better approach is the mathematical optimization technique called ordinary least squares (OLS), that was introduced by Carl Friederich Gauß in 1801. Details about the theory of OLS are given at Wikipedia [ 3 ] or MathPlanet [ 4 ]. If the function that is approximated is a polynom, theory of vector spaces, scalar products and orthonormal bases eases the calculation. Given a list with n temperature-resistance pairs ( r 0 , t 0 ), ( r 1 , t 1 ), ..., ( r n-1 , t n-1 ) (where n should be at least 3), the question for best fitting of the Steinhart-Hart coefficients leads to the minimization of n-1 sum := Σ [ t(r i ) -t i ]² i=0 where t i is the i -th temperature value and t(r i ) is the calculated temperature according to the polynom. This optimization requires a small mathematical excursion. For a positive integer n and given abscissa values x 0 , x 1 , ..., x n-1 the polynoms with degree ≤ n build a vector space V . On V a scalar product can be defined by n [ p, q ] := Σ p(x i ) · q(x i ) i=0 By |p| := [ p, p ] ½ the scalar products makes V a normed space. For a function f with ordinate values y i := f(x i ) for i = 0, .., n-1 there is exactly one...