最小二乘配置的SVD分解解法.docxVIP

最小二乘配置的SVD分解解法 最小二乘配置的SVD分解解法 第３３卷第３期 ２００８年５月 ＳｃｉｅｎｃｅｏｆＳｕｒｖｅｙｉｎｇａｎｄＭａｐｐｉｎｇ Ｖ０１．３３Ｎｏ．３ 最小二乘配置的ＳＶＤ分解解法 鲁铁定①②，宁津生①，周世健鳓，臧德彦② （①武汉大学测绘学院，武汉４３００７９；②东华理工大学地球科学与测绘工程学院， 江西抚州３４４０００；③江西省科学院，南昌３３００２９） 【摘要】最小二乘配置最初是在组合各种资料来研究地球形状与重力场的一种数学方法，目前最小二乘配置已经在测绘数据处理中得到广泛应用。本文首先分析了目前采用的最小二乘配置法解算方法，在讨论了矩阵的奇异值分解（ＳｉｎｇｕｌａｒＶａｌｕｅＤｅｃｏｍｐｏｓｉｔｉｏｎ，ＳＶＤ）方法的基础上，推导得出了矩阵ＳＶＤ分解与广义逆矩阵的关系，得出了可以直接利用ＳＶＤ分解求解矩阵的Ｍｏｏｒｅ．Ｐｅｎｒｏｓｅ广义逆，并推导了应用ＳＶＤ分解求解最小二乘配置的估值计算公式和精度估算公式，最后通过重力异常实例进行了计算，得出矩阵的ＳＶＤ分解用于最小二乘配置解算的正确性和可行性，为最小二乘配置的求解提供了一种新方法。 【关键词】最小二乘配置；奇异值分解（ＳＶＤ）；重力异常 【文章编号】１００９－２３０７（２００８）０３－００４７－０５【中图分类号】Ｐ２２【文献标识码】Ａ ＤＯＩ：１０．３７７１／ｊ．ｉｓｓｎ．１００９－２３０７．２００８．０３．０１６ 最小二乘配置起源于根据最小二乘推估来内插和外推重 力异常的课题¨．２Ｊ。在地球形状重力场的研究中，配置的普遍形式是其函数模型中除包含随机部分外，还包含非随机的系统部分（也称为倾向）。１９６９年，克拉鲁普（Ｔ．Ｋｒａｒｕｐ）把推估重力异常的方法，发展为用不同类型的数据，例如重力异常，垂线偏差等，去估计异常引力场中的任一元素，例如扰动位。大地水准面差距等，提出了最小二乘配置法。莫里兹（｝ＬＭｏｒｉｔｚ）对最小二乘配置进行了系统深入的研究，提出了带系统参数的最ｉｂ－乘配置，并概述了这种方法在大地测量其他方面的应用，进而导致几何位置和重力场的最小二乘联合求定，为整体大地测量奠定了理论基础。从２０世纪６０年代以来，国内外测绘专家学者对最小二乘配置法进行了广泛深入研究，汪洪海、罗志才、罗佳等对不同分辨率的重力数据融合问题，推导了多分辨率最小二乘配置法的具体公式，并通过算例与传统最小二乘配置法的结果进行了比较‘３１；罗志才、宁津生、晁定波研究了频域最ｔＪ、Ｚ．乘配置法，并与空域最小二乘配置法进行了比较，其方法可有效提高稳定性【４ｏ；姚宜斌、陶本藻从最／ｊ．，－－乘配置出发，推导了附加先验约束的高精度ＧＰＳ网平差的参数估计和精度估计公式一１；张继贤等讨论了顾及线性化残差相关性的最小二乘配置法地形迭代线性化技术№１；篮悦明、陶本藻讨论了数字化曲线拟合的最小二乘配置法＂ｏ；张希、江在森研究了最／ｊ，－－ 乘配置法用于地壳形变的问题，并将空间域最小二乘配置内插推广到时间域‘副；张献州根据高斯马尔可夫模型两个备选假设理论，给出了最ｔＪ、－－乘配置模型中观测粗差的可区分且可发现的表达式一１；潘雄、孙海燕导出了最小二乘配置模型中正规化矩阵正定时参数平差的计算方法，得到了参数和信 号的估计茸ｍ１。 矩阵的分解（ｄｅｃｏｍｐｏｓｉｔｉｏｎ）是通过线性变换，将某个给定或者已知的矩阵分解为二个或三个矩阵标准型的乘积（个别情况下分解为两个矩阵标准型之和）。根据矩阵分解后的矩阵的标准型，其分解大致分为对角化分解、三角化分解、三角．对角化分解和三对角化分解四大类型。矩阵的奇异值分解（ＳｉｎｇｕｌａｒＶａｌｕｅＤｅｃｏｍｐｏｓｉｔｉｏｎ，ＳＶＤ）是矩阵对角化分解的一种类型。矩阵奇异值分解自从Ｂｅｈｒａｍｉ和Ｊｏｒｄａｎ两位学者创立以来…。…，已经成为数值线性代数的最有用和最有效的工具之一，它已在最优化问题、统计分析、信号与图像处理、系统理论和控制等很多领域被广泛地应用¨“。本文在分析最小二乘配置法解算方法的基础上，推出一种基于ＳＶＤ矩阵分解的解算算法，并通过实例进行了分析，得出方法的可行性。 ２最小二乘配置的估值公式 ２．１最小二乘配置的数学模型 最小二乘配置的函数模型¨二１一般为 Ｌ＝ＢＸ４－ＧＹ＋４ 式中Ｌ为观测向量，４为观测噪声，】，为倾向参数，ｘ为滤波信号，另外还有推估信号∥，其随机模型是： 作者简介：鲁铁定（１９７４一），男，陕西富平人，副教授，博士生，主要研究方向：测绘数据处理理论和物理大地测量。 Ｅ．ｍａｉｌ：ｔｄｌｕ＠ｅｃｉｔ．ｅｄｕ．ｃｎ Ｘ和ｘ’的先验期望：