一种计算补偿最小二乘正则化参数的最优化方法.pdfVIP

第 35卷 第 1期 大 地 测 量 与 地 球 动 力 学 Vo1．35NO．1 2015年 2月 Journa1ofGeodesyandGeodynamics Feb．，2015 文章编号 2015)01—0115—03 一 种计算补偿最小二乘正则化参数的最优化方法 丁士俊 朱 留洋 姜卫平。 1 武汉大学测绘学 院，武汉市珞喻路 129号 ，430079 2 国家测绘地理信息局精密工程测量与工业测量重点实验室 ，武汉市珞喻路 129号 ，430079 3 武汉大学卫星导航定位系统技术研究 中心 ，武汉市珞喻路 129号 ，430079 摘 要 ：在半参数模型补偿最小二乘估计中，如何有效地计算正则化参数是求解的关键。在建立一种新的补 偿最小二乘模型的基础上，基于广义交叉核实与L曲线准则，采用最优化方法计算正则化参数。仿真计算表 明，本文方法是一种计算正则化参数更为有效的方法。 关键词 ：半参数模型；补偿最小二乘 ；广义交叉核实；I曲线 ；正则化参数 ；最优化方法 中图分类号 ：P207 文献标识码 ：A 在半参数模型估计方法 中，补偿最小二乘估 式中，L为 72维观测 向量 ，x为 t维参数 向量，B为 计是一种比较典型 的求解方法 ，如何快速有效地 列满秩设计矩阵，s为 维非参数向量 (用来描述函 计算正则化参数是求解的关键之一。在以往的研 数模型误差)，A为n维误差 向量 (噪声)，△～N(O， 究中，对于正则化参数求解没有给 出具体而有效 P )， 是方差因子 (测量 中又称为单位权中误 的算法_1]。研究表明，确定正则化参数的方法不 差)，P为厶的权矩 阵。模型的任务是 以观测数据 同，所得到的正则化参数也不相同。另一方面 ，正 L去估计未知参数 x、未知函数 与方差因子 。 则化参数的大小往往和噪声及信号大小有直接关 由式 (1)可得误差方程： 系 。理论上 ，正则化参数 的取值范 围在 (0，1)区 V—B + ～L (2) 间，但实际情况并非如此 ，在计算过程 中其有效取 半参数补偿最ebZ．乘是在最小二乘法的目标函数 值范围难以控制 ，给计算带来不便 。基于上述原 上增加一个非参数部分 的补偿项，即求解如下条 因，本文尝试改变补偿最小二乘估计 的准则形式， 件极值问题_l】]： 将原极值准则V PV+ Rs—min改化为 dV PV V PV+ R 一min (3) +(1--a)sRs—min，这样正则化参数的取值范围 式中，R为正定或半正定矩阵，a为正则化参数 (又 将被控制在 (0，1]，不会 因为其他因素而有所改 称光滑参数)，在极小化过程中，在二次型VTPV与 变 ，参数有了明确 的上限值和下限值。另外，在原 船 之间起平衡作用。由式(2)、(3)可得法方程： 极值准则下求解正则化参数，通常是采用人工干预 广B PB B P ]厂 ] rB PL] 的方式 ，通过迭代计算寻找极值点，计算量大且效 lPBP+Jlj—lPLJ 率低。在新极值准则下 ，由于正则化参数有明确 的 对式(4)求解 ，可得到半参数模型补偿最小二乘估 范围，采用相关 的最优化算法理论 (如黄金分割搜 计计算公式 ]： 索法等)计算正则化参数 ，计算效率得到明显提高。 f爻一(B