正则化解的单位权方差无偏估计公式.PDFVIP

维普资讯 第 27卷 第 6期 武 汉 大 学 学 报 ·信 息 科 学 版 VO1．27No．6 2002年 12月 Geomatic~andInformationScienceofW uhanUniversity Dec．2002 文章编号：1000—050X(2002)06—0604—03 文献标识码 ：A 正则化解的单位权方差无偏估计公式 沈云中 刘大杰 (1 同济大学测量与国土信息工程系，上海市四平路 1239号，200092) 摘 要 ：导 出了．tY-~4化解的单位权方差无偏估计公式，并用第一类Fredholm积分方程这一典型的病态方程作 算例验证 了导出单位权方差无偏估计公式是正确的，可以用于正则化解的单位权方差估计。 关键词：正则化解；单位权方差；无偏估计 中图法分类号：P207．2 20世纪 60年代后，人们逐步认识到大量具 和观测向量及其误差 向量 e∈R 。它们具有如 有数学物理意义的求解 问题属于不适定的模型， 下统计性质 ： 特别是不稳定模型。所谓不稳定模型就是观测数 E(L)=AX，E(e)=0，D(e)= 3P (2) 据很小的误差引起待估参数很大的误差。有关不 式中， 是x的真值；3为单位权方差；P为观测 适定方程的解法 30多年来发表了数以千计的论 值权阵；E和D分别表示期望值算子和方差算 文 ，其基本思想就是通过对不适定方程进行正则 子。观测方程式(1)的最小二乘解可表示如下： 化处理，使待估未知参数有稳定且误差较小的解。 轧 = A～ ． ] 主要的正则化处理方法为Tikhonov正则化方法， dE=AXE—L= (AN A P— )LI 在统计学中称为岭估计，在测量数据处理中常称 广 j r——— ；———一 { 为阻尼最小二乘法。此外还有主元分解法等 。 √ ，D(，L) j Tikhonov正则化方法是 Tikhonov于 1963年提出 式中，下标 L表示最 ／j,,--乘估值；N=ATPA 为 来的不适定方程解法 【1．2J。岭估计 由Hoerl于 法方程系数阵；I，是打z×打z的单位阵。如果观 1962年提出，并于 1970年作 了详细的阐述 3【J， 测值的单位权方差 8未知，可用其无偏估值 它是Tikhonov正则化方法的特例。但岭估计方 代替。 法在有偏估计意义下给出了衡量未知参数估值精 如果观测方程 (1)是病态的，法方程系数矩阵 度的指标，即均方误差。徐培亮引入了广义岭回 N 的特征值单调地趋向于0，使得很小的观测误差 归，其特点是法方程各主对角元上加不同的小参 导致参数估值较大的误差，因此，最小二乘解不能 量。有关正则化解一般只讨论参数估值的有偏 获得满意的结果。为得到稳定的解必须采用某个 性，没有见到讨论正则化引起单位权方差估值有