Advances in Applied Mathematics
Vol.
11
No.
01
(
2022
), Article ID:
48386
,
11
pages
10.12677/AAM.2022.111044
基于无相位近场的非线性积分方程方法重构粗糙曲面
李莉莉
长沙理工大学,湖南 长沙
收稿日期:2021年12月24日;录用日期:2022年1月14日;发布日期:2022年1月27日
摘要
本文研究第三类边界条件无穷粗糙曲面反散射问题。入射场 选取不同位置点源叠加的形式,点源打到散射曲面上产生散射场 ,因此测量得到无相位近场 ,于是我们利用无相位近场的非线性积分方程方法重构无穷粗糙曲面的位置和形状。该算法的有效性通过数值算例实现。
关键词
反散射问题,第三类边界条件无穷粗糙曲面,无相位近场,非线性积分方程方法
Reconstructions of a Rough Surface Based on a Nonlinear Integral Equation Method from Phaseless Near-Field Data
Lili Li
Changsha University of Science and Technology, Changsha Hunan
Received: Dec. 24th, 2021; accepted: Jan. 14th, 2022; published: Jan. 27th, 2022
ABSTRACT
This paper studies the inverse scattering problem of infinite rough surfaces with boundary condition of the third kind. The incoming field is taken as a sum of different position source points shooting at the scattering surface, which induces the scattered field . Then we reconstruct infinite rough surfaces by a nonlinear integral equation method from the phaseless near field . Numerical examples are presented to show the effectiveness of the inverse algorithm.
Keywords:Inverse Scattering Problem, Infinite Rough Surfaces with Boundary Condition of the Third Kind, Phaseless Near-Field, Nonlinear Integral Equation Method
Copyright © 2022 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
1. 绪论
在本篇文章中,我们研究无相位近场的第三类边界条件无穷粗糙曲面反散射问题。此类问题在诸如石油勘探、医学成像、信号处理、无损探测、工业控制等领域都有着广泛的应用前景和实际意义。
我们研究的第三类边界条件无穷粗糙曲面散射问题由两部分组成。第一部分是正散射问题,即已知入射场和散射曲面决定散射场。目前已有大量的研究工作比如积分方程方法 [1] - [7] 和变分法 [8] 说明了正散射问题的适定性,并且实现了正散射问题的数值计算,例如基于积分方程方法的Nystöm方法 [9] 和基于变分法的有限元方法 [10]。
第二部分则是反散射问题,即是通过测量远场或者近场来确定未知粗糙曲面的形状、位置以及物理属性的一类问题。反散射问题的数值算法主要有两大类,第一类为非迭代算法,例如线性采样方法 [11]、时域点源法 [12] 和时域探测法 [13]。第二类是迭代算法,例如非线性积分方程方法 [14] [15]、Kirsch-Kress方法 [16] [17]、牛顿类型的迭代算法 [7] [18] 和基于变换场展开的逆算法 [19] [20],此类方法的优点在于如果选取适当的初始猜测值,往往能得到非常精确的重构效果。本篇文章所用的方法即迭代算法中的非线性积分方程方法。
上述所说的反散射问题数值算法是基于我们能够获得近场或者远场的强度和相位信息。在本篇文章中,我们是利用无相位近场重构第三类边界条件粗糙曲面的形状和位置,这是因为在实际应用问题中我们很难甚至测不到数据的相位信息。于是利用无相位的测量数据来重构散射体吸引了众多学者的关注,并取得了大量的研究进展 [21] - [39]。但是利用无相位测量数据重构散射体面临的挑战是只用一个平面波作为入射场所产生的无相位远场满足平移不变性 [26] [30],即散射体任意平移前后所对应的无相位远场相同,这就导致利用无相位远场只能重构散射体的形状而无法重构其位置。在文献 [21] [22] [23] [24] [25] [27] [28] 中,我们得知利用无相位远场可以很好地重构散射体的形状。
目前,广泛应用于打破平移不变性的方法有两种,一种是在文献 [38] [39] 中,证明了由两个平面波叠加作为入射场产生的无相位远场不满足平移不变性,并利用此方法在无相位远场反散射问题取得了进一步突破,例如在文献 [34] 中证明了在单波数情形下由无穷多组两个平面波叠加作为入射场产生的无相位远场可以唯一确定散射体的形状和位置。另一种打破平移不变性的方法是在反散射问题模型中加入一个参考球组成一个整体由此唯一决定散射体 [29] [31] [32] [35] [36]。
虽然在文献 [33] [37] 中证明了以两个不同位置的点源入射所产生的无相位远场可以唯一确定散射体的位置,并且在文献 [36] 中表明点源叠加入射产生的无相位近场可以唯一决定局部扰动的粗糙曲面,但是我们没有获得多个点源叠加入射产生的无相位近场能否唯一决定全局扰动的粗糙曲面的理论证明。在本篇文章中,我们实现了在多波数情形下,利用多个点源叠加入射所产生的无相位近场重构无穷粗糙曲面的形状和位置的数值算法。通过文献 [6] 我们得知散射场由一个积分方程表示,于是由近场为入射场和散射场之和可知近场的表达形式,对其关于粗糙曲面和阻尼函数求Fréchet导数可获得关于粗糙曲面的近似Fréchet导数,由此我们提出了无相位近场的非线性积分方程方法。为了解决反散射问题的不适定性,我们利用多个波数,由此得到更精确的重构曲面,其迭代算法可以简述为给定在第一个波数下粗糙曲面和阻尼函数的初始猜测,通过无相位近场的非线性积分方程方法重构出第一个波数下的粗糙曲面,并作为第二个波数的初始猜测,由此进行下去直到迭代完所有波数以及满足停止准则。
这篇文章的第二部分构建了无相位近场的第三类边界条件无穷粗糙曲面散射问题的模型,并在第三部分提出无相位近场的非线性积分方程方法。第四部分则展示数值算例说明该算法的有效性。
2. 模型建立
我们将在这部分建立第三类边界条件无穷粗糙曲面反散射问题的数学模型,并给出散射场的积分表达。
Figure 1. The geometry of the scattering problem
图1. 散射问题的几何表示
设函数空间 ,散射曲面 ,散射曲面上方的空间记为 ,分别用 和 表示 和 依赖函数f, 表示散射曲面 的法向量,指向 的内部,函数空间 ,在 内定义法向导数为
,。
如图1所示,我们将第三类边界条件粗糙曲面散射问题表述为已知入射场 和散射曲面 ,求解在 内满足Helmholtz方程的散射场 且 ,即
,, (2.1)
这里 表示波数。散射场 在散射曲面 满足第三类边界条件,即
, (2.2)
其中 ,且 , 表示阻尼函数,是一个连续的复值函数。为了保证该散射问题的适定性,散射场 满足向上传播条件:存在 ,,使得
,, (2.3)
其中 是Helmholtz方程的基本解, 表示第一类Hankel函数。此外,散射场 需满足以下两个条件:
1) 存在 ,使得
; (2.4)
2) 存在 ,,使得
,, (2.5)
这里 ,。于是(2.1)~(2.5)表示散射场满足第三类边界条件粗糙曲面散射问题。在本篇文章中,入射场 选用点源叠加的形式,即
,。 (2.6)
在介绍接下来的结论之前,我们先设 ,,M为常数, ; 为大于0的常数; , ; ; 表示在上半平面 满足算子 的阻尼格林函数。
定理2.1设 ,, :散射问题(2.1)-(2.5)存在唯一解 ,其表达形式为
,
且
,
存在唯一解 。
为了方便书写,我们引进积分算子的形式,记
,
则 ,。
3. 无相位近场的非线性积分方程方法
在这部分我们将通过无相位近场 ( 且 )的非线性积分方程方法重构粗糙曲面的形状和位置。测量无相位近场的位置为 。
由定理2.1,我们定义密度函数 ,则
, (3.1)
利用Nystöm方法求解得到密度函数 ,详细步骤参考 [15]。由(3.1)得到已知的密度函数 ,则无相位近场
,。 (3.2)
综上,我们得出求解无相位近场的反散射问题等价于求解(3.1)和(3.2)。由于阻尼格林函数 计算极其复杂,所以我们选用Dirichlet格林函数 代替阻尼格林函数 ,其中 ,,其可行性参考 [14]。于是(3.1)和(3.2)化为
,, (3.3)
我们记
则
,, (3.4)
其中
, (3.5)
对(3.4)两边同时取平方,则有
,。 (3.6)
由于(3.6)关于f和 都是非线性的,求解f和 我们需要将其线性化,即对(3.6)关于f, 分别在 , 方向求Fréchet导数,于是我们有
,。 (3.7)
由文献 [31] [40] 关于Fréchet导数的计算,我们得出对积分算子求Fréchet导数只需对积分算子的积分核求Fréchet导数,于是有
和
,
其中
类似的,
。
在这篇文章中,表示散射曲面 的函数f,阻尼函数 以及f的更新数据 , 的更新数据 由三次样条函数的线性组合表出。设 ,,
, (3.8)
, (3.9)
其中
这里 , 为均匀分布在x轴上−d到d的点, 为三次样条函数,即
为了方便起见,我们定义列向量
,
和矩阵 ,。
我们用 表示在 上 个测量点,将 和(3.8),(3.9)式代入(3.7)式得到
, (3.10)
其中r是一个列向量,其元素为
。
设矩阵
,
,
于是有 。由于(3.5)式的积分核是光滑的,所以(3.5)式严重不适定,则(3.7)式和(3.10)式也具有同样的性质,所以我们需对(3.10)式进行Tikhonov正则化,得到正则化方程
, (3.11)
为J的伴随矩阵, 表示正则化参数。
接下来,我们说明无相位近场的非线性积分方程方法重构无穷粗糙曲面的迭代算法。为了使算法简洁明了,首先给出单波数情形无相位近场的非线性积分方程方法。
迭代算法1给定 上测量得到关于粗糙曲面的无相位近场 以及截断误差 :
1) 选定关于粗糙曲面的初始猜测 和阻尼函数的初始猜测 ;
2) 已知 ,,,在(3.3)式中运用Nystöm方法求解 ;
3) 通过以上两步,我们已经知道 ,,于是求得(3.11)的更新数据 和 ,则 , ;
4) 重复2) 3) 步直到 ,迭代停止。这里 和 分别表示近似粗糙曲面 和精确曲面f的无相位近场。
基于算法1,我们现在给出多波数情形无相位近场的非线性积分方程方法的迭代算法。
迭代算法2对于波数 ,给定在每个波数下测量得到关于粗糙曲面的无相位近场数据 ,,初始值 ,初始截断误差 以及下降率 :
1) 设 ,当 时,令 ,,,否则 ,,,且 , 和 已经通过算法1得出;
2) 令 ,若 重复第(1)步,否则迭代停止。
4. 数值算例
在这部分,我们将展示数值算例来说明第3部分算法1和算法2的有效性。在以下算例中,我们设置共同参数如下:
· 入射场 为241个点源之和,点源所在位置为 ,。求解(3.3)式密度函数 时,设置粗糙曲面的截断间隔为 ,步长 ,其中 ,波数 ,在每个波长 放入10个等距节点。测量数据的高度为 ,取值范围从 到 ,步长 。
· 在计算(3.8)和(3.9)式时,设置 ,,,。(3.11)式中,多次数值实验表明取正则化参数 效果最好。为了说明算法的稳定性,我们在测量数据中加入噪音 ,所以加入噪音的测量数据为
,
其中 ,, 服从正态分布。在算法2中,设置下降率 。
· 在图2,图3中,黑色虚线,红色虚线以及蓝色曲线分别表示粗糙曲面 的初始猜测,重构曲面以及精确值。第一、二、三排图片分别表示没有加噪音,噪音为2%和5%的重构算例。
算例1我们首先重构一个只有局部扰动的粗糙曲面,其函数表达为
。 (4.1)
阻尼函数表达如下
,
,
则 。
重构(4.1)式时,散射曲面和阻尼函数的初始猜测分别为 和 ,测量数据没有加噪音和加入2%的噪音时,初始截断误差 ,加入5%的噪音时, 。
算例2这个算例我们重构一个无穷全局扰动的粗糙曲面,其函数表达为
(4.2)
阻尼函数表达如下
,
则 。
算例2中散射曲面和阻尼函数的初始猜测分别为 ,。测量数据没有加噪音和加入2%的噪音时,初始截断误差 ,加入5%的噪音时, 。
Figure 2. Reconstruction of (4.1) from data with no noise, 2% noise, 5% noise
图2. 局部粗糙曲面(4.1)的重构,噪音为0%,2%,5%
Figure 3. Reconstruction of (4.2) from data with no noise, 2% noise, 5% noise
图3. 局部粗糙曲面(4.2)的重构,噪音0%,2%,5%
5. 结论
本文实现了多频点源叠加入射产生的无相位近场的非线性积分方程方法重构第三类边界条件无穷粗糙曲面的数值算法,并且表明该散射问题等价于求解两个联立的积分方程组。数值实例表明点源叠加入射产生的无相位近场不满足平移不变性,在接下来的工作中,我们将致力于研究点源叠加入射产生的无相位近场不满足平移不变性的理论证明。
致谢
本篇文章是在导师李建樑副教授的悉心指导下完成,在此真诚地感谢李老师的帮助。在论文的撰写与修改过程中,李老师严谨细致,诲人不倦,我将受益终身。
文章引用
李莉莉. 基于无相位近场的非线性积分方程方法重构粗糙曲面
Reconstructions of a Rough Surface Based on a Nonlinear Integral Equation Method from Phaseless Near-Field Data[J]. 应用数学进展, 2022, 11(01): 359-369. https://doi.org/10.12677/AAM.2022.111044
参考文献
- 1. Bao, G. and Lin, J. (2011) Imaging of Local Surface Displacement on an Infinite Ground Plane: The Multiple Frequency Case. SIAM Journal on Applied Mathematics, 71, 1733-1752. https://doi.org/10.1137/110824644
- 2. Chandler-Wilde, S.N., Ross, C.R. and Zhang, B. (1999) Scattering by Infinite One-Dimensional Rough Surfaces. Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences, 455, 3767-3787. https://doi.org/10.1098/rspa.1999.0476
- 3. Chandler-Wilde, S.N. and Zhang, B. (1998) Electromagnetic Scattering by an Inhomogeneous Conducting or Dielectric Layer on a Perfectly Conducting Plate. Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences, 454, 519-542. https://doi.org/10.1098/rspa.1998.0173
- 4. Chandler-Wilde, S.N. and Zhang, B. (1998) A Uniqueness Result for Scattering by Infinite Rough Surfaces. SIAM Journal on Applied Mathematics, 58, 1774-1790. https://doi.org/10.1137/S0036139996309722
- 5. Chandler-Wilde, S.N. and Zhang, B. (1999) Scattering of Electromagnetic Waves by Rough Interfaces and Inhomogeneous Layers. SIAM Journal on Applied Mathematics, 30, 559-583. https://doi.org/10.1137/S0036141097328932
- 6. Zhang, B. and Chandler-Wilde, S.N. (2003) Integral Equation Methods for Scattering by Infinite Rough Surfaces. Mathematical Methods in the Applied Sciences, 26, 463-488. https://doi.org/10.1002/mma.361
- 7. Zhang, H. and Zhang, B. (2013) A Novel Integral Equation for Scattering by Locally Rough Surfaces and Application to the Inverse Problem. SIAM Journal on Applied Mathematics, 73, 1811-1829. https://doi.org/10.1137/130908324
- 8. Chandler-Wilde, S.N. and Elschner, J. (2010) Variational Approach in Weighted Sobolev Spaces to Scattering by Unbounded Rough Surfaces. SIAM Journal on Mathematical Analysis, 42, 2554-2580. https://doi.org/10.1137/090776111
- 9. Meier, A., Arens, T., Chandler-Wilde, S.N. and Kirsch, A. (2000) A Nyström Method for a Class of Integral Equations on the Real Line with Applications to Scattering by diffraction gratings and Rough Surfaces. Journal of Integral Equations and Applications, 12, 281-321. https://doi.org/10.1216/jiea/1020282209
- 10. Chandler-Wilde, S.N. and Monk, P. (2009) The PML for Rough Surface Scattering. Applied Numerical Mathematics, 59, 2131-2154. https://doi.org/10.1016/j.apnum.2008.12.007
- 11. Ding, M., Li, J., Liu, K. and Yang, J. (2017) Imaging of Local Rough Surfaces by the Linear Sampling Method with the Near-Field Data. SIAM Journal on Imaging Sciences, 10, 1579-1602. https://doi.org/10.1137/16M1097997
- 12. Lines, C.D. and Chandler-Wilde, S.N. (2005) A Time-Domain Point Source Method for Inverse Scattering by Rough Surfaces. Computing, 75, 157-180. https://doi.org/10.1007/s00607-004-0109-8
- 13. Burkard, C. and Potthast, R. (2009) A Time-Domain Probe Method for Three-Dimensional Rough Surface Reconstructions. Inverse Problems and Imaging, 3, 259-274. https://doi.org/10.3934/ipi.2009.3.259
- 14. Li, J (2019) Simultaneous Recovery of an Infinite Rough Surface and the Impedance from Near-Field Data. Inverse Problems in Science and Engineering, 27, 17-36. https://doi.org/10.1080/17415977.2018.1436172
- 15. Li, J. and Sun, G. (2015) A Nonlinear Integral Equation Method for the Inverse Scattering Problem by Sound-Soft Rough Surfaces. Inverse Problems in Science and Engineering, 23, 557-577. https://doi.org/10.1080/17415977.2014.922077
- 16. Burkard, C. and Potthast, R. (2010) A Multi-Section Approach for Rough Surface Reconstruction via the Kirsch-Kress Scheme. Inverse Problems, 26, Article ID: 045007. https://doi.org/10.1088/0266-5611/26/4/045007
- 17. Li, J., Sun, G. and Zhang, B. (2017) The Kirsch-Kress Method for Inverse Scattering by Infinite Locally Rough Interfaces. Applicable Analysis, 96, 85-107. https://doi.org/10.1080/00036811.2016.1192141
- 18. Bao, G. and Lin, J. (2013) Near-Field Imaging of the Surface Displacement on an Infinite Ground Plane. Inverse Problems and Imaging, 7, 377-396. https://doi.org/10.3934/ipi.2013.7.377
- 19. Bao, G. and Li, P. (2013) Near-Field Imaging of Infinite Rough Surfaces. SIAM Journal on Applied Mathematics, 73, 2162-2187. https://doi.org/10.1137/130916266
- 20. Bao, G. and Li, P. (2014) Near-Field Imaging of Infinite Rough Surfaces in Dielectric Media. SIAM Journal on Imaging Sciences, 7, 867-899. https://doi.org/10.1137/130944485
- 21. Chen, Z. and Huang, G. (2017) Phaseless Imaging by Reverse Time Migration: Acoustic Waves. Numerical Mathematics: Theory, Methods and Applications, 10, 1-21. https://doi.org/10.4208/nmtma.2017.m1617
- 22. Gao, P., Dong, H. and Ma, F. (2018) Inverse Scattering via Nonlinear Integral Equations Method for a Sound-Soft Crack from Phaseless Data. Applications of Mathematics, 63, 149-165. https://doi.org/10.21136/AM.2018.0154-17
- 23. Ivanyshyn, O. (2007) Shape Reconstruction of Acoustic Obstacles from the Modulus of the Far Field Pattern. Inverse Problems and Imaging, 1, 609-622. https://doi.org/10.3934/ipi.2007.1.609
- 24. Ivanyshyn, O. and Kress, R. (2010) Identification of a Sound-Soft 3D Obstacles from Phaseless Data. Inverse Problems and Imaging, 4, 131-149. https://doi.org/10.3934/ipi.2010.4.131
- 25. Ivanyshyn, O. and Kress, R. (2011) Inverse Scattering for Surface Impedance from Phaseless Far Field Data. Journal of Computational Physics, 230, 3443-3452. https://doi.org/10.1016/j.jcp.2011.01.038
- 26. Kress, R. and Rundell, W. (1997) Inverse Obstacle Scattering with Modulus of the Far Field Pattern as Data. In: Engl, H.W., Louis, A.K. and Rundell, W., Eds., Inverse Problems in Medical Imaging and Nondestructive Testing, Springer, Vienna, 75-92. https://doi.org/10.1007/978-3-7091-6521-8_7
- 27. Lee, K.-M. (2016) Shape Reconstructions from Phaseless Data. Engineering Analysis with Boundary Elements, 71, 174-178. https://doi.org/10.1016/j.enganabound.2016.08.001
- 28. Li, J., Liu, H. and Wang, Y. (2017) Recovering an Electromagnetic Obstacle by a Few Phaseless Backscattering Measurements. Inverse Problems, 33, Article ID: 035001. https://doi.org/10.1088/1361-6420/aa5bf3
- 29. Li, J., Liu, H. and Zou, J. (2009) Strengthened Linear Sampling Method with a Reference Ball. SIAM Journal on Scientific Computing, 31, 4013-4040. https://doi.org/10.1137/080734170
- 30. Liu, J. and Seo, J. (2004) On Stability for a Translated Obstacle with Impedance Boundary Condition. Nonlinear Analysis: Theory, Methods & Applications, 59, 731-744. https://doi.org/10.1016/S0362-546X(04)00284-6
- 31. Novikov, R.G. (2015) Formulas for Phase Recovering from Phaseless Scattering Data at Fixed Frequency. Bulletin des Sciences Mathématiques, 139, 923-936. https://doi.org/10.1016/j.bulsci.2015.04.005
- 32. Novikov, R.G. (2016) Explicit Formulas and Global Uniqueness for Phaseless Inverse Scattering in Multidimensions. The Journal of Geometric Analysis, 26, 346-359. https://doi.org/10.1007/s12220-014-9553-7
- 33. Sun, F., Zhang, D., and Guo, Y. (2019) Uniqueness in Phaseless Inverse Scattering Problems with Known Superposition of Incident Point Sources. Inverse Problems, 35, Article ID: 105007. https://doi.org/10.1088/1361-6420/ab3373
- 34. Xu, X., Zhang, B. and Zhang, H. (2018) Uniqueness in Inverse Scattering Problems with Phaseless Far-Field Data at a Fixed Frequency. SIAM Journal on Applied Mathematics, 78, 1737-1753. https://doi.org/10.1137/17M1149699
- 35. Zhang, D. and Guo, Y. (2018) Uniqueness Results on Phaseless Inverse Scattering with a Reference Ball. Inverse Problems, 34, Article ID: 085002. https://doi.org/10.1088/1361-6420/aac53c
- 36. Zhang, D., Guo, Y., Li, J. and Liu, H. (2018) Retrieval of Acoustic Sources from Multi-Frequency Phaseless Data. Inverse Problems, 34, Article ID: 094001. https://doi.org/10.1088/1361-6420/aaccda
- 37. Zhang, D., Guo, Y., Sun, F. and Liu, H. (2020) Unique Determination in Inverse Scattering Problems with Phaselessnear-Field Measurements. Inverse Problems and Imaging, 14, 569-582. https://doi.org/10.3934/ipi.2020026
- 38. Zhang, B. and Zhang, H. (2017) Recovering Scattering Obstacles by Multi-Frequency Phaseless Far-Field Data. Journal of Computational Physics, 345, 58-73. https://doi.org/10.1016/j.jcp.2017.05.022
- 39. Zhang, B. and Zhang, H. (2017) Imaging of Locally Rough Surfaces from Intensity-Only Far-Field or Near-Field Data. Inverse Problems, 33, Article ID: 055001. https://doi.org/10.1088/1361-6420/aa5fc8
- 40. Potthast. R. (1994) Fréchet Differentiability of Boundary Integral Operators in Inverse Acoustic Scattering. Inverse Problems, 10, 431-447. https://doi.org/10.1088/0266-5611/10/2/016
NOTES
Email: 13688492756@163.com