同时反演阻尼区域和系数的一种方法
2010-05-29王俊杰王连堂
王俊杰，王连堂
(1.思茅师范高等专科学校数学系，云南 普洱 665000;2.西北大学数学系，陕西 西安 710127)
0 引 言
声波的反散射理论是一个典型的数学物理反问题，它是利用外部的测量信息，探测介质的内部性质或者边界性质，而这些信息通常是不能直接测量的，由于声波反散射理论在雷达及地球物理勘探等领域的需要，对反散射理论及计算方法的研究有着广泛的应用前景.D.Colton,P.Monk和R.L.Ochs[1-3]分别用远场模式的完全与不完全数据进行了区域反演，David.Colton[4]中应用线性抽样方法进行了区域反演，在Pedro.Serranho[5]中应用杂交方法对阻尼边界进行了区域和阻尼系数的同时反演.本文研究声波阻尼反问题的区域和系数同时反演的反问题，即在阻尼系数和区域未知的情况下，利用远场模式的数据来同时反演阻尼区域和系数.本文给出了具体的计算方法，并给出了数值例子.
本文考虑在均匀介质中传播的声波，此声波碰到一个无限长的柱体.设柱体的截面D⊂R2母线平行于z轴，设入射波是平面波ui(x)=eikx·α,x∈R2其中k>0是波数，α为一单位向量，入射波碰到柱体发生散射，记总体场为u=ui+us,us表示散射场，总体场满足阻尼边界条件.正散射问题是u∈C2(R2)
Δu+k2u=0 inR2D-
(1)
(2)
(3)
其中v表示外法线方向，λ为声波阻尼系数，(3) 称为sommerfeld散射条件.
若IM(λ),∂D∈C2,则问题(1)～(3)存在唯一解[6].由文献[6]可知，散射波us有以下渐进性质：
(4)
其中(r,θ)为x的极坐标，F称为散射波us的远场模式.
本文考虑的反问题是给定F(θ;k,αn);其中θ∈[0,2π],αn,n=1,2,3,…,N是N个不同的单位向量，利用这些数据来确定阻尼边界∂D和系数λ.
1 反演方法
由Green公式与第一类零阶Hankel函数的渐进性质可得:
(5)
其中x=reiθ,y=reiφ,v表示∂D上的单位外法线方向.
本文的方法基于Herglotz波函数理论[7].一个函数v称为Herglotz波函数,如果:
Δv+k2v=0 inR2
在极坐标(ρ,φ)下,Herglotz波函数有如下积分形式:
(6)
其中g∈L2[0,2π]是Herglotz波函数的核, 由(2)和(5)得:
(7)
定义集合:
S={F(θ;k,αn)-F(θ;k,α1),n=1,2…}
其中F(θ;k,αn)是对应于入射波eikx·αn的远场模式.将S中的元素看成θ的函数,S是L2[0,2π]的子集,且
定理2.1[8]:设λ>0,v是如下问题的解：
(8)
a. 如果v是Herglotz波函数,核为g, 则S⊥=Span{g}
b. 如果v不是Herglotz波函数,核为g， 则S⊥=0
因此有以下定理成立:
(F,g)=(f1,g)
周围,其中g是逼近(8)的解的Herglotz波函数的核,支集含于[0,2π],(·,·)表示L2空间上的内积.
于是反问题可归结为一个约束最优化问题.
定义集合:
U1(M)={g∈H1[0,2π]∶|g|H1≤M}
U2={ρ∈C1[0,2π]∶ρ是周期函数.
U3={λ(φ)≤M1,|λ(φ)-λ(ψ)|≤M2|φ-ψ|}
其中M,M1,M2,a,b,c为常数.
由嵌入定理U1(M)在C[0,2π]中紧，U2(M)在C[0,2π]中紧，并且由Arzela-Ascoli定理，U2(M)在C[0,2π]中紧，再由Tikhonov定理，U(M)=U1(M)×U2×U3在C[0,2π]×C[0,2π]×C[0,2π]中紧.
以下给出(F1,g)的值，因为
且文献[2]中
再由
得到：
最终得到下面的最优化问题：
(9)
其中v为(5)所定义的Herglotz波函数，由于λ>0,所以内问题不存在任何特征值，因而在此只要求k>0即可.
2 数值计算
下面对前面的反演方法给出计算方法. 反问题的数据是一系列方向的平面波和所对应的远场模式，这些数据的产生是由解正问题所得.
本文采用Nyström方法解正问题.
对于定解问题
Δu+k2u=0 inR2D-
(10)
(11)
寻找如下单层位势的解：
从而得到散射波的远场模式：
对给定的波数k和[0,2π]上不同方向的N个入射平面波，用上面的方法得到的远场模式的逼近，将此数据作为原始数据来求反问题.
记Fα为Fα(θ;,k,α)，用有限三角级数：
逼近Fα(θ;k,α)，于是
其中：
用三角级数逼近g,ρ,λ
aj,bj∈R
于是得到：
利用矩形公式积分，则(9) 离散为：
其中：
反问题参数：n1=6,n2=6,n3=6,m1=m2=30,k=1.0.
入射波数目N=16,远场模式数目n=36.图中实线表示边界ρ或系数λ，虚线表示反演的边界ρα或系数λα.
考虑在均匀介质中传播的声波，此声波碰到一个无限长的柱体.设柱体的截面D⊂R2母线平行于z轴，设入射波是平面波ui(x)=eikx·α,x∈R2其中k>0是波数，α为一单位向量，入射波碰到柱体发生散射.在一些情况下我们需要知道柱体的截面形状和边界的性质，但是不能通过直接测量得到，可以通过测量到的散射波的信息得到即：通过远场模式反演区域形状和阻尼系数.
例3.1 精确的柱体截面形状:ρ=2+0.3×cos(3t),精确的阻尼系数λ=2+sin(t).
反演区域
aj,bj∈R
其中aj,bj见表1，
反演阻尼系数
表1 例3.1 的边界和阻尼系数同时反演的数值结果
图1 例3.1 的边界和阻尼系数同时反演的结果
例3.2 精确的柱体截面形状: 边界ρ=2+cos(t)sin(t),精确的阻尼系数λ=2+cos(t).
反演区域
aj,bj∈R,
其中aj,bj见表2.
反演阻尼系数
表2 例3.2 的边界和阻尼系数同时反演的数值结果
图2 例3.2 的边界和阻尼系数同时反演的结果
从图1和图2看出本文给出的同时反演阻尼系数和边界的反演方法的反演效果是比较好的.
参考文献：
[1]Colton D, Monk P. A novel method for solving the inverse scattering problem for time harmonic acoustic wave in the resonance region[J].SIAMJ Appl Math,1985,45:1039-1053.
[2]Colton D, Monk P. A novel method for solving the inverse scattering problem for time harmonic acoustic wave in the resonance region[J].II SIAMJ Appl Math,1986,46:506-523.
[3]Ochs R L. The Limited Aperture Problem of inverse scattering: Dirichlet Boundary Condition[J].SIAMJ Appl Math,1987,47:1320-1341.
[4]David Colton, Michele Piana, Roland Potthast. A Simple method using Morozov’sdis crepancy principle for solving inverse scattering problems[J].Inverse Problem,2001,17:1997-2015.
[5]Pedro Serranho.A Hybrid Method for Inverse Scattering for Shape and Impedance[J].Inverse Problem,2006,22:663-680.
[6]Colton D, Kress R.Integral equation methods in scattering theory[M].New York:Wileyinterscience publication,1983.
[7]Colton D, Kress R. Inverse and electromagnetic scattering theory[M].New York:Springer Verlag Berlin Heidelberg,1992.
[8]王连堂,利用远场模式反演声波阻尼系数[J].计算数学,1999,2:89-98.
