用U曲線法確定嶺估計中的嶺參數

魯洋為，王振杰

用U曲線法確定嶺估計中的嶺參數

魯洋為，王振杰

(中國石油大學(華東) 地球科學與技術學院，山東 青島 266580)

在嶺估計中，嶺參數的確定是關鍵。本文針對嶺估計中嶺參數確定困難的問題，提出采用U曲線法確定嶺參數，并將U曲線法與L曲線法、GCV法進行了比較，結果表明：用U曲線法確定嶺參數能獲得較高精度的解，與L曲線法及GCV法相比，具有一定的優勢。

嶺估計；嶺參數；U曲線法；L曲線法；GCV法

0 引言

在測量數據處理中，一般采用有偏估計來改善病態條件下的最小二乘解，其中嶺估計是常用的方法[1-8]。在嶺估計中，嶺參數的選擇是關鍵，不同的嶺參數，解算的結果會有所差別。目前，嶺參數的確定方法主要有嶺跡法[9-10]、廣義交叉核實法(generalized cross validation，GCV)[11-13]和L曲線法[14-18]。嶺跡法比較直觀，但是具有一定的主觀性；GCV法在理論上能夠獲取最優的嶺參數，但有時GCV函數的變化過于平緩，定位其最小值很困難[9-10]；L曲線法可以比較容易的獲得嶺參數，并且能夠得到比較精確的解，但L曲線法存在求解結果過于依賴曲線擬合精度[19]、求解過程可能不收斂等缺點[20]。本文研究采用U曲線法來確定嶺參數[19-22]，并通過算例，將U曲線法與L曲線法、GCV法進行了比較，分析了U曲線法的性能，得出了一些有益的結論。

1 U曲線法確定嶺參數的基本原理

1.1 嶺估計的原理

設有觀測方程

L=AX+Δ

(1)

(2)

由式(1)和式(2)可得嶺估計的解為

(3)

1.2U曲線法的基本原理

嶺參數的選擇方法很多，其中最常用的方法是GCV法[11]和L曲線法[14-15]。GCV法確定嶺參數，則要求GCV函數的最小值，由(1)、式(3)得GCV函數為：

(4)

式(4)中，H(α)=A(ATPA+αI)-1ATP，n為觀測值個數，tr為矩陣的跡。

當式(4)取最小值時的α即為GCV法所確定的嶺參數，GCV法的優點是，當式(4)的最小值存在時，可以選擇一個最優的嶺參數；它的缺點是，GCV函數有時變化過于平緩，這時定位它的最小值很困難。

U曲線法與L曲線法有所不同，U曲線法根據定義的U(α)函數，得到U(α)-α曲線，曲線左側曲率最大點的α即為U曲線法所確定的嶺參數。U曲線函數的定義為：

(5)

U曲線法和L曲線法的關鍵都是定位曲線中曲率最大的點，對于U曲線，其曲率可以表示為：

(6)

式(6)中，U(α)′、U(α)″分別為U(α)的一階導數和二階導數，并且

(7)

(8)

式(7)、式(8)中x(α)、y(α)的一階導數、二階導數x(α)′、x(α)″、y(α)′、y(α)″的計算可參考文獻[21]。

由文獻[19，21]可知，U(α)函數具有以下性質：

(1)U曲線在區間 (0，)內嚴格遞減，在區間(，+∞)內嚴格遞增；

(2)U曲線在區間(，)內存在局部極小值，并且當設計矩陣A只有一個非零奇異值時，U(α)存在唯一的極小值；

采用U曲線法確定嶺參數時，對式(6)在區間 (，)內求最大值，定位出U曲線左側曲率的最大值點，該點所對應的α即為所求的嶺參數。

由以上分析可以看出，U曲線法不需要進行曲線擬合，其計算量要小于L曲線法的計算量，也不會出現曲線發散，無法確定嶺參數的現象，與L曲線法及GCV法相比，U曲線法更適合用于確定嶺參數。

2 算例分析

圖1 用U曲線法確定嶺參數的示意圖

合；GCV法在本例中也能夠準確的定位出GCV函數的最小值點，確定嶺參數，沒有出現發散的現象。

圖2 用L曲線法確定嶺參數的示意圖

圖3 用GCV法確定嶺參數的示意圖

基于這三種方法的嶺估計結果如表1所示，可以看出，基于L曲線法和GCV法的嶺估計解算的未知數估值與真值之差的二范數分別為：0.356 1、0.373 3，即L曲線法略優于GCV法；而基于U曲線法的嶺估計解算出的未知參數的估值與真值之差的二范數為0.270 6，其精度明顯優于L曲線法和GCV法的精度。

表1 嶺參數及‖△‖

3 結束語

U曲線法具有理論嚴密，定位準確的特點，與L曲線法相比，具有計算量小，不需要進行曲線擬合的優勢；與GCV法和L曲線法不同，U曲線法不會出現無法確定嶺參數的現象；基于U曲線法的嶺估計結果的精度優于基于L曲線法、GCV法的嶺估計結果的精度，值得進一步研究與推廣應用。

