• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    多元總體最小二乘問題的牛頓解法

    2016-05-16 08:38:45王樂洋趙英文陳曉勇臧德彥
    測(cè)繪學(xué)報(bào) 2016年4期

    王樂洋,趙英文,陳曉勇,臧德彥

    1. 東華理工大學(xué)測(cè)繪工程學(xué)院,江西 南昌 330013; 2. 流域生態(tài)與地理環(huán)境監(jiān)測(cè)國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室,江西 南昌 330013; 3. 江西省數(shù)字國(guó)土重點(diǎn)實(shí)驗(yàn)室,江西 南昌 330013

    ?

    多元總體最小二乘問題的牛頓解法

    王樂洋1,2,3,趙英文1,陳曉勇1,2,3,臧德彥1,2,3

    1. 東華理工大學(xué)測(cè)繪工程學(xué)院,江西 南昌 330013; 2. 流域生態(tài)與地理環(huán)境監(jiān)測(cè)國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室,江西 南昌 330013; 3. 江西省數(shù)字國(guó)土重點(diǎn)實(shí)驗(yàn)室,江西 南昌 330013

    transformation

    Foundation support: The National Natural Science Foundation of China (Nos. 41204003; 41161069; 41304020; 41464001); The National Department Public Benefit Research Foundation (Surveying,Mapping and Geoinformation)(No. 201512026); The Natural Science Foundation of Jiangxi Province (No. 20151BAB203042); The Science and Technology Project of the Education Department of Jiangxi Province (Nos. GJJ150595; KJLD12077; KJLD14049); The Found of Key Laboratory of Watershed Ecology and Geographical Environment Monitoring (No. WE2015005); The Scientific Research Foundation of ECIT(No. DHBK201113); Innovation Fund Designated for Graduate Students of Jiangxi Province (Nos. YC2015-S266; YC2015-S267); Innovation Fund Designated for Graduate Students of ECIT(No. DHYC-2015005); The Found of the Key Laboratory of Mapping from Space, NASG (No.K201502)

    摘要:為提高多元總體最小二乘問題參數(shù)估值的解算效率,推導(dǎo)了基于牛頓法的多元加權(quán)總體最小二乘算法;分析比較了基于牛頓法的多元加權(quán)總體最小二乘解和基于拉格朗日乘數(shù)法多元加權(quán)總體最小二乘解之間的關(guān)系,根據(jù)協(xié)因數(shù)傳播律給出了多元總體最小二乘平差的16種協(xié)因數(shù)陣的近似計(jì)算公式。新算法能夠解決觀測(cè)矩陣和系數(shù)矩陣元素具有相關(guān)性的問題,并且可以把觀測(cè)矩陣和系數(shù)矩陣的隨機(jī)元素和常數(shù)元素納入到一個(gè)協(xié)因數(shù)陣中進(jìn)行處理。算例結(jié)果表明,本文提出的多元總體最小二乘問題的牛頓解法可行且收斂速度更快。

    關(guān)鍵詞:多元總體最小二乘;牛頓法;協(xié)因數(shù)傳播律;仿射變換

    總體最小二乘(total least squares,TLS)方法是近30多年來發(fā)展起來的一種能同時(shí)顧及觀測(cè)值誤差和模型系數(shù)矩陣誤差的數(shù)學(xué)方法,其理論及應(yīng)用研究是目前國(guó)內(nèi)外研究的熱點(diǎn)問題[1]。文獻(xiàn)[2]首次提出總體最小二乘概念,隨著變量誤差模型(errors-in-variables,EIV)由單變量模型擴(kuò)展到多元變量模型[3],總體最小二乘也相應(yīng)擴(kuò)展到多元總體最小二乘(multivariate TLS,MTLS)[4-5]。某些模型的系數(shù)矩陣會(huì)因未知參數(shù)寫成向量形式含有重復(fù)出現(xiàn)的隨機(jī)元素和大量的零元素;多元總體最小二乘將參數(shù)向量和觀測(cè)向量改寫成矩陣形式,避免了這種現(xiàn)象,進(jìn)而削弱了系數(shù)矩陣中行列之間的相關(guān)性,而且使得系數(shù)矩陣的構(gòu)造更加簡(jiǎn)單。

    在國(guó)外相關(guān)研究中,文獻(xiàn)[5—7]認(rèn)為觀測(cè)矩陣和系數(shù)矩陣中的隨機(jī)元素是等精度的,給出了基于奇異值分解法和拉格朗日乘數(shù)法的多元總體最小二乘算法;文獻(xiàn)[8]討論了基于拉格朗日乘數(shù)法的多元加權(quán)總體最小二乘算法,將隨機(jī)元素的方差陣分解成兩個(gè)矩陣取直積的形式,但是當(dāng)觀測(cè)數(shù)據(jù)含有兩個(gè)以上的方差分量或是觀測(cè)數(shù)據(jù)含有多類觀測(cè)值時(shí),隨機(jī)變量的方差協(xié)方差陣就很難寫成或是不能寫成兩個(gè)矩陣取直積的形式;文獻(xiàn)[9]推導(dǎo)了基于奇異值分解法的多元加權(quán)總體最小二乘算法,但是這種算法只使用了行尺度權(quán)陣,其權(quán)陣不具有一般性。以上文獻(xiàn)的所有算法都不適用于系數(shù)矩陣與觀測(cè)矩陣中隨機(jī)元素具有相關(guān)性的情況。文獻(xiàn)[10]給出了基于拉格朗日乘數(shù)法的多元加權(quán)總體最小二乘算法,其權(quán)陣具有一般性,但沒有具體給出系數(shù)矩陣含有常數(shù)元素的解法以及沒有解決算法的精度評(píng)定問題,這種算法的解算效率也還需要驗(yàn)證。在國(guó)內(nèi)相關(guān)研究中,文獻(xiàn)[11]將基于奇異值分解法的多元總體最小二乘算法應(yīng)用在三維坐標(biāo)轉(zhuǎn)換中;文獻(xiàn)[12]給出的多元總體最小二乘算法受到正交Procrustes算法制約僅適用于觀測(cè)值誤差對(duì)稱獨(dú)立的情況,隨機(jī)元素相關(guān)性問題也沒有解決;文獻(xiàn)[13]只給出了觀測(cè)值等精度同方差情況下的多元總體最小二乘算法。對(duì)于系數(shù)矩陣含有常量的情況,文獻(xiàn)[10,14—16]將系數(shù)矩陣拆分成分別包含隨機(jī)元素和常數(shù)元素的系數(shù)矩陣;文獻(xiàn)[11—13]采用消元法,分兩步計(jì)算隨機(jī)元素對(duì)應(yīng)的未知參數(shù)和常數(shù)元素對(duì)應(yīng)的未知參數(shù);文獻(xiàn)[8,17]使常數(shù)元素的方差和協(xié)方差為零,然后進(jìn)行加權(quán)解算;文獻(xiàn)[18]首次將EIV模型改寫成更具一般性的PEIV(partial EIV)模型,解決了系數(shù)矩陣任意位置含有常數(shù)元素的平差問題。但以上方法都沒有把系數(shù)矩陣和觀測(cè)向量(或矩陣)中的所有元素納入到一個(gè)協(xié)因數(shù)陣中進(jìn)行處理。

    基于以上問題,本文推導(dǎo)了顧及系數(shù)矩陣含有常數(shù)元素情況下的多元總體最小二乘的牛頓解法,并討論了多元總體最小二乘的牛頓法解和拉格朗日乘數(shù)法解之間的關(guān)系以及給出了多元總體最小二乘的精度評(píng)定公式,最后通過模擬數(shù)據(jù)和實(shí)測(cè)數(shù)據(jù)驗(yàn)證了本文算法的可行性和有效性。

    1多元總體最小二乘解法

    1.1多元變量EIV模型

    多元變量EIV模型[6-8]為

    Y=(A+EA)Ξ-EY

    (1)

    (2)

    顧及系數(shù)矩陣A中某列含有常數(shù)元素,在文獻(xiàn)[8,17—18]的基礎(chǔ)上,系數(shù)矩陣和相應(yīng)的協(xié)因數(shù)陣定義為

    (3)

    (4)

    式中,As是系數(shù)矩陣A的隨機(jī)元素部分;Ad是系數(shù)矩陣A的常數(shù)元素部分;Es為隨機(jī)元素的誤差矩陣;Qi(i=1,2,…,t)∈Rn×n(t+m)。

    此時(shí),協(xié)因數(shù)陣Q為奇異陣,但在計(jì)算時(shí)并沒有直接求取Q的逆,因此并不影響解算過程。

    1.2多元總體最小二乘問題的牛頓解法

    本文在文獻(xiàn)[10]推導(dǎo)牛頓法思路的基礎(chǔ)上,進(jìn)一步推導(dǎo)了一般情形下多元總體最小二乘問題的牛頓解法。根據(jù)總體最小二乘問題平差準(zhǔn)則VTPV=min,由求條件極值的拉格朗日乘數(shù)法構(gòu)造如下目標(biāo)函數(shù)

    Γ(V,Ξ,λ)=VTPV+2λT(vec(Y+EY)-

    (ΞT?In)vec(A+EA))

    (5)

    式中,λ∈Rnm×1為拉格朗日乘數(shù)向量;?為Kronecker-Zehfuss積。

    對(duì)式(5)中的元素分別求導(dǎo)得

    (6)

    (7)

    由式(6)和式(7),得到改正數(shù)最優(yōu)估值為[10]

    (8)

    進(jìn)而得到牛頓法求解無約束問題的目標(biāo)函數(shù)

    (9)

    式中

    式(9)對(duì)vec(Ξ)求一階導(dǎo)數(shù),得

    (10)

    經(jīng)推導(dǎo)可得式(10)第1部分

    (11)

    由逆矩陣求導(dǎo)原理,得式(10)第2部分

    (12)

    式中

    (13)

    (14)

    由式(11)和式(12)得

    (15)

    式(9)對(duì)vec(Ξ)求二階導(dǎo)數(shù),得

    (16)

    并令

    (17)

    (18)

    依據(jù)牛頓法原理得到未知參數(shù)的改正數(shù)

    (19)

    以及

    (20)

    多元總體最小二乘問題的牛頓解法的迭代步驟為(k為迭代次數(shù)):

    作為迭代初始值。

    第2步,由上一步獲得的數(shù)值,分別計(jì)算

    由式(15)和式(18),進(jìn)一步計(jì)算Gk+1、Hk+1和

    第3步,根據(jù)式(20)更新未知參數(shù)估值。

    作為未知參數(shù)的最佳估值。

    1.3牛頓法解和拉格朗日乘數(shù)法解之間的關(guān)系

    拉格朗日乘數(shù)法的多元總體最小二乘解為[10]

    (21)

    式(21)在近似值處展開,寫成式(20)的形式

    (22)

    式中,拉格朗日乘數(shù)法未知參數(shù)改正數(shù)為

    (23)

    (24)

    利用矩陣反演公式[10,19],得

    (25)

    又因?yàn)?/p>

    (26)

    聯(lián)合式(23)、式(25)和式(26),得

    (27)

    2多元總體最小二乘解法精度評(píng)定

    2.1單位權(quán)方差估值公式

    根據(jù)文獻(xiàn)[6—7],多元總體最小二乘解法的單位權(quán)方差估值公式為

    (28)

    式中,m(n-t)為自由度,mn為觀測(cè)矩陣Y元素的個(gè)數(shù),mt為參數(shù)矩陣Ξ元素的個(gè)數(shù)。

    2.2協(xié)因數(shù)陣計(jì)算公式

    多元總體最小二乘法中,數(shù)據(jù)多以矩陣形式存在。為計(jì)算協(xié)因數(shù)陣,把矩陣寫成向量形式,得到

    (29)

    表1多元總體最小二乘法協(xié)因數(shù)陣

    Tab.1Cofactor matrice of multivariate total least squares adjustment

    L^V^L^ξLQ-^X5Q-^X5-Q^ΣT^XT4^V-^X5^X50Q^ΣT^XT4^LQ-^X50Q-^X50^ξ-^X4^ΣQ^X4^ΣQ0^X3

    3算例分析

    3.1模擬算例

    根據(jù)文獻(xiàn)[17]中給出的源坐標(biāo)系的真值數(shù)據(jù),本文通過給定的二維仿射變換參數(shù)真值計(jì)算得到目標(biāo)坐標(biāo)系的真值數(shù)據(jù)。在得到兩組坐標(biāo)系的真值數(shù)據(jù)后,依據(jù)文獻(xiàn)[17],本文在模擬數(shù)據(jù)中源坐標(biāo)系所加誤差的協(xié)因數(shù)陣為

    Qxy=I2?(0.01diag([1,2,3,1,5,4,2,7,2,1,8,3,6]))

    模擬數(shù)據(jù)中目標(biāo)坐標(biāo)系所加誤差的協(xié)因數(shù)陣為

    QXY=I2?(0.01diag([1,3,6,1,1,8,4,3,6,5,4,5,2]))

    共模擬500次,分別使用本文提出的牛頓法(N法)和拉格朗日-牛頓法(L-N法),以及Fang(2011)法[10]、Schaffrin(2009)法[8]和加權(quán)最小二乘法(WLS)計(jì)算二維仿射變換參數(shù),迭代法的終止條件均為10-12,將獲得的二維仿射變換參數(shù)的平均值,以及二維仿射變換參數(shù)平均值與其真值的差值2范數(shù)和平均迭代次數(shù)列于表2中。

    參數(shù)為矩陣時(shí),文獻(xiàn)[10]并沒有給出在系數(shù)矩陣含有常數(shù)列時(shí)的具體公式算法。依據(jù)其提出的在參數(shù)為向量時(shí)系數(shù)矩陣含有常數(shù)列的加權(quán)總體最小二乘算法,本文編程實(shí)現(xiàn)了這種將系數(shù)矩陣分成兩部分的多元加權(quán)總體最小二乘算法,并簡(jiǎn)稱為Fang(2011)法。拉格朗日-牛頓法(L-N法)是指把Fang(2011)法(拉格朗日乘數(shù)法)獲得的第一次多元加權(quán)總體最小二乘參數(shù)解作為本文牛頓法迭代的初始值;在選取迭代初始值時(shí),如果使用多元最小二乘參數(shù)解作為初始值計(jì)算式(20),則為牛頓法(N法),如果使用多元加權(quán)總體最小二乘參數(shù)解作為初始值計(jì)算式(20),則為拉格朗日-牛頓法(L-N法);在能夠獲得多元加權(quán)總體最小二乘參數(shù)解的情況下,建議使用收斂速度更快的L-N法。

    以上算例沒有考慮系數(shù)矩陣元素與觀測(cè)矩陣元素的相關(guān)性。為了驗(yàn)證本文的牛頓法能夠解決更一般的多元總體最小二乘問題,在文獻(xiàn)[19]的基礎(chǔ)上,使用Matlab 7.11.0軟件,利用函數(shù)randn(100,13*4)T×randn(100,13*4)生成協(xié)因數(shù)陣Qr,利用函數(shù)mvnrnd生成均值均為0,方差-協(xié)方差陣為0.001 Qr的隨機(jī)誤差矩陣,并將其施加到真值數(shù)據(jù)中以作為計(jì)算的模擬數(shù)據(jù)。

    由于模擬數(shù)據(jù)使得兩套坐標(biāo)系具有了相關(guān)性,此時(shí)Schaffrin(2009)法不再適用,選用本文提出的牛頓法和拉格朗日-牛頓法,以及Fang(2011)法模擬計(jì)算10 000次,計(jì)算時(shí)取Q=Qr,迭代法的終止條件均為10-12,并由獲得的10 000次的參數(shù)估值繪制出的參數(shù)估值頻率直方圖及其相應(yīng)的擬合出的正態(tài)分布曲線圖如圖1所示。

    計(jì)算結(jié)果表明,牛頓法的平均迭代次數(shù)為6.972 6,拉格朗日-牛頓法的平均迭代次數(shù)為6.828 7,F(xiàn)ang(2011)法的平均迭代次數(shù)為12.188 2。牛頓法和拉格朗日-牛頓法使用的協(xié)因數(shù)陣Q為奇異矩陣,但這并不影響未知參數(shù)的解算。Fang(2011)法的迭代次數(shù)最多,拉格朗日-牛頓法迭代次數(shù)要少于牛頓法,這也與上文的結(jié)論相符合。參數(shù)估值的頻率直方圖近似于正態(tài)分布,相應(yīng)的正態(tài)分布擬合曲線大約也以參數(shù)的真值為對(duì)稱軸。以上的結(jié)論證明了多元總體最小二乘牛頓解法的可行性和有效性。

    表2 5種方法參數(shù)估計(jì)結(jié)果

    圖1 所估參數(shù)的直方圖及其正態(tài)分布擬合曲線Fig.1 Histogram and fitting curve of normal distribution of estimated parameters

    3.2實(shí)測(cè)算例

    選取文獻(xiàn)[19]中二維坐標(biāo)轉(zhuǎn)換數(shù)據(jù)作為計(jì)算數(shù)據(jù),共有5組在源坐標(biāo)系和目標(biāo)坐標(biāo)系下的公共點(diǎn)。

    源坐標(biāo)系和目標(biāo)坐標(biāo)系權(quán)陣為[19]

    Pxy=diag([18.781 7,6.377 4,12.648 9,

    17.476 9,22.272 6,23.982 3,13.680 4,

    3.465 6,3.732 4,6.437 7])

    PXY=diag([9.836 1,5.535 7,12.736 9,

    12.009 9,10.181,11.366 1,11.147,

    5.883 4,9.832 2,7.567 8])

    使用本文提出的牛頓法和拉格朗日-牛頓法,以及Fang(2011)法計(jì)算二維仿射變換參數(shù),迭代法的終止條件均為10-12。3種方法的未知參數(shù)估值列于表3中。利用表3中的未知參數(shù)估值結(jié)果和表1中的公式計(jì)算未知參數(shù)估值協(xié)因數(shù)陣,并將其列于表4中。

    表3 3種方法參數(shù)估計(jì)結(jié)果

    表4 參數(shù)估值協(xié)因數(shù)陣

    4結(jié)論

    本文推導(dǎo)了多元總體最小二乘問題的牛頓解法,并比較了基于牛頓法的多元加權(quán)總體最小二乘解和基于拉格朗日乘數(shù)法的多元加權(quán)總體最小二乘解之間的關(guān)系,給出了多元總體最小二乘算法精度評(píng)定的16種協(xié)因數(shù)陣的近似計(jì)算公式。通過算例驗(yàn)證,并與現(xiàn)有的解法作了比較,得到以下幾點(diǎn)結(jié)論:

    (1) 牛頓法有著更少的迭代次數(shù),更高的解算效率。

    (2) 在觀測(cè)矩陣和系數(shù)矩陣元素對(duì)應(yīng)的協(xié)因數(shù)陣中,本文取常數(shù)元素的協(xié)因數(shù)為0。這雖然導(dǎo)致協(xié)因數(shù)陣Q奇異,但并不影響參數(shù)的估計(jì)結(jié)果;本文的算法也能解決觀測(cè)矩陣與系數(shù)矩陣元素具有相關(guān)性的問題。

    (3) 牛頓法的解算效率會(huì)受到迭代初始值的影響,本文采用的拉格朗日-牛頓法因其迭代初始值更接近參數(shù)最優(yōu)估值,而使得解算效率比起牛頓法有所提高。

    本文給出的多元總體最小二乘解法精度評(píng)定的公式只是近似公式,如何推導(dǎo)更為嚴(yán)密的總體最小二乘精度評(píng)定公式還需進(jìn)一步研究。

    參考文獻(xiàn):

    [1]王樂洋.基于總體最小二乘的大地測(cè)量反演理論及應(yīng)用研究[J]. 測(cè)繪學(xué)報(bào), 2012, 41(4): 629.

    WANG Leyang. Research on Theory and Application of Total Least Squares in Geodetic Inversion[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(4): 629.

    [2]GOLUB G H, VAN LOAN C F. An Analysis of the Total Least Squares Problem[J]. SIAM Journal on Numerical Analysis, 1980, 17(6): 883-893.

    [3]SPRENT P. Models in Regression and Related Topics[M]. London: Methuen & Co Ltd, 1969.

    [4]王樂洋, 許才軍. 總體最小二乘研究進(jìn)展[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2013, 38(7): 850-856, 878.

    WANG Leyang, XU Caijun. Progress in Total Least Squares[J]. Geomatics and Information Science of Wuhan University, 2013, 38(7): 850-856, 878.

    [5]VAN HUFFEL S, VANDEWALLE J. The Total Least Squares Problem: Computational Aspects and Analysis[D]. Philadelphia, Pennsylvania: SIAM, 1991.

    [6]SCHAFFRIN B, FELUS Y A. Multivariate Total Least-squares Adjustment for Empirical Affine Transformations[C]∥VI Hotine-Marussi Symposium on Theoretical and Computational Geodesy: Challenge and Role of Modern Geodesy. Berlin: Springer, 2008, 132: 238-242.

    [7]SCHAFFRIN B, FELUS Y A. On the Multivariate Total Least-squares Approach to Empirical Coordinate Transformations. Three Algorithms[J]. Journal of Geodesy, 2008, 82(6): 373-383.

    [8]SCHAFFRIN B, WIESER A. Empirical Affine Reference Frame Transformations by Weighted Multivariate TLS Adjustment[C]∥Geodetic Reference Frames. Berlin: Springer, 2009, 134: 213-218.

    [9]FELUS Y A, BURTCH R C. On Symmetrical Three-dimensional Datum Conversion[J]. GPS Solutions, 2009, 13(1): 65-74.

    [10]FANG X. Weighted Total Least Squares Solutions for Applications in Geodesy[D]. Hanover: Leibniz University of Hanover, 2011.

    [11]孫大雙, 張友陽, 黃令勇, 等. 多元總體最小二乘在大旋轉(zhuǎn)角三維坐標(biāo)轉(zhuǎn)換中的應(yīng)用[J]. 測(cè)繪科學(xué)技術(shù)學(xué)報(bào), 2014, 31(5): 481-485.

    SUN Dashuang, ZHANG Youyang, HUANG Lingyong, et al. Application of Multivariate Total Least Squares Method in Three-dimension Coordinate Conversion with Large Rotation Angle[J]. Journal of Geomatics Science and Technology, 2014, 31(5): 481-485.

    [12]黃令勇, 呂志平, 任雅奇, 等. 多元總體最小二乘在三維坐標(biāo)轉(zhuǎn)換中的應(yīng)用[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2014, 39(7): 793-798.

    HUANG Lingyong, Lü Zhiping, REN Yaqi, et al. Application of Multivariate Total Least Square in Three-dimensional Coordinate Transformation[J]. Geomatics and Information Science of Wuhan University, 2014, 39(7): 793-798.

    [13]錢承軍, 陳義. 多變量總體最小二乘在點(diǎn)云拼接中的應(yīng)用[J]. 測(cè)繪與空間地理信息, 2015, 38(1): 67-69, 76.

    QIAN Chengjun, CHEN Yi. Application of Multivariable Total Least Squares in the Registration of Point Clouds[J]. Geomatics & Spatial Information Technology, 2015, 38(1): 67-69, 76.

    [14]王樂洋, 許才軍, 魯鐵定. 邊長(zhǎng)變化反演應(yīng)變參數(shù)的總體最小二乘方法[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2010, 35(2): 181-184.

    WANG Leyang, XU Caijun, LU Tieding. Inversion of Strain Parameter Using Distance Changes Based on Total Least Squares[J]. Geomatics and Information Science of Wuhan University, 2010, 35(2): 181-184.

    [15]XU Caijun, WANG Leyang, WEN Yangmao, et al. Strain Rates in the Sichuan-Yunnan Region Based upon the Total Least Squares Heterogeneous Strain Model from GPS Data

    [J]. Terrestrial Atmospheric and Oceanic Sciences, 2011, 22(2): 133-147.

    [16]王樂洋, 于冬冬, 呂開云. 復(fù)數(shù)域總體最小二乘平差[J]. 測(cè)繪學(xué)報(bào), 2015, 44(8): 866-876. DOI: 10.11947/j.AGCS.2015.20130701.

    WANG Leyang, YU Dongdong, Lü Kaiyun. Complex Total Least Squares Adjustment[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(8): 866-876. DOI: 10.11947/j.AGCS.2015.20130701.

    [17]MAHBOUB V. On Weighted Total Least-squares for Geodetic Transformations[J]. Journal of Geodesy, 2012, 86(5): 359-367.

    [18]XU Peiliang,LIU Jingnan,SHI Chuang.Total Least Squares Adjustment in Partial Errors-in-variables Models: Algorithm and Statistical Analysis[J]. Journal of Geodesy, 2012, 86(8): 661-675.

    [19]JAZAERI S, AMIRI-SIMKOOEI A R, SHARIFI M A. Iterative Algorithm for Weighted Total Least Squares Adjustment[J]. Survey Review, 2014, 46(334): 19-27.

    [20]BOYD S, VANDENBERGHE L. Convex Optimization[M]. Cambridge: Cambridge University Press, 2004.

    [21]王樂洋, 魯鐵定. 總體最小二乘平差法的誤差傳播定律[J]. 大地測(cè)量與地球動(dòng)力學(xué), 2014, 34(2): 55-59.

    WANG Leyang, LU Tieding. Propagation Law of Errors in Total Least Squares Adjustment[J]. Journal of Geodesy and Geodynamics, 2014, 34(2): 55-59.

    [22]ZHOU Yongjun, KOU Xinjian, ZHU Jianjun, et al. A Newton Algorithm for Weighted Total Least-squares Solution to a Specific Errors-in-variables Model with Correlated Measurements[J]. Studia Geophysica et Geodaetica, 2014, 58(3): 349-375.

    [23]王樂洋, 許才軍, 魯鐵定. 病態(tài)加權(quán)總體最小二乘平差的嶺估計(jì)解法[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2010, 35(11): 1346-1350.

    WANG Leyang, XU Caijun, LU Tieding. Ridge Estimation Method in Ill-posed Weighted Total Least Squares Adjustment[J]. Geomatics and Information Science of Wuhan University, 2010, 35(11): 1346-1350.

    (責(zé)任編輯:叢樹平)

    修回日期: 2015-10-12

    First author: WANG Leyang (1983—), male, PhD, associate professor, majors in geodetic inversion and geodetic data processing.

    E-mail: wleyang@163.com

    A Newton Algorithm for Multivariate Total Least Squares Problems

    WANG Leyang1,2,3,ZHAO Yingwen1,CHEN Xiaoyong1,2,3,ZANG Deyan1,2,3

    1. Faculty of Geomatics, East China University of Technology, Nanchang 330013, China; 2. Key Laboratory of Watershed Ecology and Geographical Environment Monitoring, NASG,Nanchang 330013, China; 3. Jiangxi Province Key Lab for Digital Land, Nanchang 330013, China

    Abstract:In order to improve calculation efficiency of parameter estimation, an algorithm for multivariate weighted total least squares adjustment based on Newton method is derived. The relationship between the solution of this algorithm and that of multivariate weighted total least squares adjustment based on Lagrange multipliers method is analyzed. According to propagation of cofactor, 16 computational formulae of cofactor matrices of multivariate total least squares adjustment are also listed. The new algorithm could solve adjustment problems containing correlation between observation matrix and coefficient matrix. And it can also deal with their stochastic elements and deterministic elements with only one cofactor matrix. The results illustrate that the Newton algorithm for multivariate total least squares problems could be practiced and have higher convergence rate.

    Key words:multivariate total least squares; Newton method; propagation of cofactor; affine

    第一作者簡(jiǎn)介:王樂洋(1983—),男,博士,副教授,研究方向?yàn)榇蟮販y(cè)量反演及大地測(cè)量數(shù)據(jù)處理。

    收稿日期:2015-05-11

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(41204003;41161069;41304020;41464001);測(cè)繪地理信息公益性行業(yè)科研專項(xiàng)(201512026);江西省自然科學(xué)基金(20151BAB203042);江西省教育廳科技項(xiàng)目(GJJ150595;KJLD12077;KJLD14049);流域生態(tài)與地理環(huán)境監(jiān)測(cè)國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室開放基金(WE2015005);東華理工大學(xué)博士科研啟動(dòng)金(DHBK201113);江西省研究生創(chuàng)新專項(xiàng)資金(YC2015-S266;YC2015-S267);東華理工大學(xué)研究生創(chuàng)新專項(xiàng)資金(DHYC-2015005);對(duì)地觀測(cè)技術(shù)國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室開放基金(K201502)

    中圖分類號(hào):P207

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1001-1595(2016)04-0411-07

    引文格式:王樂洋,趙英文,陳曉勇,等.多元總體最小二乘問題的牛頓解法[J].測(cè)繪學(xué)報(bào),2016,45(4):411-417,424. DOI:10.11947/j.AGCS.2016.20150246.

    WANG Leyang,ZHAO Yingwen,CHEN Xiaoyong,et al.A Newton Algorithm for Multivariate Total Least Squares Problems[J]. Acta Geodaetica et Cartographica Sinica,2016,45(4):411-417,424. DOI:10.11947/j.AGCS.2016.20150246.

    色哟哟哟哟哟哟| 婷婷六月久久综合丁香| 成年女人永久免费观看视频| 这个男人来自地球电影免费观看| 国产单亲对白刺激| 亚洲在线观看片| 岛国在线观看网站| 欧美又色又爽又黄视频| 日韩欧美三级三区| 精品免费久久久久久久清纯| 国产亚洲欧美98| 中文字幕av在线有码专区| 一二三四在线观看免费中文在| av中文乱码字幕在线| 深夜精品福利| 国产淫片久久久久久久久 | 母亲3免费完整高清在线观看| 久久久国产成人精品二区| 中文字幕高清在线视频| 18禁观看日本| 久久久精品欧美日韩精品| 精品一区二区三区视频在线 | 全区人妻精品视频| 九九久久精品国产亚洲av麻豆 | 日日夜夜操网爽| 成人特级黄色片久久久久久久| 久久久成人免费电影| 国产精品久久久久久人妻精品电影| 欧美3d第一页| 欧美绝顶高潮抽搐喷水| 热99在线观看视频| 丝袜人妻中文字幕| 精品久久久久久久久久免费视频| 人妻夜夜爽99麻豆av| 丰满人妻一区二区三区视频av | 成人性生交大片免费视频hd| 久久久精品大字幕| 伊人久久大香线蕉亚洲五| 国产黄a三级三级三级人| 两性夫妻黄色片| 亚洲av成人精品一区久久| 熟妇人妻久久中文字幕3abv| 国产又色又爽无遮挡免费看| 日本 欧美在线| 亚洲av电影不卡..在线观看| 麻豆成人av在线观看| 久久久成人免费电影| 日本 欧美在线| 国产精品永久免费网站| 午夜福利在线观看免费完整高清在 | 男女那种视频在线观看| 免费无遮挡裸体视频| 国产一区二区在线观看日韩 | 真人一进一出gif抽搐免费| 最近最新中文字幕大全电影3| 全区人妻精品视频| 欧美一区二区国产精品久久精品| 国产成人一区二区三区免费视频网站| 国产亚洲精品一区二区www| 日本五十路高清| 国产精品爽爽va在线观看网站| 国产精品美女特级片免费视频播放器 | avwww免费| 91久久精品国产一区二区成人 | 精品久久久久久久久久免费视频| 日韩免费av在线播放| 国产成人福利小说| 一进一出抽搐动态| 熟妇人妻久久中文字幕3abv| 午夜福利欧美成人| 十八禁人妻一区二区| 久久午夜综合久久蜜桃| 一区二区三区高清视频在线| 色吧在线观看| 可以在线观看的亚洲视频| 国产av麻豆久久久久久久| 国产精品综合久久久久久久免费| 嫩草影院精品99| 亚洲欧美激情综合另类| 国产激情欧美一区二区| 伊人久久大香线蕉亚洲五| 舔av片在线| 美女高潮的动态| 国产午夜精品论理片| 男女午夜视频在线观看| 老鸭窝网址在线观看| 国产综合懂色| 欧美+亚洲+日韩+国产| 中文字幕熟女人妻在线| 午夜精品久久久久久毛片777| 国产探花在线观看一区二区| 男女午夜视频在线观看| 亚洲国产欧美一区二区综合| 欧美国产日韩亚洲一区| 1000部很黄的大片| 久久久久久大精品| 欧美3d第一页| 全区人妻精品视频| 亚洲性夜色夜夜综合| xxxwww97欧美| 99视频精品全部免费 在线 | 欧美激情久久久久久爽电影| 亚洲av美国av| 久久久久性生活片| 婷婷精品国产亚洲av| 国产蜜桃级精品一区二区三区| 亚洲精品久久国产高清桃花| 成人精品一区二区免费| 最好的美女福利视频网| 18美女黄网站色大片免费观看| 最近最新中文字幕大全电影3| av在线蜜桃| 麻豆国产97在线/欧美| 免费在线观看成人毛片| 又紧又爽又黄一区二区| 中文字幕人妻丝袜一区二区| 夜夜夜夜夜久久久久| 岛国在线观看网站| 99在线人妻在线中文字幕| 午夜激情福利司机影院| 一个人看的www免费观看视频| 日韩欧美三级三区| 欧美丝袜亚洲另类 | 黄片大片在线免费观看| 久久久久久久久中文| 国产精品久久久久久亚洲av鲁大| 午夜福利在线观看吧| 亚洲av成人精品一区久久| 999久久久国产精品视频| 极品教师在线免费播放| 成人亚洲精品av一区二区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av免费在线观看| 国产淫片久久久久久久久 | 亚洲精品456在线播放app | 一进一出抽搐动态| 婷婷精品国产亚洲av| 国产高清视频在线观看网站| www国产在线视频色| 变态另类成人亚洲欧美熟女| 舔av片在线| 精品久久久久久成人av| 中文亚洲av片在线观看爽| 国产av不卡久久| 老司机午夜十八禁免费视频| 欧美日韩瑟瑟在线播放| 久久久久久久久久黄片| 亚洲精品456在线播放app | 99久久国产精品久久久| 国产精品99久久99久久久不卡| 日本精品一区二区三区蜜桃| 非洲黑人性xxxx精品又粗又长| 一个人看视频在线观看www免费 | 日本一二三区视频观看| 国产午夜精品论理片| 久久这里只有精品中国| 三级毛片av免费| 欧美激情久久久久久爽电影| 美女高潮喷水抽搐中文字幕| 国产单亲对白刺激| 99久久久亚洲精品蜜臀av| 欧美性猛交╳xxx乱大交人| 欧美另类亚洲清纯唯美| 999久久久精品免费观看国产| 欧美性猛交黑人性爽| 老司机午夜福利在线观看视频| 国产精品亚洲一级av第二区| 国产精品精品国产色婷婷| 少妇的丰满在线观看| 一二三四在线观看免费中文在| 国产黄a三级三级三级人| 国产单亲对白刺激| 999久久久国产精品视频| 日日干狠狠操夜夜爽| 国产aⅴ精品一区二区三区波| 欧美一区二区精品小视频在线| 成人特级av手机在线观看| 首页视频小说图片口味搜索| 成人国产综合亚洲| 长腿黑丝高跟| 精品午夜福利视频在线观看一区| 男人舔女人的私密视频| 久久久久九九精品影院| 不卡av一区二区三区| 色吧在线观看| 亚洲av电影不卡..在线观看| 神马国产精品三级电影在线观看| 成年人黄色毛片网站| 国产亚洲精品综合一区在线观看| 亚洲激情在线av| 午夜激情欧美在线| 夜夜看夜夜爽夜夜摸| 极品教师在线免费播放| 一本精品99久久精品77| 美女免费视频网站| 无人区码免费观看不卡| 亚洲午夜精品一区,二区,三区| 亚洲av美国av| 亚洲国产精品久久男人天堂| 亚洲自拍偷在线| 国产成人影院久久av| netflix在线观看网站| 久久这里只有精品中国| 性欧美人与动物交配| 老司机在亚洲福利影院| 亚洲性夜色夜夜综合| 一本精品99久久精品77| 国产一区在线观看成人免费| 麻豆成人午夜福利视频| 亚洲色图av天堂| 九色成人免费人妻av| 久久久久国产一级毛片高清牌| 动漫黄色视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 香蕉久久夜色| 日本精品一区二区三区蜜桃| 欧美成人免费av一区二区三区| 悠悠久久av| 日日干狠狠操夜夜爽| 欧美最黄视频在线播放免费| 亚洲自拍偷在线| 久久九九热精品免费| 欧美绝顶高潮抽搐喷水| 国产成人福利小说| 亚洲国产中文字幕在线视频| 少妇人妻一区二区三区视频| 一进一出抽搐动态| 日本黄大片高清| 日本a在线网址| 午夜精品一区二区三区免费看| 97人妻精品一区二区三区麻豆| 欧美乱妇无乱码| 免费在线观看视频国产中文字幕亚洲| 欧美性猛交黑人性爽| 一本精品99久久精品77| 最新中文字幕久久久久 | 精品不卡国产一区二区三区| 久久中文字幕人妻熟女| 可以在线观看毛片的网站| 国产精品香港三级国产av潘金莲| 亚洲一区高清亚洲精品| 999久久久国产精品视频| 中亚洲国语对白在线视频| 91九色精品人成在线观看| 欧美黑人欧美精品刺激| 欧美极品一区二区三区四区| 国产精品影院久久| 亚洲成人久久爱视频| 小说图片视频综合网站| 男女之事视频高清在线观看| 色综合婷婷激情| 久久这里只有精品中国| 亚洲在线观看片| 欧美日本亚洲视频在线播放| 青草久久国产| 黑人巨大精品欧美一区二区mp4| 91av网站免费观看| 中文字幕高清在线视频| 19禁男女啪啪无遮挡网站| 欧美中文综合在线视频| 最新中文字幕久久久久 | 亚洲自偷自拍图片 自拍| 久久久水蜜桃国产精品网| 亚洲男人的天堂狠狠| 此物有八面人人有两片| 国模一区二区三区四区视频 | 嫩草影院入口| 免费看光身美女| 色综合站精品国产| 可以在线观看毛片的网站| 久久天堂一区二区三区四区| 法律面前人人平等表现在哪些方面| 欧美又色又爽又黄视频| 久久久色成人| 亚洲 欧美 日韩 在线 免费| 在线观看午夜福利视频| 免费无遮挡裸体视频| 偷拍熟女少妇极品色| 精品国内亚洲2022精品成人| 香蕉av资源在线| 精品日产1卡2卡| 国产精品影院久久| 美女 人体艺术 gogo| 国产极品精品免费视频能看的| 一个人观看的视频www高清免费观看 | 免费看a级黄色片| 男人和女人高潮做爰伦理| 国产乱人伦免费视频| 免费看美女性在线毛片视频| 韩国av一区二区三区四区| 亚洲精品美女久久av网站| 成人18禁在线播放| 久久精品国产亚洲av香蕉五月| 亚洲一区高清亚洲精品| 午夜免费成人在线视频| 操出白浆在线播放| 日本免费一区二区三区高清不卡| 亚洲国产色片| 三级国产精品欧美在线观看 | 波多野结衣高清无吗| 美女cb高潮喷水在线观看 | 欧美色视频一区免费| 国产视频一区二区在线看| 亚洲七黄色美女视频| 在线国产一区二区在线| 露出奶头的视频| 午夜精品一区二区三区免费看| 国产精品久久久久久精品电影| 欧美成人免费av一区二区三区| 最好的美女福利视频网| 99久久国产精品久久久| 成人av一区二区三区在线看| 久久久国产成人免费| 亚洲国产高清在线一区二区三| 啪啪无遮挡十八禁网站| 欧美激情久久久久久爽电影| 亚洲七黄色美女视频| 欧美乱码精品一区二区三区| 欧美性猛交黑人性爽| 99热精品在线国产| 老汉色av国产亚洲站长工具| 亚洲人成网站高清观看| 天堂动漫精品| 精品一区二区三区av网在线观看| 欧美黑人欧美精品刺激| 国内精品久久久久精免费| 免费av毛片视频| 欧美极品一区二区三区四区| av欧美777| 成年女人毛片免费观看观看9| av中文乱码字幕在线| 99久久国产精品久久久| 神马国产精品三级电影在线观看| 欧美色欧美亚洲另类二区| 大型黄色视频在线免费观看| 操出白浆在线播放| 成人国产一区最新在线观看| 非洲黑人性xxxx精品又粗又长| 久久久久久九九精品二区国产| 午夜福利视频1000在线观看| 亚洲电影在线观看av| 丝袜人妻中文字幕| 亚洲电影在线观看av| 18禁裸乳无遮挡免费网站照片| www日本在线高清视频| 亚洲精品色激情综合| 国产成人啪精品午夜网站| 亚洲美女黄片视频| 国产日本99.免费观看| 亚洲av成人不卡在线观看播放网| 久久这里只有精品19| 国产乱人视频| 国产成人一区二区三区免费视频网站| 99国产极品粉嫩在线观看| 国产成人av教育| 午夜视频精品福利| 三级毛片av免费| 伊人久久大香线蕉亚洲五| 亚洲成人久久爱视频| 日韩免费av在线播放| 精品熟女少妇八av免费久了| 日韩精品青青久久久久久| 日韩有码中文字幕| 免费在线观看影片大全网站| 成人av在线播放网站| 色哟哟哟哟哟哟| 欧美色视频一区免费| 老汉色∧v一级毛片| 亚洲成av人片免费观看| 国产精品99久久99久久久不卡| 老鸭窝网址在线观看| 国产精品99久久99久久久不卡| 久久亚洲真实| 一个人看的www免费观看视频| 嫩草影院入口| 日本黄色视频三级网站网址| av女优亚洲男人天堂 | 国产亚洲欧美98| 村上凉子中文字幕在线| 精品午夜福利视频在线观看一区| 一级毛片高清免费大全| 一个人免费在线观看的高清视频| 国产成+人综合+亚洲专区| 极品教师在线免费播放| 亚洲欧美日韩东京热| 日本a在线网址| 岛国在线观看网站| 成人无遮挡网站| 三级男女做爰猛烈吃奶摸视频| 国模一区二区三区四区视频 | aaaaa片日本免费| 婷婷精品国产亚洲av在线| 中文字幕精品亚洲无线码一区| 欧美又色又爽又黄视频| 老司机午夜福利在线观看视频| 中文字幕人成人乱码亚洲影| 欧美日韩亚洲国产一区二区在线观看| 丰满的人妻完整版| www.精华液| 日韩欧美一区二区三区在线观看| 亚洲在线自拍视频| 巨乳人妻的诱惑在线观看| 色综合欧美亚洲国产小说| 国产又色又爽无遮挡免费看| 亚洲天堂国产精品一区在线| 久久精品影院6| 成人一区二区视频在线观看| 亚洲性夜色夜夜综合| 特级一级黄色大片| 亚洲精华国产精华精| 亚洲成人免费电影在线观看| 麻豆国产av国片精品| 好男人电影高清在线观看| 国产日本99.免费观看| 亚洲欧美日韩无卡精品| 欧美中文综合在线视频| 深夜精品福利| 欧美高清成人免费视频www| 国产精品女同一区二区软件 | 毛片女人毛片| 97碰自拍视频| 精品午夜福利视频在线观看一区| 国产欧美日韩一区二区三| 真实男女啪啪啪动态图| 国产黄a三级三级三级人| 在线国产一区二区在线| 少妇裸体淫交视频免费看高清| 757午夜福利合集在线观看| 中亚洲国语对白在线视频| 全区人妻精品视频| 狂野欧美白嫩少妇大欣赏| 国产男靠女视频免费网站| 国产单亲对白刺激| 精品国产美女av久久久久小说| 久久国产精品人妻蜜桃| 视频区欧美日本亚洲| 久久久久免费精品人妻一区二区| 51午夜福利影视在线观看| 国产精品99久久99久久久不卡| 嫁个100分男人电影在线观看| 国内精品美女久久久久久| 亚洲中文字幕日韩| 精品久久久久久久末码| 999久久久精品免费观看国产| 麻豆成人午夜福利视频| 国产av在哪里看| 少妇的丰满在线观看| 少妇丰满av| 国产高清videossex| 欧美黑人欧美精品刺激| 91在线观看av| 亚洲国产精品成人综合色| 日韩中文字幕欧美一区二区| 国产成人一区二区三区免费视频网站| 欧美在线黄色| 国产 一区 欧美 日韩| 床上黄色一级片| 国产亚洲精品一区二区www| 99久久无色码亚洲精品果冻| 色哟哟哟哟哟哟| 国产成人av教育| 精华霜和精华液先用哪个| 法律面前人人平等表现在哪些方面| 亚洲av美国av| 亚洲国产日韩欧美精品在线观看 | 91麻豆av在线| 国产精品一及| 亚洲中文日韩欧美视频| 日韩av在线大香蕉| 1000部很黄的大片| 亚洲精品456在线播放app | 国产高清三级在线| 90打野战视频偷拍视频| 亚洲国产精品合色在线| 成年人黄色毛片网站| 精品人妻1区二区| 老司机福利观看| 黄片小视频在线播放| 波多野结衣巨乳人妻| 狠狠狠狠99中文字幕| 黑人巨大精品欧美一区二区mp4| 在线观看日韩欧美| 亚洲美女黄片视频| 制服人妻中文乱码| 淫秽高清视频在线观看| 精品福利观看| 久99久视频精品免费| 午夜福利18| 99久久成人亚洲精品观看| 国产欧美日韩精品亚洲av| 亚洲成av人片在线播放无| 69av精品久久久久久| 一区二区三区激情视频| 一级作爱视频免费观看| 国产一区在线观看成人免费| 人妻丰满熟妇av一区二区三区| 免费看美女性在线毛片视频| 午夜两性在线视频| 亚洲成a人片在线一区二区| 黄色女人牲交| 日日摸夜夜添夜夜添小说| 亚洲av第一区精品v没综合| 日本黄色视频三级网站网址| 日本黄色片子视频| 午夜免费观看网址| 欧美高清成人免费视频www| 天天躁日日操中文字幕| 午夜福利高清视频| 午夜福利成人在线免费观看| 久久久水蜜桃国产精品网| 日本熟妇午夜| 香蕉国产在线看| 精品久久久久久,| xxx96com| 叶爱在线成人免费视频播放| 全区人妻精品视频| 桃色一区二区三区在线观看| 18禁国产床啪视频网站| 色综合婷婷激情| 最近最新免费中文字幕在线| 日韩免费av在线播放| 一卡2卡三卡四卡精品乱码亚洲| 中文字幕人妻丝袜一区二区| 免费看光身美女| 精品欧美国产一区二区三| 亚洲成人精品中文字幕电影| 亚洲av成人精品一区久久| 免费观看人在逋| 久久精品91蜜桃| 欧美色欧美亚洲另类二区| 久久天堂一区二区三区四区| 欧美乱码精品一区二区三区| 嫁个100分男人电影在线观看| av中文乱码字幕在线| 国产人伦9x9x在线观看| 中文字幕熟女人妻在线| 国产亚洲精品久久久久久毛片| 精品无人区乱码1区二区| 国产蜜桃级精品一区二区三区| 一区二区三区国产精品乱码| 天天躁日日操中文字幕| 波多野结衣高清无吗| 精品久久久久久久久久久久久| 舔av片在线| 亚洲成人免费电影在线观看| 国产单亲对白刺激| 天堂√8在线中文| 亚洲午夜精品一区,二区,三区| 操出白浆在线播放| 黄色女人牲交| 91九色精品人成在线观看| 美女cb高潮喷水在线观看 | 男女那种视频在线观看| 性色av乱码一区二区三区2| 成年女人看的毛片在线观看| 女人被狂操c到高潮| 日韩欧美三级三区| 在线十欧美十亚洲十日本专区| 国产淫片久久久久久久久 | 黑人巨大精品欧美一区二区mp4| 日本撒尿小便嘘嘘汇集6| 成人鲁丝片一二三区免费| 99热只有精品国产| 99国产精品一区二区三区| 别揉我奶头~嗯~啊~动态视频| 国产一区二区在线观看日韩 | 国产精品一区二区三区四区免费观看 | 欧美3d第一页| 长腿黑丝高跟| 黄色日韩在线| 美女高潮的动态| 亚洲av成人av| 免费无遮挡裸体视频| 非洲黑人性xxxx精品又粗又长| 亚洲激情在线av| 熟妇人妻久久中文字幕3abv| 别揉我奶头~嗯~啊~动态视频| 草草在线视频免费看| 一区二区三区国产精品乱码| 午夜免费激情av| 在线视频色国产色| 可以在线观看毛片的网站| 午夜福利在线观看吧| 久久久国产欧美日韩av| 啦啦啦韩国在线观看视频| 最近最新免费中文字幕在线| 午夜影院日韩av| 99热精品在线国产| 久久精品91无色码中文字幕| 亚洲乱码一区二区免费版| 脱女人内裤的视频| 一个人观看的视频www高清免费观看 | 18美女黄网站色大片免费观看| 亚洲美女黄片视频| 啦啦啦韩国在线观看视频| 欧美xxxx黑人xx丫x性爽| 精品久久久久久,| 日日摸夜夜添夜夜添小说| 99riav亚洲国产免费| 亚洲欧美日韩无卡精品| 日日摸夜夜添夜夜添小说| 国产成人影院久久av| 啦啦啦免费观看视频1| 久久久久久久久中文| 给我免费播放毛片高清在线观看| 一a级毛片在线观看| 亚洲国产日韩欧美精品在线观看 | 成人欧美大片| 国产精品久久视频播放| 国产成人系列免费观看| 日日摸夜夜添夜夜添小说| 嫁个100分男人电影在线观看| 色噜噜av男人的天堂激情| 久久精品国产综合久久久| 久9热在线精品视频|