基于正则化方法一个新型数值微分算法.pdfVIP

2009年 9月 高 等 学 校 计 算 数 学 学 报 第 3l卷第 3期 基于正则化方法的一个新型数值微分算法 王希云 (山西太原科技大学数学系，太原 030024) 黄建国 陈宇 (上海交通大学数学系，上海 200240) A NEW NUM ERICAL D IFFERENTI T10N ALGoRITHM W ITH REGULARIZATIoN WangXiyun (TaiyuanuniversityofscienceandTechnology,Taiyuan030024) HuangJianguo ChenYu (ShanghaiJiaoTongUniversity，Shanghai200240) Abstract Inthispaperanew numericaldifferentiationalgorithm withregular— ization，based On theprescribed datawith measurementerrorisproposed． It retainsthesim ilarglobalerrorestimateas fortheexistingalgorithm ，butavoids somedrawbackofthelatteronefbadapproximationneartheendpointsofthe regionofinterest1．Thenew algorithm hasauniquesolutionandthecorrespond— ingcharacterizationconditionsareestablished．Inaddition，theerrorestimate is developedinthesquareintegrablenorm．Finally，numericaltestsareprovidedto sh0w theeffectivenessofouralgorithm ．Inparticular，intheCaseofsmallsizeof samples，ouralgorithm performsbetterthantheoriginalalgorithm ． Keywords numericaldifferentiation，regularizationmethods，objectivefunc- tional，regularizationparameters． AMS(2000)subjectclassifications 65D25 中图法分类号 0242．2 本文受国家自然科学基金资助 收稿 日期：2007-01—17． 2009年 9月 高 等 学 校 计 算 数 学 学 报 ．247． 1 引 言 根据带误差的测量数据重构函数使其微分能较好的拟合精确函数微分，是一个有重 要意义的研究课题，在图像处理、计算机视觉和计算力学等领域中均有应用 f2,3，4】． 如果测量数据没有误差，上述问题可归结于一维函数的拟合问题，常规的方法有 La_ grange插值法和样条函数方法等 [5]；也有直接的数值微分计算公式 [6，7]．这些方法的计 算结果都较好．但当测量数据有偏差时，即使是样条函数方法效果也不十分理想．其原因 在于重构函数虽然拟合样本点数据好，但是振荡很厉害，会产生所谓的Gibbs现象 [8】， 因此不能得到符合实际问题需要的合理计算结果． 求解带偏差数据函数拟合的一个有效方法是正则化方法 9【，10】，在文献 1【,2]中有对 用正则化方法重构数值微分问题的深入研究，在文献 2『,31中还将其应