基于 Sebastian Endrikat https://www.mathworks.com/matlabcentral/fileexchange/52343-find-vortices-in-velocity-fields并受其启发
L Graftieaux 提出的基础数学“结合 PIV、POD 和涡流识别算法用于研究不稳定的湍流涡流” https://iopscience.iop.org/article/10.1088/0957-0233/12/9/307/pdf 非常感谢!
但是,Sebastian的代码相当慢,因为它使用了4个嵌套的for循环。 通过实现 Graftieaux 论文中描述的操作可以重新排列为两个卷积,可以在 Matlab 中更有效地生成 gamma1 映射。 这段代码基本上可以做到这一点。 它还通过拟合样条曲线找到最大伽马 1,因此不受矢量场网格分辨率的限制。
2021-07-21 14:59:21
3KB
matlab
1