一种声波照射随机起伏界面时的遮挡面元快速搜索方法技术

本发明专利技术属于水下目标声散射领域，具体涉及的是一种数值计算随机起伏界面的散射声场时，搜索未被声波照射到的遮挡面元的声波照射随机起伏界面时的遮挡面元快速搜索方法。本发明专利技术包括：导入已进行三角面元划分的随机起伏界面的N个面元节点信息、面元中心点信息以及声源位置信息；求取所有面元中心点三维坐标中z分量的最大值Zmax，得到z＝Zmax平面等。本发明专利技术由于判断面元是否被遮挡时的搜索区域减小，减小了计算量，提高了遮挡面元的搜索效率；适用于随机起伏界面的遮挡问题，且起伏界面面积越大，计算效率提高的越明显；只利用了空间几何中点、线、面的有关计算，没有复杂的处理算法，易于程序实现。

【技术实现步骤摘要】

本专利技术属于水下目标声散射领域，具体涉及的是一种数值计算随机起伏界面的散 射声场时，搜索未被声波照射到的遮挡面元的声波照射随机起伏界面时的遮挡面元快速搜 索方法。
技术介绍
在利用物理声学方法预估起伏界面或及其附近目标的散射声场时，需要对面元划 分后的起伏界面模型进行遮挡面元的搜索，进而剔除掉被遮挡的面元。传统的方法是进行 遍历搜索，判断每一个面元之间是否相互遮挡，因此计算量较大。在电磁散射领域也采用了 一些遮挡面元快速搜索的方法，如拓扑踪迹方法（杨涛，陈玉林，飞机发动机复杂腔体遮挡 算法研宄，雷达科学与技术，2013,第11卷第2期）、面元分组方法（姬金祖，刘站合，基于面 元分组的电磁遮挡算法及其优化，北京航空航天大学学报，2009,第35卷第4期）等。以上 方法都是针对规则几何体目标而提出的快速搜索方法，本专利技术针对随机起伏界面的面元遮 挡问题，专利技术了。
技术实现思路
本专利技术的目的在于提供一种解决现有技术搜索过程计算量大、计算时间过长的问 题的声波照射随机起伏界面时的遮挡面元快速搜索方法。 本专利技术的目的是这样实现的： (1)导入已进行三角面元划分的随机起伏界面的N个面元节点信息、面元中心点 信息以及声源位置信息； (2)求取所有面元中心点三维坐标中z分量的最大值Zmax，得到z = Zmax平面； (3)对于第η个面元，连接声源与其中心点，求得连线与z = Zmax平面的交点 (X'，y'，Zmax)，所求交点在xoy平面的投影点为（X'，y'），第η个面元中心点在xoy平面的投 影点为（x n，yn); (4)如果 X' < xn，令 X' * = X' - Δ 1，χη* = χη+Δ 1 ;如果 X' > χη，令 X' * = X' +Δ 1， Xn* = Χη_Δ 1 ;如果 y' 彡 yn，令 y' * = y' -Δ l，yn* = yn+A 1 ;如果 y' > yn，令 y' * = y' +Δ 1， yn* = yn_△ 1，其中△ I为面元边长，在x〇y平面上得到顶点分别为（x' *，y' *)、（χη*，y' *)、 (xn*，yn*)、（X，*，yn*)的长方形区域； (5)随机起伏界面中面元中心点在xoy平面上投影点为（X，y),取所有投影点满足 X e 且y e 条件所对应的面元，得到与步骤（4)中求得的长方形区域 对应的M个面元； (6)对于步骤（5)中的M个面元中的第m个面元，求其所在平面与步骤（3)中声源 与第η个面元中心点连线的交点P，求得【主权项】1. ，其特征在于： (1) 导入已进行三角面元划分的随机起伏界面的N个面元节点信息、面元中心点信息 以及声源位置信息； (2) 求取所有面元中心点三维坐标中z分量的最大值Zmax，得到z = Zmax平面； (3) 对于第η个面元，连接声源与其中心点，求得连线与z = Zmax平面的交点 (X'，y'，Zmax)，所求交点在xoy平面的投影点为（X'，y'），第η个面元中心点在xoy平面的投 影点为（x n，yn); ⑷如果 X' 彡 xn，令 X' * = X' - Δ 1，χη* = χη+Δ 1 ;如果 X' > χη，令 X' * = X' + Δ 1，χη* =Χη_Δ I ;如果 y' 彡 yn，令 y' * = y' - Δ 1，yn* = yn+A I ;如果 y' > yn，令 y' * = y' + Δ 1， yn* = yn_△ 1，其中△ I为面元边长，在x〇y平面上得到顶点分别为（x' *，y' *)、（χη*，y' *)、 (xn*，yn*)、（X，*，yn*)的长方形区域； (5) 随机起伏界面中面元中心点在xoy平面上投影点为（x，y),取所有投影点满足 X e 且y e 条件所对应的面元，得到与步骤（4)中求得的长方形区域 对应的M个面元； (6) 对于步骤（5)中的M个面元中的第m个面元，求其所在平面与步骤（3)中声源与第 η个面元中心点连线的交点P，求得β i其中A (i = 1，2, 3, r4= r D为第m个面元的第i个顶点位置矢量，rp为交点P的位置 矢量，Vm为第m个面元的法向矢量，β i为求得的判断参数，如果所有的β i>0,则第η个 面元被第m个面元遮挡，反之未被第m个面元遮挡； (7) 重复步骤（6)，判断第η个面元与全部M个面元的遮挡关系，如果均未被遮挡，则保 留第η个面元，反之删除第η个面元； (8) 重复步骤（3)~步骤（7)，起伏界面的N个面元中所有保留下来的面元即为入射声 波照亮的面元，删除的面元即为被遮挡的面元。【专利摘要】本专利技术属于水下目标声散射领域，具体涉及的是一种数值计算随机起伏界面的散射声场时，搜索未被声波照射到的遮挡面元的声波照射随机起伏界面时的遮挡面元快速搜索方法。本专利技术包括：导入已进行三角面元划分的随机起伏界面的N个面元节点信息、面元中心点信息以及声源位置信息；求取所有面元中心点三维坐标中z分量的最大值Zmax，得到z＝Zmax平面等。本专利技术由于判断面元是否被遮挡时的搜索区域减小，减小了计算量，提高了遮挡面元的搜索效率；适用于随机起伏界面的遮挡问题，且起伏界面面积越大，计算效率提高的越明显；只利用了空间几何中点、线、面的有关计算，没有复杂的处理算法，易于程序实现。【IPC分类】G06T7-00, G01S15-88【公开号】CN104766311【申请号】CN201510140918【专利技术人】陈文剑, 孙辉 【申请人】哈尔滨工程大学【公开日】2015年7月8日【申请日】2015年3月27日本文档来自技高网...

