[发明专利]一种翼型阻力积分区域自动检测方法

权利要求书

1.一种翼型阻力积分区域自动检测方法，其特征在于，具体过程如下：

步骤1：确定待求截距C:

以尾耙测压点的坐标和得到的尾耙测压点的压力作为样本数据；对所述样本数据进行拟合以确定所述待求截距C；

对所述样本数据通过式(1)进行拟合:

X·C＝Y (1)

式中：X是尾耙测压点的坐标，X＝[x₁，x₂，……，x_n]，n是尾耙测压点总数量；C是待求截距；Y是测量得到的尾耙的压力，Y＝[y₁，y₂，……，y_n]；

步骤2：各样本点误差计算：

确定尾耙的压力Y的估计值

式中，是y的估计值；

得到估计值后，每个样本点的误差σ_i通过式(4)得到：

步骤3：计算最小二乘解误差：

通过式(5)求解步骤1中拟合结果的最小二乘解的误差σ；

步骤4：噪点检测：

设尾迹区的总压为负且给定噪点判定阈值K＝2σ、连续点数M＝5，使用下述方法判断样本点是否为噪点；

第一：遍历所有的样本点，步骤2中计算的样本点误差σ_i大于给定的噪点判定阈值K，即σ_iK的样本点标记为初始噪点；

第二：检测初始噪点的序号是否连续；

当连续的初始噪点数＞5时，将该初始噪点标记为噪点；其余不连续的初始噪点或连续点数≤5初始噪点数重新标记为样本点；

步骤5：筛选噪点：

若有所述噪点，则筛选所述噪点，并对剩余的样本点重复步骤1～步骤4，并继续对新的迭代结果进行噪点的筛选；

若本轮中没有检测到噪点，迭代结束进入步骤6；

步骤6：阻力积分区域检测：

步骤5中的迭代结束后，筛选出的各噪点之间在所述X坐标上所构成的区域为最终的阻力积分区域，其余样本点之间在所述X坐标上所构成的区域为非阻力积分区域；

完成翼型阻力积分区域的自动检测。

2.如权利要求1所述一种翼型阻力积分区域自动检测方法，其特征在于，所述样本点是读入的尾耙测量的数据。

3.如权利要求1所述一种翼型阻力积分区域自动检测方法，其特征在于，所述式(1)中的未知量X＝[x₁，x₂，……，x_n],其中的x₁，x₂，……，x_n分别一个方程,X的个数就是方程个数。

4.如权利要求1所述一种翼型阻力积分区域自动检测方法，其特征在于，拟合时,取斜率为0、y＝C的直线形式，且方程个数与样本数据中的样本点个数相同，故式(1)中方程的个数大于未知参数的个数，通过最小二乘法进行求解，得到待求截距C：

C＝(X^TX)^-1·X^TY (2)。