三维声散射反问题的一种快速算法.pdfVIP

尤云祥缪国平刘应中 (上海交通大学，上海200030) 摘要矗立研究从声散射塌的远场分布的信息珥删。三维物体形状的臣问题。提出，求南竿这类反问题 的一种‘把速算法。在这个算毡中，既小需预先知道散射体的边羿类型．也小需预先知道关于散射体的诸 如单连：厦、对葵内某一点星彤等先验几何信息．对有限孔径近场数据的数值算倒表明r这一方法是实用 的． 关键词反问题，物形重构．有限孔径 本文研究从卢散射场的远场分布的信息再现散射体边界物形的问题。众所周知，当 频率在港振区时．这类反问题通常既是不适定的义是完全1#线性的【i。1，因此其数值求解 1F常凼难。目前，求解这种完全1r线性不适定问题的数值方法主要有两类口】：一类是直 接方法．另一类是间接方法，即前线性最优化方法。对直接．方法，在其每一步迭代过程 中，都必须对迭代所得新物形求解若干相应的正问题；对间接方法，虽然其有不需求解 正问题的优点，但通常需要求解一个规模很大的非线性最优化问题。而且还可能遭遇局 部极小解。总体上，这些方法的计算耸都是很火的。此外，在这些方法中，均刚到了散 射体的请如单近通、星形以及声软、声刚性等儿何和物理的先验信息． 最近，Col／or＆Kirsch[”就二维物形的反演问题提山了一种快速求解方法。在这个方 法中．：F需要预先知道散射体的边界属性(即声软、声刚性等)、单：旺通、星形等物理和 Kirsch的方法中，闲来再现物体形状的信息，是远场分布在绕散射体周向所有观测方向 上的全孔径信息。然而，在实际应玎j中，要获得所有这些全孔径信息通常是不现实的， 充其颦只能测得远场分布在绕散射体周向某个有限孔径范围中若干有限个观测方向上的 信息。 本文发展了Colton＆Kirsch的方法，基于将远场分布在有限个观测方向上的远场数 据直接作加权台『J豆的技术．就二维物形的反演问题，建立了一种简单快速求解方法。在 这个方法中，同样不需要预先知道敞射体的边界属性、单连通、星形等物理和几何的先 验信息。与Colton＆Kirsch方法的不同之处是：住我们的方法中，心来再现散射体物形 的信息是远场分布在有限个观测方向上的远场数据，它既可以是全孔径远场数据也可以 是有限孔径远场数据．而且求解的只是一个非常简单的线性代数方程纽。 1一个简单算法及其理论分析 设D是在同质均匀介质R3中的一个具有光滑边界国的有界开区域，其外部 ’国家自然科学基金、上海交通大学科技基金资助项目 ·94· 3＼西为单连通区域。在介质申f 2=JR 一时谐平面入射波，其空问速度势为 u’fx棚)=expf÷奴·d)．式中k0IJ、3滚数d是一个l—定的挚艟向靖．它表示入射波 的入射方向。此入射卢波碰划K体后将产生散射，记“5为散射场，则介质中声波∞总场 为“=“’+“5。lE散射问题就是求一个总场“．满足 Au+k‘“=0， 在D，中 (1) 其中，A是三维Laplacian算子：在无穷远处的Sommerfeld辐射条件 lira (2) r(au。／甜一觑’J=0，r=㈧ 对所有方向量=∥H一致地成立：以及在D上的边界条件 酬J=0在彩上 (3) 界条件。 从辐射条件(2)可以得到“5有如下的渐近性质【t1 Ⅳ5(x)=0“／r)F(fc；d)+O(r。) (4) 上式中的F称为散射场“’的远场分布，i表示远场分布F的观{910方向。没国是R1 中的单位球面，Q={iI，量2，…，i～)c船魁远场分布F的观测方向的集台， A=(每，破，…，九)c船是入射方向的集合。反问题就是：对某个给定的止波数≈．设 对每一个给定的d∈A．相应的远场分布在每一个观测方向i∈Q上的值，即F(i：d)， 均是己知的，试确定散射体D的物形。在建立求婿上述反问题的一个简单快述算法之前， 我们首先