Например если данные - это кусочнолинейная функция, то вторая производная существует только как обобщенная т.е. с уходами в +/- бесконечность...

"Рабоче-крестьянский" метод такой
Вторая производная приближенно красиво определяется как
Y"(t)= (T(t+h) + T(t-h)- 2* T(t))/(h*h)
где T(t) - например кусочно-линейная интерполяция значений температуры
h - это параметр регуляризации в НЕ РАВНЫЙ и скорее всего больший чем период измерений, он подбирается экспериментально по минимальному уровню выбросов и максимальной точности что достигается на экспериментальных примерах.
Далее имеем уравнение Y"(t)=0 решаемое ТОЛЬКО численным методом откуда имеем кандидатов на перегиб.
Естественно нужно б прогнать данные сначала в компутере посмотреть на график Y"(t) и вообще что выходит...

"Научный" смотрите по ссылке
http://www.srcc.msu.su/num-meth/zhurnal/tom_2005/pdf/v6r105.zip