CN117910279A - 基于空间坐标旋转变换的三维声场快速处理方法 - Google Patents
基于空间坐标旋转变换的三维声场快速处理方法 Download PDFInfo
- Publication number
- CN117910279A CN117910279A CN202410302314.0A CN202410302314A CN117910279A CN 117910279 A CN117910279 A CN 117910279A CN 202410302314 A CN202410302314 A CN 202410302314A CN 117910279 A CN117910279 A CN 117910279A
- Authority
- CN
- China
- Prior art keywords
- matrix
- rotation transformation
- dimensional
- array
- horizontal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000009466 transformation Effects 0.000 title claims abstract description 101
- 238000003672 processing method Methods 0.000 title claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims abstract description 180
- 238000012876 topography Methods 0.000 claims abstract description 112
- 238000012821 model calculation Methods 0.000 claims abstract description 22
- 230000007613 environmental effect Effects 0.000 claims abstract description 8
- 230000001427 coherent effect Effects 0.000 claims abstract description 6
- 230000001131 transforming effect Effects 0.000 claims abstract description 4
- 238000000034 method Methods 0.000 claims description 24
- 238000006243 chemical reaction Methods 0.000 claims description 8
- 238000004804 winding Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 description 52
- 238000004088 simulation Methods 0.000 description 4
- 239000000758 substrate Substances 0.000 description 3
- 238000011426 transformation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000013535 sea water Substances 0.000 description 2
- 241000287196 Asthenes Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明涉及基于空间坐标旋转变换的三维声场快速处理方法,该方法,包括:定义元水平阵阵元位置,并计算其轴的倾角;确定三维空间坐标旋转矩阵;设置声源位置并进行三维空间坐标旋转变换;设置水平阵接收深度并进行三维空间坐标旋转变换得到旋转变换后的水平接收阵各阵元位置;确定涵盖水平阵的最小地形区域范围;设定模型计算所需的三维海底地形范围并确定空间坐标旋转变换所需的三维海底地形范围;对坐标网格矩阵和海底地形矩阵同时进行变换;截取旋转变换后的坐标网格矩阵和海底地形矩阵,并进行插值处理;设置环境参数,并计算处理水平阵各阵元位置处的相干声压,以得到处理结果并记录耗费时间。本发明提高了三维声场的处理效率。
Description
技术领域
本发明涉及三维声场处理技术领域,尤其是指基于空间坐标旋转变换的三维声场快速处理方法。
背景技术
目前存在多种三维声场计算模型,例如基于三维射线理论的Bellhop3D模型和Traceo3D模型、基于三维简正波理论的耦合简正波模型、基于三维抛物方程理论的For3D模型等。其中,三维声场计算模型的主要输入为:声源参数、接收端参数、声速剖面、三维海底地形等,通过对所需参数进行配置,可以实现对三维海洋环境中任意位置处复声压的计算。目前三维声场模型在计算中考虑了水平折射等三维效应的存在,使得计算复杂度与计算时间急剧增加。在进行三维声场仿真及阵列信号处理时,通常需要对阵列上各阵元处信号进行计算处理,目前的三维声场计算模型只支持对平行于坐标轴的阵列位置进行配置,若阵列位置与坐标轴之间存在倾角(如图2所示),那么阵元位置不会落在计算网格点上(如图3所示),目前只能对各位置的阵元处信号进行遍历计算处理,这种方式不仅仅增加了环境配置文件的写入时间,也使得每次声场计算处理都要从起始点处重新开始,极大地降低了三维声场的处理效率。
发明内容
本发明的目的在于克服现有技术的不足,提供基于空间坐标旋转变换的三维声场快速处理方法。
为了解决上述技术问题,本发明采用如下技术方案:
本实施例提供了一种基于空间坐标旋转变换的三维声场快速处理方法,包括以下步骤:
定义元水平阵阵元位置,并计算其与/>轴的倾角/>;
根据倾角确定三维空间坐标旋转矩阵,即绕/>轴逆时针旋转倾角/>的旋转矩阵;
设置声源位置,并对声源位置使用旋转矩阵进行三维空间坐标旋转变换,以得到旋转变换后的声源位置;
设置水平阵接收深度,并对水平接收阵各阵元位置使用旋转矩阵进行三维空间坐标旋转变换,以得到旋转变换后的水平接收阵各阵元位置;
根据元水平阵阵元位置确定涵盖水平阵的最小地形区域范围;
根据最小地形区域范围设定模型计算所需的三维海底地形范围,并确定空间坐标旋转变换所需的三维海底地形范围;
根据空间坐标旋转变换所需的三维海底地形范围,对坐标网格矩阵和海底地形矩阵同时进行变换,以得到旋转变换后的坐标网格矩阵和海底地形矩阵;
截取旋转变换后的坐标网格矩阵和海底地形矩阵,以得到截取坐标网格矩阵和截取海底地形矩阵;
对截取坐标网格矩阵和截取海底地形矩阵进行插值处理,以得到插值坐标网格矩阵和插值海底地形矩阵;
设置环境参数,并根据旋转变换后的声源位置、旋转变换后的水平接收阵各阵元位置、插值坐标网格矩阵、插值海底地形矩阵及环境参数计算处理水平阵各阵元位置处的相干声压,以得到处理结果并记录耗费时间。
本发明与现有技术相比的有益效果是:通过水平阵阵元位置确定与坐标轴之间的倾角,然后根据阵列与坐标轴之间的倾角确定三维空间坐标旋转矩阵,将收发位置、海底地形、声速场等参数进行空间坐标旋转变换,得到输入参数后对环境参数进行写入,进而通过三维声场模型读取环境配置文件实现对水平阵各阵元处信号进行快速计算处理,以提高三维声场的处理效率。
下面结合附图和具体实施例对本发明作进一步描述。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的基于空间坐标旋转变换的三维声场快速处理方法的流程示意图;
图2为本发明采用的坐标系以及水平阵与坐标轴之间位置关系示意图;
图3为本发明水平阵位置在平面的投影以及与计算网格的位置关系示意图;
图4为本发明涵盖接收阵列与声源位置的最小地形区域的示意图;
图5为本发明最小地形区域进行旋转变换所得结果的示意图;
图6为本发明最终确定的海底地形区域旋转变换所得结果的示意图;
图7为本发明最终确定的未进行旋转变换的海底地形区域的示意图;
图8为本发明进行截取与插值处理之后的海底地形示意图;
图9为本发明仿真海域的海水声速剖面的示意图;
图10为本发明的方法与常规方法传播损失计算结果对比示意图。
具体实施方式
请参阅图1至图10所示的具体实施例,本发明公开了一种基于空间坐标旋转变换的三维声场快速处理方法,包括以下步骤:
S1,定义元水平阵阵元位置,并计算其与/>轴的倾角/>;
具体地,水平阵各阵元位置及其与坐标轴之间的位置关系如下:
;
其中,表示/>平面内水平阵第/>个/>阵元投影的位置,表示/>平面内水平阵第1个阵元投影的位置。
S2,根据倾角确定三维空间坐标旋转矩阵,即绕/>轴逆时针旋转倾角/>的旋转矩阵/>;
具体地,旋转矩阵与倾角/>的关系如下:
;
其中,使用旋转矩阵对所需各输入参数进行处理后,得到的各输入参数与/>轴平行,水平阵各阵元位置落在平行于/>轴的计算网格线上。
S3,设置声源位置,并对声源位置使用旋转矩阵进行三维空间坐标旋转变换,以得到旋转变换后的声源位置;
具体地,旋转变换后的声源位置与旋转变换前的声源位置关系如下:
;
其中,表示旋转变换前的声源位置,/>表示旋转变换后的声源位置。
S4,设置水平阵接收深度,并对水平接收阵各阵元位置使用旋转矩阵进行三维空间坐标旋转变换,以得到旋转变换后的水平接收阵各阵元位置;
具体地,旋转变换后的水平接收阵各阵元位置与旋转变换前的水平接收阵各阵元位置关系如下:
;
其中,表示旋转变换前的水平接收阵各阵元位置,表示旋转变换后的水平接收阵各阵元位置。
S5,根据元水平阵阵元位置确定涵盖水平阵的最小地形区域范围;
具体地,涵盖水平阵的最小地形区域范围与元水平阵阵元位置关系如下:
;
其中,为涵盖水平阵的最小地形区域的/>轴方向范围,/>为涵盖水平阵的最小地形区域的/>轴方向范围。
另外,若此时声源位置不在上述区域范围内,那么需要根据下式确定同时涵盖水平阵和声源的最小地形区域范围:
;
其中,表示最远的阵元位置,/>表示距离声源位置。
S6,根据最小地形区域范围设定模型计算所需的三维海底地形范围,并确定空间坐标旋转变换所需的三维海底地形范围;
具体地,三维声场计算模型的海底地形输入格式为的矩阵,/>和/>分别代表/>轴方向和/>轴方向划分的距离网格点数目,而且地形数据中不能出现大量的(非数值符号)。上面确定了旋转变换前涵盖水平阵的最小地形区域范围为/>,为了确保海底地形输入矩阵中包含该区域旋转变换后对应的区域,那么模型计算所需的海底地形矩形区域的范围至少为:
;
其中,为模型计算所需的海底地形的/>轴方向范围,/>为模型计算所需的海底地形的/>轴方向范围。
另外,确定空间坐标旋转变换所需的三维海底地形范围,模型计算所需的矩形海底地形区域是三维空间坐标旋转变换之后的结果,为了使得变换后的区域中包含模型计算所需的海底地形区域,那么未变换的海底地形范围至少为:
;
其中,为空间坐标旋转变换所需海底地形的/>轴方向范围,/>为空间坐标旋转变换所需海底地形的/>轴方向范围。
S7,根据空间坐标旋转变换所需的三维海底地形范围,对坐标网格矩阵和海底地形矩阵同时进行变换,以得到旋转变换后的坐标网格矩阵和海底地形矩阵;
具体地,三维空间坐标旋转变换所需的地形范围为,设对应的海底地形矩阵为/>,阶数为/>,此时,将海底地形矩阵对应区域的/>坐标进行网格化,获得阶的坐标网格矩阵/>和/>,再将海底地形矩阵与两个坐标网格矩阵进行三维空间坐标旋转变换,以得到旋转变换后的坐标网格矩阵和海底地形矩阵,其表达式为:
;
其中,表示旋转变换之前的/>轴方向坐标网格矩阵,/>表示旋转变换之前的轴方向坐标网格矩阵,/>表示旋转变换之前的海底地形矩阵,/>表示旋转变换之后的/>轴方向坐标网格矩阵,/>表示旋转变换之后的/>轴方向坐标网格矩阵,/>表示旋转变换之后的海底地形矩阵。
S8,截取旋转变换后的坐标网格矩阵和海底地形矩阵,以得到截取坐标网格矩阵和截取海底地形矩阵;
具体地,过大的地形范围与输入矩阵阶数,通常会超出模型计算能力的限制,因此需要在旋转变换之后的阶矩阵中截取模型计算所需的/>阶矩阵,以得到截取坐标网格矩阵和截取海底地形矩阵。
首先确定模型计算所需的矩阵所对应区域的四个顶点坐标:
其中,表示左下顶点坐标,/>表示左上顶点坐标,/>表示右上顶点坐标,/>表示右下顶点坐标。然后根据以下判别准则截取模型计算所需的矩阵:
;
即根据选定的索引值提取同时满足横纵坐标判别准则的矩阵中元素,其中,为截取的/>轴方向坐标网格矩阵,/>为截取的/>轴方向坐标网格矩阵,/>为截取的海底地形矩阵。
需注意的是,若是倾角为45°,那么截取的矩阵阶数也应是:
;
此时会在奇数列的偶数行与偶数列的奇数行都出现空值,使用进行补齐,但是此时的/>并不会影响到下一步的插值处理;若是倾角/>不为45°,那么此时提取出的矩阵元素不会出现矩阵式的规则排布,此时设定截取的矩阵阶数/>为,即将矩阵元素按行排列得到截取坐标网格矩阵和截取海底地形矩阵。
S9,对截取坐标网格矩阵和截取海底地形矩阵进行插值处理,以得到插值坐标网格矩阵和插值海底地形矩阵;
具体地,三维海底地形的精度直接影响声场的计算结果,所以需要在计算模型允许的范围内尽可能对地形做精细化处理,根据模型计算所需海底地形区域范围选定合适的间隔距离,以保证不超过模型计算限制的同时获得最大阶数的等间隔插值坐标网格矩阵和/>,选定二维线性插值的方式对截取海底地形矩阵/>进行插值,首先设定待插值点/>,然后找到最邻近的插值区间/>和/>,/>和都为插值坐标网格矩阵/>和/>中的元素,然后计算待插值点处的值:
;
其中,表示插值处理之后的海底地形矩阵。
在一实施例中,在S9步骤之后,还包括:对声速场进行坐标旋转变换;若声场计算环境为匀速水体或声速在水平方向上无变化(各位置点采用同一声速剖面),那么不需要对声速场进行旋转变换;若声速在水平方向上异速,那么需要将声速场进行旋转变换,声速场旋转变换公式如下:
;
其中,表示未进行旋转变换的声速场/>轴方向位置坐标,/>表示未进行旋转变换的声速场/>轴方向位置坐标,/>表示未进行旋转变换的声速场/>轴方向位置坐标,/>表示未进行旋转变换的三维声速场中对应位置声速,/>表示旋转变换之后的声速场/>轴方向位置坐标,/>表示旋转变换之后的声速场/>轴方向位置坐标,/>表示旋转变换之后的声速场/>轴方向位置坐标,/>表示旋转变换之后的三维声速场中对应位置声速。
S10,设置环境参数,并根据旋转变换后的声源位置、旋转变换后的水平接收阵各阵元位置、插值坐标网格矩阵、插值海底地形矩阵及环境参数计算处理水平阵各阵元位置处的相干声压,以得到处理结果并记录耗费时间。
具体地,环境参数包括声源频率 、海面参数及底质参数。根据旋转变换后的声源位置、旋转变换后的水平接收阵各阵元位置、插值坐标网格矩阵、插值海底地形矩阵、声源频率/>、海面参数及底质参数,然后调用三维声场计算模型对水平阵各阵元位置处的相干声压/>进行计算处理,以得到处理结果并记录耗费时间。其中,调用三维声场计算模型对水平阵各阵元位置处的相干声压/>进行计算处理,此为现有公开技术,在此不做具体阐述。
本发明针对现有三维声场计算模型不支持与坐标轴存在倾角的水平阵位置配置,导致对各阵元处信号的计算处理时间急剧增加的问题,提出基于空间坐标旋转变换的三维声场计算处理方法,通过计算水平阵与坐标轴之间的倾角,获得三维空间坐标旋转矩阵,进而对环境文件中所需的声源位置、水平接收阵阵元位置、三维海底地形等参数进行旋转变换,最终通过三维声场计算模型读取配置环境参数,实现对水平阵各阵元处信号的快速计算处理,极大提升对水平接收阵列信号的计算处理效率。
为了更便于理解本发明的技术方案,在此特提供以下具体实施例:
仿真环境选取南海北部大陆坡海域,海底地形来自ETOPO1数据集,为了使各阵元处信号之间存在明显变化,从而更好验证本方法的精确性,设置23元水平接收阵的轴坐标范围1km-45km,各阵元之间/>轴方向间隔2km,/>轴坐标范围为1km-45km,各阵元之间轴间隔2km,水平阵位置如图3所示,计算与/>轴的倾角/>:
;
根据水平阵与轴的倾角/>,确定绕/>轴逆时针旋转/>使得各参数与/>轴平行的旋转矩阵/>:
;
设置声源位置为(500m,500m,200m),如图3所示,对声源位置使用旋转矩阵进行三维空间坐标旋转变换:
;
其中,表示使用/>旋转变换后的声源位置。
设置水平阵接收深度为600m,对水平接收阵各阵元位置使用旋转矩阵进行三维空间坐标旋转变换:
;
其中,表示使用/>旋转变换后的水平接收阵阵元位置。
确定涵盖水平阵和声源的最小地形区域的范围:
;
为了方便后续计算,取该范围为,如图4所示。
确定模型计算所需的海底地形区域范围,若直接对上面确定的海底地形区域(即)进行旋转变换,那么就会得到如图5所示的结果。三维声场计算模型以矩阵的形式对海底地形进行输入,且地形矩阵中不能出现大量的/>值(图5中空白部分在矩阵中需要用/>表示),为了保证海底地形输入矩阵包含所确定的区域旋转变换后的结果(图6中外层虚线所框画的区域),同时空白部分以实际地形数值的形式表示,那么模型计算所需输入矩阵对应的海底地形区域(图6中外层虚线所框画的区域)的范围至少为:
;
其中,表示三维声场模型计算所需的地形区域范围。
进一步确定旋转变换前所需的海底地形区域范围。为了保证如图6所示旋转变换后的结果中包含模型计算所需输入矩阵对应的海底地形区域,那么旋转变换后的海底地形区域范围至少为:
;
其中,表示使用/>进行旋转变换后的海底地形范围,旋转变换并不改变地形区域的大小,进而得到如图7所示的旋转变换前所需的海底地形区域范围。
根据所确定的旋转变换前海底地形范围,从ETOPO1数据集中对该区域三维海底地形进行提取,ETOPO1数据集的精度为1分,对提取出的海底地形进行插值处理,插值之后海底地形矩阵的阶数为181x181,/>轴方向和/>轴方向的距离间隔为500m,将两方向的坐标处理为与/>对应的181x181阶坐标网格矩阵/>和/>,使用/>和/>对三个矩阵中的元素进行空间坐标旋转变换:
;
其中,表示使用/>旋转变换之后的坐标网格矩阵与海底地形矩阵中的元素。
在旋转变换之后的海底地形矩阵中截取模型计算所需的输入矩阵;其中,要截取包含上面所确定的海底地形区域的矩阵,首先确定模型计算所需的矩阵所对应区域的四个顶点坐标:
其中,表示左下顶点坐标,/>表示左上顶点坐标,/>表示右上顶点坐标,/>表示右下顶点坐标。然后根据以下判别准则截取模型计算所需的矩阵:
;
即根据选定的索引值提取同时满足横纵坐标判别准则的矩阵中元素,其中,为截取的/>轴方向坐标网格矩阵,/>为截取的/>轴方向坐标网格矩阵,/>为截取的海底地形矩阵。
需注意的是,若是倾角为45°,那么截取的矩阵阶数也应是:
;
此时会在奇数列的偶数行与偶数列的奇数行都出现空值,使用进行补齐,但是此时的/>并不会影响到下一步的插值处理。
将截取的坐标网格和海底地形进行插值处理,使模型所需输入的矩阵进行精细化处理,保证三维声场模型的计算精确度。
根据模型计算所需海底地形区域范围选定间隔距离为 ,那么可得等间隔插值坐标网格矩阵 /> 和 />的阶数为/>,然后使用二维线性插值法对截取海底地形矩阵/>进行插值,首先设定待插值点/>,其中/>,然后找到最邻近的插值区间/>和/>,/>和/>都为插值坐标网格矩阵/>和/>中的元素,然后计算待插值点处的值:
;
其中,表示插值处理之后的海底地形矩阵。其中,进行截取和插值处理之后的矩阵对应的海底地形如图8所示。
仿真所需声速数据来自Argo观测数据,根据所选区域中心的经纬度提取声速剖面,将水平方向上声速设置为匀速不变,即所有位置点共用同一声速剖面,如图9所示,此时不需要对声速场进行旋转变换处理。
将声源频率设置为1kHz,海面以上设置为真空,底质参数设置为,其中,/> 表示海底纵波声速,/>表示海底横波声速,/> 表示海底底质密度,/> 表示海底纵波衰减系数,/>表示海底横波衰减系数,将以上参数与旋转变换后的声源位置/>、水平接收阵阵元位置 />、模型所需的输入矩阵/>、/>、/>分别写入环境配置文件,再使用Traceo3D模型分别计算落于平行于/>轴计算网格线上各阵元处接收信号/>,并纪录三维声场模型的计算处理时间。
另外,为了验证本发明的快速性,对水平阵23个阵元处信号进行遍历计算,即对三维声场计算模型进行23次调用,并统计计算时间;为了进一步探索阵元数对本发明计算处理效率的影响,在水平阵每两个阵元中心点处再增加一个阵元,即阵元间隔由原来的变成 />km,阵元数从23增加到45,统计条件下本发明的计算处理时间与常规遍历计算处理方法的时间,最终将不同条件以及不同方法的计算处理时间进行比较,具体如图10所示。
另外,采用不同计算处理方法及阵元数的三维声场计算时间对比,如下表所示:
所选计算处理方法与阵元数 | 23元水平阵旋转变换方法 | 23元水平阵遍历计算处理方法 | 45元水平阵旋转变换方法 | 45元水平阵遍历计算处理方法 |
三维声场软件计算时间 | 393.86s | 9090.05s | 397.21s | 17801.53s |
由上表可知,本发明的方法会显著提升三维声场的计算处理效率,对23元水平阵的计算处理效率提升23倍左右,对45元水平阵的计算效率提升45倍左右,即水平阵阵元数越多,本发明对于计算处理效率的提升效果越明显。
本发明的基本原理和实施方案经过了计算机数值仿真验证,其结果表明:利用本发明提出的基于三维空间坐标旋转变换的声场计算方法获得的接收阵阵元处声压与遍历计算获得的阵元处声压完全一致,在保证声场计算精确性的同时可以有效降低三维声场模型的计算时间,极大提升对接收水平阵阵元处信号的计算处理效率。
为便于理解附图,在此做以下说明:图2给出了本发明采用的右手坐标系,以及展示了水平阵与坐标轴之间的位置关系;图3展示了水平阵位置在平面的投影以及与计算网格的位置关系;图4展示了涵盖接收阵列与声源位置的最小地形区域;图5展示了最小地形区域进行旋转变换的结果;图6展示了既包含最小地形区域旋转变换的结果(里层虚线框画区域)又包含了符合模型计算所需地形输入矩阵形式的区域(外层虚线框画区域)的海底地形旋转变换之后的结果;图7展示了既包含最小地形区域旋转变换的结果(里层虚线框画区域)又包含了符合模型计算所需地形输入矩阵形式的区域(外层虚线框画区域)的海底地形未进行旋转变换之后的结果;图8展示了在海底地形进行旋转变换之后,对模型计算所需部分进行截取及插值处理之后所得的结果;图9展示了所选海域的海水声速剖面;图10展示了使用三维空间坐标旋转变换方法进行声场计算与未使用该方法而进行逐点遍历声场计算的结果比较。
上述实施例为本发明较佳的实现方案,除此之外,本发明还可以其它方式实现,在不脱离本技术方案构思的前提下任何显而易见的替换均在本发明的保护范围之内。
Claims (10)
1.基于空间坐标旋转变换的三维声场快速处理方法,其特征在于,包括以下步骤:
定义元水平阵阵元位置,并计算其与/>轴的倾角/>;
根据倾角确定三维空间坐标旋转矩阵,即绕/>轴逆时针旋转倾角/>的旋转矩阵/>;
设置声源位置,并对声源位置使用旋转矩阵进行三维空间坐标旋转变换,以得到旋转变换后的声源位置;
设置水平阵接收深度,并对水平接收阵各阵元位置使用旋转矩阵进行三维空间坐标旋转变换,以得到旋转变换后的水平接收阵各阵元位置;
根据元水平阵阵元位置确定涵盖水平阵的最小地形区域范围;
根据最小地形区域范围设定模型计算所需的三维海底地形范围,并确定空间坐标旋转变换所需的三维海底地形范围;
根据空间坐标旋转变换所需的三维海底地形范围,对坐标网格矩阵和海底地形矩阵同时进行变换,以得到旋转变换后的坐标网格矩阵和海底地形矩阵;
截取旋转变换后的坐标网格矩阵和海底地形矩阵,以得到截取坐标网格矩阵和截取海底地形矩阵;
对截取坐标网格矩阵和截取海底地形矩阵进行插值处理,以得到插值坐标网格矩阵和插值海底地形矩阵;
设置环境参数,并根据旋转变换后的声源位置、旋转变换后的水平接收阵各阵元位置、插值坐标网格矩阵、插值海底地形矩阵及环境参数计算处理水平阵各阵元位置处的相干声压,以得到处理结果并记录耗费时间。
2.根据权利要求1所述的基于空间坐标旋转变换的三维声场快速处理方法,其特征在于,所述定义元水平阵阵元位置,并计算其与/>轴的倾角/>步骤中,水平阵各阵元位置及其与坐标轴之间的位置关系如下:
;
其中,表示/>平面内水平阵第/>个/>阵元投影的位置,/>表示/>平面内水平阵第1个阵元投影的位置。
3.根据权利要求2所述的基于空间坐标旋转变换的三维声场快速处理方法,其特征在于,所述根据倾角确定三维空间坐标旋转矩阵,即绕/>轴逆时针旋转倾角/>的旋转矩阵步骤中,旋转矩/>与倾角/>的关系如下:
。
4.根据权利要求3所述的基于空间坐标旋转变换的三维声场快速处理方法,其特征在于,所述设置声源位置,并对声源位置使用旋转矩阵进行三维空间坐标旋转变换,以得到旋转变换后的声源位置步骤中,旋转变换后的声源位置与旋转变换前的声源位置关系如下:
;
其中,表示旋转变换前的声源位置,/>表示旋转变换后的声源位置。
5.根据权利要求4所述的基于空间坐标旋转变换的三维声场快速处理方法,其特征在于,所述设置水平阵接收深度,并对水平接收阵各阵元位置使用旋转矩阵进行三维空间坐标旋转变换,以得到旋转变换后的水平接收阵各阵元位置步骤中,旋转变换后的水平接收阵各阵元位置与旋转变换前的水平接收阵各阵元位置关系如下:
;
其中,表示旋转变换前的水平接收阵各阵元位置,表示旋转变换后的水平接收阵各阵元位置。
6.根据权利要求5所述的基于空间坐标旋转变换的三维声场快速处理方法,其特征在于,所述根据元水平阵阵元位置确定涵盖水平阵的最小地形区域范围步骤中,涵盖水平阵的最小地形区域范围与/>元水平阵阵元位置关系如下:
;
其中,为涵盖水平阵的最小地形区域的/>轴方向范围,/>为涵盖水平阵的最小地形区域的/>轴方向范围。
7.根据权利要求6所述的基于空间坐标旋转变换的三维声场快速处理方法,其特征在于,所述根据最小地形区域范围设定模型计算所需的三维海底地形范围,并确定空间坐标旋转变换所需的三维海底地形范围步骤中,模型计算所需的三维海底地形范围与最小地形区域范围关系如下:
;
其中,为模型计算所需的海底地形的/>轴方向范围,/>为模型计算所需的海底地形的/>轴方向范围。
8.根据权利要求7所述的基于空间坐标旋转变换的三维声场快速处理方法,其特征在于,所述根据空间坐标旋转变换所需的三维海底地形范围,对坐标网格矩阵和海底地形矩阵同时进行变换,以得到旋转变换后的坐标网格矩阵和海底地形矩阵步骤中,三维空间坐标旋转变换所需的地形范围为,设对应的海底地形矩阵为/>,阶数为/>,此时,将海底地形矩阵对应区域的/>坐标进行网格化,获得/>阶的坐标网格矩阵和/>,再将海底地形矩阵与两个坐标网格矩阵进行三维空间坐标旋转变换,以得到旋转变换后的坐标网格矩阵和海底地形矩阵,其表达式为:
;
其中,表示旋转变换之前的/>轴方向坐标网格矩阵,/>表示旋转变换之前的/>轴方向坐标网格矩阵,/>表示旋转变换之前的海底地形矩阵,/>表示旋转变换之后的/>轴方向坐标网格矩阵,/>表示旋转变换之后的/>轴方向坐标网格矩阵,/>表示旋转变换之后的海底地形矩阵。
9.根据权利要求8所述的基于空间坐标旋转变换的三维声场快速处理方法,其特征在于,所述截取旋转变换后的坐标网格矩阵和海底地形矩阵,以得到截取坐标网格矩阵和截取海底地形矩阵步骤中,在旋转变换之后的阶矩阵中截取模型计算所需的阶矩阵,以得到截取坐标网格矩阵和截取海底地形矩阵。
10.根据权利要求9所述的基于空间坐标旋转变换的三维声场快速处理方法,其特征在于,所述对截取坐标网格矩阵和截取海底地形矩阵进行插值处理,以得到插值坐标网格矩阵和插值海底地形矩阵步骤中,根据模型计算所需海底地形区域范围选定合适的间隔距离,以保证不超过模型计算限制的同时获得最大阶数的等间隔插值坐标网格矩阵和,选定二维线性插值的方式对截取海底地形矩阵/>进行插值,首先设定待插值点,然后找到最邻近的插值区间/>和/>,/>和都为插值坐标网格矩阵/>和/>中的元素,然后计算待插值点处的值:
;
其中,表示插值处理之后的海底地形矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410302314.0A CN117910279B (zh) | 2024-03-18 | 2024-03-18 | 基于空间坐标旋转变换的三维声场快速处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410302314.0A CN117910279B (zh) | 2024-03-18 | 2024-03-18 | 基于空间坐标旋转变换的三维声场快速处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117910279A true CN117910279A (zh) | 2024-04-19 |
CN117910279B CN117910279B (zh) | 2024-05-28 |
Family
ID=90692767
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410302314.0A Active CN117910279B (zh) | 2024-03-18 | 2024-03-18 | 基于空间坐标旋转变换的三维声场快速处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117910279B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH02112891A (ja) * | 1988-10-24 | 1990-04-25 | Fanuc Ltd | 3次元レーザ加工方法 |
CN101872490A (zh) * | 2010-07-02 | 2010-10-27 | 北京理工大学 | 一种改进的坐标轴方向突变下的三维模型变形方法 |
CN106441543A (zh) * | 2016-12-09 | 2017-02-22 | 华南理工大学 | 基于三维正交阵的水下探测路径声速测量方法及装置 |
JP2018063200A (ja) * | 2016-10-14 | 2018-04-19 | 日本電信電話株式会社 | 音源位置推定装置、音源位置推定方法、及びプログラム |
CN112198476A (zh) * | 2020-10-16 | 2021-01-08 | 昆明理工大学 | 一种基于立体视觉和波束形成的移动声源的三维定位方法 |
CA3149297A1 (en) * | 2019-07-29 | 2021-02-04 | Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. | Apparatus, method or computer program for processing a sound field representation in a spatial transform domain |
WO2021031781A1 (zh) * | 2019-08-19 | 2021-02-25 | 深圳光峰科技股份有限公司 | 投影图像校准方法、装置及投影设备 |
CN113641954A (zh) * | 2021-07-20 | 2021-11-12 | 中国科学院声学研究所 | 一种复杂海洋环境中的三维声场快速预报方法及系统 |
CN114286279A (zh) * | 2022-01-18 | 2022-04-05 | 中国科学院声学研究所 | 基于全局矩阵耦合简正波的三维声场模型的声场获取方法 |
JP7116852B1 (ja) * | 2021-07-08 | 2022-08-10 | 自然資源部第一海洋研究所 | 並列結合海洋音響予測システム及び操作方法 |
CN115327626A (zh) * | 2022-08-11 | 2022-11-11 | 中国石油大学(北京) | 一种用于三维ati介质的弹性波场矢量分解方法及系统 |
-
2024
- 2024-03-18 CN CN202410302314.0A patent/CN117910279B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH02112891A (ja) * | 1988-10-24 | 1990-04-25 | Fanuc Ltd | 3次元レーザ加工方法 |
CN101872490A (zh) * | 2010-07-02 | 2010-10-27 | 北京理工大学 | 一种改进的坐标轴方向突变下的三维模型变形方法 |
JP2018063200A (ja) * | 2016-10-14 | 2018-04-19 | 日本電信電話株式会社 | 音源位置推定装置、音源位置推定方法、及びプログラム |
CN106441543A (zh) * | 2016-12-09 | 2017-02-22 | 华南理工大学 | 基于三维正交阵的水下探测路径声速测量方法及装置 |
CA3149297A1 (en) * | 2019-07-29 | 2021-02-04 | Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. | Apparatus, method or computer program for processing a sound field representation in a spatial transform domain |
WO2021031781A1 (zh) * | 2019-08-19 | 2021-02-25 | 深圳光峰科技股份有限公司 | 投影图像校准方法、装置及投影设备 |
CN112198476A (zh) * | 2020-10-16 | 2021-01-08 | 昆明理工大学 | 一种基于立体视觉和波束形成的移动声源的三维定位方法 |
JP7116852B1 (ja) * | 2021-07-08 | 2022-08-10 | 自然資源部第一海洋研究所 | 並列結合海洋音響予測システム及び操作方法 |
CN113641954A (zh) * | 2021-07-20 | 2021-11-12 | 中国科学院声学研究所 | 一种复杂海洋环境中的三维声场快速预报方法及系统 |
CN114286279A (zh) * | 2022-01-18 | 2022-04-05 | 中国科学院声学研究所 | 基于全局矩阵耦合简正波的三维声场模型的声场获取方法 |
CN115327626A (zh) * | 2022-08-11 | 2022-11-11 | 中国石油大学(北京) | 一种用于三维ati介质的弹性波场矢量分解方法及系统 |
Non-Patent Citations (4)
Title |
---|
BOXUAN ZHANG; XIANGHAO HOU; YIXIN YANG; LONG YANG; YONG WANG: "A Fast Variational Bayesian Adaptive Extended Kalman Filter for Robust Underwater Direction-of-Arrival Tracking", IEEE, 16 May 2023 (2023-05-16) * |
JIANBO ZHOU; SHENGCHUN PIAO; YIWANG HUANG; XIAOMAN LI; HANHAO ZHU: "The effect of sea-surface roughness on ambient noise vertical directionality in shallow water waveguide", IEEE, 8 August 2016 (2016-08-08) * |
胡成博: "基于声学超构材料的空间与频域声波操控研究", 工程科技辑, 15 May 2022 (2022-05-15) * |
马晶鑫: "海底地形地貌与浅地层剖面一体化声学探测关键技术研究", 基础科学, 15 January 2023 (2023-01-15) * |
Also Published As
Publication number | Publication date |
---|---|
CN117910279B (zh) | 2024-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113259034B (zh) | 一种并行的耦合海洋声学预报系统及运行方法 | |
CN109636912B (zh) | 应用于三维声呐图像重构的四面体剖分有限元插值方法 | |
CN104637089A (zh) | 三维模型数据处理方法和装置 | |
CN110456417A (zh) | 一种地震数据多次波压制方法 | |
CN113255230B (zh) | 基于mq径向基函数的重力模型正演方法及系统 | |
CN106199519A (zh) | 一种超短基线五基元立体空间基阵及其水声定位方法 | |
CN113267754A (zh) | 一种基于三维网格的地形遮挡雷达探测范围计算方法 | |
CN111596319B (zh) | 一种地形遮挡对gnss干扰源作用区域影响的高效仿真算法 | |
CN113267822A (zh) | 基于地形约束因子权重优化提高海底地形反演精度的方法 | |
CN114925496A (zh) | 一种海洋环境噪声预报方法 | |
CN115630446A (zh) | 用于潜航器结构低频辐射噪声实时仿真的快速推演方法 | |
CN117910279B (zh) | 基于空间坐标旋转变换的三维声场快速处理方法 | |
CN112698337B (zh) | 一种宽带三维成像声呐稀疏布阵方法 | |
CN112632818B (zh) | 高阶梯度光滑无网格方法、系统、介质、计算机设备 | |
CN111369604B (zh) | 一种基于网格的地形特征点提取方法及处理终端 | |
CN112926231A (zh) | 一种基于等效源法的有限空间中近场声全息测量方法 | |
CN110287620B (zh) | 适用于地表观测面的球坐标系密度界面正演方法及系统 | |
CN116797747A (zh) | 水下检测数据可视化方法、装置、计算机设备和存储介质 | |
CN115061089A (zh) | 一种声源定位方法、系统、介质、设备及装置 | |
Lv et al. | The application of a complex composite fractal interpolation algorithm in the seabed terrain simulation | |
CN105869209A (zh) | 三维地质表面模型中的畸形三角形数据处理方法 | |
CN112230277A (zh) | 基于柱坐标系的隧道地震波传播数值模拟方法及系统 | |
CN116167303B (zh) | 一种用于直升机旋翼流固耦合模拟的曲面网格插值方法 | |
Chen et al. | Fast extraction of local underwater terrain features for underwater terrain-aided navigation | |
CN117095134B (zh) | 一种三维海洋环境数据插值处理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |