[发明专利]一种去除风洞风扇干扰的数据处理方法

权利要求书

1.一种去除风洞风扇干扰的数据处理方法，其特征在于，包括如下步骤：

步骤S10：将各个原始时域脉动压力信号均划分为原始时域脉动压力静态信号和原始时域脉动压力动态信号；

步骤S20：将各个原始时域脉动压力动态信号转换为原始频域脉动压力动态信号，计算各个原始频域脉动压力动态信号对应的风扇干扰频率，并依据所有的原始频域脉动压力动态信号对应的风扇干扰频率计算风洞风扇干扰基频；

步骤S30：判断各个原始频域脉动压力动态信号是否需要去除风洞风扇干扰，当原始频域脉动压力动态信号需要去除风洞风扇干扰时，将原始频域脉动压力动态信号中的风洞风扇干扰去除，得到去除风洞风扇干扰的频域脉动压力动态信号；

步骤S40：将步骤S30中的去除风洞风扇干扰的频域脉动压力动态信号转换为去除风洞风扇干扰的时域脉动压力动态信号，并将去除风洞风扇干扰的时域脉动压力动态信号与步骤S10中的原始时域脉动压力静态信号叠加，获得去除风洞风扇干扰的时域脉动压力信号；

所述步骤S30包括如下步骤：

步骤S31：提取第i个原始频域脉动压力动态信号中位于区间[k×f₀-d_f，k×f₀+d_f]内的频域幅值向量b_i、相位向量ϕ_i，计算频域幅值向量b_i中的元素的均值B_i-1、标准差B_i-2、最大幅值B_i-max；

步骤S32：计算第i个原始频域脉动压力动态信号的残值向量D_i，其中，D_i=b_i-B_i-1-d×B_i-2，d为风扇干扰峰判别限制阀；

步骤S33：根据残值向量D_i，判断第i个原始频域脉动压力动态信号在区间[k×f₀-d_f，k×f₀+d_f]内是否存在显著干扰峰；

步骤S34：当第i个原始频域脉动压力动态信号在区间[k×f₀-d_f，k×f₀+d_f]内不存在显著干扰峰时，不对第i个原始频域脉动压力动态信号继续进行数据处理；当第i个原始频域脉动压力动态信号在区间[k×f₀-d_f，k×f₀+d_f]内存在显著干扰峰时，计算残值向量D_i中元素的最大幅值对应的频率f_i-max，判断f_i-max是否偏离中心超过允许值；

步骤S35：当f_i-max没有偏离中心超过允许值时，对第i个原始频域脉动压力动态信号继续进行数据处理；当f_i-max偏离中心超过允许值时，不对第i个原始频域脉动压力动态信号继续进行数据处理；

其中，k为大于等于1的整数；

所述步骤S33包括如下步骤：

步骤S331：计算残值向量D_i中的各元素之和；

步骤S332：当残值向量D_i中的各元素之和小于等于0时，判断第i个原始频域脉动压力动态信号在区间[k×f₀-d_f，k×f₀+d_f]内不存在显著干扰峰；当残值向量D_i中的各元素之和大于0时，判断第i个原始频域脉动压力动态信号在区间[f₀-d_f，f₀+d_f]内存在显著干扰峰。