使用二分法查找数据进行线性插值的方法

           最近高温测试的项目接近了尾声，标定工作成了现阶段的主要任务，对光学高温传感器的标定一般是使用热电偶来进行标定。

       前期标定使用的多项式拟合的方法，这种方法的原理就是使用热电偶的数据和光学传感器的腔长数据进行来找到一个高阶的多项式，而后在程序中使用这个多项式来根据腔长计算出实际的温度，但是这种方法在程序中的效果并不是很好，计算出来的温度和热电偶的温度相差很大，跟BOSS汇报后，指示使用实验数据进行线性插值方法。

       首先来说下线性插值的原理，首先有规律的记录一组腔长和热电偶的数据，而后使用得到的腔长数据进行搜索，在最接近的两点之间使用正比例函数进行插值。公式为：

           y = y₀  + (x-x₀)*(y₁-y₀)/(x₁-x₀)

它是由直线解析函数的差分式推导而来：

           (y-y₀)/(x-x₀)=(y₁-y₀)/(x₁-x₀)

是最简单的线性插值方法，那么现在这个算法的关键部分就是在于从数据表中搜索到两个和目标数据最接近的两个数据，由于我们的数据表示有序的，我们可以使用二分查找的方法轻松的找到这两个数据。

           二分查找法的基本原理就是建立low和high这个索引量来进行二分查找，但它最初的目的是要查找到目标数据，这两个量在查找结束后本身是没什么作用的，但是在我们这里，这两个量就是最接近目标数据的两个量。但是二分法查找结束的条件是low>high，我们要进行线性插值，则要将查找结束后的两个量进行交换。算法代码如下：

 /*函数创建日期：2010-10-25*******************/ /*函数功能：该函数作用是利用腔长值进行搜索并使*/ /********* 用相应的两个热电偶温度进行线性插值*/ /********* ，从而求得当前温度***************/ /*函数作者：王晓军***************************/ /*函数参数：fIn: getGap 得到的腔长值*******/ /* fInThermo 热电偶数据数组 fInGap 腔长数据数组 */ /*返回值： fOut:tempertaure 得到的温度值****/ int __fastcall TForm1::GetTemperature(const int *fInThermo,const double *fInGap,const double getGap) { /*二分法进行查找*/ int low = 0; int high = TAB_NUM - 1; int mid = 0.0; int temperResult = 0; while(low <= high) { mid = (low + high)/2; if(fInGap[mid] == getGap) { temperResult = fInThermo[mid]; //中间值即为所要查找的值 break; } else //处理反序排列的腔长数据 if(fInGap[mid] > getGap) low = mid + 1; else high = mid - 1; } if(low <= high) return ((int)temperResult); else { /*没有查找到时找到两个相邻值进行线性插值，高位和地位复位*/ int temp = 0; temp = high; high = low; low = temp; /*线性插值*/ temperResult = fInThermo[low] + (fInThermo[high] - fInThermo[low]) * (getGap - fInGap[low]) / (fInGap[high] - fInGap[low]); return ((int)temperResult); } }

       目前读取的标定数据还是从txt文本当中读取，为了更方便的读取和修改数据，可以使用本地Acsess数据库文件来代替txt文件。