【技术保护点】
一种声波照射随机起伏界面时的遮挡面元快速搜索方法，其特征在于：(1)导入已进行三角面元划分的随机起伏界面的N个面元节点信息、面元中心点信息以及声源位置信息；(2)求取所有面元中心点三维坐标中z分量的最大值Zmax，得到z＝Zmax平面；(3)对于第n个面元，连接声源与其中心点，求得连线与z＝Zmax平面的交点(x',y',Zmax)，所求交点在xoy平面的投影点为(x',y')，第n个面元中心点在xoy平面的投影点为(xn,yn)；(4)如果x'≤xn，令x'*＝x'‑Δl，xn*＝xn+Δl；如果x'＞xn，令x'*＝x'+Δl，xn*＝xn‑Δl；如果y'≤yn，令y'*＝y'‑Δl，yn*＝yn+Δl；如果y'＞yn，令y'*＝y'+Δl，yn*＝yn‑Δl，其中Δl为面元边长，在xoy平面上得到顶点分别为(x'*,y'*)、(xn*,y'*)、(xn*,yn*)、(x'*,yn*)的长方形区域；(5)随机起伏界面中面元中心点在xoy平面上投影点为(x,y)，取所有投影点满足x∈[x'*,xn*]且y∈[y'*,yn*]条件所对应的面元，得到与步骤(4)中求得的长方形区域对应的M个面元；(6)对于步骤(5)中的M个面元中的第m个面元，求其所在平面与步骤(3)中声源与第n个面元中心点连线的交点P，求得βi其中ri(i＝1,2,3,r4＝r1)为第m个面元的第i个顶点位置矢量，rp为交点P的位置矢量，vm为第m个面元的法向矢量，βi为求得的判断参数，如果所有的βi＞0，则第n个面元被第m个面元遮挡，反之未被第m个面元遮挡；(7)重复步骤(6)，判断第n个面元与全部M个面元的遮挡关系，如果均未被遮挡，则保留第n个面元，反之删除第n个面元；(8)重复步骤(3)～步骤(7)，起伏界面的N个面元中所有保留下来的面元即为入射声波照亮的面元，删除的面元即为被遮挡的面元。...

【技术特征摘要】

【专利技术属性】
技术研发人员：陈文剑，孙辉，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：黑龙江;23