Решили мы таки попробовать меньшим количеством датчиков обойтись. Вот так выглядит график изменения сопротивления трех последовательно включенных датчиков (с интервалом в 5см):
Здесь синяя кривая — перемещение вниз, красная — вверх. Уровень жидкого азота на H=21см, по горизонтальной шкале откладывается H от горловины до самого нижнего датчика (т.е. на 31см все три уже погружены в азот). А потом я оставил один-единственный датчик (и H теперь — уровень от края горловины до датчика).
При движении вверх видно заметное отличие. Просадка на 9см скорей всего из-за того, что я не дождался стабилизации данных. Красная кривая как раз и интересна, т.к. по мере вытекания уровень азота будет уменьшаться. А вот так выглядит температура, измеренная по формуле для pt1000.
И код для октавы:
function Tout = pt1000(Rin)
% Converts resistance of pt1000 into T (degrC)
T = [-200:0.1:50];
_A = 3.9083e-03;
_B = -5.7750e-07;
_C = zeros(size(T));
_C(find(T<0.)) = -4.1830e-12;
rT = 1000.*(1 + _A*T + _B*T.^2 - _C.*T.^3*100. + _C.*T.^4);
Tout = interp1(rT, T, Rin, 'spline');
endfunction
У нас температура кипения азота примерно 197℃, а сопротивлению 192 Ом соответствует -198.4℃, т.е. врет датчик почти на полтора градуса (что вполне нормально, учитывая что он дешевый, а не HEL705 какой-нибудь за тысячу с хвостиком). В пяти сантиметрах от горловины температура около -50℃, так что силиконовым прокладкам будет вполне нормально.