Sensor Temperatura Oleo + Arduino + Matlab

Bem, preciso ler a temperatura do óleo do motor do carro.
Um motor antigo, que utiliza um sensor genérico de chavete/fusca/monza.

Neste caso o modelo é um Marflex 8011, cujo fabricante não fornece especificações técnicas do sensor, e sim apenas a aplicação em veículo/modelo/ano.

Para este sensor, quando medido sua resistência em temperatura ambiente : 25C, o mesmo apresenta resistência de 2200ohm (2k2).
Por ser um sensor tipo NTC, quanto maior a temperatura, menor a resistência apresentada pelo sensor.

Visando mapear a relação temperatura/resistência, conectamos o mesmo em um multimetro e colocamos o sensor imerso em agua, tomando o cuidado para não curto-circuitar os polos positivo (terminal) e negativo (rosca)  do sensor. Assim apenas a parte da rosca do sensor ficou submersa.

Juntamente com o sensor, foi imerso um termômetro o qual indicará a temperatura da água em um dado instante, podendo relacionar a temperatura indicada, com o valor de resistência medido.

Em verde o multimetro indicando a resistência em kOhm e em vermelho o termômetro indicando a temo em C

Iniciou-se o processo com agua "gelada" a 10C, e elevou-se até a fervura 96C.
A aproximadamente cada 5C de elevação de temperatura, foi medido o valor da resistência, fornecendo um total de 19 pontos medidos.

Temperatura C	Resistencia kOhm
10	4
15	3,3
18	3
20	2,8
25	2,2
30	1,7
35	1,35
40	1,15
45	0,9
50	0,75
55	0,65
60	0,5
65	0,45
70	0,35
75	0,3
80	0,25
85	0,2
90	0,18
95	0,15

Ao traçarmos o gráfico que relaciona a temperatura C com a resistência do sensor em kOhm, obtemos uma curva aproximada que não caracteriza bem o comportamento do sensor.

Visando aproximar esta cursa, e com o uso do Matlab, foi calculado um regressão exponencial de 2 ordem, o qual estimou os melhores coeficientes da curva para o sensor.

M = load('sensor.txt');
%Computa regressao exponecial de 2 ordem
%ft=fittype('exp2');
%cf=fit(M(:,2),M(:,1),ft);
% aplica valores lidos:
dados = M(:,2);
%
% General model Exp2:
% cf(x) = a*exp(b*x) + c*exp(d*x)
% Coefficients (with 95% confidence bounds):
% a =       73.89  (67.5, 80.28)
% b =      -5.436  (-6.305, -4.566)
% c =       69.19  (63.08, 75.3)
% d =     -0.4682  (-0.5291, -0.4074)
ptr=1;    
for i=1:22
    result(ptr) = 73.89*exp(-5.436*dados(i)) + 69.19*exp(-0.4682*dados(i));
    ptr=ptr+1;
end
result

Assim, o grafico abaixo ilustra em vermelho, a curva medida, e em azul, a curva aproximada.
(Clique na imagem para aumentar)

A curva em azul informará melhor o valor da temperatura com base na resistência medida, e mapeado conforme equação abaixo:

Temperatura = 73.89*exp(-5.436*Resistencia) + 69.19*exp(-0.4682*Resistencia);

Novos testes ainda são necessários, e quem sabe o fabricante não queira facilitar minha vida :-)

[]'s