Perhaps obvious, but I thought I'd mention it. A general matching process:
collect M measurement points for all N devices
compute sum(k) | Pik - Pjk | for i != j for the N**2-N combinations and the M measurements
sort on increasing sum, i and j and keep the first i in the list
This can be refined by including measurements at multiple temperatures, measuring noise figure and by restricting the measurements used for the sum to the desired operating conditions. Collecting a full set of curves for each device will let you match or select the best device from inventory without remeasuring.
And yeah, I'm afraid I'm a bit NF obsessed as I missed out on getting some NOS very low noise figure JFETS and many of the best are no longer made.
Edit: Sorry about that. I added the bounds N & M at the last minute. It's just the calculation of the least summed absolute error followed by brute force selection of the optimal matches. I was actually trying to work out the equations for heat transfer and temperature for TEC modules when this strayed into my head. The TEC problem involves an exponential with a negative exponent, but I still haven't figure out how you get there from the information in the datasheet. I think it's time to pull my Diff Eq book off the shelf and give myself a refresher. I'd hoped that I would recognize the solution, but it's just been too long.