测绘学报 ›› 2017, Vol. 46 ›› Issue (6): 679-688.doi: 10.11947/j.AGCS.2017.20160268

• 大地测量学与导航 • 上一篇    下一篇

顾及截断偏差影响的TSVD截断参数确定方法

林东方1,2, 朱建军1,2, 宋迎春1,2   

  1. 1. 中南大学地球科学与信息物理学院, 湖南 长沙 410083;
    2. 中南大学有色金属成矿预测与地质环境监测教育部重点实验室, 湖南 长沙 410083
  • 收稿日期:2016-06-12 修回日期:2017-05-11 出版日期:2017-06-20 发布日期:2017-06-28
  • 通讯作者: 朱建军 E-mail:zjj@csu.edu.cn
  • 作者简介:林东方(1986—),男,博士生,研究方向为测量平差数据处理及应用。E-mail:Lindongfang223@163.com
  • 基金资助:
    国家自然科学基金(41531068;41474008;41574006;41674012)

Truncation Method for TSVD with Account of Truncation Bias

LIN Dongfang1,2, ZHU Jianjun1,2, SONG Yingchun1,2   

  1. 1. School of Geosciences and Info-physics, Central South University, Changsha 410083, China;
    2. Key Laboratory of Metallogenic Prediction of Nonferrous Metals and Geological Environment Monitoring, Central South University, Changsha 410083, China
  • Received:2016-06-12 Revised:2017-05-11 Online:2017-06-20 Published:2017-06-28
  • Supported by:
    The National Natural Science Foundation of China (Nos. 41531068;41474008;41574006;41674012)

摘要: TSVD通过截断参数截掉较小的奇异值来改善病态性对估计的影响,其本质是通过引入少量偏差来降低方差,以提高估值的稳定性和可靠性。截断参数是影响TSVD解算效果的关键因素,常用的广义交叉核实法(GCV法)和L曲线法未从TSVD改善模型参数估值质量的角度确定截断参数,稳定性和可靠性不足,而最小MSE法理论依据充分但受限于MSE计算的准确性。通过分析TSVD由小到大截掉奇异值后,相应的估值方差与偏差变化,本文提出了引入偏差量小于降低方差量来确定截断参数的思想,并通过估计出较大奇异值截掉后的偏差引入量建立偏差估值可信区间,利用可信区间内偏差估值与方差下降量进行比较,避免较小奇异值截掉后的方差下降量与偏差引入量的直接比较,从而解决参数真值未知截掉较小奇异值引入偏差量难以准确计算的问题。最后通过试验验证了新方法的可行性和有效性,相比于GCV法和L曲线法,新方法确定的截断参数稳定性和可靠性更高,可有效提高TSVD的解算效果。

关键词: TSVD, 截断参数, 有偏估计, 方差, 偏差

Abstract: TSVD truncate small singular values by truncation parameter to improve the parameter estimation of ill-posed model. From the perspective of MSE (mean squared error), TSVD introduce biases to reduce variances, therefore the stability and reliability of the solution can be improved. Truncation parameter is key factor of TSVD, but it is difficult to determine in case of the gently declined singular values. The parameter determined by GCV (generalized cross-validation) and L-curve often unstable and unreliable. And the minimum MSE method is limited by the accuracy of the estimated MSE. This paper compares the changes of variance and bias produced by truncating the singular values in turn and determines the truncation parameter when the reduced variance is less than the introduced bias. In order to avoid the comparison between reduced variance and introduced bias of truncating small singular values, the confidence domain of bias is established through estimating the introduced bias of truncating big singular values that are proved to be reliable. The comparisons are replaced by comparing the reduced variance with the bias in the confidence domain. Therefore, the issue of introduced bias of truncating small singular values cannot be calculated without true values of unknown parameters is solved. Numerical examples proves the feasibility and effectiveness of the new method. Truncation parameters determined by new method are more stable and reliable than GCV and L-curve and improve the TSVD solution effectively.

Key words: TSVD, truncation parameter, biased estimation, variance, bias

中图分类号: