[发明专利]一种声波照射随机起伏界面时的遮挡面元快速搜索方法

摘要

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

主权项

一种声波照射随机起伏界面时的遮挡面元快速搜索方法，其特征在于：(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为第m个面元的第i个顶点位置矢量，i＝1,2,3,4，r4＝r1，rp为交点P的位置矢量，vm为第m个面元的法向矢量，βi为求得的判断参数，如果所有的βi＞0，则第n个面元被第m个面元遮挡，反之未被第m个面元遮挡；(7)重复步骤(6)，判断第n个面元与全部M个面元的遮挡关系，如果均未被遮挡，则保留第n个面元，反之删除第n个面元；(8)重复步骤(3)～步骤(7)，起伏界面的N个面元中所有保留下来的面元即为入射声波照亮的面元，删除的面元即为被遮挡的面元。