[1] HOERL A E，KENNARD R W.Ridge Regression：Biased Estimation for Non-orthogonal Problems[J].Technometrics，1970，12(1)：55-67.

[2] 崔希璋，於宗儔，陶本藻，等.廣義測量平差[M].武漢：武漢大學出版社，2000.

[3] 歸慶明，李國重.嶺壓縮組合估計及其在測量平差中的應用[J].大地測量與地球動力學，2002，22(1)：16-21.

[4] 王樂洋，許才軍，魯鐵定.病態加權總體最小二乘平差的嶺估計解法[J].武漢大學學報：信息科學版，2010，35(11)：1346-1350.

[5] 丁海勇，楊延珍，王新志，等.基于移動曲面擬合方法的GPS高程轉換[J].大地測量與地球動力學，2010，30(6)：1-4.

[6] 何峰，王剛，劉利等.地球靜止軌道衛星廣播星歷參數擬合與試驗分析[J].測繪學報，2011，40(增刊)：52-58.

[7] 韓松輝，杜蘭，歸慶明，等.診斷復共線性的特征分析法及其在GEO定軌中的應用[J].測繪學報，2013，42(1)：19-26.

[8] 趙坷，付政慶，劉國林.加權半參數模型及其應用效果分析[J].大地測量與地球動力學，2014，34(2)：100-104.

[9] 王振杰.大地測量中不適定問題的正則化解法研究[D].武漢：中國科學院測量與地球物理研究所，2003.

[10]王振杰，歐吉坤.用L-曲線確定嶺估計中的嶺參數[J].武漢大學學報：信息科學版，2004，29(3)：235-238.

[11]GOLUB G H，HEATH M，WAHBA G.Generalized Cross-validation as a Method for Choosing a Good Ridge Parameter[J].Technometrics.1979，21(2)：215-223.

[12]王樂洋，許才軍.附有病態約束反演問題的嶺估計法[J].測繪學報，2011，36(5)：612-616.

[13]王劭逸.基于GCV方法的線性回歸模型嶺參數估計[D].合肥：中國科學技術大學，2012.

[14]HANSEN P C.Analysis of Discrete Ill-posed Problems by Means of the L-curve[J].SIAM Review，1992，34(4)：561-580.

[15]HANSEN P C，O’LEARY D P.The Use of the L-curve in the Regularization of Discrete Ill-posed Problems[J].SIAM J.Sci.Comput，1993，14(6)：1487-1503.

[16]袁修孝，林先勇.基于嶺估計的有理多項式參數求解方法[J].測繪學報，2008，33(11)：1130-1133.

[17]王彬，高井祥，劉超，等.L曲線法在等價權抗差嶺估計模型中的應用[J].大地測量與地球動力學，2012，32(3)：97-101.

[18]劉斌，龔健雅，江萬壽，等.基于嶺參數的譜修正迭代法及其在有理多項式參數求解中的應用[J].測繪學報，2012，37(4)：399-402.

[19]袁強強，沈煥鋒，李平湘，等.自適應正則化多幅影像超分辨率重建[J].中國圖像圖形學報，2010，15(12)：1720-1727.

[20]XU Pei-liang.Truncated SVD Methods for Discrete Linear Ill-posed Problems[J].Geophysics Journal International，1998(135)：505-514.

[21]KRAWCZYK-STADO D，RUDNICKI M.Regularization Parameter Selection in Discrete-ill Posed Problems- the U-curve[J].Int.J.Appl.Match.Comput.Sci，2007，17(2)：157-164.

[22]CHAMORRO-SERVENT J，AGUIRRE J，RIPOLL J，et al.Feasibility of U-curve Method to Select the Regularization Parameter for Fluorescence Diffuse Optical Tomography in Phantom and Small Animal Studies[J].Optics Express，2011，19(12)：11490-11506.

Determining the Ridge Parameter in Ridge Estimation Using U-curve Method

LUYang-wei，WANGZhen-jie

(School of Geosciences，China University of Petroleum，Qingdao 266580，China)

The determination of ridge parameter is a key problem in ridge estimation，and the U-curve method is used to determine ridge parameter in this paper.We also compare the U-curve method with L-curve method and GCV method using an example and the results show that the U-curve method can get higher accuracy and have an advantage over L-curve method and GCV method.

ridge estimation；ridge parameter；U-curve method；L-curve method；GCV method

魯洋為，王振杰.用U曲線法確定嶺估計中的嶺參數[J].導航定位學報,2015,3(3):132-134+138.(LU Yang-wei,WANG Zhen-jie.Determining the Ridge Parameter in Ridge Estimation Using U-curve Method[J].Journal of Navigation and Positioning,2015,3(3):132-134+138.)

10.16547/j.cnki.10-1096.20150326.

2015-05-18

國家自然科學基金(41374008)。

魯洋為(1990—)，男，陜西安康人，碩士生，研究方向為GNSS定位與測量數據處理。

P228

A

2095-4999(2015)-03-0132-03