带Neumann边界条件的Helmholtz方程柯西问题的一种新的正则化方法.pdf
《带Neumann边界条件的Helmholtz方程柯西问题的一种新的正则化方法.pdf》由会员分享,可在线阅读,更多相关《带Neumann边界条件的Helmholtz方程柯西问题的一种新的正则化方法.pdf(6页珍藏版)》请在咨信网上搜索。
1、第49 卷第4期2023年8 月文章编号:16 7 3-519 6(2 0 2 3)0 4-0 151-0 6带 Neumann 边界条件的 Helmholtz 方程柯西问题的一种新的正则化方法余亚辉*1,李振平1.2(1.洛阳理工学院数学与物理教学部,河南洛阳47 10 2 3;2.西北师范大学数学与统计学院,甘肃兰州7 30 0 7 0)摘要:考虑矩形域上带Neumann边界条件的Helmholtz方程的柯西问题,该问题是一类严重不适定的偏微分方程反问题,即它的解不连续依赖于输入数据.基于经典的Tikhonov正则化方法利用自设计过滤化子修改核函数的思想,提出一种新的正则化求解方法,给出该
2、问题基于分离变量的近似解,对正则化参数的先验和后验两种选取规则下精确解与近似解进行误差分析,得到满足收敛性和稳定性的Holder型误差估计.关键词:Helmholtz方程柯西问题;Neumann边界条件;不适定问题;正则化;后验参数选取中图分类号:0 17 5A new regularization method for the Cauchy problem of the Helmholtz(1.Department of Mathematics and Physics,Luoyang Institute of Science and Technology,Luoyang 471023,Chi
3、na;2.College of Mathe-matics and Statistics,Northwest Normal University,Lanzhou 730070,China)Abstract:A Cauchy problem for the Helmholtz equation with Neumann boundary conditions in a rectan-gle domain is discussed in this paper.This problem is a serious ill-posed inverse problem of partial differ-e
4、ntial equations,that is,the solution does not depend continuously on the input data.Based on the classi-cal Tikhonov regularization method using the idea of modifying the kernel function by self-designed filterssubsets,a new regularization method is presented for approximating the solution to this p
5、roblem based onthe method of separating variables.After performing error analysis on the exact and approximate solutionsunder both a priori and a posteriori selection rules of the regularization parameters,the Holder-type errorestimates satisfying convergence and stability are obtained.Key words:Cau
6、chy problem for Helmholtz equation;Neumann boundary condition;ill-posed problem;regularization;a posteriori parameter choice rule本文考虑矩形域上带Neumann边界条件的Helmholtz方程的柯西问题1u(,y)+ku(,y)=00a元,0 y0表示波数,=/a/a是La-place算子.这里已知y=0的数据u(,0)=f()),求满足式(1)的解u(,y)(0 y 1).实际问题中,数据通常是测量得到的带有误差,这里假设精确数据f()和带噪数据f()都属于L(O
7、,元),且满足:f(.)-f(.)0 表示噪声水平.即这里已知f(),求满足问题(1,2)的解。Helmholtz方程被广泛应用于散射,波的传播、结构的振动、电磁场等众多实际领域2-3.目前关于Helmholtz方程的正问题已经得到了广泛的研究.(2)152然而,Helmholtz方程的柯西问题是一类严重不适定的反问题4,即解不连续依赖于给定的柯西数据,给定数据中的任何小的扰动可能导致解的巨大的误差因此通常的数值方法难以求解,为了恢复解的稳定性,需要一些正则化方法.目前人们已经提出了许多正则化方法来求解Helmholtz方程柯西问题如采用Fourier截断方法1.5、拟逆方法6、Tikhono
8、v型正则化方法1.7、软化方法s-9、分数次Tikhonov方法10、映射正则化方法11等.虽然已有大量文献采用各种正则化方法求解Helmholtz方程柯西问题,但很多文献是针对无限区域进行研究的,有界区域上的研究大多是考虑的Dirichlet边界条件,并且已有结果主要集中于正则化参数的先验选取,而正则化参数的先验选取必须事先知道精确解的先验界信息,这在实际问题中往往是难以估计的.本文研究矩形域上带Neumann边界条件的Helmholtz方程的柯西问题,基于经典的Tikhonov正则化利用自设计过滤化子修改核函数的思想,提出一种新的正则化方法求解问题(1,2),给出形式更为简洁的近似解的表达
9、式,对正则化参数的先验和后验两种选取规则下精确解与近似解进行误差分析,得到满足收敛性和稳定性的Holder 型误差估计.1问题的解及正则化方法利用分离变量法可求得问题(1)的解为u(,y)=Zc,cos(nc)cos(y Vk-n)+n=0.c.cos(na)cosh(y/n-k)8n=kJ+1式中Co元2元Cnf()cos(n)dx,(n 1)元J0也可表示为u(,y)=(f,Z,)Z,cosh(yVn?-k)n=0(Z,)%=。是L(0,元)中的标准正交基,(,)表示L(O,元)中的内积.注意到,这里当整数nk时,兰州理工大学学报cosh(y/n-k)=c o s(y/k n),k 表示实
10、数k的整数部分.从解的表达式可以看出,对固定的y,由于当noo时,因式cosh(yn一k)呈指数增长,输人数据f()的微小扰动将可能引起解的巨大误差,所以该问题是严重不适定的通常的数值方法难以解决,必须采用一些特殊的方法,比如正则化方法本文基于 Tikhonov正则化方法给出一种新的正则化方法.定义算子K():L(0,元)L(0,元),问题(1)可写成算子方程的形式K(y)u(,y)=f(),0y0cosh(y/n?-k)改为一1+cosh(y/n-k2)(3)作用.Qin等1通过修改过滤化子把核函数修改为cosh(y/n?-k),提出一种修正的Tikhonovf()d.c1+cosh(/n?
11、-k)(4)正则化方法,并得到正则化参数先验选取的收敛性误差估计.考虑到核函数是指数增长,本文用更为简洁的1过滤化子1+exp(2/n?-k)2化子,过滤掉不适定的高频部分,得到一种新的正则cos(nc)(n1),元第49 卷,从而起到正则化的代替Tikhonov正则化方法,即带有噪声数据的正则化近似解定义为u(,y)=b,cos(ns)cos(y/k-n)+n=0第4期2.b.cos(na)n=KJ+1式中一元下面分别在正则化参数先验和后验两种选取规则下对精确解与近似解进行误差分析.2先验参数选取及误差估计对于不适定问题,常常需要在精确解上给出一些先验假设,否则正则化近似解的收敛速度可能任意
12、的慢.这里作先验假设为Ilu(.,1)I0为固定的常数.引 理 11371)如果0 s1,0,0w1则有余亚辉等:带Neumann边界条件的Helmholtz方程柯西问题的一种新的正则化方法cosh(y Vn?-k)1+exp(2/n-k)()d.元J0f()cos(nc)dc,(n 1)1535BI2k0+s?2)如果0 s1,0,1w2则有(5)定理1设u(,y)和u(,y)分别是问题(1)的精确解和正则化近似解,假设噪声水平(2)和先验条件(6)成立,若取正则化参数=(/E),则有如下误差估计式成立:u(.,y)-u(.,y)3E-(1+o(1),0证明记式(5)中=0 时带有精确数据的
13、正则解为u。(,y),即u.(,y)=Zccos(nc)cos(y/k-n)+n=0(6)8c,cos(n)n=kj+11+exp(2n-k)(9)由于(7)1(8)cosh(y Vn?-k)c.cos(na)cosh(y Vn-k)(1Ilu(.,y)-ua(.,y)=n=k+1C,cos(nc)cosh(Vn-k)cosh(y Vn-k)n=kJ+1+exp(2Vnk)EsPa+exp(-2-r)exp(y-1)Vn-k)1(1-1+exp(2/n-k)cosh(n?-k)cosh(yVn?-k)cosh(/n?-k)令s=exp(一Vn-k),则0 s+1+exp(2Vn?-k)1+ex
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Neumann 边界条件 Helmholtz 方程 问题 一种 正则 方法
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。