【摘 要】
:
声波照射到水中目标会产生散射,不同方向的散射波声压分布与入射波性质以及目标的大小、形状、声学参数密切相关。已知入射声波、目标几何形状和周围介质的声学参数,计算散射声场分布特性是声学正问题。求解散射问题可采用解析法或数值法,但仅有少数形状规则的物体可用解析法获得场的精确解。不规则形状以及复杂材料目标的散射声场可用时域有限差分法(finite difference timedomain, FDTD)、
论文部分内容阅读
声波照射到水中目标会产生散射,不同方向的散射波声压分布与入射波性质以及目标的大小、形状、声学参数密切相关。已知入射声波、目标几何形状和周围介质的声学参数,计算散射声场分布特性是声学正问题。求解散射问题可采用解析法或数值法,但仅有少数形状规则的物体可用解析法获得场的精确解。不规则形状以及复杂材料目标的散射声场可用时域有限差分法(finite difference timedomain, FDTD)、有限元法(finite element method, FEM)、边界元(boundaryelement method, BEM)等基于网格剖分的数值方法求解,这些方法计算复杂度高,耗费时间长。逆向运用Fourier衍射定理解决声学正问题,无需网格剖分和迭代运算即可实现目标散射特性的快速预报。考虑浸没水中的无限长柱体目标,根据目标几何形状和周围介质的声学参数,建立图像模型用以描述二维声场,在图像频谱的适当位置取样可实现散射远场特性的快速预报。本论文研究基于Fourier衍射定理预报水中目标的声散射问题,围绕预报精度和计算效率的提高,取得了以下研究成果:1.通过极坐标系矩形网格表示中的平移实现目标旋转为获得声波以任意方向入射时目标产生的前后向散射声场,可保持入射方向不变而旋转目标,在不同角度下求得声波的广义投影。要得到足够精细的散射波方向特性,必须减小旋转角度间隔或提高图像采样率,从而使计算量剧增。通过将描述二维声场图像的直角坐标形式改变为极坐标下的矩形网格表示,把目标旋转中运算量较大的三角函数乘法转换成平移,大大提高了计算速度,并保证旋转精度不变。算法为散射声场快速、准确预报研究创造了条件。2.图像频谱高精度数值计算高精度频谱是逆向运用Fourier衍射定理预报散射声场的关键。受采样网格限制,二值图像在目标边缘存在锯齿,用2D-FFT得到的离散频谱与连续Fourier变换结果相比误差很大,高频部分甚至完全不可用,因此难以直接用于散射声场预报。根据Fourier切片定理,空间图像某一方向Radon投影的1D-DFT等于目标图像二维频谱对应方向切片上的样本值。由此提出频谱修正方法:对图像空间的Radon投影重复进行“插值-滤波-抽取”操作,经过多次循环获得平滑的投影,同时利用极坐标插值法将投影的1D-DFT结果内插到二维直角坐标系,得到的2D离散频谱精度大大提高。修正的图像频谱保证了精确可靠的散射声场预报。3.改进一阶Born近似提高散射声场指向特性预报精度目标受平面波照射时,在其断面图像频谱的圆形轨迹上取样,可获得散射声压在0°360°的指向图,圆半径等于入射波波数k0。该方法建立在一阶Born近似的基础上,仅适用于目标和水介质的声阻抗之比接近于1,即弱散射的情况。当阻抗比偏离1较远时误差显著增大。针对这一问题,对基于一阶Born近似的方法进行改进。考虑散射体内外波数差异引起的幅度和相位变化,引入修正项对频域采样圆弧的半径和圆心位置进行调整。根据改进的Born近似计算无限长圆柱目标的散射远场指向图,对于目标与水的声阻抗比在0.761.32的情况得到了与解析解相符的分布。对其它形状截面的柱状目标也得到了合理的计算结果。4.引入二阶Born近似,拓展Fourier衍射定理适用范围为进一步提高预报精度,并拓展散射声场快速预报方法的适用范围,考虑散射波在积分方程中的贡献,引入二阶Born近似。先用入射平面声波代替总场求解积分方程,得到散射波的一阶近似。将一阶散射波作为二次入射源引入散射积分方程,求得散射声场的增量(修正项),一阶散射与增量声场之和即为散射波的二阶近似解。散射增量的数值计算通过两步在频域一系列圆弧上采样实现,因而仍可沿用基于Fourier衍射定理的快速预报框架。计算结果表明,二阶Born近似预报精度明显提高,并进一步将散射声场快速预报方法的适用范围扩展到声阻抗比为0.721.44。
其他文献
本文主要研究了四元代数上的一些矩阵方程组。论文中给出了矩阵方程组有解的充分必要条件。当方程组有解的时候,精确解也被构造出来。矩阵方程已经在矩阵理论及其应用占有骨干地位。矩阵方程不仅在现代数学有强大的影响力,而且他们必然会影响到很多领域,例如模型简化,系统理论,最优控制,稳定性分析,控制理论,滤波技术,最优的极点配置,H-最优控制,轨道理论,线性系统,神经网络,敏感性分析,信号处理,反馈,微扰理论,
网络化Lagrange系统(Lagrange网络)的同步与控制已经成为非线性动力学与控制领域中关注的热点之一。这主要是因为经典的Lagrange方程可以用来描述诸如机械臂、飞行器及行走机器人等大量的力学系统,而Lagrange网络具有更为广泛的工程应用,尤其是在复杂与集成化的生产过程中,如多机械臂协调、飞行器编队及移动传感网络等。本文从非线性动力学与控制的角度,研究了Lagrange网络的实用同步
本文首先研究基于Kac-Moody代数s1(2,C[λ-1,λ)与so(3,R)所得新的谱问题.从so(3,R)出发,借助半单Lie代数的理论,从扩展零曲率方程得到相关矩阵谱问题的Dirac孤子方程族.在一类非半单圈代数中,结合变分恒等式构造一族可积耦合的Dirac孤子方程族,给出相应的Hamilton结构及Liouville可积性.新的谱问题可以视为AKNS谱问题的推广,由此引出耦合Burger
线性锥规划是目前优化领域中最热门的研究课题之一. Nemirovskii在2006年国际数学家大会一小时报告[72]中指出,锥规划是近20年凸优化研究工作中的一项重大突破性进展,主要包括线性规划,二阶锥规划和半正定锥规划这三类常见的对称锥规划.这三类规划可以转化成统一的线性锥规划模型.线性锥规划不仅具有良好的结构特征和对偶理论,而且具备很强的数学表达能力和实用性.许多应用领域中的实际问题都可以描述
本论文的研究内容属于凸几何分析,内容涉及仿射等周不等式与OrliczBrunn-Minkowski理论Brunn-Minkowski理论是凸几何分析的核心内容,它的发展起源于对两个凸体与它们的Minkowski和之间的体积关系的研究,由此产生了著名的Brunn-Minkowski不等式,整个Brunn-Minkowski理论的奠基石.仿射等周不等式主要研究凸体的一些仿射量的极值问题,是凸几何分析理
本文主要研究离散可积方程的Lax对在研究可积特征方面的应用.建立利用方程的Lax对获得方程的无穷多守恒律和对称的方法.这些离散的可积方程包括:·Adler-Bobenko-Suris方程,Nijhoff-Quispel-Capel方程;·离散Boussinesq方程,离散非线性Schrodinger方程等多分量方程.论文首先对Adler-Bobenko-Suris链方程、Nijhoff-Quisp
无网格方法是继有限元法之后发展起来的一种新型数值方法,因其形函数构造仅需要问题所在区域或边界的节点信息,不需要形成区域或边界网格,因而具有前处理简单、计算精度高、可消除体积闭锁现象等特点,是目前科学和工程计算的重要数值方法之一.移动最小二乘法是无网格方法中构造形函数的最重要方法之一,很多重要的无网格方法都是基于移动最小二乘法建立的.但是移动最小二乘法的形函数不满足Kronecker函数的性质,使得
本文使用非线性流固耦合偏微分方程组:Cauchy守恒方程和P-T/T方程,描述复杂接触表面的形变和应力分布变化,其中方程的初边值条件根据实验数据分析结果而设置。该方程组采用有限元和有限差分相结合的方法求解,即,使用有限元方法对问题的空间变量进行离散,得到常微分方程组的初值问题后,进一步使用有限差分方法对问题的时间变量进行离散,从而得到原问题的全离散格式。空间变量的离散格式取决于所使用的有限元形函数
在自然科学和工程技术的实际计算中,最基本的方法是用函数的Taylor展开的部分和作为该函数的近似,而Pade逼近则是一种特殊的有理逼近,它是Taylor多项式逼近的延伸.Pade逼近的研究和发展与数学中的解析函数论,逼近论,矩量理论,连分式以及差分方程等分支有着紧密的联系,并且它在数值分析,量子场论,临界现象和控制论等自然科学领域中已有若干功效卓著的应用.本文主要给出了基于广义逆多元Newton-