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

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

ID:32800481

大小:58.87 KB

页数:3页

时间:2019-02-15

结合稀疏逼近正则化方法求解非齐次双调与方程cauchy问题_第1页
结合稀疏逼近正则化方法求解非齐次双调与方程cauchy问题_第2页
结合稀疏逼近正则化方法求解非齐次双调与方程cauchy问题_第3页
资源描述:

《结合稀疏逼近正则化方法求解非齐次双调与方程cauchy问题》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库

1、结合稀疏逼近正则化方法求解非齐次双调与方程Cauchy问题在利用边界结点法(BKM)通过径向基函数和Laplace算子、重调和算子的基本解的线性组合来表示问题的解时,需利用已知的一部分边界上的边界条件来推导该线性组合中的待定系数,该过程涉及求解超定线性方程组,由于边界条件给得不充分,Cauchy问题的解不唯一,故需要用正则化方法来得到逼近解析解的结果,我们选择使用稀疏逼近正则化方法(又称lp约束正则化方法),该方法能很好地解决Cauchy问题的不适定性,能很好地反演出跳跃较大的参数部分。本文针对若干具有光滑边界或分段光滑边界的数值算例,验证了该方法的有效性,而且所得的数值计算结果关于噪

2、声是准确的并随已知数据噪声的减小而收敛。稀疏逼近正则化Cauchy问题一、引言二维双调和方程在很多实际问题中需要用到。但是实际的工程应用中,已知的边界条件可能不完整或者不准确,这样的问题就是反问题,一般来说反问题是不适定的。本文研究的Cauchy问题就是一种反问题,因而在求解过程中,用基本解方法得到的最小二乘问题的解是不唯一的,需要通过使用正则化方法提高原问题数值求解的准确度。用基本解方法求解齐次双调和方程会导致离散的Cauchy问题的线性方程不是满秩的或是超定的,其解的适定性存在问题,故本文将使用稀疏逼近的正则化方法来求解离散方程组,以避免直接求解非齐次方程时会出现的不确定性。传统的

3、Tikhonov正则化方法是将线性反问题并选择二次罚项,以使得近似解具有光滑性。本文选择的稀疏逼近正则化方法,即①(x)=xl,同时借助稀疏逼近的优势,在减少结点数目的同时,仍得到较好的结果。若该问题有解,则在求解过程中会出现病态的线性方程组,方程组不满秩或者是超定的且条件数很大。本文使用结合稀疏逼近的正则化方法的边界结点法求解满足边界条件(5)的双调和方程(3)或(4)o五、总结本文使用稀疏逼近正则化方法结合边界结点法求解非齐次双调和方程的Cauchy问题,从数值试验的结果可以看出,这是可行的。并且从以上算例中可以看出,在对源点数、结点数和内点数的关系,以及源点到边界的距离作[3]中

4、的要求的情况下,计算结果与[3]类似,但是精度更高。结果是较准确的,且在该范围内,误差非常小,即使所取点数较多,计算时间仍然较少。当边界条件有较小噪声时,计算结果仍然是稳定的,且随噪声的减小而收敛。参考文献:[l]StephenJ.Wright,RobertD,SparseReconstmctionbySeparableApproximation,IEEETransactionsonSignalProcessing,2009・2479・[2]祝家麟•边界元分析•科学出版社,2009.[3]刘晓宇,王小军,杜亚楠•用边界结点法(BKM)求解非齐次双调和方程的Cauchy问题,2009,(

5、11):154-166.[4]W.ChenandM.Tanaka,NewInsightsinBoundary-onlyandDomain-typeRBFMethods,IntJNonlinearSciNumerSimul,2000,1(3):145-52.[5]XiaolinLiandJialinZhu,Themethodoffundamentalsolutionsfornonlinearellipticproblems,EngineeringAnalysiswithBoundaryElements,2008,(10):1016.[6]RobertTibshirani・Regressi

6、onshrinkageandselectionviathelasso,JournalRoyalStatosticalSocietyB,1996(58):267-28&

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。