结合稀疏逼近的正则化方法求解非齐次双调和方程的Cauchy问题

来源 :中国校外教育·高教(下旬) | 被引量 : 0次 | 上传用户:lollipop1910
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
  在利用边界结点法(BKM)通过径向基函数和Laplace算子、重调和算子的基本解的线性组合来表示问题的解时,需利用已知的一部分边界上的边界条件来推导该线性组合中的待定系数,该过程涉及求解超定线性方程组,由于边界条件给得不充分,Cauchy问题的解不唯一,故需要用正则化方法来得到逼近解析解的结果,我们选择使用稀疏逼近正则化方法(又称lp约束正则化方法),该方法能很好地解决Cauchy问题的不适定性,能很好地反演出跳跃较大的参数部分。本文针对若干具有光滑边界或分段光滑边界的数值算例,验证了该方法的有效性,而且所得的数值计算结果关于噪声是准确的并随已知数据噪声的减小而收敛。
  稀疏逼近正则化Cauchy问题一、引言
  二维双调和方程在很多实际问题中需要用到。但是实际的工程应用中,已知的边界条件可能不完整或者不准确,这样的问题就是反问题,一般来说反问题是不适定的。本文研究的Cauchy问题就是一种反问题,因而在求解过程中,用基本解方法得到的最小二乘问题的解是不唯一的,需要通过使用正则化方法提高原问题数值求解的准确度。
  用基本解方法求解齐次双调和方程会导致离散的Cauchy问题的线性方程不是满秩的或是超定的,其解的适定性存在问题,故本文将使用稀疏逼近的正则化方法来求解离散方程组,以避免直接求解非齐次方程时会出现的不确定性。
  传统的Tikhonov正则化方法是将线性反问题
  并选择二次罚项,以使得近似解具有光滑性。本文选择的稀疏逼近正则化方法,即Φ(x)=x1,同时借助稀疏逼近的优势,在减少结点数目的同时,仍得到较好的结果。
  若该问题有解,则在求解过程中会出现病态的线性方程组,方程组不满秩或者是超定的且条件数很大。本文使用结合稀疏逼近的正则化方法的边界结点法求解满足边界条件(5)的双调和方程(3)或(4)。
  五、总结
  本文使用稀疏逼近正则化方法结合边界结点法求解非齐次双调和方程的Cauchy问题,从数值试验的结果可以看出,这是可行的。并且从以上算例中可以看出,在对源点数、结点数和内点数的关系,以及源点到边界的距离作[3]中的要求的情况下,计算结果与[3]类似,但是精度更高。结果是较准确的,且在该范围内,误差非常小,即使所取点数较多,计算时间仍然较少。当边界条件有较小噪声时,计算结果仍然是稳定的,且随噪声的减小而收敛。
  参考文献:
  [1]Stephen J. Wright, Robert D,Sparse Reconstruction by Separable Approximation,IEEE Transactions on Signal Processing,2009.2479.
  [2]祝家麟.边界元分析.科学出版社,2009.
  [3]刘晓宇,王小军,杜亚楠.用边界结点法(BKM)求解非齐次双调和方程的Cauchy问题,2009,(11):154-166.
  [4]W.Chen and M.Tanaka, New Insights in Boundary-only and Domain-type RBF Methods, Int J Nonlinear Sci Numer Simul,2000,1(3):145-52.
  [5]Xiaolin Li and Jialin Zhu, The method of fundamental solutions for nonlinear elliptic problems, Engineering Analysis with Boundary Elements, 2008,(10):1016.
  [6]Robert Tibshirani.Regression shrinkage and selection via the lasso, Journal Royal Statostical SocietyB,1996(58):267-288.
其他文献
报告30例自发性脑出血后脑血流改变的SPECT研究结果。结果发现全部患者SPECT显示脑血流量的下降区域均超过CT病灶范围,7例血流量下降局限于病灶区,余23例见广泛缺血,其中小脑失联络9例,基底节区出血
急性发病的Creutzfeldt-Jakob病一例报告赵节绪,林世和,江新梅,宋晓南,于志刚,马 君Creutzfeldt-Jakob病(CJD)多隐袭、缓慢起病,病程进展相对较快。我们接诊1例以急性起病的病人,从发病到死亡共4个月,经脑活体组织检查...
通过放免分析法测定28例脑外伤患者血清T3,T4,r-T3及TSH的变化。结果表明:脑外伤患者血清T3,T4,TSH明显低于对照组,血清r-T3明显高于对照组,且各项指标与病情轻重及血肿部位有关。同时死亡组与存活组之间各项
通過对上古“欺骗”义场中11个词项的研究,本文归纳总结了上古“欺骗”类动词的在不同时期的交替演变过程。  上古“欺骗”类动词历时演变上古汉语中表示“欺骗”义的词很多,有欺(諆)、诈、诳(诓、迋)、诬、谩、诡、谲、诞、诒(绐)、罔(誷)、谖等。这些词项具有共同的意义,即通过不真实的言语或行动给对方以假象,哄骗对方,使人上当受骗。它们并非在一个共时的平面存在,它们的分布也是不均匀的,各自出场的时间也不
在对2135柴油机的配气凸轮进行了优化设计的基础上,对其配气系统的振动进行了计算,并且针对凸轮失效的主要原因之一-胶合现象,从摩擦功的角度出发进行了分析计算,并与实验结果进行了对
报告一例克隆病所致急性脑病患者,男,37岁,临床表现头昏、头痛、记忆力减退,发作性意识丧失、抽搐,右半身瘫痪。头部CT见交替出现两侧基底节区低密度影。因肠穿孔行腹部手术,病理检查诊
通过对德国Flensbur大学社会学专业研究生课程(Interkulture)的课程性质、课时数以及从课程的准备、实习和考核总结等各个阶段的描述,指出大学教育要着眼于创新,要走国际化、社会
对55例经 CT 或 MRI 证实的早期脑梗塞患者的脑循环动力学(CVA)和血液粘度进行检测,其中多发性脑梗塞39例,脑干梗塞9例,底节区腔梗7例,结果发现其运动学和动力学的某些指标均
本文测定了60例缺血性脑血管病患者血浆内皮素水平,并与14例正常人对照。结果发现:脑血栓形成急性期患者血浆内皮素升高(P<0.01),而且其浓度与梗塞灶体积大小呈正相关;脑动脉硬化症组血浆内皮