[发明专利]一种确定混波室内搅拌器独立搅拌位置数的方法

权利要求书

1.一种确定混波室内搅拌器独立搅拌位置数的方法，其特征在于：包括以下步骤：

第一部分：获取统计样本

(a)在混波室工作空间内部确定测试场点N_w，空间相关函数验证空间场点是相互独立的，即

其中，r为场点间距，k_N为波数；当空间场点间距不小于四分之一波长，即认为场点之间是相互独立的；

(b)将独立的空间场点数记为N_w、搅拌位置总数记为N_s、独立搅拌位置数记为N_β，对于每个频点，E_1≤k≤Nw＝(E_k1,E_k2,…,E_kNs)为每个搅拌位置不同空间场点构成的行向量，经N_s次搅拌，得到N_w×N_s维矩阵样本；

第二部分：确定独立搅拌位置数

(c)任取N_w×N_s维矩阵样本中的第i列，并将其记为Eⁱ＝(E_1i,E_2i,…,E_Nwi)；

(d)计算任意两个搅拌位置的皮尔逊相关系数r_ij：

并构成N_s×N_s维相关矩阵Rs，并取该矩阵的上三角矩阵重新赋值给相关矩阵Rs；

(e)从第一行的第二列至第N_s列依次检验r_ij，若小于相关系数阈值r_s，跳过检查下一列，反之删除该列，此时独立搅拌位置数N_β由初始值0加上1，即N_β＝1；相关系数阈值r_s由空间采样点和期望的置信水平共同决定，对于给定95％的置信水平下，其阈值在N_w＝8时为0.707，在N_w＝9时为0.666，即

(f)删除第一行，并将删除后的相关矩阵重新赋值给相关矩阵Rs，将赋值后相关矩阵Rs的秩大小赋值给N_s；

(g)对赋值后的相关矩阵Rs，重复步骤(e)和(f)，且每重复1次，独立搅拌位置数N_β加1；

(h)当赋值后的为1时，重复过程停止，输出N_β即为独立搅拌位置数。