一种实用的线性插值方法.docxVIP

一种实用的线性插值方法 1 孔探针测量三维流场 5个孔检测用于测量流量数据的性能参数，如流量大小、方向、总压和静压。由于这种方法具有测量原理简单、使用方便和探针不易损坏及便于维修等优点,目前仍是测量三维流场的主要手段之一。应用五孔探针测量三维流场有3种方法,(1)对向测量法;(2)半对向测量法;(3)非对向测量法。TreasterYocum对五孔探针探针的分类、使用方法以及数据处理作了详细论述。五孔探针采用非对向测量三维流场时,首先根据五孔压力可计算出两个方向校准系数,然后将两个方向校准系数在校准曲线中插值得到偏转角和俯仰角以及总压校准系数和静压校准系数。插值方法是否适当直接影响着实验结果的准确性。在文献中对探针校准曲线使用了样条插值,而在文献中探针校准曲线采用最小二乘法拟合。本文对这两种方法进行了理论研究,并在此基础上提出了一种实用性更强的线性插值法。 2 测量三维流场 五孔探针由位于两个正交平面上的5个测压孔组成,两个平面(偏转角和俯仰角)的交线通过中心孔,偏转角α、俯仰角β空间3个分量的速度如图1所示。应用五孔探针测量三维流场有3种方法,(1)对向测量法;(2)半对向测量法;(3)非对向测量法。这3种方法在文献中均有较详细的介绍,本文仅对非对向测量的数据处理进行论述。非对向测量校准系数如下: α方向校准系数: Κα=p4-p5p2-ˉp(1) β方向校准系数: Κβ=p3-p1p2-ˉp(2) 总压校准系数: Cpt=p2-ptp2-ˉp(3) 静压校准系数: Cps=ˉp-psp2-ˉp(4)ˉp=(p1+p3+p4+p5)/4(5) 其中:pt—来流总压,ps—测点静压。 3 实验数据中的插值 应用非对向测量法测量流场时,可测得五孔探针1～5孔所对应的压力p1、p2、p3、p4、p5和来流总压p*0,这些值都是相对于参考压力值Pref(大气压力)的。将这些数据按式(1)～式(5)求出α方向校准系数KSα和β方向校准系数KSβ(为了将实验所测的数据和校准系数曲线中的参数区别开,对实验数据使用了上标S)。再由KSα和KSβ通过在校准系数中插值可以得到实际气流向对于探针坐标系的偏转角αS和俯仰角βS以及总压系数CSpt和静压系数CSps。在插值过程中国内外学者使用了不同的方法。不同插值方法所使用的方向校准系数曲线如图2所示。图2(a)为使用式(1)～式(5)计算出的校准曲线,图2(b)为样条插值后的曲线,图2(c)为最小二乘法拟合后的曲线。 3.1 探针插值效果 文献中使用样条插值法进行数据处理。这种插值法在原理上精度高,但同时对探针质量要求也高。作者对探针校准系数用第一种边界条件的3次样条插值后的曲线来表示,如图2(b)所示。从图中可以看出,就本探针而言在某些区域如插值效果并不理想,误差较大。对于性能不好的探针,两条相邻曲线还有可能相交,进而插出错误数据,严重时会出现数据溢出。所以样条插值法的使用范围就局限在了特性良好的五孔探针上。对于特性较差的探针,这种方法并不适用。 3.2 方向校准系数拟合的残差 在文献中采用五次多项式拟合Kα-Kβ、Cpt-α和Cps-α校准曲线,拟合后的方向校准系数曲线如图2(c)。校准多项式为: f(x)=α0+α1(x-x0)+α2(x-x0)2+ α3(x-x0)3+α4(x-x0)4+α5(x-x0)5(6) 其中:x0—已知数据的几何平均值;α0～α5—多项式系数。 在使用拟合曲线g(x)取代离散点(x,y)时要尽量保证这些离散点尽可能的在g(x)附近。附近的程度取决于各种逼近的标准,用残差作为拟合的标准,其残差表示为: ri=yi-g(xi) (7) 残差向量的3种范数记作: ∥r∥∞=max1sism|ri|(8)∥r∥1=mΣi=1|ri|(9)∥r∥2=(mΣi=1r2i)1/2(10) 残差可以表示拟合的误差,误差愈小拟合就愈好。最小二乘法取2—范数‖r‖2最小,同时希望1—范数‖r‖1或∞—范数‖r‖∞也较小。使用最小二乘法对方向校准系数拟合后的各曲线残差在表1中给出。 从图2(c)中可以看出,α或β每增加3°,Kα或Kβ平均增加约为0.3。如果认为角度误差在0.1°以内,那么Kα或Kβ的误差应该保证在0.01以内,这就要求残差的∞-范数‖r‖∞≤0.01。从表中可以看出,没有哪一条拟合曲线能够达到该要求。这说明最小二乘法对探针特性要求也较高。另外,对于探针特性较差的位置,最小二乘法在拟合曲线的过程中大幅度修改了这些点的参数,如α=-12、β=-12点,使得拟合后的校准曲线不再能够真实反映探针的特性。 3.3 凹四边形插值标准 该方法由两部分组成,首先根据求出KSα和KSβ在方向校准系数网格中的位置,然后再用线性插值法计算出方向角αS和βS。第一步由图3可以看出