I had an interesting Temp. control design once.
During the mid 90's I worked in Winnipeg, MB Canada as Operations Mgr / Tes t Eng for an Automated 2 way Meter Reading company for 4 years call Iris S ystems Ltd that later sold all their patents to Itron in Seattle. It used o ur custom 928 MHZ ISM band cell network from the power meter directly into any utility's Co. Sun database on 6KHz shared bandwidth MUX'd in time, dive rsity with a complex protocol developed in-house. The transmitter had to c reate a 1 ppm TCXO out of a 25 cent AT-cut Xtal.
I developed 3 order polynomials curve fit for every angle cut of AT curves then another Polynominial to generate any curve from 2 test points on the c urve at 40'C, 70'C This linear differential value could be used to predict the entire curve from -40 to 70'C which was our outdoor environmental desi gn spec.
As Hobbs said "An 0603 thermistor with one side soldered to a copper pour . .. has a thermal TC of around 100 milliseconds. "
I made a cheap OCXO design and used styrafoam around a 2oz copper FPC Lexan film with a 0603 thermistor to external Rref values and 1206 heater resist ors to heat up one side of the crystal to match resistance and measure the Xtal connected to a header with a CMOS XO while the XTAL was in a mini cop per foil around the case temp so the frequency shift of 10T was about 4 sec onds from room temp to +40'C, took a reading error in PPM error on the f co unter then switched to 70'C and did the same. Removed the crystal from the header in 10 seconds and sorted the Xtals in bins for prototype (100) level test. The bins then matched with sorted Varicaps were used by the MC6805 t o control the voltage and tune the TCO design within 1ppm over a 110'C rang e for about $1.
So we used an SMT thermistor and heater like an OCXO with a removal Xtal th en rapidly bin the parts. Using my algorithm in the final product, measure the temperature to make the XTAL < 1PPM only took 10 seconds to predict and correct the error over a 110'C range using only 2 precise points on the cu rve.
Those who understand, imagine the family of 3rd order curves going thru 25' C.
My contribution was defining the precise Polynomials for all AT crystal cur ves.
I used a DOS program I got from MJT Microwave Journal to define the coeffic ients.
While my regular job as Ops Mgr for SMT Microwave Prototype Assy was the 10 k part number Master Registry the internal part numbering system, the Foxba se Stock room Pick lists from BOMs, Purchasing, QA and inventory control sy stem and the entire SMT Prototype shop operations. We had a total staff of 15 HW designers, 15 software gurus and 15 other support staff and 3 mgrs/f und raisers. We used optical detection for the power meter rotation of existing meters a nd all the electronics of Radio and power detection fit behind the metal la bel in almost every design. However, most of the industry rejected this an d went to new digital meters without 2-way direct connection to the head of fice database and still required labour of meter readers.
Now you can buy these. 1PPM TCXO's for 1~2$