使用无相位数据反散射问题的数值方法研究

来源 :吉林大学 | 被引量 : 0次 | 上传用户:leobear
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文讨论了使用无相位数据重构裂缝和障碍物反散射问题的分析和计算,我们考虑的散射问题模型均为Helmholtz方程,由于在实际的应用中,获取散射场或者远场的相位信息是极其困难,因此,利用无相位数据求解反散射问题在数学和物理领域都得到了广泛的关注.本文的第一章为绪论,主要介绍了我们所做工作的背景以及研究现状,并且简单介绍本文的结构.第二章主要介绍本文用到的一些预备知识,包括声波散射的一些相关概念以及反问题正则化的方法和数值计算的Nystrom方法.第三章研究了利用无相位数据重构声软裂缝的问题,我们验证了平移不变性(Translation invariance),这表明利用无相位远场数据仅能重构散射体形状而不能重构散射体位置,对于反散射问题,我们采用非线性积分方程方法,此方法的优点是精度准确并且能减少计算复杂度,最后通过一些数值例子验证了算法的可行性.第四章研究了利用无相位数据重构加入参考球的裂缝问题,为了解决无相位数据无法重构裂缝位置的问题,我们考虑在散射系统中加入了声软参考球,利用非线性积分方程方法求解反问题,重构出声软裂缝的形状和位置,并且通过数值实验说明了方法的可行性.第五章研究了利用无相位数据重构加入参考球的障碍问题,我们考虑的是满足Neumann边界条件的障碍体,由平移不变性可知,仅能确定障碍的形状不能确定障碍的位置,我们利用加入声硬参考球的方法,与裂缝情况类似,使用非线性积分方程方法求解,重构出障碍体的形状和位置.本文的最后一章为结论,对全文内容进行了总结.本文几个主要工作如下:1.利用无相位数据重构裂缝问题设裂缝rc=c=:s ∈[-1,1]}为光滑弧线,z1:=z(1)和z-1:=z(-1)为裂缝Γc的端点.我们考虑如下模型问题:△u + k2uu = 0,in R2\Γc(1)u = 0,on rc,(2)反散射问题:已知一个方向入射的平面波ui以及远场模态的模|u∞|,重构声软裂缝rc的形状.使用无相位远场数据重构裂缝的反散射问题的解不是唯一的,我们给出了如下结论.定理1.(Translation invariance)设u∞(x)为声软裂缝re散射的远场模态.对于裂缝Γcε={x+εh:x∈Γc}其中h∈R2,远场模态u∞ε有如下形式u∞ε(x)=eikεh·(d-x)u∞(x),x∈Ω,(4)也就是说,利用无相位远场模态重构声软裂缝的反散射问题具有平移不变性.由平移不变性可知,仅利用一个方向入射的平面波得到的无相位远场模态数据不能重构出声软裂缝的位置.对于反散射问题,我们通过非线性积分方程方法求解,首先引进单层位势算子Se:(Sc(φ)(x)=Φ(x,y)φφ(y)ds(y),x∈Γc以及远场算子Sc,∞:(Sc,∞φ)(X)=-γ∫Γc e-ikx·yφ(y)ds(y),x∈ Ω.则可得未知曲线Γc和密度函数φ满足下述方程组Scφ = ui|Γc,(5)|Sc,∞=|2=|U∞|2.(6)然后将积分算子Sc转换为参数化算子C,表示为C(z,ψ)(t)=i/4∫π0H01(k|z(cost)-z(cosτ)|)ψ(τ)dτ,t ∈[0,π],(7)其中ψ(t):= |sint||z’(cost)|φ(z(cost)),t∈[0,π].与上述相似,我们引进参数化远场算子C∞:其中z∞(t)=(cost,sint),t∈[0,2π].此外,我们也给出入射场ui以及远场u∞的参数化表示:ωc=uioz,ωc,∞=u∞oz∞.则方程组的参数化形式可以表示为C(z,ψ)= ωc(z),(9)C∞(z,ψ)C∞(z,ψ)= |ωc,∞|2.(10)C∞(z,)关于z的Frechet导数表示为C∞ C∞关于z的导数为将方程(10)线性化得到B[z,ψ]q=fz,ψ,(12)其中迭代过程的相对误差为反问题的迭代过程为:(ⅰ)发射波数为k>0,入射方向为d ∈ Ω的平面波,然后收集无相位远场数据|u∞|.(ⅱ)给裂缝Γc一个初始近似Γc0,设k=0.(ⅲ)对于曲线Γck,从方程(9)中得到密度函数ψ.(ⅳ)求解方程(12)得到q,可得裂缝曲线新的近似Γck+1=Γck+q,计算Ek.(ⅴ)若Ek≥∈则设k=+ 1,返回(ⅲ),否则我们将近似曲线Γck+1作为最后的重构曲线.值得注意的是,我们使用的迭代方法与经典Newton迭代法是有些不同的,我们的迭代法精确,可以简单的实现并且减少计算复杂度.我们给出若干数值例子证明了方法的可行性.2.无相位数据重构加入参考球的裂缝问题由于利用无相位远场模态重构裂缝的反散射问题解是不唯一的,仅能够重构出裂缝的形状而并不能确定裂缝的位置.为了解决这个问题,我们采用在散射系统中加入声软参考球的方法.我们假定裂缝r1={z(s)= ∈[-1,1]}为光滑弧线,参考球D(?)R2,并且D ∩ =(?),Γ2:=(?)D.反散射问题:已知由固定波数k和某一个入射方向d入射的平面波ui以及无相位远场数据|uD∪Γ1∞(x)|,x∈Ω,确定未知裂缝Γ1的形状以及位置.我们定义单层位势算子以及远场算子我们得到裂缝Γ1和密度函数φj满足下面的方程组S11φ1+S21φ2=ui|Γ1,(14)S12φ1+S22φ2 =ui|Γ2)(15)|s1∞φ1+S2∞φ2|2=|u∞|2(16)我们给出裂缝r1的参数化表示Γ1 = {P1(s)= c+z(s):c=(c1,c2),s∈[-1,1 },以及参考球边界Γ2的参数化表示Γ2 = yp-2(t)= b + Rx:b=(b1,b2),x∈ Ω}.将算子Sjl,Sj∞转换为参数化算子Cjl,Cj∞,右端项的参数化形式为ω1(t)=ui(p1(coSt)),ω2(t)=ui,(p2(t)),ω∞(t)=u∞(x(t)).则可以得到积分方程(14)-(16)的参数化形式将方程(19)线性化可以得到Bq=f,(20)其中我们给出了一些数值例子用以证明方法的有效性.3.无相位数据重构加入参考球的障碍体问题我们假定障碍散射体D(?)R2,考虑如下的障碍散射模型问题:△u + k2u = 0,in R2\D,(21)(?)u/(?)v=0,on(?)D,(22)我们考虑的反散射问题是已知由固定波数k和某一个入射方向d入射的平面波ui以及无相位远场数据|uD∞(x)|,x∈Ω,确定位置障碍物D的形状以及位置.由于平移不变性|uDh∞(x,d)|=|uD∞(x,d)|,利用无相位远场数据仅能重构障碍物形状而不能够确定位置.为了解决这个问题,我们加入参考球B(?)R2,D ∩ B=(?).反散射问题:已知由固定波数kk和某一个入射方向d入射的平面波ui以及无相位远场数据|uD∪B(x)|,x∈Ω确定位置障碍物D的形状以及位置.定义法向导数算子和远场算子:可以得到障碍物D和密度函数φj满足下面的方程组我们给出D的边界Γ1的参数化表示Γ1={p1(x)=c+ r(x)x:c=(c1,C2),x ∈以及参考球边界Γ2的参数化表示Γ2 = {P2(x)= b + Rx:b=(b1,b2),x ∈ Ω}.将算子Tjl,Tj∞转换为参数化算子Ajl,Aj∞,右端项的参数化形式为我们可以得到方程组的参数化形式A11(p1,ψ1)+A21(p1,ψ2)=ω1,(27)A12(P2,ψ1)+A.22(P2,ψ2)=ω2,(28)|A1∞(p1,ψ1)+ A2∞(p2,ψ2)|2 =|ω∞|2.(29)将方程(29)线性化可以得到Bq=f,(30)其中最后我们给出了一些数值算例验证了方法的有效性。
其他文献
近年来,有机无机杂化金属卤钙钛矿材料在太阳能电池、照明以及显示方面所取得了重大的进展,推动了三维(3D)铯铅卤族钙钛矿CsPbX3(X=Cl,Br,I)材料的研究。由于其较好的稳定性,较高的荧光量子产率以及基于卤素成分诱导荧光宽色域可调性,促进了CsPbX3在光电材料领域的发展,被认为未来最具商业化前景的光电材料之一。通过降低铅卤素钙钛矿维度,能够有效的提高载流子的限域效应,从而提高自由激子结合能
宇称-时间(Parity-Time:PT)对称理论最初是在量子力学中提出的,该理论一经发表就引起了人们的广泛关注,随后被引入到光学领域中。通过精确调节介质的增益和损耗,人们可以在多种光学系统中实现PT对称性,并且发现了许多“违反”常识的现象,例如:吸收诱导的透明、相干完美吸收-激光以及超布洛赫振荡等。受此启发,人们意识到控制介质的增益和耗散可以实现传统光学中难以实现的现象和功能,对非厄米光学系统的
自然界存在四种基本相互作用即电磁相互作用、引力相互作用、强相互作用和弱相互作用。不同的相互作用有不同的对称性,我们从对称性出发能对相互作用理论有一些研究。目前人们已经相信以夸克和胶子为基本自由度的量子色动力学(Quantum Chromo Dynamics简称QCD)是描述强相互作用的基本理论。QCD中存在整体的(近似的)手征对称性,局域规范对称性,(近似的)同位旋对称性,标度对称性和分立的C、P
目的:食品安全是受到全世界关注的重要问题。现今每年有上亿人因食用受污染的食品患上疾病,超过数十万的人因此死亡。由于广泛的滥用,针对食源性致病菌传统上有效的抗生素处理也逐渐失去效果,并导致越来越多耐药性菌株的出现。金黄色葡萄球菌(Staphylococcus aureus,SA)是最常见的食源性致病菌之一,且多个菌株可针对几乎所有临床可用的抗生素产生耐药性。因此,开发能有效灭活金黄色葡萄球菌的新除菌
二维材料的发现与兴起为材料科学与纳米技术带入了具有各种可能性的新时代。与大尺寸的三维材料相比,二维材料原子级厚度和片状结构所引入的电子局限性、超强面内共价键等特性赋予了材料特殊的电子结构、更好的力学性能、更大的比表面积、更高的原子活性等优势。而对于材料而言,可以根据其原子顺序和排列分成晶体,微晶和非晶三种结构形式。其中非晶材料仅具备局部短程有序而不具有长范围的原子有序排列。与晶态的二维材料相比,非
元宇宙概念一经提出就引起了各方热议,图书馆是否应该参与元宇宙建设,能够发挥哪些功用,便成了亟需探讨的问题。图书馆与元宇宙的演进受到了相同的技术体系驱动,一方面图书馆的资源与服务是元宇宙的重要构成要素,另一方面元宇宙为图书馆带来了新的发展机遇。在元宇宙中,图书馆可以作为元宇宙与现实世界之间的交流管道,充当不同元宇宙的“交通港”,通过评价服务推动元宇宙的健康发展,并提供元宇宙原生的资源与知识服务。图书
在真核细胞中,遗传信息储存在DNA中并且以核小体为基本单位装载在染色质中。H2A.Z是存在于酵母和人细胞中的H2A一个高度保守的变体,它在染色质中的存在被认为与基因转录调控、DNA损伤修复、基因组稳定性相关。将组蛋白变体置换进、出染色质是染色质重塑中非常重要的组成部分。已经有很多发表的论文中报导了一些具有H2A.Z置换功能的酶。酵母中Swr1和它在人细胞中的两个同源物SRCAP和P400/TIP6
众所周知,Si基材料已形成丰富的产业链,尤其在电子工业中十分重要。在集成电路、芯片、光学器件、金属门限晶体管(metal gate transistor)和金属半导体接触器件(欧姆/肖特基器件)等产品中,基于Si的各种纳米结构和金属化合物使用广泛。相关的Si基纳米结构设计和材料性质研究也一直是理论研究领域的热点问题。本文通过结构设计和稳态结构搜索的方式结合密度泛函计算方法,分别对Si的表面修饰纳米
2021年是中国共产党建党一百周年,也是国家“十四五”规划开局之年。站在新的历史起点上,公共图书馆事业作为国家文化和旅游工作的重要组成部分,必须准确把握新发展阶段、深入贯彻新发展理念、加快构建新发展格局,实现公共图书馆事业的高质量发展。广东省委省政府提出要打造数字文化引领地、文化创意新高地、文旅融合示范地,加快打造粤港澳大湾区世界级旅游目的地,推动文化大发展大繁荣,努力塑造与经济实力相匹配的
期刊
本文以杉木作为生物质原料,采用热解法进行预处理得到杉木热解炭,并对杉木热解炭的理化特性、杉木热解炭与煤混燃特性及混燃过程中污染气体排放特性进行了深入研究,具体如下:在不同热解温度(350℃、450℃、550℃)下制备杉木热解炭,研究所得杉木热解炭的理化特性。结果表明:热解温度的增加使得杉木热解炭的质量产率及能量产率逐渐减小,且减小趋势随之增大,样品的孔隙也变得更加丰富,样品固定碳和灰分的含量逐渐增