2.1平滑因子的选取方法在半参数模型中，平滑因子α是一个重要的待定参数，它起到拟合程度和光滑程度的平衡作用，平滑因子的选取是否得当对估计量有很大的影响，一般采用广义交叉核实法。

　　GCV（α）=■（10）

　　式中tr（H（α）代表帽子矩阵H（α）的迹。

　　2.2正则化矩阵的求法为了求解非参数量S，在重力测量中，重力异常的影响随距离的增加而减弱，R通常选取两点间的距离d相关的量：R■=d■■（11）

　　上式中d■为重力异常点之间的距离，点d■x■，y■d■x■，y■距离：d■=■（12）

　　3算例分析

　　本算例取自文献[4]，分别应用最小二乘配置，多面函数和半参数对一测区内重力异常数据进行推估。我们选取了已知点点号为1～8，推估未知点点号为9～28，图1显示了它们的坐标关系[4]。

　　应用半参数模型中的补偿最小二乘法计算已知点重力异常估值和未测点的重力异常估值，空间重力异常与地面点高程有密切的联系，在局部重力异常计算时，重力异常不仅含有系统部分，还有随机部分，系统部分可以表示为高程H的函数[5]，Ti=X1+X2Hi，观测方程为：L=AX+S+Δ。

　　其中采用广义交叉核实法选取平滑因子，计算得α=0.1，正规化矩阵采用距离选取法，R■=d■■，根据半参数模型公式可以求得：

　　■=（-70.8409，0.0908）T，

　　s=（-2.621，1.335，-1.421，1.947，-1.526，0.09435，-4.966，

　　2.54）T

　　补偿最小二乘法与最小二乘配置法拟合值比较

　　由表3可以看出，应用半参数补偿最小二乘法，已知点拟合值与真实值极为接近，因为非参数分量S合理的解释了该模型的模型误差部分，所以残差很小，说明半参数模型有较强的适应性。平差精度也有了一定的提高，中误差由±1.85提高到±0.04。说明半参数模型在提高精度的同时，可以从观测量中分离出非参数分量S，该方法是最小二乘配置法的改进。

　　补偿最小二乘法与最小二乘配置法推估值比较

　　由表4可以发现应用半参数模型推估值与最小二乘配置法计算的推估值非常接近，满足中等山区重力异常的精度要求，符合山区复杂情况下的重力异常分布，此算例说明半参数模型应用在局部重力异常的插值格网化计算中，方法是可行的。

　　4结论

　　本文结合半参数模型在数据处理中的优越性，把半参数模型应用到重力异常格网化中。研究了半参数模型的原理和解法，详细介绍了基于正则化矩阵的补偿最小二乘原理，并推导了其参数求解方程式，介绍了平滑因子和正则化矩阵的求法。结合半参数模型特点建立了基于半参数的格网化模型，最后通过算例，验证了半参数应用到重力异常格网化中的可行性。

　　参考文献：

　　[1]汪隆六.三角形线性插值在区域重力测量数据网格化中的应用[J].物化探计算技术，1986（02）：141-146.

　　[2]丁士俊，陶本藻.半参数回归与平差模型[J].大地测量与地球动力学，2003（04）：111-114.

　　[3]孙海燕，吴云.半参数回归与模型精化[J].武汉大学学报（信息科学版），2002（02）：172-174.

　　[4]陆仲连.关于在局部扰动重力场中最小二乘配置法的应用[J].解放军测绘学院学报，1984（00）：3-16.

　　[5]彭泽辉，李辉，申重阳，等.基于最小二乘配置的重力变化插值方法[J].大地测量与地球动力学，2010（03）：43-46.

---

期刊库（http://www.zgqkk.com），是一个专门从事期刊推广、投稿辅导的网站。
　　本站提供如何投稿辅导，寻求投稿辅导合作，快速投稿辅导，投稿辅导格式指导等解决方案：省级投稿辅导/国家级投稿辅导/核心期刊投稿辅导//职称投稿辅导。

---

　　【免责声明】本文仅代表作者本人观点，与投稿辅导_期刊发表_中国期刊库专业期刊网站无关。投稿辅导_期刊发表_中国期刊库专业期刊网站站对文中陈述、观点判断保持中立，不对所包含内容的准确性、可靠性或完整性提供任何明示或暗示的保证。请读者仅作参考，并请自行承担全部责任。