一种基于模糊测度的流场涡特征检测方法
【专利摘要】本发明公开了一种基于模糊测度的流场涡特征检测方法，目的是提高复杂流场中涡特征提取准确率，技术方案是根据流场特征测度的三个原则得到临界点特征向量、流线距离特征向量和欧拉距离特征向量，进而得到流场中各个顶点的的涡特征的特征向量；然后将各个顶点的涡特征的特征向量进行标准化处理，得到标准化特征向量；最后根据模糊聚类方法计算各个顶点的隶属度，隶属度大于0.5的顶点即为涡特征的特征区域。采用本发明能够准确的提取复杂流场涡特征结构，准确率较已有方法有很大的提高。
【专利说明】一种基于模糊测度的流场涡特征检测方法
【技术领域】
[0001]本发明涉及三维流场涡特征检测方法，尤其是一种基于模糊测度的流场涡特征检测方法。
【背景技术】
[0002]随着计算机性能的飞速提升，计算流体动力学的计算网格规模越来越大，计算的解包括定义于数亿个网格点上的若干物理量。面对如此海量数据，如何利用特征提取技术捕捉其蕴含的重要特征信息，一直是科学计算可视化研究领域所关注的问题之一。涡特征作为流场中的典型特征结构，并没有准确的定义，这导致了涡特征检测比较困难。而模糊集合论比较适合于描述具有非硬性边界的集合。复杂的流场涡特征结构，难以精确界定其区域边界，因此适合采用模糊集理论对其进行提取。
[0003]已有的特征提取方法如λ 2方法、Q-标准方法等经典的涡特征提取方法都存在着误检的情况，尤其是针对复杂流场，涡特征提取的准确率较低。如何提高复杂流场中的涡特征提取的准确性，是流场特征提取的重要问题。
【发明内容】
:
[0004]本发明要解决的技术问题是:针对复杂流场中涡特征提取准确率较低的问题，提出一种基于模糊测度的流场涡特征提取方法，该方法能够准确的提取复杂流场涡特征结构，准确率较已有方法有很大的提高。
[0005]技术方案是:从模糊理论出发，将流场涡特征做为模糊集合来处理，针对涡特征结构采用不同的测度方法，根据这些测度方法计算出涡特征结构的标准特征向量；最后采用基于模糊聚类的涡特征结构提取方法实现典型流场涡特征结构的有效提取。
[0006]具体的技术方案是:
[0007]第一步:根据流场特征测度的三个原则得到临界点特征向量、流线距离特征向量和欧拉距离特征向量，得到流场中各个顶点的涡特征的特征向量。
[0008]1.1:初始化顶点循环变量k为O。
[0009]1.2:由于涡核点均为临界点，因此采用临界点测度方法对序号为k的顶点进行测度。由临界点(流场中速度为O的顶点)定义可知临界点测度原则:当流场U中的顶点Xk速度为零时，将顶点Xk的隶属度(顶点的隶属度是指该顶点隶属于涡特征的程度，取值为O到
I)映射为I。其中Xk表示流场中序号为k的顶点。因此，根据临界点测度原则可得到顶点Xk的6维(即6元组)的临界点特征向量=V1(Xk) = (U，V, W，V1, V2, V3)，该特征向量用于计算特征隶属度，其中前三个元素U，V, w为流场U中顶点Xk的三维坐标，后三个元素V1, V2, V3为顶点Xk的速度值。
[0010]1.3:利用流线距离对顶点与涡核点的关系得到如下流线距离测度原则:当流场U中的两个顶点x，y相对于同一临界点的流线距离满足sd(x，c) <sd(y, c)时，x，y的隶属度值存在μ (χ) > μ (y)关系，其中sd(x，c)为点χ和点c之间的流线距离,sd(y, c)为点y和点c之间的流线距离，两点之间的流线距离是指，当两点之间存在流线经过时，流线距离为两点之间的流线段的长度；如果两点之间不存在经过该两点的流线，则流线距离为无穷大，μ (χ)表示点χ的隶属度值，μ (y)表示点y的隶属度值。也就是说当顶点χ到临界点C的流线距离sd(x, c)小于顶点y到c的流线距离sd(y, c)时，χ的隶属度μ (χ)大于y的隶属度μ (y)。根据流线距离测度原则，采用距离各临界点的最小流线距离作为特征向量，得到如下I维(即一元组)流线距离特征向量:V2 (xk) = (Smin(Xk)),其中Smin(Xk)为顶点Xk到流场中所有临界点的最小流线距离。
[0011]1.4:根据欧拉距离的测度原则，采用距离各临界点的最小欧拉距离作为特征向量，得到如下I维欧拉距离特征向量:v3(xk) = (dmin(xk)),其中dmin(xk)为顶点xk到流场中所有临界点的最小距离。欧拉距离的测度原则是指当顶点X，Y距离同一临界点c的欧拉距离存在d(x, c) < d(y, c)关系时,隶属度有μ (χ) > μ (y)。其中d(x, c)表示x、c两点之间的距离。
[0012]1.5:综合上述测度原则，可得到涡特征的特征向量:
[0013]V (Xk) = (U，V, w, V1, V2, V3, smin (xk)，dmin (Xk))。
[0014]1.6:令 k 增加 I。
[0015]1.7:如果k小于流场U中顶点数目，则跳转至1.2步，否则执行第二步。
[0016]第二步:将各个顶点的涡特征的特征向量进行标准化处理，得到标准化特征向量
2.1:初始化顶点循环变量k为O。
[0017]2.2:对顶点Xk涡特征的特征向量的各个分量即U，V, W，V1, v2, v3, smin(xk), dmin(xk)进行标准化处理。标准化是指对其进行[0，1]区间标准化处理。设V(X)内的某一分量为m,其最大值、最小值分别为mmax和mmin，则标准化处理如下公式:
[0018]
【权利要求】
1.一种基于模糊测度的流场涡特征检测方法，其特征在于包括以下步骤: 第一步:根据流场特征测度的三个原则得到临界点特征向量、流线距离特征向量和欧拉距离特征向量，得到流场中各个顶点的涡特征的特征向量； .1.1:初始化顶点循环变量k为O ; .1.2:根据临界点测度原则得到顶点Xk的6维临界点特征向量=V1(Xk)=(U，V, W，V1, V2, V3),前三个元素U，V, W为顶点Xk的三维坐标,后三个元素V1, V2, V3为顶点Xk的速度值；临界点测度原则是指:当流场U中的顶点Xk速度为零时，将顶点Xk的隶属度映射为1，顶点的隶属度是指该顶点隶属于涡特征的程度，取值为O到I ;xk表示流场中序号为k的顶点； .1.3:根据流线距离测度原则，采用距离各临界点的最小流线距离作为特征向量，得到I维流线距离特征向量:V2 (Xk) = (Smin(Xk)),其中Smin(Xk)为顶点Xk到流场中所有临界点的最小流线距离；流线距离测度原则是指:当流场U中的两个顶点X，y相对于同一临界点的流线距离满足sd(X，c) <sd(y，c)时，x，y的隶属度值存在μ (χ) > μ (y)关系,其中sd (x, c)为点χ和点c之间的流线距离，sd (y，c)为点y和点c之间的流线距离，两点之间的流线距离是指，当两点之间存在流线经过时，流线距离为两点之间的流线段的长度；如果两点之间不存在经过该两点的流线，则流线距离为无穷大，μ (χ)表示点X的隶属度值，μ (y)表示点I的隶属度值；也就是说当顶点X到临界点C的流线距离sd(x, c)小于顶点y到c的流线距离sd(y, c)时，χ的隶属度μ (χ)大于y的隶属度； .1.4:根据欧拉距离的测度原则，采用距离各临界点的最小欧拉距离作为特征向量，得到I维欧拉距离特征向量:V3(xk) = (dmin(xk)),其中dmin(xk)为顶点χ到流场中所有临界点的最小距离；欧拉距离的测度原则是指当顶点x，y距离同一临界点c的欧拉距离存在d(x, c) < d(y, c)关系时,隶属度有μ (χ) > μ (y);其中d(x, c)表示x、c两点之间的距离； .1.5:综合1.1-1.4的测度原则，可得到涡特征的特征向量:
V (Xk) = (U, V, w, V1, V2, V3, Smin(Xk)，dmin (xk))； . 1.6:令k增加I ; .1.7:如果k小于流场U中顶点数目，则跳转至1.2步，否则执行第二步； 第二步:将各个顶点的涡特征的特征向量进行标准化处理，得到标准化特征向量: .2.1:初始化顶点循环变量k为O ; .2.2:对涡特征的特征向量的各个分量即U，V, W，V1, V2, V3, smin(xk)，dmin(xk)进行标准化处理；标准化是指对其进行[O，I]区间标准化处理；SV(x)内的某一分量为m，其最大值、最小值分别为mmax和mmin,则标准化处理如公式1:
【文档编号】G06T7/00GK103440650SQ201310375637
【公开日】2013年12月11日 申请日期:2013年8月26日 优先权日:2013年8月26日
【发明者】曾亮, 徐华勋, 王怀晖, 蔡勋, 李思昆, 王文珂, 沈恩亚 申请人:中国人民解放军国防科学技术大学