@学流体的小明

这两个图分别对应方法1 和方法2 的计算结果吗，还是不同方向的结果？

对于方法1 ，这个写法是系综平均求相关函数，最后还要归一化得到相关系数。如果时间非平稳和空间非均匀，理论上统计量跟选择的参考点位置x和时间t都相关。但槽道流的在展向和顺流向是均匀+时间平稳的，所以最后的统计特性都与参考位置参考点位置x和时间t都无关，所以可以用空间和时间一起平均得到最终统计量，of自带的postChannel工具就是这么做的。

计算空间相关函数有2种方法：
（1）方法1，根据《Shin K, Hammond J. Fundamentals of signal processing for sound and vibration engineers[M]. England: John Wiley & Sons, 2008.》，已知在固定点位置，一条时间序列 u(x=x0, y=y0, z=z0, t)，针对这条时间序列求时间相关函数公式如下截图，在MATLAB中对应函数是 xcorr 。这种方法要求序列足够长，并且点数足够多，才能算的准确。

类似的，在某一固定时刻，沿着空间展向z的一条空间序列 u(x=x0, y=y0, z, t=t0)，也一样可以用xcorr 求空间相关函数，计算结果要准确，就要求这条序列的展向宽度要足够大，且点数足够多，否则计算结果就会偏差大。 如果计算域纵向足够长，把每条展向空间相关函数做空间平均，得到的结果应该也会比较准确。想要更准确就要加上每个时刻的平均，就要输出每个时刻的数据，计算量太大了，用空间平均了，其实时间平均应该就不用了，结果差别不会太大。

（2）方法2，跟下面公式一样，计算两个序列的相关系数，对应MATLAB函数corrcoef。已知展向上空间多个点A(x0, y0, zA, t)、B(x0, y0, zB, t)、C(x0, y0, zC, t)的时间序列，就可以假设选定A为参考点，依次计算出 AB、AC的相关系数，然后就可以把连成展向空间相关系数曲线。这种方法就要设置展向监测点，测量出每个点的时程数据，相对来说时程数据比空间数据（受限于计算域大小）更长，计算结果应该会比较准确。还可以设置多个不同x向位置的测点，再对计算的空间相关系数曲线平均，结果会更准确。

CFD中文网

coolhhh