【摘 要】
:
微分方程cauchy问题是不适定的,当测量的Cauchy数据带有微小的扰动,很可能会引起反演结果的巨大偏差.因而对Cauchy问题的研究,特别是提出有效的、稳定的、可实现的、快速的数值计算方法是当今的重要课题之一本论文主要研究了满足周期条件的Laplace方程初值问题和Helmholtz方程Cauchy问题的不适定性,并给出了相应的数学模型,同时提出了有效的、稳定的数值计算方法.论文分为五部分.第
论文部分内容阅读
微分方程cauchy问题是不适定的,当测量的Cauchy数据带有微小的扰动,很可能会引起反演结果的巨大偏差.因而对Cauchy问题的研究,特别是提出有效的、稳定的、可实现的、快速的数值计算方法是当今的重要课题之一本论文主要研究了满足周期条件的Laplace方程初值问题和Helmholtz方程Cauchy问题的不适定性,并给出了相应的数学模型,同时提出了有效的、稳定的数值计算方法.论文分为五部分.第一章,首先简单介绍了Helmholtz方程的物理背景,之后对Cauchy问题及Helmholtz方程的研究状况做简短的综述.第二章,主要介绍了不适定问题及求解不适定问题的数值算法.我们针对性的介绍了Helmholtz方程Cauchy问题和初值问题,总结了相应的数值算法,以及数值算法在反问题方面的应用.第三章,考虑了周期结构Cauchy问题即椭圆方程初值问题的数值求解。我们把微分方程初值问题转化为等价的算子方程,之后对算子的谱的构成做了详细的分析,将特征空间分为两个子空间,其中一个子空间的问题严重的不适定,针对该子空间问题的求解提出了正则化策略.具体结论如下:求U(y)=(u(y),ω(y))T,u∈C1([0,1];Hp0×Hp0),满足定理1算子A的谱为σ(A)={(2nπ);n=0,±1,±2,…},其对应的特征函数为令H+=span{Vn(1),Vn(2)},H-=span{Vn(1),Vn(2)}-1n=∞,H0=span{V0},其中V0=(01),则D(A)∈H+(?)H-(?)H0,并且A(H+)∈H+,A(H)∈H-,4(H0)∈H0即H+,H-和H0都是算子A的不变子空间。记方程(1)可分解为d/dyU(y)=A+U+(y),(2) d/dyU(y)=A0+U0+(y),(3) d/dyU-(y)=A-U-(y),(4)关于方程(2)-(4),有如下结果:定理2算子-A+,A0和A分别生成空间H+,H0和H-上的C0算子半群E+(y),E0(y)和E-(y).并且,方程(2),(3),和(4)的解U+(y),U0(y)和U-(y)满足算子方程组我们提出了一种带有正则化策略的数值方法,相应的正则化问题的解的表达式:该方法尚未在其它文献中出现过,且通过数值试验证明了本文方法的可行性及易实现性。第四章,考虑利用一个入射波和在光栅上部一条直线Гb上的散射场反演单周期良导体光栅的形状。先利用Dirichlet-to-Neumann映射得到散射场在Гb上的法向导数值,然后应用求解椭圆方程初值问题的谱方法求解Helmholtz方程Cauchy问题,得到Гb以下的全场,最后逐点寻找全场的零点并连接,得到的曲线即为反演的光栅形状。主要结论如下:已知全场满足方程其中uα(x+l,y)=uα(x,y),Ωlb={(x,y);b>y>f(x),x∈[0,l]}(?)u/(?)n为边界Гb上的外法同导数,T为已知的Dirichlet-to-Neumann算子,u0已知,α=ksinθ,θ是入射角,k是波数.将该问题转化为等价的初值问题:求解(u,w)满足方程对于定义算子其定义域为其中L是微分算子那么问题(8)可表述成如下抽象初值问题:求U(t)=(u(t),ω(t))T,u∈C1[0,1]; Hp0×Hp0满足定理3算子A的谱为σ(A)={λn;n=0,±1,±2,…},其中:其对应的特征函数为这里令σi(A)={λn;m∈Zi},σ+(A)={λn;λn>0,且n∈zr},σA={λn;λn<0,且n∈Zr},H+=span{Vn,λn∈σ+(A)}n∈Zr·H=span{Vn,λn∈σA∪σ(A)}n∈Zi∪Zr,则D(4)∈H+(?)H-,并且A(H+)∈H+,A(H-)∈H-,即H+和H都是算子A的不变子空间。记方程(11)可分解为d/dtU+(t)=A+U+(t)(12) d/dtU+(t)=A+U+(t)(13)定理4算子-A+和A分别生成空间H+和H上的C0算子半群E+(t)和E-(t).并且,方程(12)和(13)的解U+(t)和U-(t)满足算子方程组定理5当t>0时,算子E+(t):H+→H+和Et:H-→H皆为紧算子.故问题(7)的正则化解为通过数值试验证明了该方法是可行的。第五章,对论文的主要工作做了总结,并提出了将要考虑的问题.
其他文献
研究多项式自同构对于仿射代数几何有着非常重要的意义.而与多项式自同构密切相关的Jacobi猜测,自被提出以来就一直受到学者们的密切关注.在对这个问题进行研究的过程中,学者们不断地提出新思想,新问题,新方法,不仅促进了相关研究领域的融合,而且对多项式自同构的研究起着强有力的推动作用.令K是一个域,X=(X1….,Xn),K[X]表示关于变量X1….,Xn的多项式环.设Fi∈K[X],则F=(F1,.
本文分别针对两类带有开弧段的微分方程外问题-斜微商问题和弧的散射问题,提出两种不同的数值方法.Ⅰ.开弧外斜微商问题的基于Chebyshev多项式的方法首先我们考虑二维平面内模拟霍尔效应的开弧外斜微商问题.用一个简单开弧Γ∈C2,λ来模拟二维平面上的电极,其中λ∈(0,1].令用z1:=y(1),z-1:=y(-1)来表示r的两个端点.我们用Γ+来表示当s增长时,Γ的左手一侧,那么用Γ来表示Γ的另一
未来的量子信息网络需要人们运用精确的手段控制光信息的动态传输,在量子化媒介中对光信息进行有效存储以及利用新颖的技术对光学数据进行相干操控。光与物质相互作用导致的量子相干效应为实现这一目标提供了强大的工具。特别地,基于激光诱导原子相干的电磁感应透明(EIT)技术被非常成功的大量应用在量子信息领域。本论文回顾了原子相干效应,重点介绍了基于EIT的光存储技术和光学前驱场的相关研究工作的发展历程和前沿方向
相干控制是利用激光的相干特性,在原子分子层面对分子激发态动力学过程的有效调控。脉冲整形技术是在实验上实现分子反应相干控制的主要手段。本文主要使用脉冲整形技术研究了一些气体分子的准直、电离、解离等过程的相干控制,观察到了碘甲烷分子准直程度随着脉冲序列子脉冲间距的变化而改变,由脉冲序列对环戊酮电离解离过程的调制解释了其激发态动力学机制,实现了对一些分子电离、解离、高次电离过程的操控和对同分异构体分子的
本文主要运用线性微分(差分)Galois理论研究非线性系统的可积性与不可积性.全文共分为五章,第一、二章分别是绪论和预备知识,第三章到第五章是本文的主体部分.在第三章中,我们运用微分和差分Galois理论分别研究了一般微分系统及映射的不可积性,从而对Goriely所提的一个公开问题给出了比较完整的回答.同时,还考虑了一般解析非线性系统在共振情形下的不可积性.第四章讨论非线性微分系统的Painlev
动力系统的研究热点之一是探讨系统的混沌性态,而符号空间上的子移位在探讨各种混沌关系过程中起到了重要作用.本文研究符号空间上几类子移位的动力性质,重点考察这些子移位的混沌性态,所得结果或改进或发展了动力系统研究中的若干已有工作.本文主要结果包括:(1)构造了拓扑熵为零且分布混沌的子移位,证明了它是Wiggins混沌,Martelli混沌,弱混合且是严格遍历的.应用此结果,我们证明:一个具有正规转移不
设f(X)∈k[X],我们称f(X)的所有同构象中所依赖的变元xi的最小个数为f的外秩(outer rank),简记为Orank(f)一直以来,具体指明给定元素的外秩是许多代数系统(如自由李代数,自由群)重要的研究课题.然而,据我们所知,关于多项式环上的外秩并没有太多的研究.在本文中我们以多项式收缩作为工具具体给出了一类多项式的外秩.另外,如果自同态将每一个多项式变成一个与该多项式具有相同外秩的多
套代数是一类非常重要的非自伴非交换算子代数,它与不变子空间问题密切相关,是Robust控制一个非常重要的应用背景.本文主要讨论套代数的相关问题.第一章介绍了套代数和斜Toeplitz算子相关背景,以及套代数可逆元群连通性问题现有结果;第二章基于前人的工作,讨论了下三角代数一类可逆元的连通性问题;第三章对k∈Z+我们给出了以“长度”小于等于k的多项式为符号的k阶斜Toeplitz算子的谱,考虑了它们
扩散方程是一类非常重要的偏微分方程,自然界中来源于物理、化学、经济和生物等领域的大量现象都可以用扩散方程数学模型来刻画.近年来,国内外愈来愈多的数学家、化学家、物理学家和生物学家关注于扩散方程领域的研究,对数学提出了许多挑战性的问题.对于扩散方程(组)解的整体存在、爆破、熄灭、爆破的临界指标、爆破时间、熄灭时间以及熄灭速率等问题的研究已成为偏微分方程理论研究中的一个重要方向.本文主要研究具有正的初
在研究油-水-表面活性剂三相系统动力学性质时得到了如下模型方程(1)是一个典型的高阶抛物方程,其中γ>0,这里u与油和水的局部浓度差成正比,f(u)=F’(u),而F(u)和a(u)由如下六次和二次多项式给出在过去的几年里,六阶抛物方程得到了广泛的关注,很多学者研究了六阶抛物方程,得到了一些结果,例如解的存在性,惟一性和正则性[4]-[8].然而,据我们所知,对于方程(1)的研究结果却很少,Paw